The slow phase of chlorophyll a fluorescence induction in silico: Origin of the S-M fluorescence rise

From leaf to multiscale models of photosynthesis: applications and challenges for crop improvement

To keep up with the growth of human population and to circumvent deleterious effects of global climate change, it is essential to enhance crop yield to achieve higher production. Here we review mathematical models of oxygenic photosynthesis that are extensively used, and discuss in depth a subset that accounts for diverse approaches providing solutions to our objective. These include models (1) to study different ways to enhance photosynthesis, such as fine-tuning antenna size, photoprotection and electron transport; (2) to bioengineer carbon metabolism; and (3) to evaluate the interactions between the process of photosynthesis and the seasonal crop dynamics, or those that have included statistical whole-genome prediction methods to quantify the impact of photosynthesis traits on the improvement of crop yield. We conclude by emphasizing that the results obtained in these studies clearly demonstrate that mathematical modelling is a key tool to examine different approaches to improve photosynthesis for better productivity, while effective multiscale crop models, especially those that also include remote sensing data, are indispensable to verify different strategies to obtain maximized crop yields.

Relationship between photosynthetic performance and yield loss in winter oilseed rape (Brassica napus L.) under frost conditions

Winter oilseed rape (Brassica napus L.), the principal oilseed crop in Europe, is notably vulnerable to spring frosts that can drastically reduce yields in ways that are challenging to predict with standard techniques. Our research focused on evaluating the efficacy of photosynthetic efficiency analysis in this crop and identifying specific chlorophyll fluorescence parameters severely impacted by frost, which could serve as noninvasive biomarkers for yield decline. The experiments were carried out in semi-controlled conditions with several treatments: a control, one day at -3°C, three days at -3°C, one day at -6°C, and three days at -6°C. We employed continuous-excitation and pulse-amplitude-modulation chlorophyll fluorescence measurements to assess plant sensitivity to frost. Also, plant gas exchange and chlorophyll content index measurements were performed. Certain parameters strongly correlated with final yield losses, thereby establishing a basis for developing new agricultural protocols to predict and mitigate frost damage in rapeseed crops accurately.

Photosynthetic Performance and Yield Losses of Winter Rapeseed (Brassica napus L. var. napus) Caused by Simulated Hail

Winter oilseed rape (Brassica napus L.), Europe's foremost oilseed crop, is significantly impacted by hailstorms, leading to substantial yield reductions that are difficult to predict and measure using conventional methods. This research aimed to assess the effectiveness of photosynthetic efficiency analysis for predicting yield loss in winter rapeseed subjected to hail exposure. The aim was to pinpoint the chlorophyll fluorescence parameters most affected by hail stress and identify those that could act as non-invasive biomarkers of yield loss. The study was conducted in partially controlled conditions (greenhouse). Stress was induced in the plants by firing plastic balls with a 6 mm diameter at them using a pneumatic device, which launched the projectiles at speeds of several tens of meters per second. Measurements of both continuous-excitation and pulse-modulated-amplitude chlorophyll fluorescence were engaged to highlight the sensitivity of the induction curve and related parameters to hail stress. Our research uncovered that some parameters such as Fs, Fm', ΦPSII, ETR, Fo, Fv/Fm, and Fv/Fo measured eight days after the application of stress had a strong correlation with final yield, thus laying the groundwork for the creation of new practical protocols in agriculture and the insurance industry to accurately forecast damage to rapeseed crops due to hail stress.

Comparative modeling of fluorescence and P700 induction kinetics for alga Scenedesmus sp. obliques and cyanobacterium Synechocystis sp. PCC 6803. Role of state 2-state 1 transitions and redox state of plastoquinone pool

The model of thylakoid membrane system (T-M model) (Belyaeva et al. Photosynth Res 2019, 140:1-19) has been improved in order to analyze the induction data for dark-adapted samples of algal (Scenedesmus obliques) and cyanobacterial (Synechocystis sp. PCC 6803) cells. The fluorescence induction (FI) curves of Scenedesmus were measured at light exposures of 5 min, while FI and P700 redox transformations of Synechocystis were recorded in parallel for 100 s intervals. Kinetic data comprising the OJIP-SMT fluorescence induction and OABCDEF P700+ absorbance changes were used to study the processes underlying state transitions qT2→1 and qT1→2 associated with the increase/decrease in Chl fluorescence emission. A formula with the Hill kinetics (Ebenhöh et al. Philos Trans R Soc B 2014, 369:20130223) was introduced into the T-M model, with a new variable to imitate the flexible size of antenna AntM(t) associated with PSII. Simulations revealed that the light-harvesting capacity of PSII increases with a corresponding decrease for that of PSI upon the qT2→1 transition induced by plastoquinone (PQ) pool oxidation. The complete T-M model fittings were attained on Scenedesmus or Synechocystis fast waves OJIPS of FI, while SMT wave of FI was reproduced at intervals shorter than 5 min. Also the fast P700 redox transitions (OABC) for Synechocystis were fitted exactly. Reasonable sets of algal and cyanobacterial electron/proton transfer (ET/PT) parameters were found. In the case of Scenedesmus, ET/PT traits remained the same irrespective of modeling with or without qT2→1 transitions. Simulations indicated a high extent (20%) of the PQ pool reduction under dark conditions in Synechocystis compared to 2% in Scenedesmus.

Revelations on photosystem II, thermoluminescence, and artificial photosynthesis: a retrospective of Govindjee from fundamentals to applications

