- 1College of Horticultural Science and Engineering, Shandong Agricultural University, Tai'an, China
- 2State Key Laboratory of Crop Biology, Tai'an, China
- 3Scientific Observing and Experimental Station of Environment Controlled Agricultural Engineering in Huang-Huai-Hai Region, Ministry of Agriculture, Beijing, China
- 4College of Horticulture, Henan Agricultural University, Zhengzhou, China
The atmospheric vapour pressure deficit (VPD) has been demonstrated to be a significant environmental factor inducing plant water stress and affecting plant photosynthetic productivity. Despite this, the rate-limiting step for photosynthesis under varying VPD is still unclear. In the present study, tomato plants were cultivated under two contrasting VPD levels: high VPD (3–5 kPa) and low VPD (0.5–1.5 kPa). The effect of long-term acclimation on the short-term rapid VPD response was examined across VPD ranging from 0.5 to 4.5 kPa. Quantitative photosynthetic limitation analysis across the VPD range was performed by combining gas exchange and chlorophyll fluorescence. The potential role of abscisic acid (ABA) in mediating photosynthetic carbon dioxide (CO2) uptake across a series of VPD was evaluated by physiological and transcriptomic analyses. The rate-limiting step for photosynthetic CO2 utilisation varied with VPD elevation in tomato plants. Under low VPD conditions, stomatal and mesophyll conductance was sufficiently high for CO2 transport. With VPD elevation, plant water stress was gradually pronounced and triggered rapid ABA biosynthesis. The contribution of stomatal and mesophyll limitation to photosynthesis gradually increased with an increase in the VPD. Consequently, the low CO2 availability inside chloroplasts substantially constrained photosynthesis under high VPD conditions. The foliar ABA content was negatively correlated with stomatal and mesophyll conductance for CO2 diffusion. Transcriptomic and physiological analyses revealed that ABA was potentially involved in mediating water transport and photosynthetic CO2 uptake in response to VPD variation. The present study provided new insights into the underlying mechanism of photosynthetic depression under high VPD stress.
Introduction
Carbon dioxide (CO2) is significant for plant photosynthesis, growth, and yield production. Although CO2 fertilisation and globally elevated trends are expected to improve crop photosynthesis and yield, large evidence has shown that the magnitude of such enhancement is constrained by other climate change-derived phenomena, such as more extreme and frequent environmental stress (Norby, 2002). The bottlenecks constraining the CO2 utilisation efficiency are limited CO2 acquisition and assimilation. It has been recognised that CO2 movement and carbon fixation are regulated by environmental factors. There is increasing evidence from physiology and crop production that high vapour pressure deficit (VPD) induces plant water stress and inhibits photosynthetic productivity (Lu et al., 2015; Zhang et al., 2015). Few previous studies have quantitatively addressed the components of photosynthetic limitation across a series of VPD. The rate-limiting step for photosynthetic CO2 transport and utilisation under different VPD conditions was highly uncertain. A quantitative limitation analysis consisting of stomatal, mesophyll, and biochemical limitations is essential to reveal the underlying mechanism by which the VPD affects the photosynthetic process.
Photosynthetic CO2 uptake and transport are constrained by a series of resistances, which have been simplified into stomatal and mesophyll resistance (Tholen and Zhu, 2011). Guard cells of stomata are the first barrier for gas exchange and modulate photosynthetic CO2 uptake and transpiration (Lawson and Blatt, 2014). Large evidence has shown that CO2 movement from the substomatal cavity to the carbon fixation site is constrained by great mesophyll resistance (Niinemets et al., 2009; von Caemmerer and Evans, 2010; Flexas et al., 2012; Kaldenhoff, 2012; Sharkey, 2012; Li et al., 2019b). In addition to stomatal resistance, mesophyll resistance also substantially constrains the photosynthetic rate, especially for C3 plants. Environmental fluctuations are thought to profoundly affect CO2 uptake and transport. Leaf anatomical properties determine the maximum potential conductance for gas or liquid phase diffusion. Some studies attributed photosynthetic limitations to anatomical adaptations under high VPD stress, such as reduced stomatal size, stomatal density, vein density, and mesophyll surface area (Fanourakis et al., 2016, 2020; Du et al., 2019). CO2 uptake and water loss share some common pathways, such as stomatal and intercellular spaces. Anatomical adaptations prevent excessive water loss, which simultaneously increases diffusion resistance for CO2 uptake. In addition to the anatomical determination over long-term adaptations, much evidence has shown that stomatal and mesophyll conductance respond rapidly and sensitively to external environmental variation (Xiong et al., 2015; Li et al., 2019a,b). The field and greenhouse VPD fluctuate dramatically over the diurnal course, which significantly affects the photosynthetic process. However, less attention has been given to reveal the mechanism of the rapid response of CO2 diffusion conductance across a series of VPD.
It has been widely reported that the plant hormone abscisic acid (ABA) is involved in various abiotic stresses and acts as a signalling molecule in response to drought, salinization, heat, and so on (Fang et al., 2019). Cellular ABA accumulation is an important dehydration-sensing and water balance-maintaining mechanism, which has special implications in stomatal closure and the decline of hydraulic conductance (Sack et al., 2018). ABA prevents excessive water loss and enhances crop drought tolerance by signalling pathways. As a C3 plant species, the photosynthetic and yield potential of tomato plants is greatly limited by the low CO2 availability inside chloroplasts. The excessive evaporative demand under high VPD exceeds root water uptake capacity and triggers plant water deficit in tomato plants, which contributes to photosynthetic depression and yield loss (Zhang et al., 2017, 2018; Li et al., 2019b). We hypothesised that ABA plays a significant role in preventing transpiration under high VPD-induced plant water deficit, which simultaneously constrains photosynthetic CO2 uptake and acquisition. Identifying the rate-limiting step for photosynthetic CO2 acquisition under contrasting VPD and revealing the mechanism has significant implications for both basic plant sciences and crop production.
To investigate the effect of long-term acclimation on the short-term rapid VPD response, the implications of leaf anatomical properties and ABA in modulating CO2 transport across a series of VPD ranges were addressed by physiological and transcriptomic analyses. Three questions were addressed in the present study: (1) How did stomatal and mesophyll conductance tune with the VPD? (2) How did the contribution of stomatal and mesophyll limitation to photosynthesis vary with the VPD? (3) How did ABA tune with the VPD and correlate with stomatal and mesophyll conductance?
Materials and Methods
Plant Materials and Growth Conditions
The experiment was conducted in two environmentally controlled greenhouses with the same characteristics (15 m in length, 10 m in width, and 3.5 m in height, north-south oriented) under spring-summer climatic conditions from May to August 2018. Two widely grown tomato cultivars (JinPeng NO.1, CV1, JinPeng & Co., Ltd., China; FenGuan, CV2, ZhongYa & Co., Ltd., China) with relatively distinct VPD responses were examined (Du et al., 2020). Seeds were sown in plugs for germination and transplanted at the four-leaf stage to 4.5 L plastic pots containing the same amount of organic substrate and perlite mixture in a 3:1 proportion (v/v). Soil moisture was maintained at ~90% container capacity according to a previous method (Zhang et al., 2017). Plants were periodically trimmed to maintain rapid vegetative growth throughout experiments. Plants were grown in two environmentally controlled greenhouses and maintained under the same growth conditions but contrasting VPD. A high VPD was achieved in a natural greenhouse environment, with a VPD of ~3–5 kPa around midday, while low VPD was maintained in the range of.5–1.5 kPa by humidification. A high-pressure micro-fog system was activated when the VPD exceeded the target values, and the characteristics of the system were described in detail in a previous study by Zhang et al. (2018). The average daily meteorological data inside the greenhouse during the growth period were ~maintained at a temperature of 20–32°C, relative humidity of 50–75%, and photosynthetically active radiation of 45–65 Wm−2.
The effects of VPD perturbations on leaf photosynthetic performance and plant water status were investigated ~50 days after treatments. Afterward, 15 uniform plants from each treatment were selected as samples and transferred to growth cabinets in the evening prior to photosynthetic measurements. The light and temperature of the growth cabinets were controlled steadily at normal levels throughout the experiment.
Leaf Gas Exchange and Chlorophyll Fluorescence
Leaf gas exchange and chlorophyll fluorescence were measured simultaneously on healthy and expanded leaflets at the same nodes by portable gas exchange systems equipped with a leaf chamber fluorometer (LI-6400, Li-Cor, Inc., Lincoln, NE, USA). All portable gas exchange systems were enclosed in growth cabinets. The VPD inside cabinets and the leaf chamber was simultaneously controlled across a series gradient of 0.5, 1.5, 2.5, 3.5, and 4.5 kPa. The temperature, light, and CO2 concentrations were controlled at the following constant and steady conditions throughout the experiment: temperature of 28 ± 1°C; saturating photosynthetic photon flux density (PPFD) of 1,100 μmol m−2 s−1; CO2 concentration of 400 μmol mol−1. The VPD was increased stepwise across the gradients for at least 60 min until photosynthesis and the plant water status achieved a new steady state.
The curve of the photosynthetic rate (Pn) vs. intercellular CO2 concentration (Ci) was determined using a previous procedure (Li et al., 2019b), across the VPD range of 0.5–4.5 kPa. Briefly, a Pn-Ci curve was generated by controlling the ambient CO2 concentration (Ca) from 400 to 300, 200, 150, 100, and 50 μmol mol−1 and then increased to 400 μmol mol−1. After re-achieving a steady-state at 400 μmol mol−1, Ca was increased gradually from 400 μmol mol−1 to 1,200 μmol mol−1. The carboxylation efficiency (CE) was estimated according to linear regression of the Pn-Ci curve in the range of Ca ≤ 200 μmol mol−1 (Sun et al., 2016). The maximum rate of Rubisco carboxylation capacity (Vcmax) and maximal rate of electron transport (Jmax) were determined according to the FvCB model (Farquhar et al., 1980).
Estimation of Photosynthetic CO2 Diffusion Conductance
Carbon dioxide diffuses via stomatal and mesophyll barriers in a series circuit, which was driven by the CO2 partial pressure gradient (Li et al., 2019b). Stomatal conductance for CO2 diffusion (gsc) was determined according to the water diffusion conductance (gsw) and the ratio between molecular diffusivities of water and CO2 in gas (Giuliani et al., 2013). The mesophyll conductance (gm) was estimated by the variable J method (Harley et al., 1992):
where Pn is the net photosynthetic rate and Ci is the intercellular CO2 concentration. Pn and Ci were measured by steady-state gas exchange. Rd is the mitochondrial respiration rate in the light, and Γ* is the CO2 compensation point inside the chloroplast. Rd and Γ* were calculated according to a previous study by Laisk and Oja (1998). Briefly, Pn-Ci curves were measured at two light intensities (75 and 500 μmol m−2 s−1) at CO2 concentrations of 30–120 μmol mol−1. Γ* (x-axis) and Rd (y-axis) were derived according to the intersection point of the Pn-Ci curves. J is the electron transport rate, which was calculated as described by a previous study (Tomas et al., 2013).
According to the series circuit, the total CO2 diffusion resistance (1/gtot) can be determined as 1/gtot = 1/gs + 1/gm (Niinemets et al., 2009). Therefore, gtot can be determined as:
Partitioning of the Photosynthetic Limitation
The photosynthetic limitation was divided into the components of stomatal limitation (Ls), mesophyll limitation (Lm), and biochemical limitation (Lm). The proportions of individual components imposed on photosynthesis were determined as follows (Muir et al., 2014; Li et al., 2019b):
By definition, Ls + Lm + Lb =1; ∂A/∂Cc was determined as the slope of the Pn-Cc curves at CO2 concentrations of 40–110 μmol mol−1.
Determination of the Plant Water Status
Once photosynthetic measurements were completed, the adjacent leaflets were harvested for the determination of the water status of the plant. The leaf water potential (Ψleaf) was measured by a pressure chamber (PMS-1000, PMS Instruments Inc., Corvallis, OR, USA). An extra test was performed where Ψleaf of adjacent leaflets was compared, and no differences in Ψleaf were detected between two adjacent leaflets. Some plants were kept in dark conditions for ~8–10 h for the determination of the soil water potential (Ψsoil) (Tsuda and Tyree, 2000). Since water movement was ~zero under dark conditions, Ψsoil remained relatively constant and can be assumed to equal the xylem pressure potential of leaves under dark conditions.
Leaf Morphology
After determination of the plant water status, the leaflet area was measured by a leaf area metre. The leaflet samples were dried at 80°C in an oven to a constant dry mass and weighed. The leaf mass area (LMA) was determined as the ratio of leaf dry mass to leaf area.
Leaf ABA Concentration
After reaching the steady-state of photosynthesis, leaflets were harvested for ABA detection and transcriptome sequencing. Phytohormone contents were determined by a liquid chromatography electrospray ionisation tandem mass spectrometry (LC-ESI-MS/MS) system (HPLC, Shim-pack UFLC SHIMADZU CBM30A system, www.shimadzu.com.cn/; MS, Applied Biosystems 6500 Triple Quadrupole, www.appliedbiosystems.com.cn/). Briefly, the leaflets for photosynthetic measurements were harvested and frozen in liquid nitrogen. The samples were extracted with methanol/water/formic acid and filtered before LC-MS/MS analysis. The detailed protocol was described on MetWare (http://www.metware.cn/) based on the AB Sciex QTRAP 6500 LC-MS/MS platform. Samples were detected with three biological replicates.
RNA Extraction and Transcriptome Sequencing
Total RNA was extracted from leaflet samples for transcriptome sequencing, according to a previous study (Zhang et al., 2019b). Sequencing libraries were constructed using the UltraTM RNA Library Prep Kit for Illumina (NEB, USA) according to the instructions of the manufacturer. The detailed protocol was described in a previous study, which was briefly described in a simplified diagram.
Sequencing Data Analysis
The clean data were obtained by processing the raw data through in-house Perl scripts. The low-quality data and sequencing adapters were trimmed. The fragments per kilobase of transcript per million fragments mapped (FPKM) was calculated based on gene length and read counts. Differentially expressed genes (DEGs) were assigned according to the adjusted P < 0.05. Gene ontology (GO) enrichment was determined by submitting DEGs to the GO database to classify the genes. Kyoto Encyclopaedia of Genes and Genomes (KEGG; https://www.genome.jp/kegg) was used to perform pathway enrichment analysis. Terms with corrected P < 0.05 were identified as significantly enriched by DEGs. To confirm the reliability of transcriptome sequencing, 10 candidates expressed genes in RNA-seq were simultaneously evaluated by qRT–PC analysis. The qRT–PCR values were linearly correlated with the RNA-seq FPKM values (P < 0.001).
Statistical Analyses
All statistical analyses were performed using SPSS 19. One-way ANOVA was used to determine the significant difference of average values according to Tukey's test (P < 0.05). Regression analysis was performed by Microsoft Excel.
Results
Effect of VPD on Water Transport Forces Along the Soil-Plant-Atmospheric Continuum
Vapour pressure deficit significantly affected the distribution of water potential along the soil-plant-atmospheric pathway (Figure 1). Atmospheric evaporative demand increased with VPD elevation, which triggered plant water stress and a linear decline in the leaf water potential (Figure 1A). With VPD elevation, the drawdown of Ψleaf in high-VPD-grown plants was less than that in low-VPD-treated plants according to the slope of linear regression (Figure 1A). The driving force for passive water flow between the soil and leaf (ΔΨ soil−leaf) increased with the VPD, and the magnitude of the increase was greater in low-VPD-grown plants than high-VPD-grown plants (Figure 1B). Since Ψleaf was negligible compared with the large negative air potential, the water driving force at the leaf-air boundary (ΔΨleaf−air) increased dramatically with an increase in the VPD (Figure 1C). The magnitude of the increase in ΔΨleaf−air was considerably greater than that of ΔΨsoil−leaf, and the difference between ΔΨleaf−air and ΔΨ soil−leaf was enlarged with the VPD: the ratio of ΔΨleaf−air to ΔΨsoil−leaf increased logarithmically from ~50 at 0.5 kPa to 150 at 1.5 kPa and then maintained at a steady level (Figure 1D). The statistical analyses of the plant water status are shown in Supplementary Table 1.
Figure 1. Effect of the vapour pressure deficit (VPD) on the spatial distribution of the water potential and driving force (ΔΨ) between two spatial positions. Values are the mean ± SE (n = 4~6 replicates). The regression lines shown are: (A) HVPD, Ψleaf = −0.242 VPD −0.358, R2 = 0.92; LVPD, Ψleaf = −0.258 VPD −0.347, R2 = 0.93. (B) HVPD, ΔΨsoil−leaf = 0.242 VPD + 0.118, R2 = 0.94; LVPD, ΔΨsoil−leaf = 0.251 VPD + 0.13, R2 = 0.93. (C) HVPD, ΔΨleaf−air = 49.8 VPD −6.86, R2 = 0.98; LVPD, ΔΨleaf−air = 49.8 VPD−6.86, R2 = 0.98. (D) HVPD, ΔΨleaf−air/ΔΨsoil−leaf = 50.82 ln (VPD) + 109.95, R2 = 0.85; LVPD, ΔΨleaf−air/ΔΨsoil−leaf = 48.63 ln (VPD) + 110.41, R2 = 0.8.
Effect of VPD on the Photosynthetic Parameters of Tomato Plants
The photosynthetic rate responded to CO2 elevation in similar patterns regardless of cultivar and VPD growth conditions: the photosynthetic rate rose rapidly across low CO2 concentrations and then reached a steady state (Supplementary Figure 1). The maximum steady-state photosynthetic rate declined as the VPD increased from 0.5 to 4.5 kPa (Supplementary Figure 1). The maximum carboxylation rate (Vcmax), maximum electron transport rate (Jmax), and CE declined linearly with VPD elevation (Supplementary Figure 2). The drawdown of Vcmax, Jmax, and CE with VPD elevation was moderated in high-VDP-grown plants compared with low-VPD-grown plants according to the slope of linear regression (Supplementary Figure 2). The statistical analyses of photosynthetic parameters across VPD ranges are shown in Supplementary Table 2.
Effect of VPD on the Photosynthetic CO2 Uptake and Transport
The stomatal, mesophyll, and total conductance for CO2 diffusion decreased linearly with VPD elevation, regardless of the cultivar and VPD growth conditions (Figure 2). The magnitudes of drawdown in the stomatal, mesophyll, and total conductance were lower in high-VPD-grown plants than in low-VPD-grown plants for two cultivars (Figure 2).
Figure 2. Effect of the VPD on the stomatal conductance (gs), mesophyll conductance (gm), and total conductance (gtot) for photosynthetic carbon dioxide (CO2) diffusion. Values are the mean ± SE (n = 4 replicates). The regression lines shown are: (A) HVPD, gs = −0.187 VPD + 0.895, R2 = 0.94; LVPD, gs = −0.200 VPD + 0.928, R2 = 0.92. (B) HVPD, gm = −0.126 VPD + 0.661, R2 = 0.95; LVPD, gm = −0.129 VPD + v0.658, R2 = 0.94. (C) HVPD, gtot = −0.0719 VPD +0.365, R2 = 0.95; LVPD, gtot = −0.0798 VPD + 0.384, R2 = 0.94.
The CO2 concentration along the “source-path-sink” was reduced to different extents with VPD elevation (Figures 3A,B). The drawdowns of Ci and Cc caused by VPD elevation were relatively lower in high-VPD-grown plants than in low-VPD-grown plants (Figures 3A,B). Consequently, the CO2 transport efficiency of Ci/Ca, Cc/Ca and Cc/Ci decreased linearly with VPD elevation. The declining slopes of Ci/Ca, Cc/Ca, and Cc/Ci vs. VPD were lower in high-VPD-grown plants than in low-VPD-grown plants (Figures 3C–E). The statistical analyses of CO2 concentrations along the “source-path-sink” across VPD ranges are shown in Supplementary Table 3.
Figure 3. Effect of the VPD on the intracellular CO2 concentration [(A); Ci], carboxylation sites inside chloroplasts CO2 concentration [(B); CC], the ratio of the intercellular to ambient CO2 concentration [(C); Ci/Ca], the ratio of the chloroplast to ambient CO2 concentration [(D); Cc/Ca] and the ratio of the chloroplast to intercellular CO2 concentration [(E); Cc/Ci]. The regression lines shown are: (A) HVPD, Ci = −17.2 VPD + 371.6, R2 = 0.86; LVPD, Ci = −22 VPD + 377.6, R2 = 0.87. (B) HVPD, CC = −31.5 VPD + 331.6, R2 = 0.91; LVPD, CC = −33.5 VPD + 332.2, R2 = 0.88. (C) HVPD, Ci/Ca = −0.0429 VPD + 0.93, R2 =0.86; LVPD, Ci/Ca= −0.055 VPD + 0.94, R2 = 0.87. (D) HVPD, Cc/Ca = −0.0788VPD + 0.83, R2 = 0.92; LVPD, Cc/Ca = −0.0837VPD + 0.83, R2 = 0.87. (E) HVPD, Cc/Ci = −0.0537 VPD + 0.89, R2 = 0.82; LVPD, Cc/Ci = −0.058VPD + 0.91, R2 = 0.89.
Partial Photosynthetic Limitation
The fractions of stomatal, mesophyll, and biochemical limitations imposed on photosynthesis varied with VPD elevation (Figure 4). Under low VPD conditions, the stomatal and mesophyll conductance for CO2 diffusion were high and imposed relatively minor limitations on photosynthesis. The stomatal and mesophyll conductance accounted for a low proportion of photosynthetic limitation, while biochemical carboxylation for carbon fixation was the most significant limitation for photosynthetic processes under low VPD conditions (Figure 4). The fraction of stomatal limitation increased linearly with the VPD, from ~15% at 0.5 kPa to 35% at 4.5 kPa (Figure 4A). A similar pattern was observed in the mesophyll limitation: the fraction of mesophyll limitation also increased linearly with VPD elevation, from ~23% at 0.5 kPa to 33% at 4.5 kPa (Figure 4B). The increments in the fractions of stomatal and mesophyll limitations tended to be less marked in high-VPD-grown plants. In contrast, the fraction of total limitations attributed to the biochemical limitation of carbon fixation gradually decreased linearly with VPD elevation, from ~65% at.5 kPa to 35% at 4.5 kPa (Figure 4C). The statistical analyses of stomatal, mesophyll, and biochemical limitation fractions across VPD ranges are shown in Supplementary Table 4.
Figure 4. Quantitative limitation analysis comparing stomatal [(A); Ls], mesophyll [(B); Lm], and biochemical [(C); Lb] limitations imposed on the photosynthetic rate under varying VPD. The regression lines shown are: (A) HVPD, Ls = 0.0472 VPD + 0.122, R2 = 0.86; LVPD, Ls = 0.0535 VPD + 0.115, R2 = 0.89. (B) HVPD, Lm = 0.0266 VPD + 0.208, R2 = 0.88; LVPD, Lm = 0.0305 VPD + 0.197, R2 = 0.91.
Biochemical limitation accounted for the greatest limitation on photosynthesis under low VPD conditions, regardless of the cultivar and growth conditions (Figure 5). The limitations that stomatal and mesophyll conductance imposed on photosynthesis gradually increased and predominated under high VPD stress (Figure 5). Diffusion limitations, i.e., the sum of the stomatal and mesophyll resistance, were the rate-limiting step for the photosynthetic process under high VPD conditions, which imposed the greatest limitation on photosynthesis in tomato plants (Figure 5).
Figure 5. The dynamic changes in the relative proportions of individual components of photosynthetic limitations across VPD ranges: (A) CV1 grown under high VPD condition; (B) CV1 grown under low VPD condition; (C) CV2 grown under high VPD condition; and (D) CV2 grown under low VPD condition.
Correlations Among gm, gs, Leaf Water Status, and LMA
The mesophyll conductance was significantly and positively correlated with the stomatal conductance (Supplementary Figure 3A). Meanwhile, the stomatal and mesophyll conductance for CO2 diffusion were closely linked to the leaf water status, wherein significant and positive correlations were found in the leaf water potential vs. the stomatal and mesophyll conductance (Supplementary Figures 3B,C). Acclimation to VPD modified leaf structural traits, wherein LMA tended to be slightly greater in high-VPD-grown plants than in low-VPD-grown plants (Supplementary Figure 4A). A significant and negative correlation between gm and LMA was observed (Supplementary Figure 4B).
Leaf ABA Concentration and Correlation With CO2 Diffusion Conductance
With VPD elevation, the foliar ABA content increased exponentially (Figure 6A). The leaf ABA content was linearly and negatively correlated with the CO2 diffusion conductance of gs and gm (Figure 6B). The slope of linear regression in gs was more negative than that in gm, indicating that gs was more sensitive to ABA in response to VPD stress (Figure 6B).
Figure 6. Leaf ABA content in response to the VPD (A) and its correlation with the CO2 diffusion conductance of gs and gm (B). The regression lines shown are ABA = 301.3e0.15VPD, R2 = 0.96, P < 0.01; gs = −0.0029 ABA + 1.76, R2 = 0.85, P < 0.01; gm = −0.0019 ABA + 1.19, R2 = 0.88, P < 0.01.
Transcriptomic Analysis of Plant Response Across Series of VPD Ranges
Kyoto Encyclopaedia of Genes and Genomes analysis showed that physiological processes of “metabolic pathway” and “plant hormone signal transduction” were involved in the response to VPD and were potentially associated with ABA biosynthesis and signal transduction (Figure 7). “Metabolic pathway” was the dominant pathway in response to VPD elevation. “Plant hormone signal transduction” also exhibited a significant pathway in response to VPD stress in the range from 1.5 to 3.5 kPa (Figure 7).
Figure 7. KEGG classification analysis of differentially expressed genes under different VPD conditions: (A) 0.5 kPa versus 1.5 kPa; (B) 0.5 kPa versus 2.5 kPa; and (C) 0.5 kPa versus 3.5 kPa.
The enriched genes in the comparison between different VPD treatments were annotated in three main GO categories: biological process, cellular component, and molecular function. A great difference in the top 50 GO enrichment terms was observed between different VPD conditions (Figure 8). In the comparison of 0 with 1.5 kPa and 0.5 kPa with 3.5 kPa, the “ABA-activated signalling pathway” was significantly enriched and was involved in the VPD response (Figures 9A,C). In the comparison of 0.5 kPa with. 2.5 kPa, ABA was potentially involved according to the enriched terms of “cellular hormone metabolic process,” “hormone biosynthetic process,” and “hormone metabolic process”. Gene expression with rising VPD can be classified into 10 patterns according to K-means analysis (Figure 9). Gene expression patterns were mostly classified into the pattern of “Subclass 6” with 768 genes, wherein gene expression remained relatively stable under mild VPD stress and increased dramatically under high VPD water stress (Figure 9). The genes associated with ABA biosynthesis and signal transduction followed different patterns.
Figure 8. Top 50 enriched GO terms of the differentially expressed genes under different VPD conditions: (A) 0.5 kPa versus 1.5 kPa; (B) 0.5 kPa versus 2.5 kPa; and (C) 0.5 kPa versus 3.5 kPa.
Figure 9. Expression model of DEGs across a series of VPD ranges. A, B, C, and D in the X-axis in the figures represent 0.5, 1.5, 2.5, and 3.5 kPa, respectively.
Discussion
The present study assessed the rate-limiting step for tomato plant photosynthesis across a series of VPD ranges and evaluated ABA-mediated regulatory mechanisms according to physiological and transcriptomic analyses. The key rate-limiting step for photosynthetic performance varied with the VPD: under low VPD conditions, stomatal, and mesophyll conductance was high for efficient CO2 transport, which facilitated sufficiently high CO2 availability inside chloroplasts for carbon fixation (Figures 2, 3). With VPD elevation, the stomatal and mesophyll conductance for CO2 transport declined gradually. Consequently, photosynthesis was substantially constrained by the low chloroplast CO2 concentration under high VPD conditions (Figure 3). Therefore, the CO2 diffusion limitation in a series of stomatal and mesophyll resistances was the key rate-limiting step for photosynthesis under high VPD conditions (Figure 4). In addition to anatomical determination, ABA accumulation and signal transduction were involved in maintaining the water balance in response to VPD. ABA accumulation was negatively correlated with CO2 diffusion conductance (Figure 6). Three steps were involved in the potential mechanism accounting for the increased limitation of stomatal and mesophyll conductance imposed on tomato plant photosynthesis with VPD elevation (Supplementary Figure 5): (I) VPD elevation caused plant water stress by disrupting the mass balance between the soil water supply and atmospheric evaporative demand; (II) plants maintained the water balance by regulating ABA accumulation and signal transduction in response to high VPD stress; (III) ABA in combination with leaf anatomical adaptation modulated CO2 uptake and transport.
VPD Elevation Triggers Plant Water Stress by Disrupting the Mass Balance Between Soil Water Supply and Atmospheric Evaporative Demand
Passive water movement was driven by the gradient of free energy along the soil-plant-atmospheric continuum, which could be quantified as the gradient in water potential in the liquid phase. Water movement at the leaf-air boundary in the gas phase was driven by the difference in the VPD. Based on physical principles, excessive air desiccation triggered a high VPD and great negative air-water potential. Δψleaf−air was substantially >Δψsoil−leaf, which drove transpiration. The substantial difference between Δψleaf−air and Δψsoil−leaf was logarithmically enlarged with an increase in the VPD (Figure 1). Quantitatively, the atmospheric driving force at the leaf-air boundary could be >100-fold larger than the soil-leaf component under high VPD conditions (Figure 1). The great asymmetry between the atmospheric evaporative demand and soil water supply triggered disruption in the water balance despite plants being well irrigated. Root water uptake and supply were inadequate to keep pace with the great atmospheric driving force under high VPD conditions, which consequently triggered leaf dehydration and decline in water potential. Therefore, the VPD is a crucial external stimulus moving water through a soil-plant-environment continuum. VPD fluctuates dramatically over the diurnal course in crop production, especially for greenhouse cultivation. Soil moisture is relatively stable over the short term compared with the VPD (Caldeira et al., 2014). Plant-water relations are regulated to a greater extent by the VPD and to a lesser extent by soil moisture. Similar to soil drought, VPD-induced plant water stress is also an important factor triggering photosynthetic depression.
Plants Maintain the Water Balance by Regulating ABA Accumulation and Signal Transduction in Response to High VPD Stress
The stoma is the “gatekeeper” for the exchange of water vapour and CO2. Guard cells surrounding the stomatal pore respond to perturbations of the soil-plant-atmospheric hydraulic continuum, which is putatively transduced into stomatal movements by feedback and feedforward mechanisms (Buckley, 2005, 2017, 2019). Stomatal control of transpired water loss is critical for sustaining physiological processes, such as leaf water status and photosynthetic CO2 uptake. It has been recognised that plants respond to drought by closing guard cells to prevent the development of water dehydration in plant tissues (Novick et al., 2016). In the present study, the atmospheric driving force was an order of magnitude greater than the water supply, which led to a great dissymmetry between the water supply and evaporative demand. The dissymmetry between the water supply and evaporative demand triggered declines in the leaf water potential and stomatal closure. However, the mechanism of VPD-triggered stomatal closure is still uncertain and is a “black box” (Buckley, 2016). Some hypotheses hold that stomatal closure in angiosperms under high VPD conditions is an active process that is regulated by hormonal and hydraulic signals (Merilo et al., 2018; Pantin and Blatt, 2018). The plant stress hormone ABA is continuously produced and delivered with a transpiration stream to guard cells (Qiu et al., 2017; Merilo et al., 2018). In the present study, leaf ABA rapidly accumulated with the rise in the VPD. Transcriptome analysis suggested that ABA biosynthesis and signal transduction were potentially involved in the response to the VPD. Based on the theory of ion channel-mediated guard cell signal transduction (Julian et al., 2001), hypothetical mechanisms of the ABA-mediated stomatal closure response to high VPD-induced water stress are proposed in Supplementary Figure 6. Ions and water flowed into guard cells under low VPD conditions and sustained turgor for stomatal openness. Under high VPD-induced water stress, ABA rapidly accumulated and promoted stomatal closure by altering ion channel activities.
Although stomatal closure prevented excess water loss to maintain physiological processes by passive or active mechanisms, closed “gatekeepers” simultaneously increased stomatal resistance for photosynthetic CO2 uptake from air to intercellular. The intercellular CO2 concentration was gradually reduced with VPD elevation (Figure 3). Consequently, the stomatal limitation imposed on photosynthesis gradually became pronounced with VPD elevation (Figure 4). The declines in the leaf water potential and stomatal conductance with VPD elevation were less marked in high-VPD-grown plants in the present research. The distinct response to the VPD is potentially modulated by physiological acclimation to growth conditions (Fanourakis et al., 2019). In contrary to the previous study, no significant differences were observed in the photosynthetic parameters across VPD ranges between the two examined cultivars. The distinct responses of two cultivars between the previous and present study were potentially caused by the examined VPD conditions. The dynamic VPD response of the present and previous studies were performed under cabinets and greenhouse conditions, respectively.
Anatomical Properties and ABA Modulate Mesophyll Conductance Under Contrasting VPD Conditions
In addition to the first barrier of stomata, CO2 movement from intercellular to carboxylation sites is constrained by mesophyll resistance. The present study demonstrated that mesophyll resistance was a significant component of diffusion resistance from air to Rubisco in tomato plants. A strong positive correlation between the mesophyll and stomatal conductance was observed among treatments (Figure 6). Similar to the stomatal conductance, the mesophyll conductance of tomato plants was also linearly reduced with VPD elevation (Figure 2). Under low VPD conditions, stomatal conductance coupled with mesophyll conductance was high for efficient CO2 transport to carboxylation sites within chloroplasts. High diffusion conductance facilitated high chloroplast CO2 concentrations for carbon fixation (Figure 3). With VPD elevation, the CO2 concentration inside chloroplasts was substantially reduced under high VPD conditions. The limitation of mesophyll conductance imposed on photosynthesis gradually dominated with VPD elevation (Figure 5). Leaf anatomical traits from the substomatal cavity to the carbon fixation site determine the maximum potential of mesophyll conductance (Muir et al., 2014; Xiong et al., 2017; Earles et al., 2018; Han et al., 2018; Carriqui et al., 2019). The LMA is a composite of underlying traits such as the lamina thickness, mesophyll thickness, cell wall thickness, cell shape, and bulk leaf density, which anatomically regulate the mesophyll conductance (Muir et al., 2014). The LMA determines the upper limit on mesophyll conductance. Meanwhile, the LMA is closely linked to abiotic stress tolerance (Xiong and Flexas, 2018; Xiong et al., 2018). Generally, a higher LMA is a good indicator of greater stress tolerance. In the present study, the LMA of high-VPD-grown plants was lower but higher than that of low-VPD-grown plants (Supplementary Figure 4). Long-term acclimation to a high VPD facilitates enhanced drought tolerance to prevent dehydration by regulating the leaf thickness, cuticular permeability, stomatal morphology, and other anatomical features (Fanourakis et al., 2016, 2020). Long-term exposure to a VPD also affects stomatal sensitivity and morphological features such as the stomatal size, density, index, and spacing, which consequently modulate transpired water loss (Fanourakis et al., 2016, 2020). As mentioned above, root water uptake and supply are inadequate to keep pacing with the great atmospheric driving force under high VPD conditions. A higher LMA indicated dense structural traits, which buffered cellular transpired water loss and prevented leaf tissue dehydration under high VPD conditions. However, CO2 and water transport share pathways through the mesophyll cell walls and perhaps plasma membranes within leaves (Barbour, 2017; Groszmann et al., 2017; Zhao et al., 2017; Drake et al., 2019). Although dense structural traits improved drought tolerance, the resistance of CO2 diffusion was simultaneously increased. The LMA was negatively correlated with mesophyll conductance in the present study, which is consistent with previous studies (Hassiotou et al., 2009).
In addition to anatomical determinations, biochemical regulations such as ABA, carbonic anhydrase, and aquaporin facilitate rapid mesophyll conductance responses to short-term changing environmental factors (Momayyezi et al., 2020). The mesophyll conductance is negatively correlated with the leaf ABA content in tomato plants, which is in accordance with a previous study (Sorrentino et al., 2016; Qiu et al., 2017). The foliar ABA content rapidly increased upon long-term and short-term exposure to a high VPD, which is in accordance with a previous study (McAdam and Brodribb, 2016). However, the ABA-mediated regulatory mechanism has rarely been reported. CO2 entering from intercellular to carboxylation sites inside chloroplasts must pass through plasma membranes. The resistance of transport across the membrane accounts for a great proportion of the mesophyll resistance. It is now well established that aquaporins function as water pores for water transport across membranes and play significant roles in maintaining water homeostasis in response to drought and salinity (Qian et al., 2015; Zhang et al., 2019a). There is increasing evidence that some specific aquaporins (which localise to the plasma membrane and chloroplast inner envelope membrane) are permeable to CO2 and contribute to the mesophyll conductance (Uehlein et al., 2012; Groszmann et al., 2017; Zhao et al., 2017). Similar to guard cells, some PIPs pores mediate CO2 uptake and water transport across the plasma membrane (Zhang et al., 2021). Therefore, the specific PIPs potentially reconciled the trade-off between carbon gain and water loss in response to VPD-induced water stress. PIPs were sensitive to drought signals and responded rapidly to enclose gating and inhibit activity.
Gating is a general mechanism of membrane-mediated channels for controlling the permeability of water and CO2. Although the inhibition of PIPs channels prevents water loss under high VPD stress, CO2 uptake across the membrane is also restricted by the gating enclosure. ABA has been reported as a signal-inducing variation in the aquaporin content and activity (Fang et al., 2019). Therefore, VPD potentially modulated PIPs gating for CO2 and water permeability via ABA signalling, which contributed to mesophyll conductance and the photosynthetic rate.
Conclusions
The present study revealed the rate-limiting step for photosynthetic CO2 utilisation under contrasting VPD conditions and proposed ABA-mediated regulatory mechanisms according to transcriptomic and physiological evidence. The photosynthetic performance of tomato plants was gradually constrained with VPD elevation. The key rate-limiting steps for photosynthetic performance varied with the rise in the VPD. With VPD elevation, plant water stress was gradually pronounced and triggered linear declines in the stomatal and mesophyll conductance. The contributions of stomatal and mesophyll limitations to photosynthesis increased gradually with VPD elevation. Consequently, the low CO2 availability inside chloroplasts substantially constrained photosynthesis under high VPD conditions. Leaf ABA accumulated rapidly with pronounced water stress under a high VPD and negatively correlated with the stomatal and mesophyll conductance for CO2 diffusion. Transcriptomic combined with physiological analyses revealed that ABA biosynthesis and signal transduction were potentially involved in mediating CO2 transport in response to the VPD.
Data Availability Statement
The original contributions presented in the study are publicly available. This data can be found here: National Center for Biotechnology Information (NCBI) BioProject database under accession number PRJNA762604.
Author Contributions
DZ, QL, and MW conceived and designed the experiments. PS, JL, and XL conducted the experiments. QD analysed the data and wrote the draft. All the authors reviewed and approved the manuscript.
Funding
This project was supported by the National Natural Science Foundation of China (32102466), the Natural Science Foundation of Shandong Province (ZR2019BC035), and the Major Scientific Innovation Project of Shandong Province (2019JZZY010715).
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.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.745110/full#supplementary-material
Abbreviations
VPD, vapour pressure deficit; Ψleaf, leaf water potential; Ψsoil, soil water potential; Ψair, air-water potential; ΔΨ leaf−air, the drawdown of water potential between leaf and air; ΔΨ soil−leaf, the drawdown of water potential between soil and leaf; Vcmax, Maximum carboxylation rate; Jmax, maximum electron transport rate; CE, carboxylation efficiency; gs, stomatal conductance; gm, mesophyll conductance; gtot, total conductance; Ca, ambient CO2 concentration; Ci, intracellular CO2 concentration; CC, CO2 concentration of carboxylation sites inside chloroplast; Ls, stomatal limitations imposed on the photosynthetic rate; Lm, mesophyll limitations imposed on the photosynthetic rate; Lb, biochemical limitations imposed on the photosynthetic rate; LMA, leaf mass area; Pn, net photosynthetic rate; Rd, the rate of mitochondrial respiration in the light; Γ, chloroplastic CO2 compensation point.
References
Barbour, M. M. (2017). Understanding regulation of leaf internal carbon and water transport using online stable isotope techniques. New phytol. 213, 83–88. doi: 10.1111/nph.14171
Buckley, T. N. (2005). The control of stomata by water balance. New phytol. 168, 275–292. doi: 10.1111/j.1469-8137.2005.01543.x
Buckley, T. N. (2016). Stomatal responses to humidity: has the 'black box' finally been opened? Plant Cell Environ. 39, 482–484. doi: 10.1111/pce.12651
Buckley, T. N. (2017). Modeling stomatal conductance. Plant physiol. 174, 572–582. doi: 10.1104/pp.16.01772
Buckley, T. N. (2019). How do stomata respond to water status? New phytol. 224, 21–36 doi: 10.1111/nph.15899
Caldeira, C. F., Bosio, M., Parent, B., Jeanguenin, L., Chaumont, F., and Tardieu, F. (2014). A hydraulic model is compatible with rapid changes in leaf elongation under fluctuating evaporative demand and soil water status. Plant Physiol. 164, 1718–1730. doi: 10.1104/pp.113.228379
Carriqui, M., Roig-Oliver, M., Brodribb, T. J., Coopman, R., Gill, W., Mark, K., et al. (2019). Anatomical constraints to nonstomatal diffusion conductance and photosynthesis in lycophytes and bryophytes. New phytol. 222, 1256–1270 doi: 10.1111/nph.15675
Drake, P. L., Boer, H. J. D., Schymanski, S. J., and Veneklaas, E. J. (2019). Two sides to every leaf: water and CO2 transport in hypostomatous and amphistomatous leaves. New Phytolol. 222, 1179–1187. doi: 10.1111/nph.15652
Du, Q., Song, X., Bai, P., Jiao, X., Ding, J., Zhang, J., et al. (2020). Effect of different vapor pressure deficits on gas exchange parameters and growth of tomatoes and comprehensive evaluation. Acta Agriculturae Boreli-occidentalis Sinica. 29, 66–74.
Du, Q., Liu, T., Jiao, X., Song, X., Zhang, J., and Li, J. (2019). Leaf anatomical adaptations have central roles in photosynthetic acclimation to humidity. J. Exp. Bot. 70, 4949–4962. doi: 10.1093/jxb/erz238
Earles, J. M., Theroux-Rancourt, G., Roddy, A. B., Gilbert, M. E., McElrone, A. J., and Brodersen, C. R. (2018). Beyond porosity: 3D leaf intercellular airspace traits that impact mesophyll conductance. Plant Physiol. 178, 148–162. doi: 10.1104/pp.18.00550
Fang, L., Abdelhakim, L. O. A., Hegelund, J. N., Li, S., Liu, J., Peng, X., et al. (2019). ABA-mediated regulation of leaf and root hydraulic conductance in tomato grown at elevated CO2 is associated with altered gene expression of aquaporins. Hortic. Res. 6, 1–10. doi: 10.1038/s41438-019-0187-6
Fanourakis, D., Aliniaeifard, S., Sellin, A., Giday, H., Korner, O., Rezaei Nejad, A., et al. (2020). Stomatal behavior following mid- or long-term exposure to high relative air humidity: A review. Plant Physiol. Biochem. 153, 92–105. doi: 10.1016/j.plaphy.2020.05.024
Fanourakis, D., Bouranis, D., Giday, H., Carvalho, D. R., Rezaei Nejad, A., and Ottosen, C. O. (2016). Improving stomatal functioning at elevated growth air humidity: A review. J. Plant Physiol. 207, 51–60. doi: 10.1016/j.jplph.2016.10.003
Fanourakis, D., Giday, H., Hyldgaard, B., Bouranis, D., Körner, O., and Ottosen, C. O. (2019). Low air humidity during cultivation promotes stomatal closure ability in rose. Eur. J. Hortic. Sci. 84, 245–252. doi: 10.17660/eJHS.2019/84.4.7
Farquhar, G. D., von Caemmerer, S., and Berry, J. A. (1980). A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species. Planta. 149, 78–90. https://doi.org/10.1007/BF00386231.
Flexas, J., Barbour, M. M., Brendel, O., Cabrera, H. M., Carriqui, M., Diaz-Espejo, A., et al. (2012). Mesophyll diffusion conductance to CO2: an unappreciated central player in photosynthesis. Plant Sci. 193–194, 70–84. doi: 10.1016/j.plantsci.2012.05.009
Giuliani, R., Koteyeva, N., Voznesenskaya, E., Evans, M. A., Cousins, A. B., and Edwards, G. E. (2013). Coordination of leaf photosynthesis, transpiration, and structural traits in rice and wild relatives (Genus Oryza). Plant Physiol. 162, 1632–1651. doi: 10.1104/pp.113.217497
Groszmann, M., Osborn, H. L., and Evans, J. R. (2017). Carbon dioxide and water transport through plant aquaporins. Plant Cell Environ. 40, 938–961. doi: 10.1111/pce.12844
Han, J., Lei, Z., Flexas, J., Zhang, Y., Carriqui, M., Zhang, W., et al. (2018). Mesophyll conductance in cotton bracts: anatomically-determined internal CO2 diffusion constraints on photosynthesis. J. Exp. Bot. 69, 5433–5443. doi: 10.1093/jxb/ery296
Harley, P. C., Loreto, F., Di Marco, G., and Sharkey, T. D. (1992). Theoretical considerations when estimating the mesophyll conductance to CO2 flux by analysis of the response of photosynthesis to CO2. Plant Physiol. 98, 1429–1436. doi: 10.1104/pp.98.4.1429
Hassiotou, F., Ludwig, M., Renton, M., Veneklaas, E. J., and Evans, J. R. (2009). Influence of leaf dry mass per area, CO2, and irradiance on mesophyll conductance in sclerophylls. J. Exp. Bot. 60, 2303–2314. doi: 10.1093/jxb/erp021
Julian, I., Schroeder, Gethyn, J., Allen, Veronique Hugouvieux, June, M., Kwak, a., and Waner, D. (2001). Guard cell signal transduction. Annu. Rev. Plant Physiol. Plant Mol.Biol. 52, 627–658. doi: 10.1146/annurev.arplant.52.1.627
Kaldenhoff, R. (2012). Mechanisms underlying CO2 diffusion in leaves. Curr Opin Plant Biol.15, 276–281. doi: 10.1016/j.pbi.2012.01.011
Lawson, T., and Blatt, M. R. (2014). Stomatal size, speed, and responsiveness impact on photosynthesis and water use efficiency. Plant Physiol. 164, 1556–1570. doi: 10.1104/pp.114.237107
Li, Q., Liu, Y., Tian, S., Liang, Z., Li, S., Li, Y., et al. (2019a). Effect of supplemental lighting on water transport, photosynthetic carbon gain and water use efficiency in greenhouse tomato. Sci. Hortic. 256, 108630. doi: 10.1016/j.scienta.2019.108630
Li, Q., Wei, M., Li, Y., Feng, G., Wang, Y., Li, S., et al. (2019b). Effects of soil moisture on water transport, photosynthetic carbon gain and water use efficiency in tomato are influenced by evaporative demand. Agric. Water Manage. 226, 105818. doi: 10.1016/j.agwat.2019.105818
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
McAdam, S. A., and Brodribb, T. J. (2016). Linking turgor with ABA biosynthesis: implications for stomatal responses to vapor pressure deficit across land plants. Plant Physiol. 171, 2008–2016. doi: 10.1104/pp.16.00380
Merilo, E., Yarmolinsky, D., Jalakas, P., Parik, H., Tulva, I., Rasulov, B., et al. (2018). Stomatal VPD response: there is more to the story than ABA. Plant Physiol. 176, 851–864. doi: 10.1104/pp.17.00912
Momayyezi, M., McKown, A. D., Bell, S. C. S., and Guy, R. D. (2020). Emerging roles for carbonic anhydrase in mesophyll conductance and photosynthesis. Plant J. 101, 831–844. doi: 10.1111/tpj.14638
Muir, C. D., Hangarter, R. P., Moyle, L. C., and Davis, P. A. (2014). Morphological and anatomical determinants of mesophyll conductance in wild relatives of tomato (Solanum sect. Lycopersicon, sect. Lycopersicoides; Solanaceae). Plant Cell Environ. 37, 1415–1426. doi: 10.1111/pce.12245
Niinemets, U., Diaz-Espejo, A., Flexas, J., Galmes, J., and Warren, C. R. (2009). Role of mesophyll diffusion conductance in constraining potential photosynthetic productivity in the field. J. Exp. Bot. 60, 2249–2270. doi: 10.1093/jxb/erp036
Norby, R. J. (2002). Plant water relations at elevated CO2-implications for water-limited environments. Plant Cell Environ.25, 319–331. doi: 10.1046/j.1365-3040.2002.00796.x
Novick, K. A., Miniat, C. F., and Vose, J. M. (2016). Drought limitations to leaf-level gas exchange: results from a model linking stomatal optimization and cohesion-tension theory. Plant Cell Environ. 39, 583–596. doi: 10.1111/pce.12657
Pantin, F., and Blatt, M. R. (2018). Stomatal response to humidity: blurring the boundary between active and passive movement. Plant Physiol. 176, 485–488, doi: 10.1104/pp.17.01699
Qian, Z. J., Song, J. J., Chaumont, F., and Ye, Q. (2015). Differential responses of plasma membrane aquaporins in mediating water transport of cucumber seedlings under osmotic and salt stresses. Plant Cell Environ. 38, 461–473, doi: 10.1111/pce.12319
Qiu, C., Ethier, G., Pepin, S., Dube, P., Desjardins, Y., and Gosselin, A. (2017). Persistent negative temperature response of mesophyll conductance in red raspberry (Rubus idaeus L.) leaves under both high and low vapour pressure deficits: a role for abscisic acid? Plant Cell Environ. 40, 1940–1959, doi: 10.1111/pce.12997
Sack, L., John, G. P., and Buckley, T. N. (2018). ABA accumulation in dehydrating leaves is associated with decline in cell volume, not turgor pressure. Plant Physiol 176, 489–495. doi: 10.1104/pp.17.01097
Sharkey, T. D. (2012). Mesophyll conductance: constraint on carbon acquisition by C3 plants. Plant Cell Environ.35, 1881–1883. doi: 10.1111/pce.12012
Sorrentino, G., Haworth, M., Wahbi, S., Mahmood, T., Zuomin, S., and Centritto, M. (2016). Abscisic acid induces rapid reductions in mesophyll conductance to carbon dioxide. PloS ONE. 11, e0148554. doi: 10.1371/journal.pone.0148554
Sun, J., Ye, M., Peng, S., and Li, Y. (2016). Nitrogen can improve the rapid response of photosynthesis to changing irradiance in rice (Oryza sativa L.) plants. Sci Rep. 6, 31305. doi: 10.1038/srep31305
Tholen, D., and Zhu, X. G. (2011). The mechanistic basis of internal conductance: a theoretical analysis of mesophyll cell photosynthesis and CO2 diffusion. Plant Physiol.156, 90–105. doi: 10.1104/pp.111.172346
Tomas, M., Flexas, J., Copolovici, L., Galmes, J., Hallik, L., Medrano, H., et al. (2013). Importance of leaf anatomy in determining mesophyll diffusion conductance to CO2 across species: quantitative limitations and scaling up by models. J. Exp. Bot. 64, 2269–2281. doi: 10.1093/jxb/ert086
Tsuda, M., and Tyree, M. T. (2000). Hydraulic conductance measured by the high pressure flow meter in crop plants. J. Exp. Bot. 51, 823–828. doi: 10.1093/jexbot/51.345.823
Uehlein, N., Sperling, H., Heckwolf, M., and Kaldenhoff, R. (2012). The Arabidopsis aquaporin PIP1;2 rules cellular CO2 uptake. Plant Cell Environ. 35, 1077–1083. doi: 10.1111/j.1365-3040.2011.02473.x
von Caemmerer, S., and Evans, J. R. (2010). Enhancing C3 photosynthesis. Plant Physiol.154, 589–592. doi: 10.1104/pp.110.160952
Xiong, D., Douthe, C., and Flexas, J. (2018). Differential coordination of stomatal conductance, mesophyll conductance, and leaf hydraulic conductance in response to changing light across species. Plant Cell Environ. 41, 436–450. doi: 10.1111/pce.13111
Xiong, D., and Flexas, J. (2018). Leaf economics spectrum in rice: leaf anatomical, biochemical, and physiological trait trade-offs. J. Exp. Bot. 69, 5599–5609. doi: 10.1093/jxb/ery322
Xiong, D., Flexas, J., Yu, T., Peng, S., and Huang, J. (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
Xiong, D., Liu, X., Liu, L., Douthe, C., Li, Y., et al. (2015). Rapid responses of mesophyll conductance to changes of CO2 concentration, temperature and irradiance are affected by N supplements in rice. Plant Cell Environ. 38, 2541–2550. doi: 10.1111/pce.12558
Zhang, D., Du, Q., Zhang, Z., Jiao, X., Song, X., and Li, J. (2017). Vapour pressure deficit control in relation to water transport and water productivity in greenhouse tomato production during summer. Sci. Rep. 7, 43461. doi: 10.1038/srep43461
Zhang, D., Jiao, X., Du, Q., Song, X., and Li, J. (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
Zhang, D., Li, Y., and Li, Y. (2021). The potential implications of a plasma membrane aquaporin in improving CO2 transport capacity, photosynthetic potential and water use efficiency under contrasting CO2 source in Solanum lycopersicum (tomato). Sci. Hortic. 283: 110122. doi: 10.1016/j.scienta.2021.110122
Zhang, D., Zhang, Z., Li, J., Chang, Y., Du, Q., and Pan, T. (2015). Regulation of vapor pressure deficit by greenhouse micro-fog systems improved growth and productivity of tomato via enhancing photosynthesis during summer season. PLoS ONE. 10, e0133919. doi: 10.1371/journal.pone.0133919
Zhang, S., Feng, M., Chen, W., Zhou, X., Lu, J., Wang, Y., et al. (2019a). In rose, transcription factor PTM balances growth and drought survival via PIP2;1 aquaporin. Nat Plants. 5, 290–299. doi: 10.1038/s41477-019-0376-1
Zhang, Z., Cao, B., Li, N., Chen, Z., and Xu, K. (2019b). Comparative transcriptome analysis of the regulation of ABA signaling genes in different rootstock grafted tomato seedlings under drought stress. Environ Exp Bot. 166, 103814. doi: 10.1016/j.envexpbot.2019.103814
Keywords: abscisic acid, evaporative demand, mesophyll conductance, plant water status, stomatal conductance
Citation: Zhang D, Du Q, Sun P, Lou J, Li X, Li Q and Wei M (2021) Physiological and Transcriptomic Analyses Revealed the Implications of Abscisic Acid in Mediating the Rate-Limiting Step for Photosynthetic Carbon Dioxide Utilisation in Response to Vapour Pressure Deficit in Solanum Lycopersicum (Tomato). Front. Plant Sci. 12:745110. doi: 10.3389/fpls.2021.745110
Received: 21 July 2021; Accepted: 24 September 2021;
Published: 10 November 2021.
Edited by:
Thorsten M. Knipfer, University of British Columbia, CanadaReviewed by:
Wencheng Wang, Huazhong Agricultural University, ChinaDimitrios Fanourakis, Technological Educational Institute of Crete, Greece
Copyright © 2021 Zhang, Du, Sun, Lou, Li, Li and Wei. 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: Min Wei, bWlud2VpJiN4MDAwNDA7c2RhdS5lZHUuY24=