Next Article in Journal
The Mammary Gland: Basic Structure and Molecular Signaling during Development
Previous Article in Journal
Asthma and Allergy: Unravelling a Tangled Relationship with a Focus on New Biomarkers and Treatment
Previous Article in Special Issue
Integration of Light and Auxin Signaling in Shade Plants: From Mechanisms to Opportunities in Urban Agriculture
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Leaf Photosynthesis and Its Temperature Response Are Different between Growth Stages and N Supplies in Rice Plants

1
Ministry of Agriculture and Rural Affairs Key Laboratory of Crop Ecophysiology and Farming System in the Middle Reaches of the Yangtze River, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, China
2
Key Laboratory of Crop Genetics and Physiology of Jiangsu Province, Jiangsu Key Laboratory of Crop Cultivation and Physiology, Co-Innovation Center for Modern Production Technology of Grain Crops, Yangzhou University, Yangzhou 225009, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(7), 3885; https://doi.org/10.3390/ijms23073885
Submission received: 26 February 2022 / Revised: 18 March 2022 / Accepted: 29 March 2022 / Published: 31 March 2022
(This article belongs to the Special Issue Response of Photosynthesis to Abiotic Stresses)

Abstract

:
Leaf photosynthesis is highly correlated with CO2-diffusion capacities, which are determined by both leaf anatomical traits and environmental stimuli. In the present study, leaf photosynthetic rate (A), stomatal conductance (gs), mesophyll conductance (gm) and the related leaf anatomical traits were studied on rice plants at two growth stages and with two different N supplies, and the response of photosynthesis to temperature (T) was also studied. We found that gm was significantly higher at mid-tillering stage and at high N treatment. The larger gm was related to a larger chloroplast surface area facing intercellular air spaces and a thinner cell wall in comparison with booting stage and zero N treatment. At mid-tillering stage and at high N treatment, gm showed a stronger temperature response. The modelling of the gm-T relationships suggested that, in comparison with booting stage and zero N treatment, the stronger temperature response of gm was related to the higher activation energy of the membrane at mid-tillering stage and at high N treatment. The findings in the present study can enhance our knowledge on the physiological and environmental determinants of photosynthesis.

1. Introduction

Rice is one of the most important cereal crops in the world, feeding more than half of the world’s population [1]. Studies relating to the physiological and anatomical determinants on rice photosynthesis, and to the response of photosynthesis to environmental stimuli, are of great importance to further improve rice yield. Photosynthesis in C3 plants, including rice plants, is limited by both CO2-diffusion capacities and biochemical functions [2,3,4,5,6,7]. Before being fixed by the key Calvin cycle enzyme of Rubisco, CO2 in the air should firstly diffuse across stomata to reach the substomatal cavity; thereafter, it will further diffuse across the cell wall, plasma membrane, cytoplasm and chloroplast envelope to reach the carboxylation sites [8,9]. The CO2-diffusion capacities through stomata and mesophyll cells are called stomatal conductance (gs) and mesophyll conductance (gm), respectively. Leaf photosynthesis has been frequently found to positively correlate with stomatal conductance and mesophyll conductance [10,11,12].
Both gs and gm are closely related to leaf anatomical traits. Previous studies have shown that stomatal conductance in rice plants is not correlated with either stomatal size or stomatal number [13,14], but it is significantly correlated with leaf hydraulic conductance (Kleaf) [15]. The two major determinants of leaf hydraulic conductance are leaf vein density and xylem size [16,17,18], because more leaf veins can provide more parallel water flow paths through the vein system [19] and hydraulic conductance through leaf xylems is positively correlated with the xylem conduits diameter [18]. In rice plants, Kleaf is positively related to the area of xylem conduits within the bundle sheath, but it is not related to leaf vein density [15]. Both leaf vein density and xylem size can vary largely across different growth stages, which may in turn have a great effect on gs. Therefore, the first objective of this study was to investigate whether the variation in gs across different growth stages is related to the variations in leaf vein density and xylem size.
Cell-wall thickness (Tcw) and chloroplast surface area facing intercellular air spaces (Sc) are two important leaf anatomical traits determining gm. It has been frequently found that gm is negatively correlated with cell-wall thickness (Tcw), and it is positively correlated with chloroplast surface area facing intercellular air spaces (Sc) [7,20]. Leaf N content has a significant effect on Tcw and Sc. In comparison with low N supply, high N supply can both decrease Tcw and increase Sc, which are the major reasons for the increased gm under high N supply [3,21]. There are large variations in Tcw and Sc, which can subsequently have a large effect on gm. Therefore, the second objective of this study was to investigate the responses of gm, Tcw and Sc to growth stages under two different N supplies, which can improve our understanding on the determinants of gm.
Photosynthesis is sensitive to environmental variations, and temperature (T) is one of the major environmental stimuli that have great impacts on photosynthesis and crop production. However, it is not known whether the impact of temperature on photosynthesis varies across different growth stages, although it is known that N supply has a significant effect on the temperature response of photosynthesis in rice plants [22]. The variation in photosynthesis in response to temperature is largely related to mesophyll conductance [23], and the mechanisms relating to temperature response of gm have been intensively studied [24,25,26,27]. CO2 diffusion through mesophyll cells can be divided into two processes, which are the liquid phase and membrane phase [25]. The liquid phase refers to CO2 diffusion through cell walls, cytoplasm and chloroplast stroma; the membrane phase refers to CO2 diffusion through plasma membrane and chloroplast envelope. It has been hypothesized that the sensitivity of mesophyll conductance to temperature (Ea,gm) is determined by the activation energy of membrane (Ea,mem) and by the ratio of CO2-diffusion conductance through liquid phase to that through membrane phase (gliq/gmem) [25]. The variations in leaf anatomical traits in response to growth stage and N supply may change the ratio of gliq/gmem, and thus the sensitivity of mesophyll conductance to temperature. Therefore, the third objective of this study was to investigate the differential responses of gm to temperature at different growth stages and N supplies.
In addition to Ea,mem and gliq/gmem, the response of leaf water potential (Ψleaf) to temperature is also an important determinant in the response of gm to temperature [27]. The increment of gm with temperature can be substantially inhibited if Ψleaf decreases with temperature. It is not known whether Ψleaf can significantly decrease in response to increasing temperature at the late growth stage in rice plants, although Ψleaf is insensitive to temperature at tillering stage [23]. Therefore, the fourth objective was to investigate the variation in Ψleaf in response to the temperature, and its influence on the response of gm to temperature.
To this end, a rice cultivar of Fengliangyouxiang 1 was grown in pots under two N supplies. At both tillering and booting stages, the gas-exchange parameters, the leaf anatomical traits and the temperature response of photosynthesis were investigated. We hypothesized that Tcw is significantly higher at booting stage than that at tillering stage, which in turn leads to a lower gm and to a lower sensitivity of gm to temperature at booting stage. The results will potentially improve our knowledge on rice photosynthesis, which will be beneficial to improving leaf photosynthesis and crop yields.