Photosynthesis, as one of the most important chemical reactions, has powered our planet for over four billion years on a massive scale. This review summarizes and highlights the major contributions of Govindjee from fundamentals to applications in photosynthesis. His research included primary photochemistry measurements, in the picosecond time scale, in both Photosystem I and II and electron transport leading to NADP reduction, using two light reactions. He was the first to suggest the existence of P680, the reaction center of PSII, and to prove that it was not an artefact of Chlorophyll a fluorescence. For most photobiologists, Govindjee is best known for successfully exploiting Chlorophyll a fluorescence to understand the various steps in photosynthesis as well as to predict plant productivity. His contribution in resolving the controversy on minimum number of quanta in favor of 8-12 vs 3-4, needed for the evolution of one molecule of oxygen, is a milestone in the area of photosynthesis research. Furthermore, together with Don DeVault, he is the first to provide the correct theory of thermoluminescence in photosynthetic systems. His research productivity is very high: ~ 600 published articles and total citations above 27,000 with an h-index of 82. He is a recipient of numerous awards and honors including a 2022: Lifetime Achievement Award of the International Society of Photosynthesis Research. We hope that the retrospective of Govindjee described in this work will inspire and stimulate the readers to continue probing the photosynthetic apparatuses with new discoveries and breakthroughs.

A Synthetic Review of Various Dimensions of Non-Destructive Plant Stress Phenotyping

Non-destructive plant stress phenotyping begins with traditional one-dimensional (1D) spectroscopy, followed by two-dimensional (2D) imaging, three-dimensional (3D) or even temporal-three-dimensional (T-3D), spectral-three-dimensional (S-3D), and temporal-spectral-three-dimensional (TS-3D) phenotyping, all of which are aimed at observing subtle changes in plants under stress. However, a comprehensive review that covers all these dimensional types of phenotyping, ordered in a spatial arrangement from 1D to 3D, as well as temporal and spectral dimensions, is lacking. In this review, we look back to the development of data-acquiring techniques for various dimensions of plant stress phenotyping (1D spectroscopy, 2D imaging, 3D phenotyping), as well as their corresponding data-analyzing pipelines (mathematical analysis, machine learning, or deep learning), and look forward to the trends and challenges of high-performance multi-dimension (integrated spatial, temporal, and spectral) phenotyping demands. We hope this article can serve as a reference for implementing various dimensions of non-destructive plant stress phenotyping.

The pattern of photosynthetic response and adaptation to changing light conditions in lichens is linked to their ecological range

Epiphytic lichens constitute an important component of biodiversity in both deforested and forest ecosystems. Widespread occurrence is the domain of generalist lichens or those that prefer open areas. While, many stenoecious lichens find shelter only in a shaded interior of forests. Light is one of the factors known to be responsible for lichen distribution. Nevertheless, the effect of light intensity on photosynthesis of lichen photobionts remain largely unknown. We investigated photosynthesis in lichens with different ecological properties in relation to light as the only parameter modified during the experiments. The aim was to find links between this parameter and habitat requirements of a given lichen. We applied the methods based on a saturating light pulse and modulated light to perform comprehensive analyses of fast and slow chlorophyll fluorescence transient (OJIP and PSMT) combined with quenching analysis. We also examined the rate of CO₂ assimilation. Common or generalist lichens, i.e. Hypogymnia physodes, Flavoparmelia caperata and Parmelia sulcata, are able to adapt to a wide range of light intensity. Moreover, the latter species, which prefers open areas, dissipates the excess energy most efficiently. Conversely, Cetrelia cetrarioides considered an old-growth forest indicator, demonstrates definitely lower range of energy dissipation than other species, although it assimilates CO₂ efficiently both at low and high light. We conclude that functional plasticity of the thylakoid membranes of photobionts largely determines the dispersal abilities of lichens and light intensity is one of the most important factors determining the specificity of a species to a given habitat.

Detection of antibiotic and microplastic pollutants in Chrysanthemum coronarium L. based on chlorophyll fluorescence

Large amounts of antibiotics and microplastics are used in daily life and agricultural production, which affects not only plant growth but also potentially the food safety of vegetables and other plant products. Fast detection of the presence of antibiotics and microplastics in leafy vegetables is of great interest to the public. In this work, a method was developed to detect sulfadiazine and polystyrene, commonly used antibiotics and microplastics, in vegetables by measuring and modeling photosystem II chlorophyll a fluorescence (ChlF) emission from leaves. Chrysanthemum coronarium L., a common beverage and medicinal plant, was used to verify the developed method. Scanning electron microscopy, transmission electron microscopy, and liquid chromatograph-mass spectrometer analysis were used to show the presence of the two pollutants in the samples. The developed kinetic model could describe measured ChlF variations with an average relative error of 0.6%. The model parameters estimated for the chlorophyll a fluorescence induction kinetics curve (OJIP) induction can differentiate the two types of stresses while the commonly used ChlF OJIP induction characteristics cannot. This work provides a concept to detect antibiotic pollutants and microplastic pollutants in vegetables based on ChlF.

Molecular, Brownian, kinetic and stochastic models of the processes in photosynthetic membrane of green plants and microalgae

