Mesophyll conductance and reaction-diffusion models for CO2 transport in C3 leaves; needs, opportunities and challenges

Elevated CO2 and ammonium nitrogen promoted the plasticity of two maple in great lakes region by adjusting photosynthetic adaptation

Introduction

Climate change-related CO2 increases and different forms of nitrogen deposition are thought to affect the performance of plants, but their interactions have been poorly studied.

Methods

This study investigated the responses of photosynthesis and growth in two invasive maple species, amur maple (Acer ginnala Maxim.) and boxelder maple (Acer negundo L.), to elevated CO2 (400 µmol mol-1 (aCO2) vs. 800 µmol mol-1 (eCO2) and different forms of nitrogen fertilization (100% nitrate, 100% ammonium, and an equal mix of the two) with pot experiment under controlled conditions.

Results and discussion

The results showed that eCO2 significantly promoted photosynthesis, biomass, and stomatal conductance in both species. The biochemical limitation of photosynthesis was switched to RuBP regeneration (related to Jmax) under eCO2 from the Rubisco carboxylation limitation (related to Vcmax) under aCO2. Both species maximized carbon gain by lower specific leaf area and higher N concentration than control treatment, indicating robust morphological plasticity. Ammonium was not conducive to growth under aCO2, but it significantly promoted biomass and photosynthesis under eCO2. When nitrate was the sole nitrogen source, eCO2 significantly reduced N assimilation and growth. The total leaf N per tree was significantly higher in boxelder maple than in amur maple, while the carbon and nitrogen ratio was significantly lower in boxelder maple than in amur maple, suggesting that boxelder maple leaf litter may be more favorable for faster nutrient cycling. The results suggest that increases in ammonium under future elevated CO2 will enhance the plasticity and adaptation of the two maple species.

Simulated Analysis of Influence of Changes in H+-ATPase Activity and Membrane CO2 Conductance on Parameters of Photosynthetic Assimilation in Leaves

Photosynthesis is an important process in plants which influences their development and productivity. Many factors can control the efficiency of photosynthesis, including CO₂ conductance of leaf mesophyll, which affects the CO₂ availability for Rubisco. It is known that electrical stress signals can decrease this conductance, and the response is probably caused by inactivation of H⁺-ATPase in the plasma membrane. In the current work, we analyzed the influence of both CO₂ conductance in the plasma membrane, and chloroplast envelopes and H⁺-ATPase activity on photosynthetic CO₂ assimilation, using a two-dimensional mathematical model of photosynthesis in leaves. The model included a description of assimilation on the basis of the Farquhar-von Caemmerer-Berry model, ion transport through the plasma membrane, diffusion of CO₂ in the apoplast, and transport of CO₂ through the plasma membrane and chloroplast envelope. The model showed that the photosynthetic CO₂ assimilation rate was mainly dependent on the plasma membrane and chloroplast envelope conductance; direct influence of the H⁺-ATPase activity (through changes in pH and CO₂/HCO₃^- concentration ratio) on this rate was weak. In contrast, both changes in CO₂ conductance of the plasma membrane and chloroplast envelopes and changes in the H⁺-ATPase activity influenced spatial heterogeneity of the CO₂ assimilation on the leaf surface in the simulated two-dimensional system. These effects were also observed under simultaneous changes in the CO₂ conductance of the plasma membrane and H⁺-ATPase activity. Qualitatively similar influence of changes in the CO₂ conductance of the plasma membrane and chloroplast envelopes, and changes in the H⁺-ATPase activity on photosynthesis were shown for two different densities of stomata in the simulated leaf; however, lowering the density of stomata decreased the assimilation rate and increased the heterogeneity of assimilation. The results of the model analysis clarify the potential influence of H⁺-ATPase inactivation on photosynthesis, and can be the basis for development of new methods for remote sensing of the influence of electrical signals.

Development of Two-Dimensional Model of Photosynthesis in Plant Leaves and Analysis of Induction of Spatial Heterogeneity of CO2 Assimilation Rate under Action of Excess Light and Drought

Photosynthesis is a key process in plants that can be strongly affected by the actions of environmental stressors. The stressor-induced photosynthetic responses are based on numerous and interacted processes that can restrict their experimental investigation. The development of mathematical models of photosynthetic processes is an important way of investigating these responses. Our work was devoted to the development of a two-dimensional model of photosynthesis in plant leaves that was based on the Farquhar-von Caemmerer-Berry model of CO₂ assimilation and descriptions of other processes including the stomatal and transmembrane CO₂ fluxes, lateral CO₂ and HCO₃^- fluxes, transmembrane and lateral transport of H⁺ and K⁺, interaction of these ions with buffers in the apoplast and cytoplasm, light-dependent regulation of H⁺-ATPase in the plasma membrane, etc. Verification of the model showed that the simulated light dependences of the CO₂ assimilation rate were similar to the experimental ones and dependences of the CO₂ assimilation rate of an average leaf CO₂ conductance were also similar to the experimental dependences. An analysis of the model showed that a spatial heterogeneity of the CO₂ assimilation rate on a leaf surface should be stimulated under an increase in light intensity and a decrease in the stomatal CO₂ conductance or quantity of the open stomata; this prediction was supported by the experimental verification. Results of the work can be the basis of the development of new methods of the remote sensing of the influence of abiotic stressors (at least, excess light and drought) on plants.

Photosynthetic physiological response of water-saving and drought-resistant rice to severe drought under wetting-drying alternation irrigation

Water-saving and drought-resistant rice (WDR) is widely grown in central China in recent years. However, studies have not explored the interaction effect of WDR and irrigation regimes on drought-resistance capacities under severe drought at sensitive growth periods. A pot experiment was conducted using a WDR cultivar Hanyou73 (HY73) and traditional high-yielding and drought-sensitive cultivar Huiliangyou 898 (HLY898). Three irrigation regimes, including flooding irrigation (W1), mild wetting-drying alternation irrigation (W2), and severe wetting-drying alternation irrigation (W3), were applied before heading. At heading, severe drought with -50 KPa soil water potential was established for all treatments and cultivars. The findings showed that cultivar HY73 under W2 treatment had the highest yield, 1000-grain yield, filled grain, relative water content, and photosynthesis potential compared with the other combinations. The higher net photosynthetic rate (P_n ) was attributed to larger mesophyll conductance (g_m ) in drought for cultivar HY73 under W2 treatment compared with that for cultivar HLY898 and the other water treatments. Enhanced photo-respiration rate may be an important photoprotection mechanism for achieving high P_n for cultivar HY73 coupled with W2 treatment than for other combinations in drought. The relative expression level of OsPIP1;1 gene was significantly down-regulated during drought in all cultivars and water regimes. But OsPIP1;2, OsPIP2;3, OsTIP2;2, and OsTIP3;1 genes were upregulated to alleviate the significant decrease in g_s and g_m under drought. These results suggest that WDR and mild wetting-drying alternation irrigation (W2) have significant interaction effects in improving photosynthetic production potential by maintaining higher g_m under severe drought.

Anatomical variation of mesophyll conductance due to salt stress in Populus cathayana females and males growing under different inorganic nitrogen sources

Synergistic regulation in leaf architecture and photosynthesis is essential for salt tolerance. However, how plant sex and inorganic nitrogen sources alter salt stress-dependent photosynthesis remains unknown. Leaf anatomical characteristics and photosynthesis of Populus cathayana Rehder females and males were investigated under salt stress conditions combined with nitrate NO3- and ammonium NH4+ supplies to clarify the underlying mechanisms. In salt-stressed females, we observed an increased mesophyll spongy cell density, a reduced chloroplast density, a decreased surface area of chloroplasts adjacent to the intercellular air space (Sc/S) and an increased mesophyll cell area per transverse section width (S/W), consequently causing mesophyll conductance (gm) and photosynthesis inhibition, especially under NH4+ supply. Conversely, males with a greater mesophyll palisade tissue thickness and chloroplast density, but a lower spongy cell density had lower S/W and higher Sc/S, and higher gm and photosynthesis. NH4+-fed females had a lower CO2 conductance through cell wall and stromal conductance perpendicular to the cell wall, but a higher chloroplast conductance from the cell wall (gcyt1) than females supplied with NO3-, whereas males had a higher chloroplast conductance and lower CO2 conductance through cell wall when supplied with NO3- instead of NH4+ under salt stress. These findings indicate sex-specific strategies in coping with salt stress related to leaf anatomy and gm under both types of nitrogen supplies, which may contribute to sex-specific CO2 capture and niche segregation.

Constant ratio of Cc to Ci under various CO2 concentrations and light intensities, and during progressive drought, in seedlings of Japanese white birch

Constant mesophyll conductance (g_m), and two-resistance g_m model (involved in resistances of cell wall and chloroplast), where g_m reaches maximum under higher CO₂ concentrations, cannot describe the phenomenon that g_m decreases with increasing intercellular CO₂ concentration (C_i) under relatively higher CO₂ concentrations. Yin et al. (2020) proposed a g_m model, according to which the ratio of chloroplastic CO₂ concentration (C_c) to C_i is constant in the two-resistance g_m model, which can describe the decreasing g_m with increasing C_i. In the present study, we investigated the relationship between C_c and C_i in leaves of Japanese white birch by using simultaneous measurements of gas exchange and chlorophyll fluorescence under various CO₂ concentrations, light intensities, and during progressive drought. Across the range of ambient CO₂ from 50 to 1000 μmol mol^-1, and light intensities of 50 to 2000 μmol m^-2 s^-1, measured under well irrigation, the ratio of C_c to C_i kept constant. During the progressive drought, overestimated C_i due to stomatal patchiness and/or cuticular transpiration was empirically corrected (threshold: stomatal conductance < 0.08 mol H₂O m^-2 s^-1) from the A/C_i response measured under adequate irrigation. The ratio of C_c to C_i during progressive drought (predawn leaf potential reached ≈ - 2 MPa) also remained constant irrespective of soil drying rate in various pot sizes. The present study suggests the involvement of some physiologically regulative mechanisms to keep C_c:C_i ratio constant, which might act on g_m in addition to the physical interaction of diffusive resistances in the cell components.

Mesophyll conductance: the leaf corridors for photosynthesis

Besides stomata, the photosynthetic CO2 pathway also involves the transport of CO2 from the sub-stomatal air spaces inside to the carboxylation sites in the chloroplast stroma, where Rubisco is located. This pathway is far to be a simple and direct way, formed by series of consecutive barriers that the CO2 should cross to be finally assimilated in photosynthesis, known as the mesophyll conductance (gm). Therefore, the gm reflects the pathway through different air, water and biophysical barriers within the leaf tissues and cell structures. Currently, it is known that gm can impose the same level of limitation (or even higher depending of the conditions) to photosynthesis than the wider known stomata or biochemistry. In this mini-review, we are focused on each of the gm determinants to summarize the current knowledge on the mechanisms driving gm from anatomical to metabolic and biochemical perspectives. Special attention deserve the latest studies demonstrating the importance of the molecular mechanisms driving anatomical traits as cell wall and the chloroplast surface exposed to the mesophyll airspaces (Sc/S) that significantly constrain gm. However, even considering these recent discoveries, still is poorly understood the mechanisms about signaling pathways linking the environment a/biotic stressors with gm responses. Thus, considering the main role of gm as a major driver of the CO2 availability at the carboxylation sites, future studies into these aspects will help us to understand photosynthesis responses in a global change framework.

Recent developments in mesophyll conductance in C3, C4, and crassulacean acid metabolism plants

The conductance of carbon dioxide (CO₂ ) from the substomatal cavities to the initial sites of CO₂ fixation (g_m ) can significantly reduce the availability of CO₂ for photosynthesis. There have been many recent reviews on: (i) the importance of g_m for accurately modelling net rates of CO₂ assimilation, (ii) on how leaf biochemical and anatomical factors influence g_m , (iii) the technical limitation of estimating g_m , which cannot be directly measured, and (iv) how g_m responds to long- and short-term changes in growth and measurement environmental conditions. Therefore, this review will highlight these previous publications but will attempt not to repeat what has already been published. We will instead initially focus on the recent developments on the two-resistance model of g_m that describe the potential of photorespiratory and respiratory CO₂ released within the mitochondria to diffuse directly into both the chloroplast and the cytosol. Subsequently, we summarize recent developments in the three-dimensional (3-D) reaction-diffusion models and 3-D image analysis that are providing new insights into how the complex structure and organization of the leaf influences g_m . Finally, because most of the reviews and literature on g_m have traditionally focused on C₃ plants we review in the final sections some of the recent developments, current understanding and measurement techniques of g_m in C₄ and crassulacean acid metabolism (CAM) plants. These plants have both specialized leaf anatomy and either a spatially or temporally separated CO₂ concentrating mechanisms (C₄ and CAM, respectively) that influence how we interpret and estimate g_m compared with a C₃ plants.

In silico study of the role of cell growth factors in photosynthesis using a virtual leaf tissue generator coupled to a microscale photosynthesis gas exchange model

Computational tools that allow in silico analysis of the role of cell growth and division on photosynthesis are scarce. We present a freely available tool that combines a virtual leaf tissue generator and a two-dimensional microscale model of gas transport during C3 photosynthesis. A total of 270 mesophyll geometries were generated with varying degrees of growth anisotropy, growth extent, and extent of schizogenous airspace formation in the palisade mesophyll. The anatomical properties of the virtual leaf tissue and microscopic cross-sections of actual leaf tissue of tomato (Solanum lycopersicum L.) were statistically compared. Model equations for transport of CO2 in the liquid phase of the leaf tissue were discretized over the geometries. The virtual leaf tissue generator produced a leaf anatomy of tomato that was statistically similar to real tomato leaf tissue. The response of photosynthesis to intercellular CO2 predicted by a model that used the virtual leaf tissue geometry compared well with measured values. The results indicate that the light-saturated rate of photosynthesis was influenced by interactive effects of extent and directionality of cell growth and degree of airspace formation through the exposed surface of mesophyll per leaf area. The tool could be used further in investigations of improving photosynthesis and gas exchange in relation to cell growth and leaf anatomy.

Using a reaction-diffusion model to estimate day respiration and reassimilation of (photo)respired CO2 in leaves

Methods using gas exchange measurements to estimate respiration in the light (day respiration R d ) make implicit assumptions about reassimilation of (photo)respired CO2 ; however, this reassimilation depends on the positions of mitochondria. We used a reaction-diffusion model without making these assumptions to analyse datasets on gas exchange, chlorophyll fluorescence and anatomy for tomato leaves. We investigated how R d values obtained by the Kok and the Yin methods are affected by these assumptions and how those by the Laisk method are affected by the positions of mitochondria. The Kok method always underestimated R d . Estimates of R d by the Yin method and by the reaction-diffusion model agreed only for nonphotorespiratory conditions. Both the Yin and Kok methods ignore reassimilation of (photo)respired CO2 , and thus underestimated R d for photorespiratory conditions, but this was less so in the Yin than in the Kok method. Estimates by the Laisk method were affected by assumed positions of mitochondria. It did not work if mitochondria were in the cytosol between the plasmamembrane and the chloroplast envelope. However, mitochondria were found to be most likely between the tonoplast and chloroplasts. Our reaction-diffusion model effectively estimates R d , enlightens the dependence of R d estimates on reassimilation and clarifies (dis)advantages of existing methods.

Critical review: incorporating the arrangement of mitochondria and chloroplasts into models of photosynthesis and carbon isotope discrimination

The arrangement of mitochondria and chloroplasts, together with the relative resistances of cell wall and chloroplast, determine the path of diffusion out of the leaf for (photo)respired CO₂. Traditional photosynthesis models have assumed a tight arrangement of chloroplasts packed together against the cell wall with mitochondria located behind the chloroplasts, deep inside the cytosol. Accordingly, all (photo)respired CO₂ must cross the chloroplast before diffusing out of the leaf. Different arrangements have recently been considered, where all or part of the (photo)respired CO₂ diffuses through the cytosol without ever entering the chloroplast. Assumptions about the path for the (photo)respiratory flux are particularly relevant for the calculation of mesophyll conductance (g_m). If (photo)respired CO₂ can diffuse elsewhere besides the chloroplast, apparent g_m is no longer a mere physical resistance but a flux-weighted variable sensitive to the ratio of (photo)respiration to net CO₂ assimilation. We discuss existing photosynthesis models in conjunction with their treatment of the (photo)respiratory flux and present new equations applicable to the generalized case where (photo)respired CO₂ can diffuse both into the chloroplast and through the cytosol. Additionally, we present a new generalized Δ¹³C model that incorporates this dual diffusion pathway. We assess how assumptions about the fate of (photo)respired CO₂ affect the interpretation of photosynthetic data and the challenges it poses for the application of different models.

Embracing 3D Complexity in Leaf Carbon-Water Exchange

Leaves are a nexus for the exchange of water, carbon, and energy between terrestrial plants and the atmosphere. Research in recent decades has highlighted the critical importance of the underlying biophysical and anatomical determinants of CO₂ and H₂O transport, but a quantitative understanding of how detailed 3D leaf anatomy mediates within-leaf transport has been hindered by the lack of a consensus framework for analyzing or simulating transport and its spatial and temporal dynamics realistically, and by the difficulty of measuring within-leaf transport at the appropriate scales. We discuss how recent technological advancements now make a spatially explicit 3D leaf analysis possible, through new imaging and modeling tools that will allow us to address long-standing questions related to plant carbon-water exchange.

Functioning of potassium and magnesium in photosynthesis, photosynthate translocation and photoprotection

Potassium (K) and magnesium (Mg) are mineral nutrients that are required in large quantities by plants. Both elements critically contribute to the process of photosynthesis and the subsequent long-distance transport of photoassimilates. If K or Mg is not present in sufficient quantities in photosynthetic tissues, complex interactions of anatomical, physiological and biochemical responses result in a reduction of photosynthetic carbon assimilation. As a consequence, excessive production of reactive oxygen species causes photo-oxidation of the photosynthetic apparatus and causes an up-regulation of photoprotective mechanisms. In this article, we review the functioning of K and Mg in processes directly or indirectly associated with photosynthesis. Focus is given to chloroplast ultrastructure, light-dependent and -independent reactions of photosynthesis and the diffusion of CO2  - a major substrate for photosynthesis - into chloroplasts. We further emphasize their contribution to phloem-loading and long-distance transport of photoassimilates and to the photoprotection of the photosynthetic apparatus.

Localization of (photo)respiration and CO2 re-assimilation in tomato leaves investigated with a reaction-diffusion model

The rate of photosynthesis depends on the CO₂ partial pressure near Rubisco, C_c, which is commonly calculated by models using the overall mesophyll resistance. Such models do not explain the difference between the CO₂ level in the intercellular air space and C_c mechanistically. This problem can be overcome by reaction-diffusion models for CO₂ transport, production and fixation in leaves. However, most reaction-diffusion models are complex and unattractive for procedures that require a large number of runs, like parameter optimisation. This study provides a simpler reaction-diffusion model. It is parameterized by both leaf physiological and leaf anatomical data. The anatomical data consisted of the thickness of the cell wall, cytosol and stroma, and the area ratios of mesophyll exposed to the intercellular air space to leaf surfaces and exposed chloroplast to exposed mesophyll surfaces. The model was used directly to estimate photosynthetic parameters from a subset of the measured light and CO₂ response curves; the remaining data were used for validation. The model predicted light and CO₂ response curves reasonably well for 15 days old tomato (cv. Admiro) leaves, if (photo)respiratory CO₂ release was assumed to take place in the inner cytosol or in the gaps between the chloroplasts. The model was also used to calculate the fraction of CO₂ produced by (photo)respiration that is re-assimilated in the stroma, and this fraction ranged from 56 to 76%. In future research, the model should be further validated to better understand how the re-assimilation of (photo)respired CO₂ is affected by environmental conditions and physiological parameters.

Simple generalisation of a mesophyll resistance model for various intracellular arrangements of chloroplasts and mitochondria in C3 leaves

The classical definition of mesophyll conductance (g m) represents an apparent parameter (g m,app) as it places (photo)respired CO2 at the same compartment where the carboxylation by Rubisco takes place. Recently, Tholen and co-workers developed a framework, in which g m better describes a physical diffusional parameter (g m,dif). They partitioned mesophyll resistance (r m,dif = 1/g m,dif) into two components, cell wall and plasmalemma resistance (r wp) and chloroplast resistance (r ch), and showed that g m,app is sensitive to the ratio of photorespiratory (F) and respiratory (R d) CO2 release to net CO2 uptake (A): g m,app = g m,dif/[1 + ω(F + R d)/A], where ω is the fraction of r ch in r m,dif. We herein extend the framework further by considering various scenarios for the intracellular arrangement of chloroplasts and mitochondria. We show that the formula of Tholen et al. implies either that mitochondria, where (photo)respired CO2 is released, locate between the plasmalemma and the chloroplast continuum or that CO2 in the cytosol is completely mixed. However, the model of Tholen et al. is still valid if ω is replaced by ω(1-σ), where σ is the fraction of (photo)respired CO2 that experiences r ch (in addition to r wp and stomatal resistance) if this CO2 is to escape from being refixed. Therefore, responses of g m,app to (F + R d)/A lie somewhere between no sensitivity in the classical method (σ =1) and high sensitivity in the model of Tholen et al. (σ =0).