2. Results

2.1. Response of Leaf N Content to N Supplies and Growth Stages

Regardless of growth stages, area-based leaf N content (Narea) under high N (HN) treatment was significantly higher than that under zero N (N0) treatment (Figure 1a). In comparison with N0 treatment, Narea under HN treatment was increased by 54.7% and 21.1% at mid-tillering and booting stages, respectively. Under HN treatment, Narea showed no obvious difference between the two growth stages; under N0 treatment, however, Narea at booting stage was 29.8% higher than that at mid-tillering stage.

2.2. The Responses of Gas-Exchange Parameters to Temperature under Different N Supplies and at Different Growth Stages

Net photosynthetic rate, stomatal conductance and mesophyll conductance at 25 °C were represented by A25, gs,25 and gm,25, respectively. At mid-tillering stage, A25, gs,25 and gm,25 under HN treatment were significantly higher than those under N0 treatment (Table 1). In comparison with N0 treatment, A25, gs,25 and gm,25 under HN treatment were increased by 55.8%, 74.6% and 101.5%, respectively. At booting stage, however, they were not significantly different between the two N treatments, although Narea under HN treatment was significantly higher than that under N0 treatment (Figure 1a).
Under HN treatment, A25, gs,25 and gm,25 at mid-tillering stage were 67.4%, 113.9% and 135.7%, respectively, higher than those at booting stage (Table 1), although Narea was similar between the two growth stages under HN treatment (Figure 1a). Under N0 treatment, A25, gs,25 and gm,25 at mid-tillering stage were not significantly different with those at booting stage (Table 1), although Narea at mid-tillering stage was lower than that at booting stage (Figure 1a). Consequently, photosynthetic nitrogen-use efficiency (PNUE), which was calculated as the ratio of A25/Narea, was significantly higher at mid-tillering stage than that at booting stage, although it was not significantly different between the two N treatments (Figure 1b).
Both A and gm increased dramatically with the increasing temperature; in contrast, temperature had no significant effect on gs (Table 1). Across different temperatures, A was significantly correlated with gm, but it was not correlated with gs (Figure 2). At mid-tillering stage, the modelled values of Ea,gm were 38.4 and 28.2 kJ mol−1 under HN and N0 treatments, respectively, and the modelled value at booting stage was 29.9 kJ mol−1 under HN treatment (Table 2). This suggested that gm was less sensitive to temperature at N0 treatment and at booting stage.
Using the two-component model, the fitted membrane permeability to CO2 at 25 °C (Pmem,25) varied from 0.722 mm s−1 to 1.272 mm s−1, and Ea,mem varied from 60.4 kJ mol−1 to 81.0 kJ mol−1 (Table 3), which fell well within the reported ranges [2,26]. The modelled Pmem,25 and Ea,mem under HN treatment were larger than those under N0 treatment, and they were also larger at mid-tillering stage than those at booting stage.

2.3. Effects of N Supply and Growth Stage on Leaf Anatomical Traits

N supply had significant effects on Sc and Tcw, but mesophyll cell surface area facing intercellular air spaces (Sm) was not significantly different between N supplies (Table 4). In comparison with N0 treatment, Sc under HN treatment was increased by 25.0% and 21.8% at mid-tillering and booting stages, respectively. At mid-tillering stage, Tcw under HN treatment was significantly lower than that under N0 treatment; at booting stage, however, N supply had no significant effect on Tcw. In comparison with mid-tillering stage, Sc was significantly lower, while Sm and Tcw were significantly higher at booting stage. Across different N supplies and growth stages, gm,25 was negatively correlated with Tcw (Figure 3), but it was not significantly related to either Sm or Sc (Figure S1).

2.4. The Responses of Leaf Hydraulic Traits to Temperature under Different N Supplies and at Different Growth Stages

Regardless of growth stage, temperature had no significant effect on Ψleaf under HN treatment; under N0 treatment, however, Ψleaf was more negative at 35 °C than that at 15 °C (Table 5). Ψleaf was not significantly different across different N supplies or growth stages.
Regardless of growth stage, temperature had no significant effect on Kleaf (Table 5). In contrast, Kleaf under HN was significantly higher than that under N0 treatment. In comparison with N0 treatment, Kleaf under HN treatment was increased by 37.3% and 43.5% at mid-tillering and booting stages, respectively. Under both N supplies, Kleaf at mid-tillering stage was significantly higher than that at booting stage.
Leaf venation traits were significantly affected by N supplies and growth stages (Table S1). Regardless of growth stage, the inter-vein distance between major veins (IVDmajor) and the inter-vein distance between minor veins (IVDminor) under HN treatment were significantly higher than those under N0 treatment, although IVDminor was not significantly different between different N supplies at mid-tillering stage. At mid-tillering stage, the area of xylem conduits in veins per leaf width (Sx) under HN treatment was significantly lower than that under N0 treatment; at booting stage, however, Sx was not significantly different between two N supplies. In comparison with mid-tillering stage, IVDmajor, IVDminor and Sx were all larger at booting stage.

3. Discussion

The results obtained in the present study support our hypotheses that, in comparison with mid-tillering stage, gm was lower and was less sensitive to temperature at booting stage. Moreover, gm was significantly larger and was more sensitive to temperature under HN treatment than that under N0 treatment.

3.1. The Variation in Photosynthesis between Growth Stage and N Supply Is Related to CO2-Diffusion Capacity