The paper presents the results of recent work at the Department of Biophysics of the Biological Faculty, Lomonosov Moscow State University on the kinetic and multiparticle modeling of processes in the photosynthetic membrane. The detailed kinetic models and the rule-based kinetic Monte Carlo models allow to reproduce the fluorescence induction curves and redox transformations of the photoactive pigment P700 in the time range from 100 ns to dozens of seconds and make it possible to reveal the role of individual carriers in their formation for different types of photosynthetic organisms under different illumination regimes, in the presence of inhibitors, under stress conditions. The fitting of the model curves to the experimental data quantifies the reaction rate constants that cannot be directly measured experimentally, including the non-radiative thermal relaxation reactions. We use the direct multiparticle models to explicitly describe the interactions of mobile photosynthetic carrier proteins with multienzyme complexes both in solution and in the biomembrane interior. An analysis of these models reveals the role of diffusion and electrostatic factors in the regulation of electron transport, the influence of ionic strength and pH of the cellular environment on the rate of electron transport reactions between carrier proteins. To describe the conformational intramolecular processes of formation of the final complex, in which the actual electron transfer occurs, we use the methods of molecular dynamics. The results obtained using kinetic and molecular models supplement our knowledge of the mechanisms of organization of the photosynthetic electron transport processes at the cellular and molecular levels.

Natural variation in the fast phase of chlorophyll a fluorescence induction curve (OJIP) in a global rice minicore panel

Photosynthesis can be probed through Chlorophyll a fluorescence induction (FI), which provides detailed insight into the electron transfer process in Photosystem II, and beyond. Here, we have systematically studied the natural variation of the fast phase of the FI, i.e. the OJIP phase, in rice. The OJIP phase of the Chl a fluorescence induction curve is referred to as "fast transient" lasting for less than a second; it is obtained after a dark-adapted sample is exposed to saturating light. In the OJIP curve, "O" stands for "origin" (minimal fluorescence), "P" for "peak" (maximum fluorescence), and J and I for inflection points between the O and P levels. Further, F_o is the fluorescence intensity at the "O" level, whereas F_m is the intensity at the P level, and F_v (= F_m - F_o) is the variable fluorescence. We surveyed a set of quantitative parameters derived from the FI curves of 199 rice accessions, grown under both field condition (FC) and growth room condition (GC). Our results show a significant variation between Japonica (JAP) and Indica (IND) subgroups, under both the growth conditions, in almost all the parameters derived from the OJIP curves. The ratio of the variable to the maximum (F_v/F_m) and of the variable to the minimum (F_v/F_o) fluorescence, the performance index (PI_abs), as well as the amplitude of the I-P phase (A_I-P) show higher values in JAP compared to that in the IND subpopulation. In contrast, the amplitude of the O-J phase (A_O-J) and the normalized area above the OJIP curve (S_m) show an opposite trend. The performed genetic analysis shows that plants grown under GC appear much more affected by environmental factors than those grown in the field. We further conducted a genome-wide association study (GWAS) using 11 parameters derived from plants grown in the field. In total, 596 non-unique significant loci based on these parameters were identified by GWAS. Several photosynthesis-related proteins were identified to be associated with different OJIP parameters. We found that traits with high correlation are usually associated with similar genomic regions. Specifically, the thermal phase of FI, which includes the amplitudes of the J-I and I-P subphases (A_J-I and A_I-P) of the OJIP curve, is, in turn, associated with certain common genomic regions. Our study is the first one dealing with the natural variations in rice, with the aim to characterize potential candidate genes controlling the magnitude and half-time of each of the phases in the OJIP FI curve.

Inhibition of Primary Photosynthesis in Desiccating Antarctic Lichens Differing in Their Photobionts, Thallus Morphology, and Spectral Properties

Five macrolichens of different thallus morphology from Antarctica (King George Island) were used for this ecophysiological study. The effect of thallus desiccation on primary photosynthetic processes was examined. We investigated the lichens' responses to the relative water content (RWC) in their thalli during the transition from a wet (RWC of 100%) to a dry state (RWC of 0%). The slow Kautsky kinetics of chlorophyll fluorescence (ChlF) that was recorded during controlled dehydration (RWC decreased from 100 to 0%) and supplemented with a quenching analysis revealed a polyphasic species-specific response of variable fluorescence. The changes in ChlF at a steady state (Fs), potential and effective quantum yields of photosystem II (F_V/F_M, Φ_PSII), and nonphotochemical quenching (NPQ) reflected a desiccation-induced inhibition of the photosynthetic processes. The dehydration-dependent fall in F_V/F_M and Φ_PSII was species-specific, starting at an RWC range of 22-32%. The critical RWC for Φ_PSII was below 5%. The changes indicated the involvement of protective mechanisms in the chloroplastic apparatus of lichen photobionts at RWCs of below 20%. In both the wet and dry states, the spectral reflectance curves (SRC) (wavelength 400-800 nm) and indices (NDVI, PRI) of the studied lichen species were measured. Black Himantormia lugubris showed no difference in the SRCs between wet and dry state. Other lichens showed a higher reflectance in the dry state compared to the wet state. The lichen morphology and anatomy data, together with the ChlF and spectral reflectance data, are discussed in relation to its potential for ecophysiological studies in Antarctic lichens.

Electron transport in Tradescantia leaves acclimated to high and low light: thermoluminescence, PAM-fluorometry, and EPR studies

