- 1Key Laboratory of Agricultural Water Resources, Hebei Key Laboratory of Soil Ecology, Center for Agricultural Resources Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Shijiazhuang, China
- 2School of Life Sciences, Anqing Normal University, Anqing, China
- 3Environmental Sciences Division and Climate Change Science Institute, Oak Ridge National Laboratory, Oak Ridge, TN, United States
Introduction: The photosynthetic electron transport chain (ETC) is the bridge that links energy harvesting during the photophysical reactions at one end and energy consumption during the biochemical reactions at the other. Its functioning is thus fundamental for the proper balance between energy supply and demand in photosynthesis. Currently, there is a lack of understanding regarding how the structural properties of the ETC are affected by nutrient availability and plant developmental stages, which is a major roadblock to comprehensive modeling of photosynthesis.
Methods: Redox parameters reflect the structural controls of ETC on the photochemical reactions and electron transport. We conducted joint measurements of chlorophyll fluorescence (ChlF) and gas exchange under systematically varying environmental conditions and growth stages of maize and sampled foliar nutrient contents. We utilized the recently developed steady-state photochemical model to infer redox parameters of electron transport from these measurements.
Results and discussion: We found that the inferred values of these photochemical redox parameters varied with leaf macronutrient content. These variations may be caused either directly by these nutrients being components of protein complexes on the ETC or indirectly by their impacts on the structural integrity of the thylakoid and feedback from the biochemical reactions. Also, the redox parameters varied with plant morphology and developmental stage, reflecting seasonal changes in the structural properties of the ETC. Our findings will facilitate the parameterization and simulation of complete models of photosynthesis.
1 Introduction
The photosynthetic electron transport chain (ETC) between photosystem II (PSII) and photosystem I (PSI) is a key bridge in photosynthesis. It links photophysical reactions at one end and biochemical reactions at the other, while photochemical reactions occur on this bridge (Kamen, 1963; Gu et al., 2023a). The photophysical reactions are responsible for photon harvesting, exciton transfer to reaction centers, and the dissipation of excess energy as heat and fluorescence. The photochemical reactions convert the excitons into free electrons and transport them down the ETC, resulting in water splitting in the lumen, proton translocation from the stroma to the lumen, and the synthesis of NADPH and ATP. The biochemical reactions use the NADPH and ATP produced to reduce CO2 to sugars and to oxygenate RuBP (Kamen, 1963; Gu et al., 2023a). Because the three stages of the reactions are sequentially connected, the balance between them is crucial for the photosynthetic machinery to operate productively and safely in fluctuating environments.
Among the three reaction stages of photosynthesis at the leaf scale, the photochemical reactions at the ETC are the least studied. The photophysical reactions have long been probed with Pulse Amplitude-Modulated (PAM) fluorometry (Baker, 2008), whereas extensive research has been conducted on the biochemical reactions with the gas exchange technique (von Caemmerer and Farquhar, 1981). In contrast, there has been a lack of a technique that can be used to directly and quantitatively monitor the ETC and the photochemical reactions that occur on it. This lack of quantitative study of the ETC and the photochemical reactions has hindered progress on multiple fronts. For example, many components of the photophysical and biochemical reactions have been targeted for bioengineered modification to improve the efficiency of photosynthesis, including the size of the light-harvesting antenna complex (Kirst et al., 2017; Cardona et al., 2018), nonphotochemical quenching (Kromdijk et al., 2016; De Souza et al., 2022), and Rubisco kinetics (Lin et al., 2014). Currently, it is not clear whether the ETC will be capable of supporting these proposed modifications in such a way that the modified photosynthetic machinery will operate as expected under field conditions. This issue is of relevance for crop bioengineering because the ETC is likely suboptimal under current environmental conditions (Chida et al., 2007; Simkin et al., 2017; Ermakova et al., 2019) and may also need to be modified to support increased electron transport demand without suffering from excessive photooxidative stress (Sonoike, 2011; Vass, 2012). Resolving this issue will require a better understanding of the ETC and the factors that control photochemical reactions (Gu, 2023b).
Another example of areas that call for increased research on the ETC is applications of sun-induced chlorophyll fluorescence (SIF) to monitor CO2 assimilation in real time under field conditions (Sun et al., 2023a; Sun et al., 2023b). SIF is emitted during the photophysical reactions as one of the dissipation pathways of photon energy harvested by antenna complexes, while CO2 assimilation is part of the biochemical reactions. Using SIF to infer CO2 assimilation requires a complete model of photosynthesis that integrates the photophysical, photochemical, and biochemical reactions. To make this complete model of photosynthesis a reality, knowledge of electron transport along the ETC is essential (Gu et al., 2019).
Recently, we developed a mechanistic photochemical model of electron transport that relates the rate of linear electron transport to the redox state of PSII (Gu et al., 2023a). This model can be coupled with the photophysical model of Gu et al. (2019) and the biochemical model of Farquhar et al. (1980) to form a complete model of photosynthesis. However, before such a complete model of photosynthesis can be applied, it is necessary to determine the variabilities of the parameters of the photochemical model and which biotic and abiotic factors may affect these variabilities. This is analogous to the study of parameter variabilities in the biochemical model of photosynthesis (e.g., Walker et al., 2014). The parameters of the photochemical model reflect the characteristics of the redox reactions, which are controlled by the structural properties of the ETC. In this study, we conducted PAM fluorometry and gas exchange measurements on maize leaves at different developmental stages and canopy depths. After these measurements, leaf samples were taken for laboratory analyses of nutrient content and morphology. We inferred photochemical parameters from PAM fluorometry measurements and related them to foliar nutrient content and morphological indices to answer the following questions: How variable are the photochemical parameters of the ETC? Are photochemical parameters affected by leaf nutrient contents and morphology? How are variations in photochemical parameters related to each other? As the first study in this area, the answers to these questions should guide future research on other species.
2 Materials and methods
2.1 Site description and plant material
The study was conducted at the Luancheng Agro-ecological Experimental Station (37°53′ N, 114°41′ E, elevation 50.1 m). The dominant cropping system in this region is winter wheat (Triticum aestivum L.) and summer maize (Zea mays L.) rotations without fallow. Summer maize is planted from mid-June to early October, and winter wheat is cultivated from late October to early June of the following year. The climate is a warm temperate semi-humid monsoon, with cold winters and warm summers. The mean annual air temperature and precipitation for the period 1984 - 2016 were 12.6°C and 438.3 mm, and in 2022, 25.3°C and 492 mm occurred in the summer maize season (June-September), respectively. The soils are mainly cinnamon soils, containing 23.4 g kg-1 organic matter, 1.42 g kg-1 total N, 135 mg kg-1 available N, 18.9 mg kg-1 available P, and 98.4 mg kg-1 available K in the topsoil (0-20 cm).
Summer maize, Zhengdan 958, was sown on 18 June 2022 at 60 cm row spacing and 25 cm plant spacing, and harvested on 4 October 2022. Before sowing, winter wheat was harvested, and the straw was chopped and returned to the field. Subsequently, summer maize was planted with a no-till seeder equipped with a fertilizer distributor, and 600 kg ha-1 compound fertilizer (N: P2O5: K2O = 30:5:5) was applied as basal fertilizer. On the fifth day after sowing, the plants were fully irrigated to ensure even germination. With the exception of a herbicide application in early July, no other management measures were applied during plant growth.
2.2 Gas exchange and chlorophyll fluorescence measurements
Using fully expanded leaves from the bottom to the top of the canopy, combined gas exchange and chlorophyll fluorescence (ChlF) measurements of summer maize were conducted at the seedling (8 July), jointing (17 July), flowering (29 July, 11 August), filling (20 August), and maturity (7 September, 21 September) stages. Measurements were made between 9:00 and 16:00 using a LI-6400XT portable photosynthesis system (Li-Cor, Lincoln, NE, USA) fitted with a LI-6400-40 fluorescence leaf chamber. It is important to note that photosynthetic light response and CO2 response curves were measured together with ChlF on one leaf without moving the leaf chamber in between. Minimum fluorescence (Fo) and maximum fluorescence (Fm) in the dark-adapted state were measured pre-dawn. To produce light response curves and ChlF, the portable photosynthesis system was set to 400 μmol mol-1 CO2, airflow at 500 μmol s-1, photosynthetically active radiation (PAR) at 1800 μmol m-2 s-1, and leaf temperature and relative humidity at ambient levels. Once net photosynthesis rate (A) and stomatal conductance (gs) were stable, PAR was lowered sequentially to 1600, 1400, 1200, 1000, 800, 600, 400, 200, 150, 100, 50, 20, and 0 μmol m-2 s-1, and measurements were taken at each PAR as A reached steady state. Following completion of the light response curves, PAR was switched to 1600 μmol m-2 s-1, and CO2 response curves and ChlF were measured using the following sequence of CO2 concentrations: 400, 300, 200, 150, 100, 50, 400, 400, 600, 800, 1000, and 1200 μmol mol-1. Before each measurement, the CO2 and water vapor concentrations were automatically matched between the leaf and reference chambers. Finally, a total of 22 ChlF and gas exchange data were obtained for parameter calculation.
2.3 Leaf characteristics
After gas exchange measurements, leaves were harvested for the determination of leaf area with an electronic area meter (LI-3000A, Li-Cor, Lincoln, NE), leaf thickness with a digital micrometer (Mitutoyo, Japan), and dry mass after drying at 60°C to a constant mass. Specific leaf weight was then calculated as the ratio of leaf dry weight to leaf area. Subsequently, dried leaf samples were ground, nitrogen content (N) was determined by the Kjeldahl method (Kjeltec 8400, Foss, Sweden), phosphorus content (P) was measured by the molybdate colorimetric method (UV-2450, Shimadzu, Japan), and potassium (K) and calcium (Ca) contents were determined by a flame atomic absorption spectrophotometer (Analytik Jena, Germany). In total, six leaf characteristics were quantified: leaf thickness, specific leaf weight, foliar N, P, K, and Ca contents.
2.4 Inference of redox reaction parameters
Gu et al. (2023a) derived a steady-state redox model to study the relationship between electron transport and the redox state of the ETC. According to this model, the linear electron transport (LET) rate from PSII to PSI (JPSII) is photochemically related to the fraction of open PSII reaction centers (q) via the following equations:
where U is the maximum oxidation potential of the combined mobile plastoquinone/plastoquinol (PQ/PQH2) pool by the cytochrome b6f complex (Cyt). R1 and R2 are the first and second electron transport resistances, respectively. u is the second-order rate constant for the oxidation of PQH2 by the RieskeFeS protein of Cyt. rd and rr are the second-order rate constants for the electron transfer from the reduced acceptor to PQ to form PQH2 and for the reverse reaction, respectively. , , and are the total foliar concentrations of PSII, the combined PQ and PQH2 pool, and Cyt for linear electron transport, respectively. qr is the fraction of reversible PSII reaction centers, which may be less than unity due to the presence of inhibited and QB-nonreducing PSII reaction centers and the two-electron gate.
Equation 1 also contains three function modifiers, fT, fs, and fq. fT (Equation 5) is the standardized temperature (T) response function for modifying redox reactions derived from the Marcus theory of electron transfer in proteins. ET is a composite temperature sensitivity parameter related to the Gibbs free energy of activation. T0 = 298.15 K is the reference temperature. fs (Equation 6) is the light-induced thylakoid ultrastructure dynamic function, which quantifies the degree of thylakoid ultrastructural control over electron transport. This ultrastructural control is achieved by regulating the effect of macromolecular crowding on the diffusion of mobile electron carriers and the effective availability of Cyt for linear electron transport (LET) (Gu et al., 2022). v is the total thylakoid volume at a given PAR level and swells/shrinks in response to osmotic water fluxes into and out of the lumen, similar to the guard cell turgor pressure dynamics. vmax is the maximum thylakoid volume when it is fully swollen. bs controls the speed of light-induced swelling/shrinking, while cs inversely determines the maximum net impact of macromolecular crowding on the effective availability of Cyt for LET. fs varies between a value determined by cs (thylakoid minimally shrunk in the dark) and 1 (thylakoid maximally expanded in full light). fq is the photosynthetically controlled redox poise balance function between Cyt and PSII, with aq as the redox poise stoichiometry parameter. This function relates the fraction of Cyt available for LET, denoted by hCyt, to the fraction of PSII open reaction centers (i.e., q) via hCyt = fq × q. aq = 0 gives the redox isocline between Cyt and PSII (hCyt = q, Gu et al., 2023a), and is a special case that has been assumed in previous studies (e.g., Johnson and Berry, 2021). If aq > 0, hCyt > q, it indicates that PSII is more strained than Cyt for LET. If aq< 0, hCyt< q, it indicates that Cyt is more strained than PSII for LET.
It should be noted that the photochemical model of the JPSII – q relationship of Gu et al. (2023a) is complementary to, but fundamentally different from, the corresponding photophysical JPSII- q relationship of Gu et al. (2019). The former is a consequence of redox reactions at the ETC, while the latter reflects the partitioning of the absorbed energy into different dissipation pathways according to the principle of energy conservation. These two models can be used together to couple the photophysics with the photochemistry of photosynthesis. A detailed discussion of this topic is available in Gu et al. (2023a).
Using the steady-state photochemical model in conjunction with ChlF and gas exchange data, we can infer the redox parameters of summer maize at different growth stages. To do so, we calculated JPSII (μmol m-2 s-1) with . Here, is the photochemical yield of PSII, PAR (μmol m-2 s-1) is the incident photosynthetically active radiation, = 0.85 is the leaf absorptance in , and = 0.5 is the fraction of absorbed PAR allocated to PSII. We assumed a lake model of photosynthetic unit connectivity and calculated q as (Kramer et al., 2004), (Oxborough and Baker, 1997) and are minimum, maximum fluorescence yield in the light-adapted state, and Fo and Fm are minimum, maximum fluorescence yield in the dark-adapted state, respectively, and Fs is the steady-state fluorescence yield. As explained in Gu et al. (2023a), the form of Equation 1 is independent of any assumption regarding the connectivity of the photosynthetic unit, but its parameter values may depend on whether the lake or puddle model is assumed, as the value of q (but not and therefore JPSII, Kramer et al., 2004) will be different.
Equation 1 contains eight independent composite redox parameters in total (U, R1, R2, qr, ET, cs, bs, and aq). We estimated these eight redox reaction parameters by fitting the modeled JPSII as closely as possible to the measured values with the Tool for Optimizing the Open-Closed Redox Model (TOOCRM) of photosynthetic electron transport developed by and provided in Gu et al. (2023a). TOOCRM is a convenient Excel spreadsheet-based tool that uses the evolutionary method for parameter optimization.
2.5 Statistical analyses
We evaluated model performance by comparing predicted and measured JPSII – q relationships at different growth stages and canopy heights and examining their correlation. The relationships between redox parameters (U, R1, R2, qr, aq, ET, bs, and cs) and leaf characteristics (specific leaf weight, leaf thickness, N, P, K, and Ca), in addition to the relationships among these redox parameters, were examined by multiple regression. The goodness of fit of the model was evaluated by the coefficient of determination (R2) and analysis of variance (p-value). Multivariate partial least squares regression, including leave-one-out cross-validation and jackknife estimation of regression coefficients, was performed to examine the associations between redox parameters and specific leaf weight, leaf thickness, N, P, K, and Ca. Additionally, we analyzed the temporal variation of redox parameters and leaf characteristics at different growth stages.
3 Results
3.1 Model performance
As shown in Figures 1, 2, we found close agreement between the observed and modeled JPSII – q relationships at different growth stages and canopy heights and almost perfect correlations between modeled and measured JPSII. Table 1 summarizes the test statistics for model performance. Our evaluation is consistent with that reported by Gu et al. (2023a), and indicated that the steady-state photochemical model performed well in reproducing the relationship between the linear electron transport rate and the fraction of open PSII reaction centers and that our use of the JPSII – q equation to infer redox parameters was warranted.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g001.jpg)
Figure 1 Examples demonstrating the performance of the steady-state photochemical model for predicting the linear electron transport rate (JPSII) as a function of the fraction of open PSII reaction centers (ql) for top canopy leaves at different growth stages. Inset: Comparison of measured vs. modeled JPSII. (A): Jointing (17 July); (B): Flowering (29 July); (C): Filling (20 August); (D): Maturity (7 September).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g002.jpg)
Figure 2 Examples demonstrating the performance of the steady-state photochemical model for predicting the linear electron transport rate (JPSII) as a function of the fraction of open PSII reaction centers () at the flowering stage. Inset: Comparison of measured vs. modeled JPSII. (A): Bottom; (B): Lower middle; (C): Upper middle; (D): Top.
3.2 Relationships between redox parameters and leaf characteristics
We examined all 48 relationships between the eight redox parameters and six leaf characteristic parameters (8 × 6 = 48). 18 of the 48 relationships were statistically significant and are shown in Figures 3–13. U decreased with leaf thickness and Ca (Figures 3A, B) but increased with leaf N and K (Figures 3C, D). Power functions can reflect the relationships between U and leaf thickness (Figure 3A, R2 = 0.489, p = 0.000) and between U and N (Figure 3C, R2 = 0.023, p = 0.014), while a linear function and a quadratic polynomial can fit the relationship between U and Ca (Figure 3B, R2 = 0.201, p = 0.036), and that between U and K (Figure 3D, R2 = 0.745, p = 0.000), respectively. Partial least squares regression revealed that specific leaf weight (negative, p=0.031), N (positive, p = 0.015), P (negative, p=0.013), and K (positive, p = 0.041) significantly influenced U (R2 = 0.566, Figures 4A, 5A). R1 decreased with leaf K, and an exponential function can describe the relationship between them (Figure 6, R2 = 0.122, p = 0.033). R2 decreased with specific leaf weight, leaf thickness, leaf N, P, and Ca (Figures 7A–E). Linear relationships were found between R2 and specific leaf weight (Figure 7A, R2 = 0.293, p = 0.009), leaf N (Figure 7C, R2 = 0.183, p = 0.047), and P (Figure 7D, R2 = 0.262, p = 0.015), while exponential relationships were found between R2 and leaf thickness (Figure 7B, R2 = 0.113, p = 0.019) and leaf Ca (Figure 7E, R2 = 0.081, p = 0.027). Partial least squares regression revealed that specific leaf weight (p = 0.047) and N (p = 0.025) were negatively related to R2 (R2 = 0.073, Figures 4B, 5B). qr appeared to be a peaked function of specific leaf weight, leaf N, and P, with both low and high values of these leaf characteristics lowering qr (Figures 8A–C). As a result, quadratic polynomial functions can describe the relationships between qr and specific leaf weight (Figure 8A, R2 = 0.512, p = 0.008), between qr and N (Figure 8B, R2 = 0.417, p = 0.006), and between qr and P (Figure 8C, R2 = 0.328, p = 0.023). aq increased with specific leaf weight (Figure 9), and ET and bs increased with leaf P (Figures 10, 11). Quadratic polynomial functions can model the relationships between aq and specific leaf weight (Figure 9, R2 = 0.342, p = 0.023), between ET and P (Figure 10, R2 = 0.396, p = 0.008), and between bs and P (Figure 11, R2 = 0.299, p = 0.041). cs decreased with leaf thickness and leaf Ca (Figures 12A, B). Linear functions can fit the relationships between cs and leaf thickness (Figure 12A, R2 = 0.318, p = 0.006) and between cs and Ca (Figure 12B, R2 = 0.243, p = 0.020). Partial least squares regression revealed that Ca was negatively correlated with cs (R2 = 0.079, p = 0.017, Figures 4C, 5C). As shown in Figure S1, the remaining 30 of the 48 relationships are highly scattered and are not statistically significant; these plots will not be discussed here.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g003.jpg)
Figure 3 The relationships between the redox parameter U and leaf characteristics during summer maize growth. (A): U-leaf thickness; (B): U-Ca; (C): U-N; (D): U-K.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g004.jpg)
Figure 4 Loading plot of partial least squares regression for redox parameters and leaf characteristics. (A): U; (B): R2; (C): cs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g005.jpg)
Figure 5 Histogram of regression coefficients of leaf characteristics affecting redox parameters. (A): U; (B): R2; (C): cs. S: Specific leaf weight; L: Leaf thickness. ∗ Significant impact at p< 0.05.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g006.jpg)
Figure 6 The relationships between redox parameter R1 and leaf K content during summer maize growth.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g007.jpg)
Figure 7 The relationships between redox parameter R2 and leaf characteristics during summer maize growth. (A): R2-specific leaf weight; (B): R2-leaf thickness; (C): R2-N; (D): R2-P; (E): R2-Ca.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g008.jpg)
Figure 8 The relationships between redox parameter qr and leaf characteristics during summer maize growth. (A): qr-specific leaf weight; (B): qr-N; (C): qr-P.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g009.jpg)
Figure 9 The relationships between redox parameter aq and specific leaf weight during summer maize growth.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g010.jpg)
Figure 10 The relationships between redox parameter ET and leaf P content during summer maize growth.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g011.jpg)
Figure 11 The relationships between redox parameter bs and leaf P content during summer maize growth.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g012.jpg)
Figure 12 The relationships between redox parameter cs and leaf characteristics during summer maize growth. (A): cs- leaf thickness; (B): cs-Ca.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g013.jpg)
Figure 13 The relationships among redox parameters during summer maize growth. (A): aq-U; (B): aq-R2; (C): aq-bs; (D): aq-cs; (E): bs-qr; (F): bs-cs; (G): ET-R1; (H): U-R2.
3.3 Relationships between redox parameters
We examined all 28 relationships between any two of the eight redox parameters (8×7/2 = 28). Eight of the 28 relationships were statistically significant and are shown in Figure 12. aq decreased with U, R2, bs, and cs (Figures 13A–D). Logarithmic functions can describe the relationships between aq and U (Figure 13A, R2 = 0.260, p = 0.015) and between aq and R2 (Figure 13B, R2 = 0.211, p = 0.031), while quadratic polynomial functions can describe the relationships between aq and bs (Figure 13C, R2 = 0.658, p = 0.000) and between aq and cs (Figure 13D, R2 = 0.464, p = 0.003). bs decreased with qr but increased with cs (Figures 13E, F). A power function can fit the relationship between bs and qr (Figure 13E, R2 = 0.168, p = 0.000), while a quadratic polynomial function can model the relationship between bs and cs (Figure 13F, R2 = 0.285, p = 0.049). ET increased with R1 and a linear function fit the relationship between them well (Figure 13G, R2 = 0.456, p = 0.000). U increased with R2 and a logarithmic function can describe the relationship between them (Figure 13H, R2 = 0.236, p = 0.022). As shown in Figure S2, the remaining 20 of the 28 relationships did not reach statistical significance.
3.4 Temporal variations of redox parameters
Redox parameters varied significantly from seedling to maturity. U increased significantly from the seedling to the jointing stage and then decreased rapidly to the filling stage but remained almost constant at the maturity stage (Figure 14A). R1 increased from the seedling to the flowering stage and then decreased with grain maturation (Figure 14B). R2 fluctuated at high levels until the flowering stage and then remained nearly constant during the maturity stages (Figure 14C). qr increased steeply at the jointing stage and then fluctuated at relatively high levels until maturity (Figure 14D). ET fluctuated at low levels until the flowering period and increased steeply at the filling stage but decreased with grain maturation (Figure 14E). aq remained nearly constant until the flowering period and then gradually increased during the maturity stages (Figure 14F). bs was highest in seedlings, and remained relatively constant throughout the growth period, but decreased at the end of late maturation (Figure 14G). cs remained nearly constant until the flowering period and then decreased in a fluctuating way until maturity (Figure 14H).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g014.jpg)
Figure 14 The temporal variations of redox parameters during summer maize growth. (A): U; (B): R1; (C): R2; (D): qr; (E): ET; (F): aq; (G): bs; (H): cs.
3.5 Temporal variations of leaf characteristics
Leaf physical traits and nutrient contents varied dynamically between different growth stages (Figure 15). Specific leaf weight increased gradually with the growth of summer maize (Figure 15A). Leaf thickness remained relatively stable before flowering and increased thereafter (Figure 15B). Leaf N content increased significantly from seedling to jointing stages and then remained almost constant throughout the growth period (Figure 15C). Leaf P and Ca contents increased with the growth of summer maize and decreased at the late stages of grain maturation (Figures 15D, E). Leaf K content increased rapidly at the jointing stage, decreased at the flowering stage, and then remained relatively constant at the filling stage until maturity (Figure 15F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279963/fpls-14-1279963-HTML/image_m/fpls-14-1279963-g015.jpg)
Figure 15 The temporal variations of leaf characteristics during summer maize growth. (A): specific leaf weight; (B): leaf thickness; (C): N; (D): P; (E): Ca; (F): K.
4 Discussion and conclusion
Our study found that there is considerable variation in photochemical redox parameters inferred from PAM fluorometry measurements in maize. Part of these variations may be due to measurement noise and a lack of data constraints on parameter optimization (i.e., overfitting), while statistically significant relationships between these parameters and independently measured leaf macronutrient content and morphology, as well as their systematic variations with developmental stage, suggest that at least part of these variations may be due to changes in the structure of the ETC. All parameters in the photochemical model of electron transport are composite parameters that are influenced by multiple components of the ETC. For example, U is the main parameter that determines the capacity of the ETC. It is the product of the rate constant for the oxidation of PQH2 by the RieskeFeS protein of Cyt and the abundances of the PQ and Cyt pools. Although PQ does not contain nitrogen, Cyt contains nitrogen in its hemes (Tikhonov, 2023). Therefore, the positive correlation between U and leaf nitrogen probably reflects the higher abundance of Cyt at higher leaf nitrogen contents. The positive correlation between U and leaf nitrogen contributes to the balance between electron transport and CO2 assimilation, which is catalyzed by Rubisco, a dominant sink of leaf nitrogen (Evans and Clarke, 2019). Potassium is not involved in the structure of organic macromolecules. However, it is the most important inorganic osmotic ion in plant cells and plays an important role in controlling stomatal conductance and in maintaining the structural integrity of granal thylakoids, where photosynthetic macromolecules are located (Tränkner et al, 2018), which likely explains the positive correlation between U and potassium.
The resistance parameters R1 and R2 are given by R1 = and , respectively. rd and rr are the second-order rate constants for the forward and reverse reactions between the reduced acceptor of PSII and PQ, respectively, while u is the second-order rate constant for the oxidation of PQH2 by the RieskeFeS protein of Cyt. Thus, R1 is inversely related to the efficiency of PQ to take away electrons from PSII, and R2 is inversely related to the capacity of PSII to supply electrons relative to the capacity of Cyt to oxidize PQH2. The smaller the values of R1 and R2, the higher the efficiency of the ETC. Thus, the decrease of R1 with potassium and R2 with nitrogen and phosphorus is consistent with the expected impact of these nutrients on electron transport.
A significant fraction of PSII reaction centers may be QB-nonreducing (Tomek et al., 2003; Schansker and Strasser, 2005; Vredenberg et al., 2006; Vredenberg, 2011), which can result in qr being significantly less than 1. Currently, little is known about why such reaction centers exist at all and what their photochemical functions may be. The significant correlation between qr and leaf nitrogen may suggest that the fraction of QB-nonreducing centers is affected by leaf nitrogen content.
Any factors that affect the relative abundances of PSII and Cyt can affect the value of the parameter aq. A negative aq indicates that there is a greater constraint on Cyt than on PSII for electron transport, whereas a positive aq indicates the opposite (Gu et al., 2023a). The increase in aq with specific leaf weight likely indicates that the relative abundance of PSII vs. Cyt in a leaf depends on leaf-specific weight. Variations in bs and cs may reflect systematic changes in thylakoid structure. Such changes may include the number of grana per thylakoid, the height of a granum, and the density of ion channels, which affect ion exchange between the lumen and stroma and, therefore, osmotic water fluxes. These changes can affect the degree and speed of thylakoid swelling/shrinking and, therefore, the values of bs and cs. Variations in ET can be caused by changes in the Gibbs free energy of activation (Gu et al., 2023a), which in turn can be caused by changes in the reduction potentials of redox reactions (Silverstein, 2012; Bostick et al., 2018). Thus, the variation of ET with leaf phosphorus content may indicate that the latter affects the reduction potentials of redox reactions along the ETC.
Currently, it is difficult to provide an in-depth discussion of all the relationships reported in this study. There are very few previous studies that have investigated how the structures and redox reactions of the ETC may be affected by nutrient contents and plant developmental stages, and thus may provide guidance for our analysis. Nevertheless, our study provides conclusive evidence that the structural properties of the electron transport chain are not static and may vary with plant growth conditions and developmental stages, allowing the photophysical, photochemical, and biochemical reactions to be balanced on seasonal time scales. Our findings also show that it may be possible to parameterize photochemical parameters of electron transport with leaf nutrient content and morphological properties for large-scale modeling applications.
Author’s note
This manuscript was co-authored by UT-Battelle, LLC under contract number DE-AC05-00OR22725 with the U.S. Department of Energy. The publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, worldwide license to publish or reproduce, or allow others to do so, the published form of this manuscript for the purposes of the United States Government. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan).
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
XL: Investigation, Formal analysis, Writing – original draft. YQ: Investigation, Writing – review & editing. WZ: Investigation, Writing – review & editing. WD: Supervision, Writing – review & editing. LG: Methodology, Writing – original draft.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by the National Key Research and Development Program of China (2021YFD1901104, 2022YFD1901604), the “Strategic Priority Research Program” of the Chinese Academy of Sciences (XDA28020303, XDA26040103), and the Key Research and Development Program of Hebei Province (22326412D). LG is supported by the U.S. Department of Energy (DOE), Office of Science, Biological and Environmental Research Program. ORNL is managed by UT-Battelle, LLC, for DOE under contract DE-AC05-00OR22725.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1279963/full#supplementary-material
References
Baker, N. R. (2008). Chlorophyll fluorescence: a probe of photosynthesis in vivo. Annu. Rev. Plant Biol. 59, 89–113. doi: 10.1146/annurev.arplant.59.032607.092759
Bostick, C. D., Mukhopadhyay, S., Pecht, I., Sheves, M., Cahen, D., Lederman, D. (2018). Protein bioelectronics: a review of what we do and do not know. Rep. Prog. Phys. 81, 026601. doi: 10.1088/1361-6633/aa85f2
Cardona, T., Shao, S., Nixon, P. J. (2018). Enhancing photosynthesis in plants: the light reactions. Essays Biochem. 62, 85–94. doi: 10.1042/EBC20170015
Chida, H., Nakazawa, A., Akazaki, H., Hirano, T., Suruga, K., Ogawa, M., et al. (2007). Expression of the algal cytochrome c6 gene in Arabidopsis enhances photosynthesis and growth. Plant Cell Physiol. 48, 948–957. doi: 10.1093/pcp/pcm064
De Souza, A. P., Burgess, S. J., Doran, L., Hansen, J., Manukyan, L., Maryn, N., et al. (2022). Soybean photosynthesis and crop yield are improved by accelerating recovery from photoprotection. Science 377, 851–854. doi: 10.1126/science.adc9831
Ermakova, M., Lopez-Calcagno, P. E., Raines, C. A., Furbank, R. T., von Caemmerer, S. (2019). Overexpression of the RieskeFeS protein of the Cytochrome b6f complex increases C4 photosynthesis in Setaria viridis. Commun. Biol. 2, 314. doi: 10.1038/s42003-019-0561-9
Evans, J. R., Clarke, V. C. (2019). The nitrogen cost of photosynthesis. J. Exp. Bot. 70, 7–15doi: 10.1093/jxb/ery366
Farquhar, G. D., von Caemmerer, S., Berry, J. A. (1980). A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species. Planta 149, 78–90. doi: 10.1007/BF00386231
Gu, L. (2023b). Optimizing the electron transport chain to sustainably improve photosynthesis. Plant Physiol., kiad490. doi: 10.1093/plphys/kiad490
Gu, L., Grodzinski, B., Han, J., Marie, T., Zhang, Y.-J., Song, Y. C., et al. (2022). Granal thylakoid structure and function: Explaining an enduring mystery of higher plants. New Phytol. 236, 319–329. doi: 10.1111/nph.18371
Gu, L., Grodzinski, B., Han, J., Marie, T., Zhang, Y.-J., Song, Y. C., et al. (2023a). An exploratory steady-state redox model of photosynthetic linear electron transport for use in complete modelling of photosynthesis for broad applications. Plant Cell Environ. 46, 1540–1561. doi: 10.1111/pce.14563
Gu, L., Han, J., Wood, J. D., Chang, C. Y. Y., Sun, Y. (2019). Sun-induced Chl fluorescence and its importance for biophysical modeling of photosynthesis based on light reactions. New Phytol. 223, 1179–1191. doi: 10.1111/nph.15796
Johnson, J. E., Berry, J. A. (2021). The role of Cytochrome b6f in the control of steady-state photosynthesis: a conceptual and quantitative model. Photosynth. Res. 148, 101–136. doi: 10.1007/s11120-021-00840-4
Kirst, H., Gabilly, S. T., Niyogi, K. K., Lemaux, P. G., Melis, A. (2017). Photosynthetic antenna engineering to improve crop yields. Planta 245, 1009–1020. doi: 10.1007/s00425-017-2659-y
Kramer, D. M., Johnson, G., Kiirats, O., Edwards, G. E. (2004). New fluorescence parameters for the determination of QA redox state and excitation energy fluxes. Photosynth. Res. 79, 209–218. doi: 10.1023/B:PRES.0000015391.99477.0d
Kromdijk, J., Głowacka, K., Leonelli, L., Gabilly, S. T., Iwai, M., Niyogi, K. K., et al. (2016). Improving photosynthesis and crop productivity by accelerating recovery from photoprotection. Science 354, 857–861. doi: 10.1126/science.aai8878
Lin, M. T., Occhialini, A., Andralojc, P. J., Parry, M. A. J., Hanson, M. R. (2014). A faster Rubisco with potential to increase photosynthesis in crops. Nature 513, 547–550. doi: 10.1038/nature13776
Oxborough, K., Baker, N. R. (1997). Resolving chlorophyll a fluorescence images of photosynthetic efficiency into photochemical and non-photochemical components – calculation of qP and Fv‘/Fm‘without measuring Fo‘. Photosynth. Res. 54, 135–142. doi: 10.1023/A:1005936823310
Schansker, G., Strasser, R. J. (2005). Quantification of non-QB-reducing centers in leaves using a far-red pre-illumination. Photosynth. Res. 84, 145–151. doi: 10.1007/s11120-004-7156-z
Silverstein, T. P. (2012). Marcus theory: thermodynamics can control the kinetics of electron transfer reactions. J. Chem. Educ. 89, 1159–1167. doi: 10.1021/ed1007712
Simkin, A. J., McAusland, L., Lawson, T., Raines, C. A. (2017). Overexpression of the RieskeFeS protein increases electron transport rates and biomass yield. Plant Physiol. 175, 134–145. doi: 10.1104/pp.17.00622
Sonoike, K. (2011). Photoinhibition of photosystem I. Physiol. Plant 142, 56–64. doi: 10.1111/j.1399-3054.2010.01437.x
Sun, Y., Gu, L., Wen, J., van der Tol, C., Porcar-Castell, A., Joiner, J., et al. (2023a). From remotely sensed solar-induced chlorophyll fluorescence to ecosystem structure, function, and service: Part I—Harnessing theory. Glob. Change Biol. 29, 2926–2952. doi: 10.1111/gcb.16634
Sun, Y., Wen, J., Gu, L., Joiner, J., Chang, C. Y., van der Tol, C., et al. (2023b). From remotely-sensed solar-induced chlorophyll fluorescence to ecosystem structure, function, and service: Part II—Harnessing data. Glob. Change Biol. 29, 2893–2925. doi: 10.1111/gcb.16646
Tikhonov, A. N. (2023). The cytochrome b6f complex: plastoquinol oxidation and regulation of electron transport in chloroplasts. Photosynth. Res. doi: 10.1007/s11120-023-01034-w
Tomek, P., Ilík, P., Lazár, D., Štroch, M., Nauš, J. (2003). On the determination of QB-non-reducing photosystem II centers from chlorophyll a fluorescence induction. Plant Sci. 164, 665–670. doi: 10.1016/S0168-9452(03)00029-3
Tränkner, M., Tavakol, E., Jákli, B. (2018). Functioning of potassium and magnesium in photosynthesis, photosynthate translocation and photoprotection. Physiol. Plant 163, 414–431. doi: 10.1111/ppl.12747
Vass, I. (2012). Molecular mechanisms of photodamage in the photosystem II complex. Biochim. Biophys. Acta 1817, 209–217. doi: 10.1016/j.bbabio.2011.04.014
von Caemmerer, S., Farquhar, G. D. (1981). Some relationships between the biochemistry of photosynthesis and the gas exchange of leaves. Planta 153, 376–387. doi: 10.1007/BF00384257
Vredenberg, W. (2011). Kinetic analyses and mathematical modeling of primary photochemical and photoelectrochemical processes in plant photosystems. Biosystems 103, 138–151. doi: 10.1016/j.biosystems.2010.10.016
Vredenberg, W., Kasalicky, V., Durchan, M., Prasil, O. (2006). The chlorophyll a fluorescence induction pattern in chloroplasts upon repetitive single turnover excitations: Accumulation and function of QB-nonreducing centers. Biochim. Biophys. Acta 1757, 173–181. doi: 10.1016/j.bbabio.2006.02.006
Walker, A. P., Beckerman, A. P., Gu, L., Kattge, J., Cernusak, L. A., Domingues, T. F., et al. (2014). The relationship of leaf photosynthetic traits – Vcmax and Jmax – to leaf nitrogen, leaf phosphorus, and specific leaf area: a meta-analysis and modeling study. Ecol. Evol. 4, 3218–3235. doi: 10.1002/ece3.1173
Keywords: photosynthesis, photosynthetic electron transport, redox parameters, leaf characteristics, plant growth stages, maize
Citation: Liu X, Qiao Y, Zhou W, Dong W and Gu L (2023) Determinants of photochemical characteristics of the photosynthetic electron transport chain of maize. Front. Plant Sci. 14:1279963. doi: 10.3389/fpls.2023.1279963
Received: 19 August 2023; Accepted: 25 October 2023;
Published: 20 November 2023.
Edited by:
Qingpeng Yang, Institute of Applied Ecology, Chinese Academy of Sciences (CAS), ChinaReviewed by:
Songzhu Zhang, Shenyang Agricultural University, ChinaWenhua Xu, Jilin Academy of Agricultural Sciences (CAAS), China
Copyright © 2023 Liu, Qiao, Zhou, Dong and Gu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wenxu Dong, dongwx@sjziam.ac.cn; Lianhong Gu, lianhong-gu@ornl.gov