The photosynthetic rate in C3 plants at an ambient CO2 concentration of ~400 μmol mol−1 is suggested to be limited by Rubisco carboxylation capacity, which is related to both the content and the specific activity of Rubisco [28,29], and the latter is dependent on CO2 partial pressure inside chloroplasts [3,21]. However, due to the significant resistance during the CO2-diffusion pathways through stomata and mesophyll cells, CO2 partial pressure inside chloroplasts is usually not saturated for C3 plants [24,30,31]. Therefore, CO2-diffusion capacities, namely gs and gm, are major limitations to C3 photosynthesis [7,32]. In the present study, across different N treatments and growth stages, the variation trends of gs and gm were identical to that of net photosynthetic rate (Table 1). This suggested that the variation in photosynthesis between different growth stages and between different N supplies is also related to CO2-diffusion capacity.
The anatomical trait of Sc is an important determinant to gm [33], and the value of Sc is related to chloroplast development and N supply [3,21]. The Sc value has been found to positively relate to Narea [34]. In the present study, however, the variations of Sc and Narea were uncoupled between two growth stages (Figure 1a and Table 4). The values of Sc at booting stage were significantly lower than those at mid-tillering stage (Table 4), although Narea at booting stage was similar to, or was even larger than, that at mid-tillering stage (Figure 1). The unparallel changes of Sc and Narea between different growth stages may be caused by differential N partitions. It is suggested that 10% of leaf N is distributed in cell walls [35], but the proportion can be increased to 30% in leaves with large leaf mass per area (LMA) and thick cell walls [36]. In the present study, Tcw was dramatically higher at booting stage than that at mid-tillering stage (Table 4). This suggested that N distribution to cell walls may be larger at booting stage than that at mid-tillering stage, which may in turn resulted in the unparallel changes of Sc and Narea between different growth stages.
In addition to Sc, Tcw is also an important anatomical determinant in gm, which is negatively correlated with Tcw [2,20]. In line with previous studies [3,21], at mid-tillering stage, the larger gm at HN treatment in comparison with N0 treatment can be explained by the increased Sc and the decreased Tcw (Table 1 and Table 4). In contrast, with a similar Tcw between two N supplies at booting stage, the larger Sc at HN treatment did not lead to an increased gm in comparison to N0 treatment (Table 1 and Table 4). Moreover, gm was significantly correlated with Tcw across different growth stages and N supplies (Figure 3), while it was not significantly correlated with Sc (Figure S1). This suggested that Tcw is a more important anatomical trait than Sc in determining gm, and the manipulation of cell walls is suggested to be an efficient approach to improve leaf photosynthesis [37]. In comparison with booting stage, the larger gm at mid-tillering stage was related to the higher Sc and the lower Tcw.
Stomatal conductance has been frequently found to be positively correlated with Kleaf. In the present study, however, gs was not correlated with either Kleaf or leaf venation traits (Table 1, Table 5 and Table S1). In the soil-plant-atmosphere continuum, the water diffusion resistance through leaves contributes ~30% of the whole plant hydraulic resistance [38]. In contrast, recent studies have suggested that root hydraulic resistance, including the resistance through radial pathway from root surface to the xylem and that through the root-soil interface, is the major resistance for water diffusion [39,40]. The relative resistance through leaves and roots in rice plants is not known, but we speculate that root hydraulic conductance may vary with growth stages and N supplies, which in turn determines the variation in gs.

3.2. Temperature Response of gm Varies with Growth Stage and N Supply

Photosynthesis is sensitive to temperature, increasing with increasing temperature, but decreasing dramatically at supra-optimal temperature [27]. The increase in A with temperature is suggested to at least partially correlate with the increase in gm [23]. In line with previous studies, in the present study, the response of A to temperature was similar to that of gm, while the response of gs to temperature was different from both A and gm (Table 1). This suggested that the response of photosynthesis to temperature is mainly driven by the gm-T relationship.
The gm-T relationship is correlated with the variation in Ψleaf in response to temperature [27], because leaf dehydration can severely depress gm, probably through the deactivation of aquaporins [5,41]. Therefore, the decrease in Ψleaf in response to temperature can lead to a lower increment of gm with temperature [27]. In the present study, however, there was little variation in Ψleaf in response to temperature, although Ψleaf was significantly decreased at N0 treatment at both growth stages (Table 5). This suggested that the differential sensitivity of gm to temperature between growth stages and N treatments is not related to the variation in Ψleaf.
The two-component model hypothesized that gm would show a strong temperature response if Ea,mem were large and the ratio of gliq,25’/gmem,25’ were high [25], where gliq,25’ and gmem,25’ represent the CO2 conductance through the liquid phase per Sc and the CO2 conductance through the membrane phase per Sc, respectively. The modelling of the gm-T relationships suggested that, in comparison with N0 treatment, Ea,mem was larger, while gliq,25’/gmem,25’ was lower at HN treatment; in comparison with booting stage, Ea,mem was larger while gliq,25’/gmem,25’ was comparable at mid-tillering stage (Table 3). This suggested that the larger Ea,mem was accounted for by the stronger temperature response of gm at mid-tillering stage and at HN treatment. However, the factors that determine Ea,mem are not known, but the membrane compositions, such as cholesterol and aquaporins, are suggested to affect Ea,mem [26]. More research is needed to investigate the mechanisms underlying the differential Ea,mem.

4. Materials and Methods

4.1. Plant Materials and N Treatments

A rice cultivar of Fengliangyouxiang 1, which has been widely grown locally, was planted in pots in Huazhong Agricultural University (114.37° E, 30.48° N), Wuhan, Hubei province, China. After germination on moist filters on 15 July 2019, seeds were transferred to nursery boxes. When the seedlings had developed an average of three leaves, which usually requires 15 days, they were transplanted to 11.0 L pots with a density of three hills per pot and two seedlings per hill. Each pot was filled with 10.0 kg of soil. Phosphorus (P) and potassium (K) were applied as basal fertilizers at the rates of 1.50 and 1.89 g pot−1, respectively, in the form of KH2PO4. N was applied with urea at a rate of 1.60 g pot−1 at HN treatment, 40% of which was applied as the basal fertilizer, and another two topdressings of 30% each were applied at mid-tillering and booting stages. No N was applied for N0 treatment. The soil used in this study had the following properties: pH 7.1, 6.7 g kg−1 of organic matter, 6.27 mg kg−1 of Olsen-P, 129 mg kg−1 of exchangeable K, and 0.63‰ total N. There were 10 pots per treatment. Plants were irrigated daily with tap water, and a minimum 2 cm water layer was maintained to avoid drought stress. Measurements were conducted on the newest fully expanded leaves at mid-tillering stage and on the flag leaves at booting stage.

4.2. Gas-Exchange Measurements