Using thermoluminescence, PAM-fluorometry, and electron paramagnetic resonance (EPR) for assaying electron transport processes in chloroplasts in situ, we have compared photosynthetic characteristics in Tradescantia fluminensis leaves grown under low light (LL, 50-125 µmol photons m^-2 s^-1) or high light (HL, 875-1000 µmol photons m^-2 s^-1) condition. We found differences in the thermoluminescence (TL) spectra of LL- and HL-acclimated leaves. The LL and HL leaves show different proportions of the Q (~ 0 °C) and B (~ 25-30 °C) bands in their TL spectra; the ratios of the "light sums" of the Q and B bands being S_Q/S_B ≈ 1/1 (LL) and S_Q/S_B ≈ 1/3 (HL). This suggests the existence of different redox states of electron carriers on the acceptor side of PSII in LL and HL leaves, which may be affected, in particular, by different capacities of their photo-reducible PQ pools. Enhanced content of PQ in chloroplasts of LL leaves may be the reason for an efficient performance of photosynthesis at low irradiance. Kinetic studies of slow induction of Chl a fluorescence and measurements of P₇₀₀ photooxidation by EPR demonstrate that HL leaves have faster (about 2 times) response to switching on actinic light as compared to LL leaves grown at moderate irradiation. HL leaves also show higher non-photochemical quenching (NPQ) of Chl a fluorescence. These properties of HL leaves (faster response to light and generation of enhanced NPQ) reflect the flexibility of their photosynthetic apparatus, providing sustainability and rapid response to fluctuations of environmental light intensity and solar stress resistance. Analysis of time-courses of the EPR signals of [Formula: see text] induced by far-red (λ_max = 707 nm), exciting predominantly PSI, and white light, exciting both PSI and PSII, suggests that there is a contribution of cyclic electron flow around PSI to electron flow through PSI in HL leaves. The data obtained are discussed in terms of photosynthetic apparatus sustainability of HL and LL leaves under variable irradiation conditions.

Temperature-dependent regulation of electron transport and ATP synthesis in chloroplasts in vitro and in silico

The significance of temperature-dependent regulation of photosynthetic apparatus (PSA) is determined by the fact that plant temperature changes with environmental temperature. In this work, we present a brief overview of temperature-dependent regulation of photosynthetic processes in class B chloroplasts (thylakoids) and analyze these processes using a computer model that takes into account the key stages of electron and proton transport coupled to ATP synthesis. The rate constants of partial reactions were parametrized on the basis of experimental temperature dependences of partial photosynthetic processes: (1) photosystem II (PSII) turnover and plastoquinone (PQ) reduction, (2) the plastoquinol (PQH₂) oxidation by the cytochrome (Cyt) b₆f complex, (3) the ATP synthase activity, and (4) the proton leak from the thylakoid lumen. We consider that PQH₂ oxidation is the rate-limiting step in the intersystem electron transport. The parametrization of the rate constants of these processes is based on earlier experimental data demonstrating strong correlations between the functional and structural properties of thylakoid membranes that were probed with the lipid-soluble spin labels embedded into the membranes. Within the framework of our model, we could adequately describe a number of experimental temperature dependences of photosynthetic reactions in thylakoids. Computer modeling of electron and proton transport coupled to ATP synthesis supports the notion that PQH₂ oxidation by the Cyt b₆f complex and proton pumping into the lumen are the basic temperature-dependent processes that determine the overall electron flux from PSII to molecular oxygen and the net ATP synthesis upon variations of temperature. The model describes two branches of the temperature dependence of the post-illumination reduction of [Formula: see text] characterized by different activation energies (about 60 and ≤ 3.5 kJ mol^-1). The model predicts the bell-like temperature dependence of ATP formation, which arises from the balance of several factors: (1) the thermo-induced acceleration of electron transport through the Cyt b₆f complex, (2) deactivation of PSII photochemistry at sufficiently high temperatures, and (3) acceleration of the passive proton outflow from the thylakoid lumen bypassing the ATP synthase complex. The model describes the temperature dependence of experimentally measured parameter P/2e, determined as the ratio between the rates of ATP synthesis and pseudocyclic electron transport (H₂O → PSII → PSI → O₂).

Development of a minimized model structure and a feedback control framework for regulating photosynthetic activities

In this work, the main activities of the plant photosynthesis process are discussed to yield a minimized mathematical model structure with photosystem II (PSII) chlorophyll a fluorescence (ChlF) as a measurable output. After experimental validation of the model structure, we demonstrate that the states of the photosynthetic process may be observed by using this model and the extended Kalman filter method. We then show a feedback control framework that can be used to alter a given photosynthetic activity. The control framework is demonstrated with an example in which PSII ChlF is used as the feedback signal and light intensity is used as a controllable process input to regulate plastoquinone reduction. Although there are caveats, and further research is needed, the results lay the groundwork for further research on novel methods for optimization and regulation of photosynthetic activities, with a goal for sustainability.

Model quantification of the light-induced thylakoid membrane processes in Synechocystis sp. PCC 6803 in vivo and after exposure to radioactive irradiation

