Abstract
Classical approaches to estimate mesophyll conductance ignore differences in resistance components for CO2 from intercellular air spaces (IAS) and CO2 from photorespiration (F) and respiration (Rd). Consequently, mesophyll conductance apparently becomes sensitive to (photo)respiration relative to net photosynthesis, (F + Rd)/A. This sensitivity depends on several hard-to-measure anatomical properties of mesophyll cells. We developed a method to estimate the parameter m (0 ≤ m ≤ 1) that lumps these anatomical properties, using gas exchange and chlorophyll fluorescence measurements where (F + Rd)/A ratios vary. This method was applied to tomato and rice leaves measured at five O2 levels. The estimated m was 0.3 for tomato but 0.0 for rice, suggesting that classical approaches implying m = 0 work well for rice. The mesophyll conductance taking the m factor into account still responded to irradiance, CO2, and O2 levels, similar to response patterns of stomatal conductance to these variables. Largely due to different m values, the fraction of (photo)respired CO2 being refixed within mesophyll cells was lower in tomato than in rice. But that was compensated for by the higher fraction via IAS, making the total re-fixation similar for both species. These results, agreeing with CO2 compensation point estimates, support our method of effectively analysing mesophyll resistance.
Introduction
Quantifying the CO2 diffusion inside leaves of C3 plants is important in both physiological and ecological contexts. Physiologists assess leaf photosynthetic efficiency and capacity, and both of them depend on how CO2 from the atmosphere travel to the chloroplast stroma and how much CO2 released by respiration and photorespiration [“(photo)respired CO2” hereafter] can be refixed by Rubisco (Busch et al. 2013; von Caemmerer 2013). Ecologists often project the impact of global land CO2 fertilization (Sun et al. 2014). The model of Farquhar, von Caemmerer and Berry (1980; “the FvCB model” hereafter), which is widely used as a component for this projection, requires the CO2 level at carboxylation sites of Rubisco (Cc) as its input. The drawdown of Cc, relative to the CO2 level in the ambient air (Ca), depends not only on stomatal conductance for CO2 transfer (gsc) but also on mesophyll conductance (gm), such that (von Caemmerer and Evans 1991):
where Ci is the intercellular air space (IAS) CO2 level and A is the net photosynthesis rate.
The FvCB model calculates A as the minimum of the Rubisco activity limited rate (Ac) and electron transport-limited rate (Aj) of photosynthesis, and Sharkey (1985) added a third limitation, accounting for the rate set by triose phosphate utilization (Ap) (see Supplementary Text S1). Equation (1) has been combined with the FvCB model to estimate gm from combined data of gas exchange and chlorophyll fluorescence measurements on photosystem II (PSII) electron transport efficiency Φ2 (Harley et al. 1992; Yin and Struik 2009). The most commonly used method to estimate gm is the ‘variable J method’ (Harley et al. 1992), derived from the Aj part of the FvCB model, using measurements that have to include photorespiratory conditions (Laisk et al. 2006; see Supplementary Text S2). Equation (1) has also been used to estimate gm from online carbon isotope discrimination measurements (e.g. Evans et al. 1994; Tazoe et al. 2011; Barbour et al. 2016a) or oxygen isotope techniques (Barbour et al. 2016b).
Equation (1), as the classical gm model, treats the (photo)respired CO2 in the same way as it treats the CO2 flux that comes from the IAS. Mesophyll resistance (the inverse of mesophyll conductance) consists of components imposed by IAS, cell wall, plasmalemma, cytosol, chloroplast envelope and stroma (Evans et al. 2009; Terashima et al. 2011). Unlike the CO2 from the IAS, the (photo)respired CO2, mainly coming from the mitochondria, does not need to cross the cell wall and plasmalemma, and thus experiences a different resistance. For this reason, Tholen et al. (2012) developed an Equation for the drawdown of Cc, relative to Ci:
where F and Rd are CO2 fluxes from photorespiration and respiration, respectively, rwp is the combined cell wall and plasma membrane resistance, and rch is the chloroplast envelope and stroma resistance (rch). Combining Eqs. (1) and (2) results in gm = 1/[rwp + rch + rch(F + Rd)/A]. Tholen et al. (2012) concluded that mesophyll conductance, as defined by Eq. (1), is influenced by the ratio of (photo)respired CO2 release to net CO2 uptake, (F + Rd)/A, thereby resulting in an apparent sensitivity of mesophyll conductance to [CO2] and [O2]. As this sensitivity does not imply a change in intrinsic diffusion properties, gm as defined by Eq. (1) is an apparent parameter. We shall call it the apparent mesophyll conductance (gm,app). In developing their model, Tholen et al. (2012) assumed a negligible IAS and cytosol resistance, but Eq. (2) still holds if the IAS resistance is lumped into rwp, and part of cytosol resistance is lumped into rwp, and the remaining part is lumped into rch (Berghuijs et al. 2015). If rwp and rch both represent physical resistances, the total mesophyll diffusion resistance (rm,dif) is rwp + rch, and the model of Tholen et al. can be rewritten as
where ω is the fraction of rch in rm,dif.
However, the relative position of mitochondria and chloroplasts is underrepresented in the model of Tholen et al. (2012). Considering six scenarios of the arrangement of these organelles, Yin and Struik (2017) derived the model: gm,app = 1/{rm,dif[1 + ω(1 − λk)(F + Rd)/A]}, where λ is the fraction of mitochondria located in the inner cytosol (i.e. the cytosol area between chloroplasts and vacuole), and k is a factor allowing an increase (k > 1), no change (k = 1), and a decrease (0 ≤ k < 1) in the fraction of inner (photo)respired CO2, caused by gaps when chloroplasts are not continuously aligned. The gaps largely depend on the anatomical parameter Sc/Sm, the ratio of chloroplast area to the mesophyll area exposed to IAS (Sage and Sage 2009). As (1 − λk) is between 0 and 1, the model predicts that the sensitivity of gm,app to (F + Rd)/A is lower than Tholen et al. (2012) initially stated (Yin and Struik 2017). The model of Tholen et al. applies to an extreme case, either where mitochondria are located exclusively in the outer cytosol between plasmalemma and chloroplasts (λ = 0) or where (photo)respired CO2 are completely mixed in cytosol if cytosol resistance is negligible and there are chloroplast gaps (k → 0). In another extreme case where mitochondria are located exclusively in the inner cytosol (λ = 1) and chloroplasts cover completely the cell periphery (k = 1), the model predicts no sensitivity of gm,app to (F + Rd)/A, and Eq. (1) would work well as gm,app becomes gm,dif (= 1/rm,dif). Equation (1) also works when rch is negligible compared to rwp (ω = 0) as if (photo)respired CO2 is released in the same organelle where RuBP carboxylation occurs. Either situation (λk = 1 or ω = 0) can be approximately represented by leaves where mitochondria lies only in the inner cytosol, intimately behind chloroplasts that form a continuum.
Most likely scenarios are somewhere between the two extremes defined by Eqs. (1) and (2), that is, 0 < λk < 1 and 0 < ω < 1. All these scenarios result in different fractions of re-assimilation of (photo)respired CO2 (Yin and Struik 2017), both within mesophyll cells and via IAS (see Supplementary Text S3). It would be useful if ω, λ and k can be measured. One way to derive ω is to use individual resistances that can be calculated from microscopic measurements on leaf anatomy (Evans et al. 1994; Peguero-Pino et al. 2012; Tosen et al. 2012a, b; Tomas et al. 2013; Berghuijs et al. 2015), despite uncertainties in the value of gas diffusion coefficients. Another possible method to estimate ω is to first estimate rwp from oxygen isotope techniques assuming that the outer limit of carbonic anhydrase activity represents the cytosol immediately adjacent to the cell wall (Barbour 2017). Parameter λ can be assessed using electron microscope images for mitochondria distribution (Hatakeyama and Ueno 2016). Most difficult is to measure k, which depends on Sc/Sm. However, whether a high Sc/Sm would make k > 1 or < 1 would depend on the λ value as well as on cytosol resistance, and such a complex relationship is hard to quantify with a simple resistance model. However, because ω, λ and k lump together co-defining the sensitivity of gm,app to (F + Rd)/A, the model of Yin and Struik (2017) can be rewritten to
where m = ω (1 − λk). Although Eq. (3) looks the same as Eq. (2a), their underlying intracellular fluxes for CO2 gradient and re-assimilation differ (see Supplementary Text S3). Equation (3) may be used for estimating m from noninvasive gas exchange measurements where (F + Rd)/A varies.
Many reports (e.g. Flexas et al. 2007a; Vrábl et al. 2009; Yin et al. 2009; Tazoe et al. 2011) showed that gm,app responds to changes in [CO2] or irradiance levels. gm,app was shown in tobacco to increase when [O2] was decreased from 21 to 1% (Tholen et al. 2012). All these responses can be described using a phenomenological equation (Yin et al. 2009). Tholen et al. (2012) explained the O2 response and the commonly observed decline of gm,app with decreasing CO2 below the ambient level (e.g. Flexas et al. 2007a; Vrábl et al. 2009; Yin et al. 2009), based on the earlier introduced sensitivity of gm,app to (F + Rd)/A, because both increasing O2 and decreasing Ci increase (F + Rd)/A. However, the sensitivity of gm,app to (F + Rd)/A cannot explain the observed response of gm,app to irradiances. Moreover, it is unknown whether gm,dif would be conserved across irradiance, CO2 and O2 levels.
In this study, we described a method that explores varying (F + Rd)/A ratios to analyse mesophyll resistance from combined gas exchange and chlorophyll fluorescence measurements. The varying (F + Rd)/A ratios were mainly created using five levels of O2, on two contrasting species tomato and rice. Using these data, we assessed (i) the value of the m factor and whether it differs between species, (ii) whether gm,dif responds to [CO2], irradiance and [O2], and (iii) how the re-assimilation of (photo)respired CO2 is affected by the m factor.
Materials and methods
Experiments and growth conditions
Seeds of tomato and rice were sown, and uniform seedlings were transplanted into pots 2 weeks after sowing, in glasshouse compartments. Pots were filled with soil, and after assessing initial soil nutrient contents, extra nutrients were applied (Table 1). Tomato plants were watered regularly, while rice plants were maintained submerged.
About 60% of the radiation incident on the glasshouse was transmitted to the plant level. During daytime supplemental light from 600 W HPS Hortilux Schréder lamps (Monster, NL) was automatically switched on when the incident solar flux dropped below a threshold and off when it exceeded a threshold outside glasshouse. These threshold levels were set different for tomato and rice (Table 1), to mimic growth environments of the two species.
Simultaneous gas exchange and chlorophyll fluorescence measurements
We used the Li-Cor-6400XT open gas exchange system with an integrated fluorescence head enclosing a 2-cm2 area (Li-Cor Inc, Lincoln-NE, USA). Young but fully expanded leaves of four replicated plants from staggered sowings were measured for incident irradiance (Iinc) and Ca response curves in each species (Table 1).
Curves were measured at five O2 concentrations (Table 1). Additional light response curves were obtained at 1000 µmol mol−1Ca and 2% O2 to establish nearly nonphotorespiratory conditions for calibration (see later). Gas from a cylinder containing a mixture of O2 and N2 was humidified and supplied via an overflow tube to the air inlet of the Li-Cor where CO2 was blended with the gas, and the IRGA was adjusted for O2 composition of the gas mixture according to the manufacturer’s instructions. Based on pre-test measurements, we used 7–8 min for each step of an A − Iinc curve, and 3–4 min for each step of an A − Ci curve, to reach a steady state. All CO2 exchange data were corrected for CO2 leakage into and out of the leaf cuvette, using measurements on boiled leaves (Flexas et al. 2007b), and then Ci was re-calculated.
When A reached steady state at each light or CO2 step, steady-state fluorescence (\({F}_{\mathrm{s}}\)) was recorded. Maximum fluorescence (\({F}_{\mathrm{m}}^{^{\prime}}\)) was measured using a 0.8 s light pulse of > 8000 µmol m−2 s−1, or the multiphase flash with each phase of 300 ms and ramp depth of 40% (Loriaux et al. 2013). The PSII operating efficiency (\(\Delta F/{F}_{\mathrm{m}}^{^{\prime}}\)) was set as \(({F}_{\mathrm{m}}^{^{\prime}}-{F}_{\mathrm{s}})/{F}_{\mathrm{m}}^{^{\prime}}\) (Genty et al. 1989).
Calibration and pre-determination of Rd and Rubisco parameters
Setting that \({\Phi }_{2}=\Delta F/{F}_{\mathrm{m}}^{^{\prime}}\), Rd was estimated as the negative intercept of a linear regression of A against (IincΦ2/4) using data of A − Iinc curves within the electron transport-limited range for the nonphotorespiratory condition (Yin et al. 2009, 2011). The slope of the regression yields a calibration factor (s), which lumps (1) absorptance by leaf photosynthetic pigments, (2) the factor for excitation partitioning to PSII, (3) basal forms of alternative electron transport, (4) any difference between real efficiency of PSII electron transport (Φ2) and \(\Delta F/{F}_{\mathrm{m}}^{^{\prime}}\), and (5) possibly difference in chloroplast populations sampled by gas exchange and by chlorophyll fluorescence (van der Putten et al. 2018). The electron transport rate J can then be obtained as \(J=s{I}_{\mathrm{i}\mathrm{n}\mathrm{c}}(\Delta F/{F}_{\mathrm{m}}^{^{\prime}})\) (Yin et al. 2009). Like other calibration methods, this procedure assumes that the calibration factor is the same for photorespiratory and nonphotorespiratory conditions, for which photosynthetic rates differ by a factor of (Cc − Γ*)/(Cc + 2 Γ*) (see Eqs. S1.1 and S1.3 in Supplementary Text S1; but with cautions from recent literature, Busch et al. 2018; Tcherkez and Limami 2019).
The parameter Γ* was calculated as 0.5O2/Sc/o, where Sc/o is the relative CO2/O2 specificity of Rubisco (von Caemmerer et al. 1994). Values from in vitro measurements of Cousins et al. (2010) on Sc/o (= 3.022 mbar μbar−1) and Michaelis–Menten coefficients of Rubisco for CO2 (KmC = 291 μbar) and for O2 (KmO = 194 mbar) were taken, assuming that Rubisco kinetic constants are conserved among C3 species. This assumption was checked by in vivo estimates of Sc/o from the lower parts of A − Ci curves of five O2 levels (see “Results”).
Model method
After the above parameters were quantified, we first checked whether gm,dif was variable based on the combined data of gas exchange and chlorophyll fluorescence. Using measured A, Ci and a tentative value for m across its range (0 ≤ m ≤ 1), gm,dif was calculated as
where F and Cc can be solved from the Aj equation of the FvCB model, see Eq. (S1.6) in Supplementary Text S1 and Eq. (S2.1) in Supplementary Text S2, respectively. Equation (4) was derived by Yin and Struik (2017, see their Eq. 19), in analogy to the variable J method of Harley et al. (1992; also see Eq. S2.2 in Supplementary Text S2).
The obtained gm,dif responded to a change in both Ci and irradiance (see “Results”). Explaining these responses would need a separate study; to estimate m, here we adopted the generic phenomenological equation of Yin et al. (2009) to describe this response:
where gmo,dif and δ are parameters. If δ = 0, Eq. (5) becomes a constant gm,dif mode (= gmo,dif). Any nonzero δ would predict a variable gm,dif in response to CO2, O2 and irradiance levels, and if gmo,dif = 0, parameter δ, as discussed later, represents the carboxylation: mesophyll resistance ratio. Equation (5) was combined with the FvCB and other equations to solve for A (Supplementary Text S1, where reasons for using Eq. 5 are also explained). This results in an equation expressing A as a function of Ci and other variables:
where
where x1 = Vcmax (maximum carboxylation activity of Rubisco) and x2 = KmC(1 + O2/KmO) for the Ac-limited conditions; x1 = J/4 and x2 = 2Γ* for the Aj-limited conditions, and for the Ap-limited conditions: x1 = 3Tp (where Tp is the rate of triose phosphate export from the chloroplast) and x2 = − (1 + 3α) Γ* (where α is the fraction of glycolate carbon not returned to the chloroplast).
We found that the \(\sqrt{{b}^{2}-4ac}\) term of Eq. (6) should always take the – sign for either Ac- or Aj-limited rate, but the solution for Ap is mathematically complicated if α > 0 (see Supplementary Text S4). Our data showed that A often declined with increasing Ci within high Ci ranges (see “Results”), suggesting the limitation by triose phosphate utilization with α > 0 (Harley and Sharkey 1991). We conducted sensitivity analyses to choose a value of α although metabolic flux data (Abadie et al. 2018) suggest that its value might be small. We then used Eq. (6) to estimate four parameters: m (0 ≤ m ≤ 1), δ, Vcmax and Tp, by a nonlinear fitting to all data of A − Ci and A − Iinc curves of the five O2 levels (gmo,dif was set to zero, see “Results”). For that, J, as defined earlier as\(s{I}_{\mathrm{i}\mathrm{n}\mathrm{c}}(\Delta F/{F}_{\mathrm{m}}^{^{\prime}})\), were used as input. Our method assumed that Rd does not vary with [O2], and was based on the expectation that neither Vcmax nor Tp varies with [O2], as confirmed experimentally for Vcmax (von Caemmerer et al. 1994). The fitting minimizes the sum of squared differences between estimated and measured A values, using the GAUSS method in PROC NLIN (SAS Institute, NC, USA). SAS scripts can be obtained upon request.
Once A was calculated from Eq. (6), Cc could be solved from Eq. (S2.1) in Supplementary Text S2. Then, gm,dif was re-calculated from Eq. (4) using the estimated m and measured A and Ci, where x1 and x2 terms were chosen according to whether the modelled A was Ac-, Aj- or Ap-limited. This showed gm,dif in response to CO2, irradiance, and O2 levels.
With gm,dif and other parameters, we calculated the fraction of (photo)respired CO2 being refixed (frefix), the fraction of (photo)respired CO2 being refixed within the mesophyll cells (frefix,cell), and the fraction of (photo)respired CO2 being refixed via IAS (frefix,ias), using Eqs. (S3.4), (S3.5) and (S3.6), respectively, in Supplementary Text S3. In these equations, rsc is the stomatal resistance to CO2 diffusion (being 1.6 times measured stomatal resistance to water vapour), and rcx is the carboxylation resistance (which is \(({C}_{\mathrm{c}}+{x}_{2})/{x}_{1}\), von Caemmerer 2000). As discussed in Supplementary Text S3, these calculations need ω and λk as inputs. The estimate for m was 0.3 for tomato and 0.0 for rice (see “Results”). For tomato, we measured ω (0.65) for leaves of the same age in the same cultivar “Growdena” (see Berghuijs et al. 2015) for calculating λk, from m = ω (1 − λk). For rice, λk was set to 1.0 to agree with the estimate that m = 0. In such a case, ω is not needed as Eqs. (S3.4) and (S3.5) become simplified as Eqs. (S3.3) and (S3.9) in Supplementary Text S3, respectively.
Results
Use of the five O2 levels generated diverse shapes of photosynthetic responses to irradiance and CO2 levels (Fig. 1). Our model approach, combined with data for A (Fig. 1) and for \(\Delta F/{F}_{\mathrm{m}}^{^{\prime}}\) (Fig. S1), yielded an estimation of a set of parameters as described below.
Estimated Rd and s
Data of A − Iinc curves within the range of Iinc ≤ 200 µmol m−2 s−1 showed that the relationship between A and (IincΦ2/4) was linear for the conditions with a gas mixture of 2% O2 with 1000 μmol mol−1Ca (Fig. 2), where Φ2 was set to be \(\Delta F/{F}_{\mathrm{m}}^{^{\prime}}\). The value of Rd estimated from this linear relationship was 1.2 (standard error or s.e. 0.1) µmol m−2 s−1 for tomato and 1.1 (s.e. 0.1) µmol m−2 s−1 for rice. The slope of the A − (IincΦ2/4) linearity (i.e. calibration factor s) was 0.4570 (s.e. 0.0076) for tomato and 0.5488 (s.e. 0.0076) for rice. Values of s were also re-estimated, together with other parameters, in fitting Eq. (6) to all data; but the re-estimated s remained the same, suggesting that we reached a nonphotorespiratory condition using the gas mixture.
The first few data points of the A − Ci curves were linear, and gross leaf photosynthesis values A + Rd were plotted versus Ci within this linear range. The intercept of this line with the Ci-axis gives the estimate of the Ci-based CO2 compensation point, commonly noted as Ci*. The value of Ci* increased linearly with increasing O2 levels (Fig. 3). Half of the reciprocal of this linear slope gives an in vivo estimate of Sc/o, which was 2.71 mbar μbar−1 for tomato and 3.13 mbar μbar−1 for rice. Using the method of Yin et al. (2009) gave similar in vivo estimates of Sc/o (results not shown). These values are close to 3.02 mbar μbar−1 measured in vitro for wheat by Cousins et al. (2010), confirming that Sc/o is conserved among C3 species. We will use 3.02 mbar μbar−1 for further analysis (but see sensitivity analysis later).
Dependence of gm,dif on CO2 and irradiance level
Equation (4) assuming an electron transport limitation, was applied to check the pattern of gm,dif across a range of Iinc and Ci levels, by setting m either to 0 (equivalent to the variable J method of Harley et al. (1992) for gm,app) or to a value between 0 and 1. A similar response was obtained for various O2 levels, except for 2% O2. At that oxygen concentration, Eq. (4), like the variable J method, cannot be reliably applied due to insufficient photorespiration (see Supplementary Text S2). Although the obtained gm,dif sometimes had unrealistic values largely due to unrealistic values of Cc (as often occurs when using the variable J method, see Yin and Struik 2009), an overall trend of gm,dif in response to Iinc and to Ci was obtained. An example of the response is shown in Fig. 4 for the case of 10% O2 level for tomato. gm,dif increased monotonically with increasing Iinc (Fig. 4a), and decreased gradually with an increase in Ci (Fig. 4b). Changing m did not change the response pattern, but only the absolute value of gm,dif, and a nonzero m resulted in higher gm,dif than the value obtained from setting m = 0 (Fig. 4).
Estimates of parameters δ, m, Vcmax and Tp
Equation (6) for describing A was applied to estimate gmo,dif, δ and m, using data of both A − Iinc and A − Ci curves. The obtained gmo,dif did not differ significantly from zero (p > 0.05), which is supported by the result that the calculated gm,dif by Eq. (4) at low Iinc was close to zero (Fig. 4a). Also, model fit became worse if δ was fixed to zero than if gmo,dif was fixed to zero, supporting the variable gm,dif mode. Sensitivity analysis with respect to α suggested that a change within its relevant range had no impact on the estimates of parameters other than Tp (see below). We set gmo,dif to zero, and α to 0.3 (Busch and Sage 2017), in the subsequent analysis.
Equation (6) describes well both A − Iinc and A − Ci curves (Fig. 1), with an overall R2 being > 0.99 for either species (Table 2). Most of the data points (> 80%) were Aj-limited, indicating that chlorophyll fluorescence signals generally echoed gas exchange data since we calculated J from chlorophyll fluorescence measurements as \(s{I}_{\mathrm{i}\mathrm{n}\mathrm{c}}(\Delta F/{F}_{\mathrm{m}}^{^{\prime}})\). Only a few points at low Ci of A − Ci curves or at high Iinc of A − Iinc curves were Ac-limited, and a few points at high Ci of A − Ci curves under low O2 conditions were Ap-limited. The estimated m was ca 0.3 for tomato but was 0.0 for rice (Table 2). The estimated δ was also higher for tomato (1.4) than for rice (1.0) (Table 2). Other parameter values were similar for the two species: 113.7 and 111.0 µmol m−2 s−1 for Vcmax, and 8.3 and 7.8 µmol m−2 s−1 for Tp, for tomato and rice, respectively.
Sensitivity analysis
Given that any uncertainty in estimated s and Rd and in other parameters (Sc/o, KmC, KmO and α) may have an impact on the major estimated parameters (m and δ in this study), we carried out sensitivity analyses. The estimation of δ and m was very sensitive to s and Sc/o, and less sensitive to Rd (Fig. S2), but virtually insensitive to KmC, KmO and α (results not shown). Both δ and m decreased monotonically with increasing s (Fig. S2a). The estimate of δ decreased with increasing Sc/o, whereas that of m changed in an opposite direction (Fig. S2b). The obtained response of δ (the parameter in Eq. (5) on mesophyll conductance) to both Sc/o and s is expected in the same way as gm,app responds to these parameters (Harley et al. 1992). The opposite response of m to Sc/o and s is probably because photorespiration, i.e. the F term in Eq. (3), which is relevant to determining m, has an opposite response to Sc/o and s. As Rd has the same effect as the F term has (see Eq. 3), the estimated m decreased with increasing Rd, whereas δ changed in an opposite direction (Fig. S2c). As expected, any sensitivity to KmC and KmO occurred with the estimated Vcmax, whereas a sensitivity to α occurred with Tp (results not shown).
Calculated fractions for re-assimilation of (photo)respired CO2
The calculated fractions of (photo)respired CO2 being refixed, using Eqs. (S3.4–S3.6) in Supplementary Text S3, are shown in Fig. 5, using the result at 21% O2 as the example. The trends were similar for O2 levels above 2%. Except for very low Iinc or Ci levels, the refixed fractions were quite consistent over a wide range of conditions. frefix,cell was lower in tomato (0.25) than in rice (0.49) (Fig. 5), largely due to the fact that the estimated m was 0.3 for tomato but 0.0 for rice (Table 2). In contrast, frefix,ias was higher in tomato than in rice. As a result, the total re-fixation fraction frefix was comparable for the two species, i.e. up to ca 0.6.
Responses of stomatal and mesophyll conductance to O2
Except for a few cases, gsc generally decreased with increasing [O2], and was lower in tomato than in rice (Fig. 6). The calculated value of gm,dif also decreased with increasing [O2], except for very high CO2 conditions which lowered gm,dif to the extent that the O2 response of gm,dif was no longer significant (Fig. 6d,j). gm,dif was higher in tomato than in rice.
Discussion
Analysing mesophyll resistance
Compared with the CO2 flux coming from IAS, (photo)respired CO2 experiences different resistances. This suggests the need to dissect rm,dif into sub-components. Anatomical measurements can partition mesophyll resistance into individual sub-components (Peguero-Pina et al. 2012; Tosens et al. 2012a, b; Tomas et al. 2013; Carriquí et al. 2019). The calculation of these sub-components relies on many assumed diffusion or permeability coefficients that are uncertain (Berghuijs et al. 2015). Furthermore, this approach does not quantify the effect of the arrangement of mitochondria and chloroplasts on the intracellular CO2 diffusion.
In line with anatomical measurements, Eq. (2) dissects rm,dif into two sub-components rwp and rch (Tholen et al. 2012). Oxygen isotope techniques may estimate rwp based on certain assumptions (Barbour 2017), but so far have been explored to separate rwp and rch within the framework of the classical gm model, Eq. (1) (Barbour et al. 2016b). Equations (1) and (2) both underrepresent the intracellular arrangements of organelles. In contrast, the model of Yin and Struik (2017), Eq. (3), has a factor lumping (i) the rch:rm,dif ratio (ω), (ii) the fraction of (photo)respired CO2 that are released in the inner cytosol (λ), and (iii) k, the factor for the change in λ as a result of the chloroplast gaps. The factor k is particularly hard to assess. Since ω, λ and k lump as such that m = ω (1 − λk), Eq. (3) provides an approach by exploring nondestructive gas exchange and chlorophyll fluorescence measurements under different levels of O2 that created large variations in photorespiration. Instead of estimating individual resistances, our nonlinear fitting approach estimates rm,dif as a whole, as well as the m factor. Common nonlinear procedures, typically by fitting A − Ci curves, estimate four or even more parameters of the FvCB model, such as Vcmax, J, Tp, and gm (e.g. Sharkey et al. 2007). In our method, J was measured from chlorophyll fluorescence. Despite a wide range of O2 levels exploited, we restricted the number of estimated parameters to four from Eq. (6)-based nonlinear fitting. All four were reliably estimated for both species with small standard errors (Table 2).
Our approach is still a simplified representation of complex diffusion pathways. Some respiratory flux may originate in the chloroplasts, in cytosol, and in the heterotrophic tissues such as epidermis, vasculature, and bundle sheath (Tcherkez et al. 2017). These components of Rd could be incorporated as additional terms into the Ci–Cc gradient equation, Eq. (S1.7) in Supplementary Text S1. However, they are ignored here as fractions of these components in Rd are generally unknown. There may also be some activity of phosphoenolpyruvate carboxylase (Vpepc) in cytosol (Douthe et al. 2012; Abadie and Tcherkez 2019), which would counteract the effect of (F + Rd) on gm,app. But our procedure of estimating Rd may have accounted for this, i.e. the estimated Rd represents the net rate of true Rd minus Vpepc. Tholen et al. (2012) showed that small amounts of Vpepc have little impact on gm,app.
Variation of gm,dif with CO2, irradiance and O2 levels
Reports using chlorophyll fluorescence data consistently showed that gm,app initially increases and then decreases with increasing Ci and increases monotonically with increasing Iinc (e.g. Flexas et al. 2007a; Yin et al. 2009). Similar results for gm,app in response to Ci (Vrábl et al. 2009; Tazoe et al. 2011) and to Iinc (Douthe et al. 2012) were sometimes reported, using the carbon isotope discrimination method. No change in anatomical arrangements was observed that could explain the variable gm,app (Carriquí et al. 2019). Gu and Sun (2014) showed that the reported response of gm,app to a change in CO2 or in Iinc may be due to the artefact of errors in experimental measurements. Although resolving experimental uncertainties is urgently needed, consistent variations of gm,app cannot be ascribed only to experimental errors because responses due to random errors would be irregular and inconsistent among various reports. Théroux-Rancourt and Gilbert (2017) demonstrated that changing patterns of light penetration within the leaf 3D-structure leads to different contributions of each cell layer to bulk-leaf mesophyll conductance, resulting in an apparent response of the bulk-leaf gm,app to light intensity. However, their theory cannot explain the response of gm,app to Ci.
Most results using the variable J method of Harley et al. (1992) showed that within a low Ci range, gm,app typically decreases with decreasing Ci (e.g. Flexas et al. 2007a; Vrábl et al. 2009; Yin et al. 2009; Fig. 4b). Tholen et al. (2012) also noted a decrease of gm,app with increasing O2 level. Tholen et al. (2012) explained these responses to Ci and O2 as a consequence of the fact that gm,app as an apparent parameter decreases with an increase in the (F + Rd)/A ratio. If this is the only explanation of variable gm,app, one would expect that gm,dif would be independent of Ci because Eq. (4) for gm,dif already accounts for the (F + Rd)/A ratio. However, gm,dif still declined with decreasing Ci within its low range, albeit to a lesser extent (Fig. 4b). In fact, within the low Ci range, A is limited by Rubisco activity; so, as noted by Yin et al. (2009), using the variable J method or Eq. (4) assuming an electron transport limitation results in an artefactual decline of estimated gm because the occurrence of additional alternative electron transport is wrongly attributed to the mesophyll diffusional limitation. This assertion was supported by the result that the decrease of gm,dif with decreasing Ci within the low Ci range was no longer obtained once gm,dif was calculated from the full FvCB model (Fig. 4b). This suggests that the decrease of gm,app with decreasing Ci is explained more by the occurrence of alternative electron transport than by the theory of Tholen et al. Anyway, the theory does not explain the decreases of gm,app with increasing Ci within its high range, or with decreasing Iinc, or with lowering temperature as reported previously (Bernacchi et al. 2002; Warren and Dreyer 2006; Yamori et al. 2006; Scafaro et al. 2011; Evans and von Caemmerer 2013).
We found that gm,dif increased with Iinc (Fig. 4a). This increase continued within the high Iinc range, where some additional alternative electron transport is also expected, suggesting that the increase of gm,dif with Iinc overrode any artefactual decline caused by alternative electron fluxes. Also, gm,dif decreased with increasing O2 (Fig. 6), in the same direction as the O2 response of gm,app reported by Tholen et al. (2012). Literature on O2 responses of diffusional conductance is scarce (Farquhar and Wong 1984; Buckley et al. 2003). Our data showed that both gsc and gm,dif generally declined with increasing O2. So, gm,dif is variable, in response to Ci, Iinc, and O2, in a similar pattern as gsc responds to these variables (e.g. Morison and Gifford 1983; Farquhar and Wong 1984; Buckley et al. 2003).
The identified variable gm,dif was based on the assumption that m (= ω (1-λk)) is constant, independent of short-term changes (within 3–8 min) in irradiance or [CO2]. This is supported by Carriquí et al. (2019), who reported that anatomical parameters determining ω and k hardly vary with short-term changes in irradiance or [CO2]. Chloroplasts and mitochondria in some plants may move under varying light, but they always colocalize (Islam et al. 2009), suggesting that λ also hardly varies. We are unable to find evidences supporting quantitative changes that m or its components must have to obtain invariable gm,dif with irradiance, [CO2] and [O2].
gm,dif defined here is still a bulk-leaf trait. Like bulk-leaf gm,app, it may not represent intrinsic transport properties. Also, our result on the variable gm,dif is subject to experimental confirmation by other methods. If proven true, future studies are needed to examine if the variable gm,dif can emerge from fluxes and concentrations across the real 3D-structure of leaves, as well as in relation to membrane permeability and other properties. Here we only describe the response from bulk-leaf equations themselves. gm,dif can be formulated from Eq. (4) and A = V − F − Rd (where V is the carboxylation rate) as
When Iinc increases, only the numerator increases significantly; so Eq. (7) predicts that gm,dif increases with increasing Iinc. If the CO2 gradient from Ci and Cc is regulated such that the Cc:Ci ratio is roughly constant for a given O2 level (results not shown), Eq. (7) also predicts that gm,dif will decrease monotonically with Ci because according to the FvCB model, the V increment per Ci increment decreases with increasing Ci. Finally, the F term increases when O2 increases; as a result, Eq. (7) predicts that gm,dif decreased with increasing O2 (Fig. 6).
Interpretation of the model and estimated parameter values
Our method is based on Eq. (3), the equation summarized by Yin and Struik (2017) from considering six possible scenarios for the intracellular organelle arrangement. Recently, Ubierna et al. (2019) came up with the same model but formulated gm,app in a Michaelis–Menten-like equation, i.e. gm,app = A·gm,dif/[A + m(F + Rd)] (see their Eq. 15; note that gm,app was written as gm in their notations). The maximum value of gm,app is gm,dif, while the Michaelis–Menten constant “Km” is m(F + Rd). For the case of tomato where m = 0.3 and Rd = 1.2, the “Km” occurs at A ≈ 2.0 μmol m−2 s−1 for the ambient O2 condition. This suggests that gm,app and gm,dif only differ significantly when A is low, which our results (Fig. 4) confirmed.
In view of the variation of gm,dif shown in Fig. 4, we adopted Eq. (5), which accommodates either constant or variable gm,dif in relation to Ci, Iinc and O2 levels. Although the equation is phenomenological and has an a priori assumption that gm,dif grows with relative carboxylation and the estimates of its parameters are expectedly sensitive to the pre-input values of s, Sc/o and Rd (Fig. S2), the model generated useful insights.
Our results supported no constant gm,dif, but a variable gm,dif with parameter δ being 1.0 for rice and 1.4 for tomato (Table 2). Equation (5) with gmo,dif = 0 for our variable gm,dif mode can be rewritten to \({r}_{\mathrm{m},\mathrm{d}\mathrm{i}\mathrm{f}}=({C}_{\mathrm{c}}-{\Gamma }_{*})/[\delta \left(A+{R}_{\mathrm{d}}\right)]\). As \(\left(A+{R}_{\mathrm{d}}\right)\) can be calculated from the FvCB model as \(({C}_{\mathrm{c}}-{\Gamma }_{*}){x}_{1}/\left({C}_{\mathrm{c}}+{x}_{2}\right)\), the above equation becomes \({r}_{\mathrm{m},\mathrm{d}\mathrm{i}\mathrm{f}}=({C}_{\mathrm{c}}+{x}_{2})/\left(\delta {x}_{1}\right)\). As \(({C}_{\mathrm{c}}+{x}_{2})/{x}_{1}\) is defined as carboxylation resistance rcx (von Caemmerer 2000), it follows that
Thus, parameter δ of Eq. (5) has a meaning, representing the carboxylation: mesophyll resistance ratio. Our estimates for δ (Table 2) suggest that \({{r}_{\mathrm{c}\mathrm{x}}\, \mathrm{a}\mathrm{n}\mathrm{d}\, r}_{\mathrm{m},\mathrm{d}\mathrm{i}\mathrm{f}}\) had similar values in rice leaves, whereas rcx was ca 40% higher than rm,dif in tomato leaves.
Our estimate of the factor m was ca 0.3 for tomato and 0.0 for rice (Table 2). Thus, using Eq. (1), which is the special case of the generalized model when m = 0, actually suits for rice leaves but does not work for tomato leaves when (F + Rd)/A is high. As stated in Introduction, the classical model works well if mitochondria are located exclusively in the inner cytosol (λ = 1) and chloroplasts cover fully the mesophyll periphery that k = 1. Sage and Sage (2009) and Busch et al. (2013) showed that compared with other species, in rice leaves, there are stromules that effectively extend chloroplast coverage of the cell periphery and mitochondria locate in the cell interior and are intimately associated with chloroplasts/stromules. These features engender such a structure as if (photo)respired CO2 is released in the same compartment where RuBP carboxylation occurs. This is the case when Eq. (1) works well. Therefore, our results with curve-fitting to gas exchange data actually agree with anatomical differences between species.
Such differences are also shown in the fractions of re-fixation of (photo)respired CO2 calculated from resistance components (Fig. 5). With the distinct anatomical feature of rice leaves, (photo)respired CO2, if to exit mesophyll cells, will have to travel via the stroma, thereby maximizing the re-fixation of (photo)respired CO2 within the cell. Therefore, rice had higher values of frefix,cell than tomato (Fig. 5). For a given set of resistance values, the organelle arrangements as in rice leaves that make the highest frefix,cell can result in low frefix,ias (see Supplementary Text S3). Moreover, in line with the observation of Ouyang et al. (2017) on rice ‘IR64′, the cultivar we used, rice had high stomatal conductance, compared with tomato (Fig. 6). A low gsc would make (photo)respired CO2 more difficult to exit into the atmosphere via IAS. This also contributed to higher values of frefix,ias in tomato than in rice (Fig. 5). As a result, the two species had similar values (up to 60%) of the total re-fixation, frefix. The calculated frefix,ias and frefix varied with Iinc or Ci levels (Fig. 5), because resistance components rsc and rcx varied with these variables. The calculated frefix,cell was more constant (Fig. 5). Substituting Eq. (8) into Eq. (S3.5) in Supplementary Text S3 gives
As all terms are constant, Eq. (9) describes why frefix,cell stayed invariant. Using isotope mass spectrometry and gas exchange measurements, Busch et al. (2013) determined frefix,cell, frefix,ias and frefix, being 0.29, 0.22, and 0.51 for rice under ambient CO2 and high light conditions. Our estimates for rice somewhat differed from their values for comparable conditions (Fig. 5d).
The difference in the value of factor m between the species also has implications on values of Ci* and the relationship between Ci* and Γ*. Ci* was lower in rice than in tomato at a given O2 level (Fig. 3), and Ci* at the lowest O2 in rice was even negative (Fig. 3b). A negative Ci* could be due to measurement noises, uncertainties in assuming constant Rd, and the influence of varying amounts of Vpepc (see earlier discussions). However, for the case where m = 0, it can be seen from Eq. (1) that Ci* = Γ* – Rd/gm (von Caemmerer et al. 1994); so, Ci* is always lower than Γ*. This agrees with our linear relation for rice in Fig. 3b (where the term 0.16O2 can be considered as Γ*, given that Γ* = 0.5O2/Sc/o). Setting the intercept of this relation equal to –Rd/gm and knowing that Rd = 1.064 µmol m−2 s−1 (Fig. 2b), gm for rice can be solved as ca 0.1 mol m−2 s−1 bar−1, comparable with its value calculated in the other way for low CO2 conditions (Fig. 6l). So, a negative Ci* for low O2 conditions could actually represent biological realities, i.e. high intracellular re-fixation of both respired CO2 and photorespired CO2 sufficed to (over)compensate for photorespiratory losses. In contrast, for cases where m ≥ 0, the relation between Ci* and Γ* can be formulated from Eq. (S1.7) in Supplementary Text S1 as
Equation (10) means that Ci* is no longer necessarily lower than Γ* (see also Tholen et al. 2012), depending on relative values of (1–m)Rd versus mF. Our result in Fig. 3a suggests that Ci* is 2.134 μbar higher than Γ* for tomato. Thus, different m values estimated by curve-fitting for the two species are supported by the Ci* vs Γ* relationships in Fig. 3, suggesting that our approach is internally consistent.
References
Abadie C, Tcherkez G (2019) In vivo phosphoenolpyruvate carboxylase activity is controlled by CO2 and O2 mole fractions and represents a major flux at high photorespiration rates. New Phytol 221:1843–1852
Abadie C, Bathellier C, Tcherkez G (2018) Carbon allocation to major metabolites in illuminated leaves is not just proportional to photosynthesis when gaseous conditions (CO2 and O2) vary. New Phytol 218:94–106
Barbour MM (2017) Understanding regulation of leaf internal carbon and water transport using online stable isotope techniques. New Phytol 213:85–88
Barbour M, Bachmann S, Bansal U, Bariana H, Sharp P (2016a) Genetic control of mesophyll conductance in common wheat. New Phytol 209:461–465
Barbour MM, Evans JR, Simonin KA, von Caemmerer S (2016b) Online CO2 and H2O oxygen isotope fractionation allows estimation of mesophyll conductance in C4 plants, and reveals that mesophyll conductance decreases as leaves age in both C4 and C3 plants. New Phytol 210:875–889
Berghuijs HNC, Yin X, Ho QT, van der Putten PEL, Verboven P, Retta MA, Nicolaï BM, Struik PC (2015) Modelling the relationship between CO2 assimilation and leaf anatomical properties in tomato leaves. Plant Sci 238:297–311
Bernacchi CJ, Portis AR, Nakano H, von Caemmerer S, Long SP (2002) Temperature response of mesophyll conductance. Implication for the determination of Rubisco enzyme kinetics and for limitations to photosynthesis in vivo. Plant Physiol 130:1992–1998
Buckley TN, Mott KA, Farquhar GD (2003) A hydromechanical and biochemical model of stomatal conductance. Plant Cell Environ 26:1767–1785
Busch FA, Sage RF (2017) The sensitivity of photosynthesis to O2 and CO2 concentration identifies strong Rubisco control above the thermal optimum. New Phytol 213:1036–1051
Busch FA, Sage TL, Cousins AB, Sage RF (2013) C3 plants enhance rates of photosynthesis by reassimilating photorespired and respired CO2. Plant Cell Environ 36:200–212
Busch FA, Sage RF, Farquhar GD (2018) Plants increases CO2 uptake by assimilating nitrogen via the photorespiratory pathway. Nat Plants 4:46–54
Carriquí M, Douthe C, Molins A, Flexas J (2019) Leaf anatomy does not explain apparent short-term responses of mesophyll conductance to light and CO2 in tobacco. Physiol Plant 165:604–618
Cousins AB, Ghannoum O, von Caemmerer S, Badger MR (2010) Simultaneous determination of Rubisco carboxylase and oxygenase kinetic parameters in Triticum aestivum and Zea mays using membrane inlet mass spectrometry. Plant Cell Environ 33:444–452
Douthe C, Dreyer E, Brendel O, Warren CR (2012) Is mesophyll conductance to CO2 in leaves of three Eucalyptus species sensitive to short-term changes of irradiance under ambient as well as low O2? Funct Plant Biol 39:435–448
Evans JR, von Caemmerer S (2013) Temperature response of carbon isotope discrimination and mesophyll conductance in tobacco. Plant Cell Environ 36:745–756
Evans JR, von Caemmerer S, Setchell BA, Hudson GS (1994) The relationship between CO2 transfer conductance and leaf anatomy in transgenic tobacco with a reduced content of Rubisco. Aust J Plant Physiol 21:475–495
Evans JR, Kaldenhoff R, Genty B, Terashima I (2009) Resistances along the CO2 diffusion pathway inside leaves. J Exp Bot 60:2235–2248
Farquhar GD, Wong SC (1984) An empirical model of stomatal conductance. Aust J Plant Physiol 11:191–210
Farquhar GD, von Caemmerer S, Berry JA (1980) A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species. Planta 149:78–90
Flexas J, Diaz-Espejo A, Galmes J, Kaldenhoff R, Medrano H, Ribas-Carbό M (2007a) Rapid variation of mesophyll conductance in response to changes in CO2 concentration around leaves. Plant Cell Environ 30:1284–1298
Flexas J, Diaz-Espejo A, Berry JA, Cifre J, Galmes J, Kaldenhoff R, Medrano H, Ribas-Carbό M (2007b) Analysis of leakage in IRGA’s leaf chambers of open gas exchange systems: quantification and its effects in photosynthesis parameterization. J Exp Bot 58:1533–1543
Genty B, Briantais J-M, Baker N (1989) The relationship between the quantum yield of photosynthetic electron transport and quenching of chlorophyll fluorescence. Biochem Biophys Acta 990:87–92
Gu L, Sun Y (2014) Artefactual responses of mesophyll conductance to CO2 and irradiance estimated with the variable J and online isotope discrimination methods. Plant Cell Environ 37:1231–1249
Harley PC, Sharkey TD (1991) An improved model of C3 photosynthesis at high CO2: Reversed O2 sensitivity explained by lack of glycerate reentry into the chloroplast. Photosynth Res 27:169–178
Harley PC, Loreto F, Di Marco G, Sharkey TD (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
Hatakeyama Y, Ueno O (2016) Intracellular position of mitochondria and chloroplasts in bundle sheath and mesophyll cells of C3 grasses in relation to photorespiratory CO2 loss. Plant Prod Sci 19:540–551
Islam MS, Niwa Y, Takagi S (2009) Light-dependent intracellular positioning of mitochondria in Arabidopsis thaliana mesophyll cells. Plant Cell Physiol 50:1032–1040
Kebeish R, Niessen M, Thirshnaveni K, Bari R, Hirsch H-J, Rosenkranz R, Stäbler N, Schönfeld B, Kreuzaler F, Peterhänsel C (2007) Chloroplastic photorespiratory bypass increases photosynthesis and biomass production in Arabidopsis thaliana. Nat Biotechnol 25:593–599
Laisk A, Eichelmann H, Oja V, Rasulov B, Rämma H (2006) Photosystem II cycle and alternative electron flow in leaves. Plant Cell Physiol 47:972–983
Loriaux SD, Avenson TJ, Wells JM, McDermitt DK, Eckles RD, Riensche B, Genty B (2013) Closing in on maximum yield of chlorophyll fluorescence using a single multiphase flash of sub-saturating intensity. Plant Cell Environ 36:1755–1770
Morison JIL, Gifford RM (1983) Stomatal sensitivity to carbon dioxide and humidity. Plant Physiol 71:789–796
Ouyang W, Struik PC, Yin X, Yang J (2017) Stomatal conductance, mesophyll conductance, and transpiration efficiency in relation to leaf anatomy in rice and wheat genotypes under drought. J Exp Bot 68:5191–5205
Peguero-Pina JJ, Flexas J, Galmes J, Niinemets Ü, Sancho-Knapik D, Barredo G, Villarroya D, Gil-Pelegrin E (2012) Leaf anatomical properties in relation to differences in mesophyll conductance to CO2 and photosynthesis in two related Mediterranean Abies species. Plant Cell Environ 35:2121–2129
Sage TL, Sage RF (2009) The functional anatomy of rice leaves: Implications for refixation of photorespiratory CO2 and effects to engineer C4 photosynthesis into rice. Plant Cell Physiol 50:756–772
Scafaro AP, von Caemmerer S, Evans JR, Atwell BJ (2011) Temperature response of mesophyll conductance in cultivated and wild Oryza species with contrasting mesophyll cell wall thickness. Plant Cell Environ 34:1999–2008
Sharkey TD (1985) O2-insensitive photosynthesis in C3 plants: Its occurrence and a possible explanation. Plant Physiol 78:71–75
Sharkey TD, Bernacchi CJ, Farquhar GD, Singsaas EL (2007) Fitting photosynthetic carbon dioxide response curves for C3 leaves. Plant Cell Environ 30:1035–1040
Sun Y, Gu L, Dickinson RE, Norby RJ, Pallardy SG, Hoffman FM (2014) Impact of mesophyll diffusion on estimated global land CO2 fertilization. Proc Natl Acad Sci USA 111:15774–15779
Tazoe Y, von Caemmerer S, Badger MR, Evans JR (2009) Light and CO2 do not affect the mesophyll conductance to CO2 diffusion in wheat leaves. J Exp Bot 60:2291–2301
Tazoe Y, von Caemmerer S, Estavillo GM, Evans JR (2011) Using tunable diode laser spectroscopy to measure carbon isotope discrimination and mesophyll conductance to CO2 diffusion dynamically at different CO2 concentrations. Plant Cell Environ 34:580–591
Tcherkez G, Limami AN (2019) Net photosynthetic CO2 assimilation: more than just CO2 and O2 reduction cycles. New Phytol 223:520–529
Tcherkez G, Gauthier P, Buckley TN et al (2017) Leaf day respiration: low CO2 flux but high significance for metabolism and carbon balance. New Phytol 216:986–1001
Terashima I, Hanba YT, Tholen D, Niinemets Ü (2011) Leaf functional anatomy in relation to photosynthesis. Plant Physiol 155:108–116
Théroux-Rancourt G, Gilbert ME (2017) The light response of mesophyll conductance is controlled by structure across leaf profiles. Plant Cell Environ 40:726–740
Tholen D, Ethier G, Genty B, Pepin S, Zhu X-G (2012) Variable mesophyll conductance revisited: theoretical background and experimental implications. Plant Cell Environ 35:2087–2103
Tomas M, Flexas J, Copolovici L, Galmes J, Hallik L, Medrano H, Ribas-Carbo M, Tosens T, Vislap V, Niinemets Ü (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
Tosens T, Niinemets Ü, Vislap V, Eichelmann H, Castro Diez P (2012a) Developmental changes in mesophyll diffusion conductance and photosynthetic capacity under different light and water availabilities in Populus tremula: how structure constrains function. Plant Cell Environ 35:839–856
Tosens T, Niinemets Ü, Westoby M, Wright IJ (2012b) Anatomical basis of variation in mesophyll resistance in eastern Australian sclerophylls: news of a long and winding path. J Exp Bot 63:5105–5119
Ubierna N, Cernusak LA, Holloway-Phillips M, Busch FA, Cousins AB, Farquhar GD (2019) Critical review: incorporating the arrangement of mitochondria and chloroplasts into models of photosynthesis and carbon isotope discrimination. Photosynth Res 141:5–31
van der Putten PEL, Yin X, Struik PC (2018) Calibration matters: on the procedure of using the chlorophyll fluorescence method to estimate mesophyll conductance. J Plant Physiol 220:167–172
von Caemmerer S (2000) Biochemical models of leaf photosynthesis, vol 2. CSIRO Publishing, Collingwood
von Caemmerer S (2013) Steady-state models of photosynthesis. Plant Cell Environ 36:1617–1630
von Caemmerer S, Evans JR (1991) Determination of the average partial pressure of CO2 in chloroplasts from leaves of several C3 plants. Aust J Plant Physiol 18:287–305
von Caemmerer S, Evans JR, Hudson GS, Andrews TJ (1994) The kinetics of ribulose-1,5-bisphosphate carboxylase/oxygenase in vivo inferred from measurements of photosynthesis in leaves of transgenic tobacco. Planta 195:88–97
Vrábl D, Vašková M, Hronková M, Flexas J, Šantrůček J (2009) Mesophyll conductance to CO2 transport estimated by two independent methods: effect of variable CO2 concentration and abscisic acid. J Exp Bot 60:2315–2323
Warren CR, Dreyer E (2006) Temperature response of photosynthesis and internal conductance to CO2: results from two independent approaches. J Exp Bot 57:3057–3067
Yamori W, Noguchi K, Hanba YT, Terashima I (2006) Effects of internal conductance on the temperature dependence of the photosynthetic rate in spinach leaves from contrasting growth temperatures. Plant Cell Physiol 47:1069–1080
Yin X, Struik PC (2009) Theoretical reconsiderations when estimating the mesophyll conductance to CO2 diffusion in leaves of C3 plants by analysis of combined gas exchange and chlorophyll fluorescence measurements. Plant Cell Environ 32:1513–1524 (corrigendum in 33:1595)
Yin X, Struik PC (2017) Simple generalisation of a mesophyll resistance model for various intracellular arrangements of chloroplasts and mitochondria in C3 leaves. Photosynth Res 132:211–220
Yin X, Struik PC, Romero P, Harbinson J, Evers JB, van der Putten PEL, Vos J (2009) Using combined measurements of gas exchange and chlorophyll fluorescence to estimate parameters of a biochemical C3 photosynthesis model: a critical appraisal and a new integrated approach applied to leaves in a wheat (Triticum aestivum) canopy. Plant Cell Environ 32:448–464
Yin X, Sun Z, Struik PC, Gu J (2011) Evaluating a new method to estimate the rate of leaf respiration in the light by analysis of combined gas exchange and chlorophyll fluorescence measurements. J Exp Bot 62:3489–3499
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Yin, X., van der Putten, P.E.L., Belay, D. et al. Using photorespiratory oxygen response to analyse leaf mesophyll resistance. Photosynth Res 144, 85–99 (2020). https://doi.org/10.1007/s11120-020-00716-z
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11120-020-00716-z