To minimise the effects of environmental fluctuations and midday depression on photosynthesis, rice plants were transferred to an environmentally controlled growth chamber (Conviron GR48, Controlled Environments Ltd., Winnipeg, MB, Canada) in the afternoon before the day of measurement. The air temperatures of the growth chamber were controlled to match the desired leaf temperatures of 15, 25 or 35 °C. The CO2 concentration and light intensity in the growth chamber were controlled at 400 μmol mol−1 and 1000 μmol m−2 s−1, respectively. A portable photosynthesis system (Licor-6800; Li-Cor Inc., Lincoln, NE, USA) with an integrated fluorescence leaf chamber (6800-01A) was used to measure leaf gas exchange and chlorophyll fluorescence between 08:00 and 16:00. Before the measurement, the leaves were attached to the leaf chamber to stabilize. CO2 concentration inside the leaf chamber was controlled to 400 µmol mol−1, and photosynthetic photon flux density (PPFD) was set to 1500 µmol m−2 s−1. The vapour-pressure deficit between leaf and air (VPD) increased dramatically with temperature. When leaf photosynthetic parameters were stabilised, which usually takes 15–25 min, gas-exchange parameters and chlorophyll fluorescence were simultaneously recorded with a light saturating pulse of 8000 µmol m−2 s−1. The actual photochemical efficiency of photosystem II (ΦPSII) was calculated as follows:
Φ PSII = F m F s F m
where Fs and Fm′ are steady-state fluorescence and the maximum fluorescence, respectively. The electron transport rate (J) was calculated as follows:
J = P P F D × α × β × Φ PSII
where α is the leaf absorptance and β is the partitioning of absorbed quanta between photosystem II and photosystem I. The product α × β was determined from the slope of the relationship between ΦPSII and the quantum efficiency of CO2 uptake (ΦCO2), which was measured by varying light intensity under non-photorespiratory conditions at <2% O2 [42].
The variable J method described in Harley et al. [43] was used to calculate chloroplastic CO2 concentration (Cc) and gm:
C c = Γ * × ( J + 8 * ( A + R d ) ) J 4 * ( A + R d )
g m = A C i C c
where Ci is the intercellular CO2 concentration, Γ* is the CO2 compensation point in the absence of day respiration, and Rd is the day respiration rate. When calculating gm at different temperatures, the values of Γ* and Rd at different temperatures were calculated using the following equation:
P = e ( c E a R × ( 273 + T ) )
where P is calculated parameter, c is the scaling factor, Ea is the activation energy and R is the molar gas constant of 8.314 J K−1 mol−1. The values of c and Ea for Γ* and Rd were taken from Bernacchi et al. [24,44]. The values of c and Ea for Γ* are 13.49 and 24.46 kJ mol−1, respectively, while for Rd, they are 18.72 kJ mol−1 and 46.39 kJ mol−1, respectively. At booting stage, gas-exchange measurements under N0 treatment were conducted at 25 °C only (Table 1). The leaf photosynthetic rate, stomatal conductance and mesophyll conductance at 25 °C are represented by A25, gs,25 and gm,25, respectively.

4.3. Measurements of Leaf Hydraulic Parameters

After gas-exchange measurements were complete, the leaves were immediately detached and placed in a previously exhaled-in and sealable bag. After equilibration for at least 10 min, Ψleaf was measured using a pressure chamber (PMS Instrument Company, Albany, OR, USA).
The measurements of Kleaf at different temperatures were conducted in the same growth chamber that was used for the gas-exchange measurements. The environments during the measurement of Kleaf were similar to those used for the gas-exchange measurements. Kleaf was measured using the evaporating flux method (EFM) [45]. Briefly, the newest fully expanded leaves were excised with a fresh razor blade, and then immediately recut under water. Then, the leaf was connected to silicone tubing with a compression fitting under water to prevent air entering the system. The tubing connected the leaf to a hard tube connected to a graduated cylinder on a balance capable of reading 0.1 mg. The balance logged data every 30 s to a computer. The excised leaves were placed under LED lights for transpiration; the PPFD at the leaf level was 1500 μmol m−2 s−1. After equilibration to a steady state, which required ∼30 min after excising the leaves, leaf transpiration rate (Tr) was calculated after measuring the leaf area. Afterwards, the leaves were immediately placed in a previously exhaled-in and sealable bag. After equilibration for at least 20 min, Ψleaf was measured using a pressure chamber. The unnormalised Kleaf (Kleaf′) was calculated as follows [46]:
K leaf = T r 0 Ψ leaf
During the measurement, leaf temperature was measured using a Multi-channel Digital Thermometer (AZ88598, AZ Instrument Corp. Ltd., Taichung, China). Water viscosity has significant effects on leaf hydraulic traits [47,48]. To exclude these effects, leaf hydraulic conductance at different temperatures were normalised to the water viscosity at 25 °C:
K leaf = K leaf × μ μ 25
where μ is the water viscosity at the measured leaf temperature and μ25 is the water viscosity at 25 °C.

4.4. Leaf N Content Measurement

Immediately after the gas-exchange measurements, newly expanded leaves were detached to measure leaf area using a LI-Cor 3000C (LI-COR Inc., Lincoln, NE, USA) leaf area analyser. Then, the leaves were oven-dried at 80 °C to reach a constant weight. Afterwards, leaf dry mass was weighed, and LMA was calculated as the ratio of leaf dry mass to leaf area. Mass-based leaf N content (Nmass, %) was measured using a stable isotope ratio mass spectrometer (IsoPrime 100, IRMS, Isoprime Ldt., Cheadle, UK), and Narea was calculated as: Narea = Nmass × LMA. PNUE was calculated as: PNUE = A25/Narea.

4.5. Measurements of Leaf Anatomical Traits

Paraffin and ultrathin sections were made from three leaves per treatment to analyse leaf anatomy using light microscope (LM) and transmission electron microscope (TEM). For the paraffin sections, leaf discs of about 5.0 mm length were cut from the middle of the leaves, and they were then fixed in FAA buffer (5% formaldehyde, 5% glacial acetic acid and 63% alcohol (v/v) in pure water) at 4 °C for 24 h. Thereafter, they were vacuumed in a vacuum chamber (DZF-6050, Shanghai Hasuc Co., Ltd., Shanghai, China). The samples were embedded in paraffins, and the leaf cross-sections were made by professionals from Wuhan Google Biotechnology Co. Ltd. The paraffin sections were stained with safranin-fast green, and they were photographed at a magnification of ×300 with a Nikon Eclipse E100 LM (Nikon Optical, Tokyo, Japan). There were 6–9 LM images taken for each treatment. The Sx, IVDmajor and IVDminor were measured using the ImageJ software.
For the ultrathin sections, small leaf sections of 2.0 × 2.0 mm were cut from middle of the leaves (avoiding midrib). The leaf sections were infiltrated with fixative 2.5% (v/v) glutaric aldehyde in 0.1 M phosphate buffer (pH = 7.6) in a vacuum chamber for 2 h. Ultrathin sections were made from Wuhan Google Biotechnology Co. Ltd. Images were acquired using a transmission electron microscope (H-7650; Hitachi-Science & Technology, Tokyo, Japan). The LM and TEM images (Figures S2 and S3) were used to measure Sc and Sm following the methods described in Evans et al. [34] and Ye et al. [49]. Cell-wall thickness was measured with ×10,000 TEM images using the ImageJ software.

4.6. Quantification of the Temperature Response of gm

The response of gm to temperature was modelled using the equation:
g m = e ( c E a , gm R × T k ) 1 + e ( Δ S × T k E d , g m R × T k )
where ΔS is an entropy term and Ed,gm is a term for deactivation of gm [24].

4.7. Modelling the Temperature Response of gm