Measurements of OJIP-SMT patterns of fluorescence induction (FI) in Synechocystis sp. PCC 6803 (Synechocystis) cells on a time scale up to several minutes were mathematically treated within the framework of thylakoid membrane (T-M) model (Belyaeva et al., Photosynth Res 140:1-19, 2019) that was renewed to account for the state transitions effects. Principles of describing electron transfer in reaction centers of photosystems II and I (PSII and PSI) and cytochrome b₆f complex remained unchanged, whereas parameters for dissipative reactions of non-radiative charge recombination were altered depending on the oxidation state of Q_B-site (neutral, reduced by one electron, empty, reduced by two electrons). According to our calculations, the initial content of plastoquinol (PQH₂) in the total quinone pool of Synechocystis cells adapted to darkness for 10 min ranged between 20 and 40%. The results imply that the PQ pool mediates photosynthetic and respiratory charge flows. The redistribution of PBS antenna units responsible for the increase of Chl fluorescence in cyanobacteria (qT_{2 → 1}) upon state 2 → 1 transition or the fluorescence lowering (qT_{1 → 2}) due to state 1 → 2 transition were described in the model by exponential functions. Parameters of dynamically changed effective cross section were found by means of simulations of OJIP-SMT patterns observed on Synechocystis cells upon strong (3000 μmol photons m^-2s^-1) and moderate (1000 μmol photons m^-2s^-1) actinic light intensities. The corresponding light constant values kL_ΣAnt = 1.2 ms^-1 and 0.4 ms^-1 define the excitation of total antenna pool dynamically redistributed between PSII and PSI reaction centers. Although the OCP-induced quenching of antenna excitation is not involved in the model, the main features of the induction signals have been satisfactorily explained. In the case of strong illumination, the effective cross section decreases by approximately 33% for irradiated Synechocystis cells as compared to untreated cells. Under moderate light, the irradiated Synechocystis cells showed in simulations the same cross section as the untreated cells. The thylakoid model renewed with state transitions description allowed simulation of fluorescence induction OJIP-SMT curves detected on time scale from microseconds to minutes.

Comparative effect of tenuazonic acid, diuron, bentazone, dibromothymoquinone and methyl viologen on the kinetics of Chl a fluorescence rise OJIP and the MR820 signal

In this study, the comparative effect of TeA, DCMU, bentazone, DBMIB and MV on prompt fluorescence and the MR₈₂₀ signal was simultaneously analyzed to provide an insight into how to elucidate their precise influence on Ageratina adenophora photosystems. The herbicides that interrupt electron transport beyond Q_A, such as TeA, DCMU and bentazone, mainly increased the J-step level of fluorescence rise kinetics as a result of accumulation of Q_A^-, but showed differences in detail. The IP phase disappeared in the presence of DCMU and bentazone with a significant increase in F_O value. TeA treatment retained the IP phase with lowering F_M. As an inhibitor of plastoquinone re-oxidation, DBMIB increased the I-step (IP phase almost unnoticable) without changing F_O and F_M values. MV blocking PSI electron transfer through intercepting electrons from the FeS clusters suppressed the IP phase by decreasing the P level. Considering the W_IP kinetics, TeA and DBMIB also affected PSI activity. After DCMU and MV treatment, the major change in the MR₈₂₀ kinetics was the loss of the slow phase due to the complete prevention of electron movement from PSII to re-reduce PC⁺ and P₇₀₀⁺. TeA, bentazone and DBMIB clearly suppressed the MR₈₂₀ slow phase and decreased the re-reduction rate of PC⁺ and P₇₀₀⁺ (V_red), significantly. However, there were still parts of electrons being donated to PC⁺ and P₇₀₀⁺, showing a smaller slow phase and PC⁺ and P₇₀₀⁺ re-reduction rate. Additionally, TeA and DBMIB also somewhat declined the fast phase and PC and P₇₀₀ oxidation rate (V_ox).

Transient depletion of transported metabolites in the streaming cytoplasm of Chara upon shading the long-distance transmission pathway

Export of reducing power from chloroplasts to cytoplasm serves to balance the NADPH/ATP ratio that is optimal for CO₂ assimilation. Rapid cytoplasmic streaming in characean algae conveys the exported metabolites downstream towards the shaded plastids where envelope transporters may operate for the import of reducing power in accordance with the direction of concentration gradients. Import of reducing equivalents by chloroplasts in the analyzed area transiently enhances the pulse-modulated chlorophyll fluorescence F' controlled by the redox state of photosystem II acceptor Q_A. When the microfluidic pathway was transferred to darkness while the analyzed cell area remained in dim background light, the amplitude of cyclosis-mediated F' changes dropped sharply and then recovered within 5-10 min. The suppression of long-distance signaling indicates temporal depletion of transmitted metabolites in the streaming cytoplasm. The return to overall background illumination induced an exceptionally large F' response to the first local light pulse admitted to a remote cell region. This indicates the appearance of excess reductants in the streaming cytoplasm at a certain stage of photosynthetic induction. The results suggest highly dynamic exchange of metabolites between stationary chloroplasts lining the microfluidic pathway and the streaming cytoplasm upon light-dark and dark-light transitions. Evidence is obtained that slow stages of chlorophyll fluorescence induction in algae with rapid cytoplasmic streaming directly depend on cyclosis-mediated long-distance delivery of metabolites produced far beyond the analyzed cell area.

Photosynthesis: basics, history and modelling

BACKGROUND

With limited agricultural land and increasing human population, it is essential to enhance overall photosynthesis and thus productivity. Oxygenic photosynthesis begins with light absorption, followed by excitation energy transfer to the reaction centres, primary photochemistry, electron and proton transport, NADPH and ATP synthesis, and then CO2 fixation (Calvin-Benson cycle, as well as Hatch-Slack cycle). Here we cover some of the discoveries related to this process, such as the existence of two light reactions and two photosystems connected by an electron transport 'chain' (the Z-scheme), chemiosmotic hypothesis for ATP synthesis, water oxidation clock for oxygen evolution, steps for carbon fixation, and finally the diverse mechanisms of regulatory processes, such as 'state transitions' and 'non-photochemical quenching' of the excited state of chlorophyll a.