In order to interpret the difference in the gm-T relationships between growth stages and N supplies, we used the equations in von Caemmerer and Evans [26] to model the temperature response of gm. Generally, gm was separated into liquid phase and membrane phase:
g m = 1 1 g liq + 1 g mem
The CO2 conductance through the liquid phase per Sc ( g liq = g liq S c ) can be given by
g liq = ρ H D l
where ρ (mol m−3) is the molar density of water, H (bar−1) is the Henry coefficient for CO2, D (m2 s−1) is the diffusivity of CO2 in water and l (m) is the effective pathlength.
Solubility of CO2 in water decreases with temperature and
ρ H = 33.06 × e ( 2400 × ( 1 273 + T 1 298 ) )
Diffusivity of CO2 in water increases with temperature:
D = 1 . 81 × 10 6 × e ( 16900 R ( 273 + T ) )
l was calculated from Tcw and cell-wall porosity (Pcw) [50]:
l = T c w P c w
Pcw varied with Tcw according to Tosens et al. [51]:
P c w = 0.3733 × T c w + 0.3378
The temperature dependence of CO2-diffusion across biological membranes per Sc ( g mem = g mem S c ) is assumed to be exponential:
g mem = ρ H P m e m , 25 × e ( T 25 ) × E a , m e n R × 298 × ( 273 + T )

4.8. Statistical Analysis

One-way and two-way analyses of variance (ANOVA) were used to assess the effects of leaf temperature, N supply and growth stage, as well as their interactions on parameters using Statistix 9.0 software (Analytical Software, Tallahassee, FL, USA). Parameters were compared between treatments based on the least significant difference (LSD) test level at the 0.05 probability level. Graphs were created, and a linear regression analysis was performed to test the correlations between parameters using SigmaPlot 10.0 (Systat Software Inc., San Jose, CA, USA).

5. Conclusions

Leaf photosynthesis can be significantly affected by growth stage and N supply, and the larger photosynthetic rate at mid-tillering stage and with HN treatment was related to the higher gs and gm. In comparison with Sc, Tcw is a more important anatomical trait in determining gm at different growth stages. The response of leaf photosynthesis to temperature can also be affected by growth stage and N supply, and these effects are related to the strong temperature response of gm. The stronger response of gm to temperature at mid-tillering stage and HN treatment was related to the larger Ea,mem.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23073885/s1.

Author Contributions

Y.L. conceived and designed the research. Z.Z. and M.Y. conducted the experiments and collected the data. Y.L., M.Y. and G.H. analysed the data and wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was Funded by National Natural Science Foundation of China 31871532 and 32172103; Fundamental Research Funds for the Central Universities 2021ZKPY017; China Postdoctoral Science Foundation 2021M702769.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kathuria, H.; Giri, J.; Tyagi, H.; Tyagi, A.K. Advances in transgenic rice biotechnology. Crit. Rev. Plant Sci. 2007, 26, 65–103. [Google Scholar] [CrossRef]
  2. Evans, J.R.; Kaldenhoff, R.; Genty, B.; Terashima, I. Resistances along the CO2 diffusion pathway inside leaves. J. Exp. Bot. 2009, 60, 2235–2248. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Li, Y.; Gao, Y.; Xu, X.; Shen, Q.; Guo, S. Light-saturated photosynthetic rate in high-nitrogen rice (Oryza sativa L.) leaves is related to chloroplastic CO2 concentration. J. Exp. Bot. 2009, 60, 2351–2360. [Google Scholar] [CrossRef] [PubMed]
  4. Yamori, W.; Nagai, T.; Makino, A. The rate-limiting step for CO2 assimilation at different temperatures is influenced by the leaf nitrogen content in several C3 crop species. Plant Cell Environ. 2011, 34, 764–777. [Google Scholar] [CrossRef] [PubMed]
  5. Flexas, J.; Barbour, M.M.; Brendel, O.; Cabrera, H.M.; Carriqui, M.; Diaz-Espejo, A.; Douthe, C.; Dreyer, E.; Ferrio, J.P.; Cago, J.; et al. Mesophyll diffusion conductance to CO2: An unappreciated central player in photosynthesis. Plant Sci. 2012, 193, 70–84. [Google Scholar] [CrossRef]
  6. Adachi, S.; Nakae, T.; Uchida, M.; Soda, K.; Takai, T.; Oi, T.; Hirasawa, T. The mesophyll anatomy enhancing CO2 diffusion is a key trait for improving rice photosynthesis. J. Exp. Bot. 2013, 64, 1061–1072. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Gago, J.; Carriquí, M.; Nadal, M.; Clemente-Moreno, M.J.; Coopman, R.E.; Fernie, A.R.; Flexas, J. Photosynthesis optimized across land plant phylogeny. Trends Plant Sci. 2019, 24, 947–958. [Google Scholar] [CrossRef] [PubMed]
  8. Terashima, I.; Hanba, Y.T.; Tholen, D.; Niinemets, Ü. Leaf functional anatomy in relation to photosynthesis. Plant Physiol. 2011, 155, 108–116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Tholen, D.; Boom, C.; Zhu, X.G. Opinion: Prospects for improving photosynthesis by altering leaf anatomy. Plant Sci. 2012, 197, 92–101. [Google Scholar] [CrossRef] [PubMed]
  10. Giuliani, R.; Koteyeva, N.; Voznesenskaya, E.; Evans, M.A.; Cousins, A.B.; Edwards, G.E. Coordination of leaf photosynthesis, transpiration, and structural traits in rice and wild relatives (Genus Oryza). Plant Physiol. 2013, 162, 1632–1651. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Carriquí, M.; Cabrera, H.M.; Conesa, M.À.; Coopman, R.E.; Douthe, C.; Gago, J.; Gallé, A.; Galmés, J.; Ribas-Carbo, M.; Tomás, M.; et al. Diffusional limitations explain the lower photosynthetic capacity of ferns as compared with angiosperms in a common garden study. Plant Cell Environ. 2015, 38, 448–460. [Google Scholar] [CrossRef]
  12. Liu, X.; Li, Y. Varietal difference in the correlation between leaf nitrogen content and photosynthesis in rice (Oryza sativa L.) plants is related to specific leaf weight. J. Integr. Agr. 2016, 15, 2002–2011. [Google Scholar] [CrossRef]
  13. Zhang, Q.Q.; Peng, S.B.; Li, Y. Increase rate of light-induced stomatal conductance is related to stomatal size in the genus Oryza. J. Exp. Bot. 2019, 70, 5259–5269. [Google Scholar] [CrossRef]
  14. Zhang, Q.Q.; Yang, Y.H.; Peng, S.B.; Li, Y. Nighttime transpirational cooling enabled by circadian regulation of stomatal conductance is related to stomatal anatomy and leaf morphology in rice. Planta 2021, 254, 12. [Google Scholar] [CrossRef] [PubMed]
  15. Huang, G.J.; Shu, Y.; Peng, S.B.; Li, Y. Leaf photosynthesis is positively correlated with xylem and phloem areas in leaf veins in rice (Oryza sativa) plants. Ann. Bot. 2022, 2022, mcac020. [Google Scholar] [CrossRef]
  16. Flexas, J.; Scoffoni, C.; Gago, J.; Sack, L. Leaf mesophyll conductance and leaf hydraulic conductance: An introduction to their measurement and coordination. J. Exp. Bot. 2013, 64, 3965–3981. [Google Scholar] [CrossRef] [PubMed]
  17. Sack, L.; Scoffoni, C.; John, G.P.; Poorter, H.; Mason, C.M.; Mendez-Alonzo, R.; Donovan, L.A. How do leaf veins influence the worldwide leaf economic spectrum? Review and synthesis. J. Exp. Bot. 2013, 64, 4053–4080. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. North, G.B.; Lynch, F.H.; Maharaj, F.D.; Phillips, C.A.; Woodside, W.T. Leaf hydraulic conductance for a tank bromeliad: Axial and radial pathways for moving and conserving water. Front. Plant Sci. 2013, 4, e78. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Buckley, T.N. The contributions of apoplastic, symplastic and gas phase pathways for water transport outside the bundle sheath in leaves. Plant Cell Environ. 2015, 38, 7–22. [Google Scholar] [CrossRef] [PubMed]
  20. Evans, J.R. Mesophyll conductance: Walls, membranes and spatial complexity. New Phytol. 2020, 229, 1864–1876. [Google Scholar] [CrossRef] [PubMed]
  21. Li, Y.; Ren, B.; Ding, L.; Shen, Q.; Peng, S.; Guo, S. Does chloroplast size influence photosynthetic nitrogen use efficiency? PLoS ONE 2013, 8, e62036. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Huang, G.J.; Zhang, Q.Q.; Wei, X.H.; Peng, S.B.; Li, Y. Nitrogen can alleviate the inhibition of photosynthesis caused by high temperature stress under both steady-state and flecked irradiance. Front. Plant Sci. 2017, 8, 945. [Google Scholar] [CrossRef] [PubMed]
  23. Yang, Y.; Zhang, Q.Q.; Huang, G.J.; Peng, S.B.; Li, Y. Temperature responses of photosynthesis and leaf hydraulic conductance in rice and wheat. Plant Cell Environ. 2020, 43, 1437–1451. [Google Scholar] [CrossRef] [PubMed]
  24. Bernacchi, C.J.; Portis, A.R.; Nakano, H.; von Caemmerer, S.; Long, S.P. Temperature response of mesophyll conductance. Implications for the determination of Rubisco enzyme kinetics and for limitations to photosynthesis in vivo. Plant. Physiol. 2002, 130, 1992–1998. [Google Scholar] [CrossRef] [Green Version]
  25. Evans, J.R.; von Caemmerer, S. Temperature response of carbon isotope discrimination and mesophyll conductance in tobacco. Plant Cell Environ. 2013, 36, 745–756. [Google Scholar] [CrossRef] [PubMed]
  26. von Caemmerer, S.; Evans, J.R. Temperature responses of mesophyll conductance differ greatly between species. Plant Cell Environ. 2015, 38, 629–637. [Google Scholar] [CrossRef] [PubMed]
  27. Li, Y.; Song, X.; Li, S.; Salter, W.T.; Barbour, M.M. The role of leaf water potential in the temperature response of mesophyll conductance. New Phytol. 2020, 225, 1193–1205. [Google Scholar] [CrossRef] [PubMed]
  28. Makino, A.; Shimada, T.; Takumi, S.; Kaneko, K.; Matsuoka, M.; Shimamoto, K.; Nakano, H.; Miyao-Tokutomi, M.; Mae, T.; Yamamoto, N. Does decrease in ribulose-1,5-bisphosphate carboxylase by antisense rbcS lead to a higher N-use efficiency of photosynthesis under conditions of saturating CO2 and light in rice plants? Plant Physiol. 1997, 114, 483–491. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Spreitzer, R.J.; Salvucci, M.E. Rubisco: Structure, regulatory interactions, and possibilities for a better enzyme. Annu. Rev. Plant Biol. 2002, 53, 449–479. [Google Scholar] [CrossRef] [Green Version]
  30. Leegood, R.C.; Sharkey, T.D.; von Caemmerer, S. Photosynthesis: Physiology and Metabolism; Kluwer Academic Publishers: New York, NY, USA, 2004. [Google Scholar]
  31. Flexas, J.; Ribas-Carbó, M.; Hansom, D.T.; Bota, J.; Otto, B.; Cifre, J.; McDowell, N.; Medrano, H.; Kaldenhoff, R. Tobacco aquaporin NtAQP1 is involved in mesophyll conductance to CO2 in vivo. Plant J. 2006, 48, 427–439. [Google Scholar] [CrossRef]
  32. Flexas, J.; Carriquí, M. Photosynthesis and photosynthetic efficiencies along the terrestrial plant’s phylogeny: Lessons for improving crop photosynthesis. Plant J. 2020, 101, 964–978. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Evans, J.R.; von Caemmerer, S.; Setchell, B.A.; Hudson, S. The relationship between CO2 transfer conductance and leaf anatomy in transgenic tobacco with a reduced content of Rubisco. Funct. Plant. Biol. 1994, 21, 475–495. [Google Scholar] [CrossRef]
  34. Xiong, D.; Liu, X.; Liu, L.; Douthe, C.; Li, Y.; Peng, S.; Huang, J. Rapid responses of mesophyll conductance to changes of CO2 concentration, temperature and irradiance are affected by N supplies in rice. Plant Cell Environ. 2015, 38, 2541–2550. [Google Scholar] [CrossRef]
  35. Evans, J.R.; Clarke, V.C. The nitrogen cost of photosynthesis. J. Exp. Bot. 2019, 70, 7–15. [Google Scholar] [CrossRef]
  36. Harrison, M.T.; Edwards, E.J.; Farquhar, G.D.; Nicotra, A.B.; Evans, J.R. Nitrogen in cell walls of sclerophyllous leaves accounts for little of the variation in photosynthetic nitrogen-use efficiency. Plant Cell Environ. 2009, 32, 259–270. [Google Scholar] [CrossRef] [PubMed]
  37. Flexas, J.; Clemente-Moreno, M.J.; Bota, J.; Brodribb, T.J.; Gago, J.; Mizokami, Y.; Nadal, M.; Perera-Castro, A.V.; Roig-Oliver, M.; Sugiura, D.; et al. Cell wall thickness and composition are involved in photosynthetic limitation. J. Exp. Bot. 2021, 72, 3971–3986. [Google Scholar] [CrossRef]
  38. Sack, L.; Holbrook, N.M. Leaf hydraulics. Annu. Rev. Plant. Biol. 2006, 57, 361–381. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Rodriguez-Dominguez, C.M.; Brodribb, T.J. Declining root water transport drives stomatal closure in olive under moderate water stress. New Phytol. 2020, 225, 126–134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Cai, G.C.; Ahmed, M.A.; Abdalla, M.; Carminati, A. Root hydraulic phenotypes impacting water uptake in drying soils. Plant Cell Environ. 2022, 45, 650–663. [Google Scholar] [CrossRef]
  41. Perez-Martin, A.; Michelazzo, C.; Torres-Ruiz, J.M.; Flexas, J.; Fernandez, J.E.; Sebastiani, L.; Diaz-Espejo, A. Regulation of photosynthesis and stomatal and mesophyll conductance under water stress and recovery in olive trees: Correlation with gene expression of carbonic anhydrase and aquaporins. J. Exp. Bot. 2014, 65, 3143–3156. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Valentini, R.; Epron, D.; De Angelis, P.; Matteucci, G.; Dreyer, E. In situ estimation of net CO2 assimilation, photosynthetic electron flow and photorespiration in Turkey oak (Q. cerris L.) leaves: Diurnal cycles under different levels of water supply. Plant Cell Environ. 1995, 18, 631–640. [Google Scholar] [CrossRef]
  43. Harley, P.C.; Loreto, F.; Di Marco, G.; Sharkey, T.D. Theoretical considerations when estimating the mesophyll conductance to CO2 flux by analysis of the response of photosynthesis to CO2. Plant Physiol. 1992, 98, 1429–1436. [Google Scholar] [CrossRef] [Green Version]
  44. Bernacchi, C.J.; Singsaas, E.L.; Pimentel, C.; Portis, A.R.; Long, S.P. Improved temperature response functions for models of Rubisco-limited photosynthesis. Plant Cell Environ. 2001, 24, 253–259. [Google Scholar] [CrossRef]
  45. Sack, L.; Scoffoni, C. Measurement of leaf hydraulic conductance and stomatal conductance and their responses to irradiance and dehydration using the Evaporative Flux Method (EFM). J. Vis. Exp. 2012, 70, e4179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Taylaran, R.D.; Adachi, S.; Ookawa, T.; Usuda, H.; Hirasawa, T. Hydraulic conductance as well as nitrogen accumulation plays a role in the higher rate of leaf photosynthesis of the most productive variety of rice in Japan. J. Exp. Bot. 2011, 62, 4067–4077. [Google Scholar] [CrossRef] [Green Version]
  47. Sack, L.; Streeter, C.M.; Holbrook, N.M. Hydraulic analysis of water flow through leaves of sugar maple and red oak. Plant Physiol. 2004, 134, 1824–1833. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Sellin, A.; Kupper, P. Temperature, light and leaf hydraulic conductance of little-leaf linden (Tilia cordata) in a mixed forest canopy. Tree Physiol. 2007, 27, 679–688. [Google Scholar] [CrossRef] [PubMed]
  49. Ye, M.; Zhang, Z.C.; Huang, G.J.; Xiong, Z.; Peng, S.B.; Li, Y. High leaf mass per area Oryza genotypes invest more leaf mass to cell wall and show a low mesophyll conductance. AoB Plants 2020, 12, plaa028. [Google Scholar] [CrossRef]
  50. Tosens, T.; Niinemets, Ü.; Westoby, M.; Wright, I.J. Anatomical basis of variation in mesophyll resistance in eastern australian sclerophylls: News of a long and winding path. J. Exp. Bot. 2012, 63, 5105–5119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Tosens, T.; Nishida, K.; Gago, J.; Coopman, R.E.; Cabrera, H.M.; Carriquí, M.; Laanisto, L.; Morales, L.; Nadal, M.; Rojas, R.; et al. The photosynthetic capacity in 35 ferns and fern allies: Mesophyll CO2 diffusion as a key trait. New Phytol. 2015, 209, 1576–1590. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Area-based leaf N content (Narea, (a)) and photosynthetic nitrogen use efficiency (PNUE, (b)) of plants under high N (HN) and zero N (N0) treatments at mid-tillering and booting stages. One-way analysis of variance (ANOVA) was used to assess the effect of N supply on parameters. ***, p < 0.001; ns, non-significant at p < 0.05 level.