SCOPE

In this review, we emphasize that mathematical modelling is a highly valuable tool in understanding and making predictions regarding photosynthesis. Different mathematical models have been used to examine current theories on diverse photosynthetic processes; these have been validated through simulation(s) of available experimental data, such as chlorophyll a fluorescence induction, measured with fluorometers using continuous (or modulated) exciting light, and absorbance changes at 820 nm (ΔA820) related to redox changes in P700, the reaction centre of photosystem I.

CONCLUSIONS

We highlight here the important role of modelling in deciphering and untangling complex photosynthesis processes taking place simultaneously, as well as in predicting possible ways to obtain higher biomass and productivity in plants, algae and cyanobacteria.

Slow induction of chlorophyll a fluorescence excited by blue and red light in Tradescantia leaves acclimated to high and low light

Tradescantia is a good model for assaying induction events in higher plant leaves. Chlorophyll (Chl) fluorescence serves as a sensitive reporter of the functional state of photosynthetic apparatus in chloroplasts. The fluorescence time-course depends on the leaf growth conditions and actinic light quality. In this work, we investigated slow induction of Chl a fluorescence (SIF) excited by blue light (BL, λ_max = 455 nm) or red light (RL, λ_max = 630 nm) in dark-adapted leaves of Tradescantia fluminensis acclimated to high light (~ 1000 µmol photons m^-2 s^-1; HL) or low light (~ 100 µmol photons m^-2 s^-1; LL). Our special interest was focused on the contribution of the avoidance response to SIF kinetics. Bearing in mind that BL and RL have different impacts on photoreceptors that initiate chloroplast movements within the cell (accumulation/avoidance responses), we have compared the SIF patterns during the action of BL and RL. The time-courses of SIF and kinetics of non-photochemical quenching (NPQ) of Chl a fluorescence revealed a certain difference when leaves were illuminated by BL or RL. In both cases, the yield of fluorescence rose to the maximal level P and then, after the lag-phase P-S-M₁, the fluorescence level decreased toward the steady state T (via the intermediate phases M₁-M₂ and M₂-T). In LL-acclimated leaves, the duration of the P-S-M₁ phase was almost two times longer that in HL-grown plants. In the case of BL, the fluorescence decay included the transient phase M₁-M₂. This phase was obscure during the RL illumination. Non-photochemical quenching of Chl a fluorescence has been quantified as [Formula: see text], where [Formula: see text] and [Formula: see text] stand for the fluorescence response to saturating pulses of light applied to dark-adapted and illuminated samples, respectively. The time-courses of such a formally determined NPQ value were markedly different during the action of RL and BL. In LL-grown leaves, BL induced higher NPQ as compared to the action of RL. In HL-grown plants, the difference between the NPQ responses to BL and RL illumination was insignificant. Comparing the peculiarities of Chl a fluorescence induced by BL and RL, we conclude that the avoidance response can provide a marked contribution to SIF and NPQ generation. The dependence of NPQ on the quality of actinic light suggests that chloroplast movements within the cell have a noticeable impact on the formally determined NPQ value. Analyzing kinetics of post-illumination decay of NPQ in the context of solar stress resistance, we have found that LL-acclimated Tradescantia leaves are more vulnerable to strong light than the HL-grown leaves.

Modelling and simulation of chlorophyll fluorescence from PSII of a plant leaf as affected by both illumination light intensities and temperatures

The emission of chlorophyll fluorescence (ChlF) from photosystem II (PSII) of plant leaves the couple with photoelectron transduction cascades in photosynthetic reactions and can be used to probe photosynthetic efficiency and plant physiology. Because of population increase, food shortages, and global warming, it is becoming more and more urgent to enhance plant photosynthesis efficiency by controlling plant growth rate. An effective model structure is essential for plant control strategy development. However, there is a lack of reporting on modelling and simulation of PSII activities under the interaction of both illumination light intensities and temperatures, which are the two important controllable factors affecting, plant growth, especially for a greenhouse. In this work, the authors extended their work on modelling photosynthetic activities as affected by light and temperature to cover both the interaction effects of illumination light intensities and temperature on ChlF emission. Experiments on ChlF were performed under different light intensities and temperatures and used to validate the developed model structure. The average relative error between experimental data and model fitting is <0.3%, which shows the effectiveness of the developed model structure. Simulations were performed to show the interaction effect of light and temperature effects on photosynthetic activities.

Modelling and simulation of photosystem II chlorophyll fluorescence transition from dark-adapted state to light-adapted state

Green houses play a vital role in modern agriculture. Artificial light illumination is very important in a green house. While light is necessary for plant growth, excessive light in a green house may not bring more profit and even damages plants. Developing a plant-physiology-based light control strategy in a green house is important, which implies that a state-space model on photosynthetic activities is very useful because modern control theories and techniques are usually developed according to model structures in the state space. In this work, a simplified model structure on photosystem II activities was developed with seven state variables and chlorophyll fluorescence (ChlF) as the observable variable. Experiments on ChlF were performed. The Levenberg-Marquardt algorithm was used to estimate model parameters from experimental data. The model structure can fit experimental data with a small relative error (<2%). ChlF under different light intensities were simulated to show the effect of light intensity on ChlF emission. A simplified model structure with fewer state variables and model parameters will be more robust to perturbations and model parameter estimation. The model structure is thus expected useful in future green-house light control strategy development.

Analyzing both the fast and the slow phases of chlorophyll a fluorescence and P700 absorbance changes in dark-adapted and preilluminated pea leaves using a Thylakoid Membrane model

The dark-to-light transitions enable energization of the thylakoid membrane (TM), which is reflected in fast and slow (OJIPSMT or OABCDE) stages of fluorescence induction (FI) and P700 oxidoreduction changes (ΔA810). A Thylakoid Membrane model (T-M model), in which special emphasis has been placed on ferredoxin-NADP+-oxidoreductase (FNR) activation and energy-dependent qE quenching, was applied for quantifying the kinetics of FI and ΔA810. Pea leaves were kept in darkness for 15 min and then the FI and ΔA810 signals were measured upon actinic illumination, applied either directly or after a 10-s light pulse coupled with a subsequent 10-s dark interval. On the time scale from 40 µs to 30 s, the parallel T-M model fittings to both FI and ΔA810 signals were obtained. The parameters of FNR activation and the buildup of qE quenching were found to differ for dark-adapted and preilluminated leaves. At the onset of actinic light, photosystem II (PSII) acceptors were oxidized (neutral) after dark adaptation, while the redox states with closed and/or semiquinone QA(-)QB(-) forms were supposedly generated after preillumination, and did not relax within the 10 s dark interval. In qE simulations, a pH-dependent Hill relationship was used. The rate constant of heat losses in PSII antenna kD(t) was found to increase from the basic value kDconst, at the onset of illumination, to its maximal level kDvar due to lumenal acidification. In dark-adapted leaves, a low value of kDconst of ∼ 2 × 106 s-1 was found. Simulations on the microsecond to 30 s time scale revealed that the slow P-S-M-T phases of the fluorescence induction were sensitive to light-induced FNR activation and high-energy qE quenching. Thus, the corresponding time-dependent rate constants kD(t) and kFNR(t) change substantially upon the release of electron transport on the acceptor side of PSI and during the NPQ development. The transitions between the cyclic and linear electron transport modes have also been quantified in this paper.

A sixty-year tryst with photosynthesis and related processes: an informal personal perspective

After briefly describing my early collaborative work at the University of Allahabad, that had laid the foundation of my research life, I present here some of our research on photosynthesis at the University of Illinois at Urbana-Champaign, randomly selected from light absorption to NADP⁺ reduction in plants, algae, and cyanobacteria. These include the fact that (i) both the light reactions I and II are powered by light absorbed by chlorophyll (Chl) a of different spectral forms; (ii) light emission (fluorescence, delayed fluorescence, and thermoluminescence) by plants, algae, and cyanobacteria provides detailed information on these reactions and beyond; (iii) primary photochemistry in both the photosystems I (PS I) and II (PS II) occurs within a few picoseconds; and (iv) most importantly, bicarbonate plays a unique role on the electron acceptor side of PS II, specifically at the two-electron gate of PS II. Currently, the ongoing research around the world is, and should be, directed towards making photosynthesis better able to deal with the global issues (such as increasing population, dwindling resources, and rising temperature) particularly through genetic modification. However, basic research is necessary to continue to provide us with an understanding of the molecular mechanism of the process and to guide us in reaching our goals of increasing food production and other chemicals we need for our lives.

Low temperature induced modulation of photosynthetic induction in non-acclimated and cold-acclimated Arabidopsis thaliana: chlorophyll a fluorescence and gas-exchange measurements

Cold acclimation modifies the photosynthetic machinery and enables plants to survive at sub-zero temperatures, whereas in warm habitats, many species suffer even at non-freezing temperatures. We have measured chlorophyll a fluorescence (ChlF) and CO₂ assimilation to investigate the effects of cold acclimation, and of low temperatures, on a cold-sensitive Arabidopsis thaliana accession C24. Upon excitation with low intensity (40 µmol photons m^- 2 s^- 1) ~ 620 nm light, slow (minute range) ChlF transients, at ~ 22 °C, showed two waves in the SMT phase (S, semi steady-state; M, maximum; T, terminal steady-state), whereas CO₂ assimilation showed a linear increase with time. Low-temperature treatment (down to - 1.5 °C) strongly modulated the SMT phase and stimulated a peak in the CO₂ assimilation induction curve. We show that the SMT phase, at ~ 22 °C, was abolished when measured under high actinic irradiance, or when 3-(3, 4-dichlorophenyl)-1, 1- dimethylurea (DCMU, an inhibitor of electron flow) or methyl viologen (MV, a Photosystem I (PSI) electron acceptor) was added to the system. Our data suggest that stimulation of the SMT wave, at low temperatures, has multiple reasons, which may include changes in both photochemical and biochemical reactions leading to modulations in non-photochemical quenching (NPQ) of the excited state of Chl, "state transitions," as well as changes in the rate of cyclic electron flow through PSI. Further, we suggest that cold acclimation, in accession C24, promotes "state transition" and protects photosystems by preventing high excitation pressure during low-temperature exposure.

Simulation of chlorophyll fluorescence rise and decay kinetics, and P700-related absorbance changes by using a rule-based kinetic Monte-Carlo method

A model of primary photosynthetic reactions in the thylakoid membrane was developed and its validity was tested by simulating three types of experimental kinetic curves: (1) the light-induced chlorophyll a fluorescence rise (OJIP transients) reflecting the stepwise transition of the photosynthetic electron transport chain from the oxidized to the fully reduced state; (2) the dark relaxation of the flash-induced fluorescence yield attributed to the Q_A^- oxidation kinetics in PSII; and (3) the light-induced absorbance changes near 820 or 705 nm assigned to the redox transitions of P₇₀₀ in PSI. A model was implemented by using a rule-based kinetic Monte-Carlo method and verified by simulating experimental curves under different treatments including photosynthetic inhibitors, heat stress, anaerobic conditions, and very high light intensity.

Acclimation of shade-tolerant and light-resistant Tradescantia species to growth light: chlorophyll a fluorescence, electron transport, and xanthophyll content

In this study, we have compared the photosynthetic characteristics of two contrasting species of Tradescantia plants, T. fluminensis (shade-tolerant species), and T. sillamontana (light-resistant species), grown under the low light (LL, 50-125 µmol photons m^-2 s^-1) or high light (HL, 875-1000 µmol photons m^-2 s^-1) conditions during their entire growth period. For monitoring the functional state of photosynthetic apparatus (PSA), we measured chlorophyll (Chl) a emission fluorescence spectra and kinetics of light-induced changes in the heights of fluorescence peaks at 685 and 740 nm (F ₆₈₅ and F ₇₄₀). We also compared the light-induced oxidation of P₇₀₀ and assayed the composition of carotenoids in Tradescantia leaves grown under the LL and HL conditions. The analyses of slow induction of Chl a fluorescence (SIF) uncovered different traits in the LL- and HL-grown plants of ecologically contrasting Tradescantia species, which may have potential ecophysiological significance with respect to their tolerance to HL stress. The fluorometry and EPR studies of induction events in chloroplasts in situ demonstrated that acclimation of both Tradescantia species to HL conditions promoted faster responses of their PSA as compared to LL-grown plants. Acclimation of both species to HL also caused marked changes in the leaf anatomy and carotenoid composition (an increase in Violaxanthin + Antheraxantin + Zeaxanthin and Lutein pools), suggesting enhanced photoprotective capacity of the carotenoids in the plants grown in nature under high irradiance. Collectively, the results of the present work suggest that the mechanisms of long-term PSA photoprotection in Tradescantia are based predominantly on the light-induced remodeling of pigment-protein complexes in chloroplasts.

High spatio-temporal-resolution detection of chlorophyll fluorescence dynamics from a single chloroplast with confocal imaging fluorometer

BACKGROUND

Chlorophyll fluorescence (CF) is a key indicator to study plant physiology or photosynthesis efficiency. Conventionally, CF is characterized by fluorometers, which only allows ensemble measurement through wide-field detection. For imaging fluorometers, the typical spatial and temporal resolutions are on the order of millimeter and second, far from enough to study cellular/sub-cellular CF dynamics. In addition, due to the lack of optical sectioning capability, conventional imaging fluorometers cannot identify CF from a single cell or even a single chloroplast.

RESULTS AND DISCUSSION

Here we demonstrated a fluorometer based on confocal imaging, that not only provides high contrast images, but also allows CF measurement with spatiotemporal resolution as high as micrometer and millisecond. CF transient (the Kautsky curve) from a single chloroplast is successfully obtained, with both the temporal dynamics and the intensity dependences corresponding well to the ensemble measurement from conventional studies. The significance of confocal imaging fluorometer is to identify the variation among individual chloroplasts, e.g. the temporal position of the P-S-M phases, and the half-life period of P-T decay in the Kautsky curve, that are not possible to analyze with wide-field techniques. A linear relationship is found between excitation intensity and the temporal positions of P-S-M peaks/valleys in the Kautsky curve. Based on the CF transients, the photosynthetic quantum efficiency is derived with spatial resolution down to a single chloroplast. In addition, an interesting 6-order increase in excitation intensity is found between wide-field and confocal fluorometers, whose pixel integration time and optical sectioning may account for this substantial difference.

CONCLUSION

Confocal imaging fluorometers provide micrometer and millisecond CF characterization, opening up unprecedented possibilities toward detailed spatiotemporal analysis of CF transients and its propagation dynamics, as well as photosynthesis efficiency analysis, on the scale of organelles, in a living plant.

Frequently asked questions about chlorophyll fluorescence, the sequel

Using chlorophyll (Chl) a fluorescence many aspects of the photosynthetic apparatus can be studied, both in vitro and, noninvasively, in vivo. Complementary techniques can help to interpret changes in the Chl a fluorescence kinetics. Kalaji et al. (Photosynth Res 122:121-158, 2014a) addressed several questions about instruments, methods and applications based on Chl a fluorescence. Here, additional Chl a fluorescence-related topics are discussed again in a question and answer format. Examples are the effect of connectivity on photochemical quenching, the correction of F V /F M values for PSI fluorescence, the energy partitioning concept, the interpretation of the complementary area, probing the donor side of PSII, the assignment of bands of 77 K fluorescence emission spectra to fluorescence emitters, the relationship between prompt and delayed fluorescence, potential problems when sampling tree canopies, the use of fluorescence parameters in QTL studies, the use of Chl a fluorescence in biosensor applications and the application of neural network approaches for the analysis of fluorescence measurements. The answers draw on knowledge from different Chl a fluorescence analysis domains, yielding in several cases new insights.