Figure 1. Area-based leaf N content (Narea, (a)) and photosynthetic nitrogen use efficiency (PNUE, (b)) of plants under high N (HN) and zero N (N0) treatments at mid-tillering and booting stages. One-way analysis of variance (ANOVA) was used to assess the effect of N supply on parameters. ***, p < 0.001; ns, non-significant at p < 0.05 level.
Ijms 23 03885 g001
Figure 2. The relationships between net photosynthetic rate and (a) stomatal conductance and (b) mesophyll conductance across different temperatures under both N supplies and growth stages. A, net photosynthetic rate; gs, stomatal conductance; gm, mesophyll conductance. Filled cycles represent zero N (N0) treatment at mid-tillering stage, open cycles represent high N (HN) treatment at mid-tillering stage; filled triangles represent N0 treatment at booting stage, open triangles represent HN treatment at booting stage.
Figure 2. The relationships between net photosynthetic rate and (a) stomatal conductance and (b) mesophyll conductance across different temperatures under both N supplies and growth stages. A, net photosynthetic rate; gs, stomatal conductance; gm, mesophyll conductance. Filled cycles represent zero N (N0) treatment at mid-tillering stage, open cycles represent high N (HN) treatment at mid-tillering stage; filled triangles represent N0 treatment at booting stage, open triangles represent HN treatment at booting stage.
Ijms 23 03885 g002
Figure 3. The relationship between mesophyll conductance at 25 °C (gm,25) and cell-wall thickness (Tcw).
Figure 3. The relationship between mesophyll conductance at 25 °C (gm,25) and cell-wall thickness (Tcw).
Ijms 23 03885 g003
Table 1. Temperature responses of leaf gas-exchange parameters under different nitrogen (N) supplies at two growth stages.
Table 1. Temperature responses of leaf gas-exchange parameters under different nitrogen (N) supplies at two growth stages.
Growth StageN SupplyLeaf Temperature (°C)A (µmol m−2 s−1)gs (mol m−2 s−1)gm (mol m−2 s−1)
Mid-tilleringN01510.2 ± 1.1 c0.350 ± 0.124 a0.063 ± 0.006 c
2520.3 ± 1.0 b0.213 ± 0.020 b0.149 ± 0.027 b
3524.6 ± 2.6 a0.294 ± 0.068 ab0.207 ± 0.043 a
HN1514.0 ± 1.4 b0.452 ± 0.161 a0.096 ± 0.032 c
2531.7 ± 4.1 a0.372 ± 0.087 a0.300 ± 0.073 b
3535.1 ± 2.7 a0.430 ± 0.105 a0.388 ± 0.066 a
BootingN015---
2518.7 ± 2.50.178 ± 0.1740.132 ± 0.121
35---
HN158.6 ± 1.5 c0.188 ± 0.024 b0.084 ± 0.016 c
2518.9 ± 2.2 b0.174 ± 0.018 b0.127 ± 0.019 b
3524.0 ± 2.2 a0.257 ± 0.044 a0.188 ± 0.021 a
ANOVA
T ***ns***
N **ns*
S ******
T × N **ns***
T × S **ns***
N × S *******
Data are means ± SD of 4–7 replicates. *, p < 0.05; **, p < 0.01; ***, p < 0.001; ns, non-significant. T, temperature; S, growth stage; A, net photosynthetic rate; gs, stomatal conductance; gm, mesophyll conductance. The data followed by different letters in the same N and growth stages are significant at p < 0.05 level.
Table 2. The modelled parameters using Equation (8) for the temperature response of mesophyll conductance (gm).
Table 2. The modelled parameters using Equation (8) for the temperature response of mesophyll conductance (gm).
Growth StageN SupplycEa,gm (kJ mol−1)ΔS (kJ K−1 mol−1)Ed,gm (kJ mol−1)
Mid-tilleringN011.228.21.25437.4
HN15.338.41.25437.4
BootingN0----
HN12.129.91.25437.4
c, a scaling factor; Ea,gm, activation energy of gm; ΔS, an entropy term; Ed,gm, deactivation energy of gm.
Table 3. The modelled parameters using the two-component model for the temperature response of mesophyll conductance (gm) under two N supplies at different growth stages.
Table 3. The modelled parameters using the two-component model for the temperature response of mesophyll conductance (gm) under two N supplies at different growth stages.
Growth StageN Supplyl (µm)Pmem,25 (mm s−1)Ea,mem
(kJ mol−1)
gliq,25
(mol CO2 m−2 s−1)
gmem,25
(mol CO2 m−2 s−1)
gm,25
(mol CO2 m−2 s−1)
gliq,25’/gmem,25
Mid-tilleringN00.890.72260.40.07320.02390.01803.07
HN0.731.27281.00.08970.04210.02862.13
BootingN0-------
HN1.010.91066.40.06430.03010.02052.14
l, effective pathlength; Pmem,25, combined membrane permeability to CO2 at 25 °C; Ea,mem, activation energy of CO2-diffusion conductance through the membrane; gliq,25’, CO2-diffusion conductance through the liquid phase per surface areas of chloroplasts facing the intercellular air spaces at 25 °C; gmem,25’, CO2-diffusion conductance through the membrane per surface areas of chloroplasts facing the intercellular air spaces at 25 °C; gm,25’, mesophyll conductance at 25 °C per surface areas of chloroplasts facing the intercellular air spaces.
Table 4. Leaf anatomical traits under different N supplies at two growth stages.
Table 4. Leaf anatomical traits under different N supplies at two growth stages.
Growth StageN SupplySm
(µm2 µm−2)
Sc
(µm2 µm−2)
Tcw
(µm)
Mid-tilleringN010.95 ± 1.70 a6.92 ± 1.82 b0.226 ± 0.030 a
HN10.95 ± 1.51 a8.64 ± 1.90 a0.193 ± 0.024 b
BootingN011.37 ± 1.01 a5.25 ± 0.97 b0.245 ± 0.022 a
HN11.69 ± 0.98 a6.40 ± 1.76 a0.249 ± 0.042 a
ANOVA
N ns***
S ********
N × S nsns*
*, p < 0.05; **, p < 0.01; ***, p < 0.001; ns, non-significant. S, growth stage; Sm, the mesophyll cell surface area facing intercellular airspace per leaf area; Sc, the surface area of chloroplasts facing intercellular airspace per leaf area; Tcw, cell-wall thickness. The data followed by different letters in the same growth stages are significant at p < 0.05 level.
Table 5. Temperature responses of leaf hydraulic parameters under two N supplies at different growth stages.
Table 5. Temperature responses of leaf hydraulic parameters under two N supplies at different growth stages.
Growth StageN SupplyLeaf Temperature (°C)Ψleaf (MPa)Kleaf
(mmol m−2 s−1 MPa−1)
Mid-tilleringN015−0.293 ± 0.087 a8.57 ± 1.11 a
25−0.285 ± 0.090 a7.86 ± 1.80 a
35−0.423 ± 0.097 b7.05 ± 1.56 a
HN15−0.360 ± 0.099 a11.39 ± 2.53 a
25−0.430 ± 0.114 a10.71 ± 1.45 a
35−0.370 ± 0.106 a10.13 ± 2.58 a
BootingN015−0.155 ± 0.064 a7.48 ± 2.25 a
25−0.278 ± 0.068 b5.98 ± 1.78 a
35−0.560 ± 0.140 c6.03 ± 1.55 a
HN15−0.395 ± 0.137 a8.83 ± 1.16 a
25−0.310 ± 0.081 a9.20 ± 2.17 a
35−0.425 ± 0.044 a9.59 ± 1.38 a
ANOVA
T ***ns
N ns***
S ns**
T × N ***ns
T × S *ns
N × S nsns
Data are means ± SD of 4–7 replicates. *, p < 0.05; **, p < 0.01; ***, p < 0.001; ns, non-significant. S represents growth stage. Ψleaf, leaf water potential; Kleaf, the normalised leaf hydraulic conductance with water viscosity. The data followed by different letters in the same N and growth stages are significant at p < 0.05 level.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ye, M.; Zhang, Z.; Huang, G.; Li, Y. Leaf Photosynthesis and Its Temperature Response Are Different between Growth Stages and N Supplies in Rice Plants. Int. J. Mol. Sci. 2022, 23, 3885. https://doi.org/10.3390/ijms23073885

AMA Style

Ye M, Zhang Z, Huang G, Li Y. Leaf Photosynthesis and Its Temperature Response Are Different between Growth Stages and N Supplies in Rice Plants. International Journal of Molecular Sciences. 2022; 23(7):3885. https://doi.org/10.3390/ijms23073885

Chicago/Turabian Style

Ye, Miao, Zhengcan Zhang, Guanjun Huang, and Yong Li. 2022. "Leaf Photosynthesis and Its Temperature Response Are Different between Growth Stages and N Supplies in Rice Plants" International Journal of Molecular Sciences 23, no. 7: 3885. https://doi.org/10.3390/ijms23073885

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop