植物生态学报  2016 , 40 (11): 1208-1217 https://doi.org/10.17521/cjpe.2015.0470

Orginal Article

光系统II实际光化学量子效率对光的响应模型的比较

叶子飘1*, 胡文海2, 闫小红23

1井冈山大学数理学院, 江西吉安 343009
2井冈山大学生命科学学院, 江西吉安 343009
3南京林业大学森林资源与环境学院, 南京 210037

Comparison on light-response models of actual photochemical efficiency in photosystem II

YE Zi-Piao1*, HU Wen-Hai2, YAN Xiao-Hong23

1Maths & Physics College, Jinggangshan University, Ji’an, Jiangxi 343009, China;
2School of Life Sciences, Jinggangshan University, Ji’an, Jiangxi 343009, China;
and 3College of Forest Resources and Environment, Nanjing Forestry University, Nanjing 210037, China

通讯作者:  * E-mail: yezp@jgsu.edu.cn

责任编辑:  YE Zi-PiaoHU Wen-HaiYAN Xiao-Hong

收稿日期: 2015-12-22

接受日期:  2016-07-19

网络出版日期:  2016-11-10

版权声明:  2016 植物生态学报编辑部 本文是遵循CCAL协议的开放存取期刊,引用请务必标明出处。

基金资助:  国家自然科学基金(31560069)和江西省自然科学基金(20142BAB20402)

展开

摘要

为了比较光系统II实际光化学量子效率(ΦPSII)对光的响应机理模型(简称机理模型)、负指数模型和指数模型的优缺点, 用LI-6400-40B光合作用测定仪控制CO2浓度和温度, 测量了剑叶金鸡菊(Coreopsis lanceolata)、黄荆(Vitex negundo)和大狼杷草(Bidens frondosa)的电子传递速率(ETR)对光的响应曲线(ETR-I)和ΦPSII对光的响应曲线(ΦPSII-I), 然后用这3个模型分别拟合了这些数据。拟合结果表明: 3个模型都可以较好地拟合这3种植物的ETR-I的响应数据和ΦPSII-I的响应数据, 但由指数模型拟合ETR-IΦPSII-I的响应数据得到相应的饱和光强(PARsat)和光系统II最大光能利用效率(Fv/Fm)之间存在显著差异, 且估算的饱和光强远低于实测值。由机理模型可知, ΦPSII不仅与光强的函数有关, 还与植物的内禀特性有关, 即与天线色素分子的本征光能吸收截面、激子的传递效率、能级的简并度、光化学反应常数、热耗散常数和处于最低激发态的平均寿命等参数有关。此外, 由机理模型还可知, ΦPSII随光强的增加而下降的原因是捕光色素分子的有效光能吸收截面随光强增加而降低。

关键词: 光系统II ; 实际光化学量子效率 ; 电子传递速率 ; 机理模型 ; 天线色素分子

Abstract

Aims The objective of this study was to compare the merits and demerits of three models (i.e., a mechanistic model, a negative exponential model and an exponential model) to simulate the light-response curves of actual photochemical efficiency (ΦPSII-I). Moreover, it was to reveal the mechanism that ΦPSII decreased with light intensity increasing.
Methods The electron transport rate (ETR) and the ΦPSII of Coreopsis lanceolata, Vitex negundo and Bidens frondosa were measured by LI-6400-40B under controlled CO2 concentrations and temperatures, then light-response curves of ETR-I and ΦPSII-I were simulated by a mechanistic model, a negative exponential model and an exponential model, respectively.
Important findings The fitted results showed that ETR-I and ΦPSII-I data of the three plants fit well to the three models. However, the saturation light intensity (PARsat) and maximum efficiency of photosystem II (Fv/Fm) estimated by exponential model for ETR-I and ΦPSII-I were greatly different from the measured data. Moreover, the mechanistic model revealed that the ΦPSII for the three species decreased with increasing I as the effective light energy absorption cross-section of light-harvesting pigments decreased. At the same time, it showed that ΦPSII depended not only on I, but also on eigen-absorption cross-section of light-harvesting pigment, efficiency of exciton, photochemical constant, heat dissipation constant and average lifetime of light-harvesting pigment-protein complex in the lowest state.

Keywords: photosystem II ; actual photochemical efficiency ; electron transport rate ; a mechanistic model ; light- harvesting pigment-protein complex

0

PDF (584KB) 元数据 多维度评价 相关文章 收藏文章

本文引用格式 导出 EndNote Ris Bibtex

叶子飘, 胡文海, 闫小红. 光系统II实际光化学量子效率对光的响应模型的比较. 植物生态学报, 2016, 40(11): 1208-1217 https://doi.org/10.17521/cjpe.2015.0470

YE Zi-Piao, HU Wen-Hai, YAN Xiao-Hong. Comparison on light-response models of actual photochemical efficiency in photosystem II. Chinese Journal of Plant Ecology, 2016, 40(11): 1208-1217 https://doi.org/10.17521/cjpe.2015.0470

光系统II (PSII)的实际光化学量子效率(ΦPSII)为PSII反应中心部分关闭时的实际原初光能捕获效率, 表征作用光下PSII吸收光能用于原初醌电子受体QA还原的实际运行效率, 是一个被广泛应用的荧光参数。由于它与通过PSII反应中心的线性电子传递(LET)有直接的联系, 可快速地确定光强变化时或其他环境条件下PSII反应中心的作用效率(Baker, 2008), 所以, 它又被称为PSII的运行效率(PSII operating efficiency)。因此, 研究它的变化规律具有重要的理论价值和实际意义(Ritchie & Bunthawin, 2010; Takahashi & Badger, 2011; Silsbe & Kromkamp, 2012; Niyogi & Truong, 2013; 安东升等, 2015)。

按照Genty等(1989)关于电子传递速率(ETR)或光合电子流(J)的理论, ETR的大小可以由公式ETR = α' × β' × ΦPSII × I 得到, 式中的ΦPSII就是PSII的实际光化学量子效率, α'、β'、I分别为光能在PSII和PSI的分配系数(无量纲)、叶片的光吸收系数(无量纲)和光强。而ΦPSII = (F′m - Fs)/F′m (Genty et al., 1989; Krall & Edward, 1992), 其中F′m为光下的最大荧光, Fs为光下稳态荧光。实验结果表明: ΦPSII是一个依赖光强的函数, 它随光强的增加而下降(Genty et al., 1989; Krall & Edward, 1992; Ritchie, 2008; Ritchie & Bunthawin, 2010; Silsbe & Kromkamp, 2012; 冯汉青等, 2014; 安东升等, 2015; 林琭等, 2015; Ware et al., 2015)。同时, ΦPSII也可以由几个经验模型直接拟合(Laws et al., 2002; Smyth et al., 2004; Ritchie, 2008; Ritchie & Bunthawin, 2010; Silsbe & Kromkamp, 2012), 但这些经验模型主要应用于藻类和浮游植物(Laws et al., 2002; Smyth et al., 2004; Silsbe & Kromkamp, 2012)和少数几种陆生植物(Ritchie, 2008; Ritchie & Bunthawin, 2010), 其他陆生植物应用较少。已有的研究结果表明: 由指数模型分别拟合ETR对光的响应曲线(ETR-I)和ΦPSII对光的响应曲线(ΦPSII-I)的响应数据得到相应的饱和光强(PARsat)和植物经过暗适应的PSII最大光能利用效率Fv/Fm (Fv为可变荧光, Fm为最大荧光产量)差异较大, 且所估算的饱和光强与实际测量值差异较大(Laws et al., 2002; Smyth et al., 2004; Silsbe & Kromkamp, 2012); 负指数模型拟合ETR-IΦPSII-I的响应数据得到的饱和光强与实测值符合, 但由此模型得到光强为0时的最大光化学量子效率(ΦPSIImax)之间存在较大的差异(Ritchie, 2008)。由此可知, 由负指数和指数模型拟合ETR-IΦPSII-I的响应曲线分别得到的ΦPSIImax之间和Fv/Fm之间存在明显差异, 但尚不清楚产生这种现象的原因。此外, 这两个模型也不能解释为什么植物的ΦPSII会随光强的增加而下降。 正是由于存在这两个问题, 所以限制了这些模型的 应用。

而由ΦPSII = (F′m - Fs)/F′m的定义可知, ΦPSII与光合作用原初反应中的PSII的开闭和天线色素分子发射的荧光过程紧密相连。首先, 天线色素分子吸收光能后将由基态跃迁到高能态, 然后这些激发能将分配到光化学反应、热耗散和发射叶绿素荧光这3条主要途径且它们彼此竞争(Oxborough, 2004; Baker, 2008), 这种分配过程直接影响到随后的光合电子传递与同化力的形成。而原初反应是一个纯粹的物理过程, 涉及天线色素分子的光能吸收、量子态的改变、激子的共振传递和退激发等, 这些过程又是由天线色素分子的内禀特性所决定的(Govindjee, 2002; Baker, 2008; Richter et al., 2008; Panitchayangkoon et al., 2010; Sarovar et al., 2010)。此外, 光合作用也与光合能力的形成(如ATP和NADPH的形成) (Govindjee, 2002; Nelson & Yocum, 2006; Ye et al., 2013a, 2013b)和与碳同化相连的酶动力过程有关(Lin et al., 2015)。由此推测: (1) ΦPSII的大小应该与植物的天线色素分子、PSII的内禀特性和酶动力学过程密切相关; (2) ΦPSII随光强的增加而下降可能与捕光色素分子的有效光能有关; (3)由负指数和指数模型拟合ETR-IΦPSII-I的响应曲线分别得到的ΦPSIImax之间和Fv/Fm之间明显存在差异可能是由于LI-6400-40B荧光仪设定的叶片光能吸收系数为0.84和光能在PSII和PSI之间的分配比例为0.5并非适合所有的植物。

针对以上的3个推测, 在本研究中, 首先, 我们基于光合电子流对光的响应机理模型(Ye et al., 2013a, 2013b)推导出ΦPSII对光的响应模型, 确定其变化规律与PSII天线色素分子的内禀特性之间的关系; 其次, 利用ΦPSII对光响应的机理模型、负指数模型和指数模型分别拟合剑叶金鸡菊(Coreopsis lanceolata)、黄荆(Vitex negundo)和大狼杷草(Bidens frondosa)在正常条件下ΦPSII对光的响应曲线, 判断所推导的ΦPSII对光的响应模型是否合理和有效, 并与负指数模型和指数模型的拟合结果进行比较, 以确定这些模型的优缺点; 再次, 用所测量的ETR对光的响应数据验证我们的推测, 即植物的ΦPSII随光强的增加而下降可能与其天线色素分子的有效光能吸收截面特性密切相关; 最后, 利用所拟合的结果确定LI-6400-40B荧光仪设定的叶片光能吸收系数为0.84和光能在PSII和PSI的分配比例为0.5是否适合剑叶金鸡菊、黄荆和大狼杷草。因此, 这些问题的探讨将为研究植物的ΦPSII对光的响应以及ΦPSII随光强的增加而下降的机理提供理论基础, 同时也为研究植物的ETR-IΦPSII-I提供数学工具。

1 材料和方法

1.1 材料

供试材料为自然生长在江西省吉安市井冈山大学校园内的剑叶金鸡菊、黄荆和大狼杷草。

1.2 测定及计算方法

在2013年6月连续的几个晴天内, 在9:00-11:30和14:30-16:30, 随机选5株长势较为一致的健壮植株用LI-6400-40B (LI-COR, Lincoln, USA)测定其ETR对光的响应曲线和其他参数。在测定过程中, 气室CO2浓度控制为(390 ± 2) μmol·mol-1, 气室温度设置为(35 ± 1) ℃, 气室的空气相对湿度设置为50%-70%。光合有效辐射(PAR)设置为2000、1800、1600、1400、1200、1000、800、600、400、200、100、50和0 μmol·m-2·s-1共13个水平, 并分别于各光强下测定剑叶金鸡菊、黄荆和大狼杷草的光合电子流(J)或ETR和其他参数, 每个测定点的响应时间在120 s以上。然后分别用植物光合电子流对光响应的机理模型、负指数模型和指数模型拟合这3种植物的光合电子流对光的响应曲线, 得到其光响应曲线的初始斜率(α)或光强为0时的最大光化学量子效率(ΦPSIImax)或植物经过暗适应的PSII最大光化学量子产量(Fv/Fm)、饱和光强(PARsat)、最大光合电子流(Jmax)或最大电子传递速率(ETRmax), 并用相关的参数估算天线色素分子的有效光能吸收截面对光的响应曲线等。

1.3 机理模型的推导

植物叶片的天线色素分子吸收的光能将在光化学反应、热耗散和荧光发射等3条途径上竞争分配(Oxborough, 2004; Baker, 2008), Ye等(2013a, 2013b)在此基础上构建了一个ETR对光的响应机理模型。其数学表达式为:

式中的符号如φN0SgigkkPkDξ1ξ2ξ3σikτ等皆为描述天线色素分子物理特性的参数, 详见相关文献(Ye et al., 2013a, 2013b; 叶子飘等, 2014), α'为叶片的光吸收系数, β'为光能在PSII和PSI两个光系统中的分配比例。其中α'一般取值为0.84 (Ehleringer & Pearcy, 1983; Krall & Edwards, 1992), β'一般取值为0.5 (Ehleringer & Pearcy, 1983; Krall & Edwards, 1992; Major & Dunton, 2002)。

由于在一定的环境下植物叶片天线色素分子的物理参数是确定的, 可以假设: α = α'β'N0φσik/S (μmol electrons·(μmol photons)-1), β = ((1-gi/gk)σikτ)/ (ξ3 + (ξ1kP + ξ2kD)τ) (m2·s·(μmol photons)-1), γ = ((1+ gi/gk)σikτ)/(ξ3 + (ξ1kP + ξ2kD)τ) (m2·s·(μmol photons)-1), 那么公式(1)可以简化为:

由公式(2)可以得到植物叶片的饱和光强(PARsat)和最大电子传递速率(ETRmax)分别为:

此外, 按照Krall和Edward (1992)给出的光合电子流对光的响应公式:

ETR = α'β 'ΦPSIII (3)

对比公式(1)和公式(3), 可以得到ΦPSII与光强之间的关系式为:

由公式(4)可知, ΦPSII不仅与光强存在非线性关系, 还与天线色素分子的本征光能吸收截面(σik)、天线色素分子数(N0)、激子的传递效率(φ)、激子处于最低激发态的平均寿命(τ)、光化学反应速率(kP)和热耗散速率(kD)、天线色素分子处于基态i的简并度(gi)和激发态k的简并度(gk)等参数存在密切的关系。

假设ω = N0φσik/S (μmol electrons·(μmol photons)-1), 那么公式(4)可以简化为:

这里的ωΦPSII-I曲线的初始斜率。由公式(5)可知, ΦPSII随光强的增加而非线性下降。对比公式(2)和公式(5)可知, α/ω = α'β', 由此也可以确定Dual-PAM- 100所推荐的α'β'的值分别为0.84和0.5是否合理。

由公式(1)还我们还可以得到天线色素分子有效光能吸收截面(σ'ik)与光强之间的关系(Ye et al., 2013a, 2013b)为:

由公式(6)可知, σ'ik随光强的增加呈非线性下降。对比公式(6)和公式(4)可知, 对于一定植物在一定环境条件下, ΦPSIIσ'ik成正比, 即:

ΦPSII = N0φσ'ik/S (7)

1.4 经验模型的介绍

现有拟合ΦPSII-I曲线的经验模型有负指数模型(Ritchie, 2008; Ritchie et al., 2010)、指数模型(Webb et al., 1974; Smyth et al., 2004; Silsbe & Kromkamp, 2012)和正切函数模型(Webb et al., 1974; Smyth et al., 2004), 其中以负指数模型和指数模型应用最为广泛。负指数模型(Ritchie, 2008; Ritchie et al., 2010)的数学表达式是:

ΦPSII = ΦPSIImax e-kwI (8)

式中的ΦPSIImax是光强为0时的最大光化学量子效率, kw是一个常数, I为光强。与公式(8)相对应的ETR对光的响应模型的数学表达式为:

ETR = α'β'IΦPSIImax e-kwI (9)

由公式(9)可以得到植物叶片的饱和光强和最大电子传递速率分别为: PARsat = 1/kwETRmax =α'β'PARsatΦPSIImaxe-1

指数模型(Webb et al., 1974; Smyth et al., 2004; Silsbe & Kromkamp, 2012)的数学表达式是:

式中的PARsatFv/FmI与以上的定义相同。由公式(10)可以直接估算出植物叶片的饱和光强(PARsat)。与公式(10)相对应的光合电子流对光的响应模型的数学表达式为:

ETR = α'β'(Fv/Fm) PARsat [1-exp(-I/PARsat)] (11)

此外, 由公式(11)也可以直接估算植物叶片的PARsat, 最大电子传递速率则为: ETRmax = α'β'(Fv/Fm) PARsat [1-exp(-1)]。

1.5 数据处理

ETR-I曲线和ΦPSII-I曲线的实验数据处理用SPSS 12.5; 用Origin 7.0作图, 然后用Adobe Illustrator CS5进行图形处理, 采用SPSS 12.5作方差差异显著性分析、ANOVA检验和配对样本t检验, 用DPS 7.05软件作模型分析。本文将用AIC信息准则并结合确定系数(R2)对3个模型的拟合结果进行比较, 并判断3个模型的优缺点。其中AIC的计算参考Akaike (1973)的方法, R2则是直接由SPSS 12.5在拟合数据后给出。AIC值越小, 说明拟合的模型既精度高又简洁; R2是衡量所建立模型效果好坏的指标之一, 其值越大表明模型拟合得越好。此外, 为了比较α/ω与0.42 (α'β' = 0.42)之间是否存在显著性差异, 本文对此进行了配对样本t检验, 显著性水平设定为α = 0.05。文中的图和表中的部分数据均用平均值±标准误差表示。

1.6 叶绿素含量的测定

取测定完光合电子流和净光合速率对光的响应曲线的成熟叶片, 在不同叶片上取直径0.7 cm的叶圆片30片, 分3组, 分别浸入10 mL 80%丙酮的提取液中, 盖紧样品瓶盖并密封, 避光低温放置48 h至叶圆片无色。然后参照Arnon (1949)的方法测量叶绿素的含量, 重复5次取平均值。

2 实验结果和分析

2.1 3种植物叶片的电子传递速率对光的响应曲线(ETR-I响应曲线)

植物的光合电子流对光的响应曲线一般可以分为3个阶段(Ralph & Gademann, 2005): 第一阶段是光合电子流随光强的增加而几乎线性增加; 第二阶段是光合电子流随光强的增加缓慢增加到一个平台; 第三阶段是光合电子流随光强的增加而下降, 出现光PSII的动力学下调(每个光强点的测量时间小于60 s)或光抑制现象(每个光强点的测量时间大于180 s)。由图1可知, 剑叶金鸡菊和大狼杷草的光合电子流随光强的增加而下降, 而黄荆的光合电子流随光强的增加一直在增加, 尚未出现下降。

图1   剑叶金鸡菊(A)、黄荆(B)和大狼杷草(C)的电子传递速率对光的响应曲线(ETR-I)(平均值±标准误差)。公式(2)为ETR = α (1-β I) / (1 + γI) I, 其中αETR-I的初始斜率, β为光抑制项, γ为光饱和项, I为光强; 公式(9)为ETR =α'β'IΦPSIImax e-kwI, 其中α'为叶片光能吸收系数, β'为光能在2个光系统的分配比例, ΦPSIImax为光强为0时的最大光化学量子效率, kw是一个常数; 公式(11)为ETR = α'β'(Fv/Fm) PARsat [1-exp (-I/PARsat)], 其中Fv/Fm为光系统II的最大光量子效率, PARsat为饱和光强。

Fig. 1   Light-response curves of electron transport rate (ETR- I) for Coreopsis lanceolata (A), Vitex negundo (B) and Bidens frondosa (C) (mean ± SE). Equation (2) is ETR =α (1-β I) / (1 + γI) I. Here α is the initial slope of ETR-I, β is the photoinhibition term, γ is the saturation term, I is light intensity; Equation (9) is ETR =α'β'IΦPSIImaxe-kwI. Here α' is the light absorption coefficient of leaf, β' is the light energy distribution coefficient between photosystem II and photosystem I, ΦPSIImax is the maximum photochemical quantum efficiency while I = 0, kw is a constant; Equation (11) is ETR = α'β'(Fv/Fm) PARsat [1-exp(-I/PARsat)]. Here Fv/Fm is the maximum quantum of photosystem II, PARsat is the saturation irradiance.

图1表1可知, 由公式(2)和(9)得到剑叶金鸡菊、黄荆和大狼杷草的最大光合电子流和饱和光强与实测值高度符合, 而由公式(11)得到的ETRmax分别为156.5、155.99和86.67 μmol·m-2·s-1, 测量值则分别约为230.4、228.8和135.8 μmol·m-2·s-1, 估算值分别约为实测值的67.8%、68.2%和63.8%; 由公式(11)得到的PARsat分别为617.83、633.1和395.2 μmol·m-2·s-1, 测量值分别约为1600、1750和1400 μmol·m-2·s-1, 估算值尚未达到实测值的一半(表1)。此外, 由表1的数据可知, 公式(2)、(9)和(11)拟合3种植物的光合电子流对光的响应曲线的R2都在0.994以上。如果用确定系数的大小表示模型的拟合与实际测量曲线的符合度, 那么(2)式的拟合效果最好, (11)式的拟合效果最差。如果以AIC值的大小作为衡量标准, 则以公式(9)的拟合效果最好, 公式(11)的拟合效果最差(表1)。

表1   三种模型对三种植物光合电子流对光曲线的拟合参数

Table 1   The measured data and fitted values for Coreopsis lanceolata, Vitex negundo and Bidens frondosa using three models

光合参数
Photosynthetic
parameter
剑叶金鸡菊 C. lanceolata黄荆 V. negundo大狼杷草 B. frondosa
公式(2)
Equation (2)
公式(9)
Equation (9)
公式(11)
Equation (11)
测量值
Measured value
公式(2)
Equation (2)
公式(9)
Equation (9)
公式(11)
Equation (11)
测量值
Measured value
公式(2)
Equation (2)
公式(9)
Equation (9)
公式(11)
Equation (11)
测量值
Measured value
ETR-I的初始斜率
Initial slope of ETR-I (μmol·μmol -1)
0.347ΦPSIImax
(0.859)
Fv/Fm
(0.954)
-0.352ΦPSIImax
(0.834)
Fv/Fm
(0.928)
-0.315ΦPSIImax
(0.675)
Fv/Fm
(0.820)
-
最大电子传递速率
Maximu electron
transport rate
(μmol·m-2·s-1)
232.3231.3156.5230.4230.2229.8155.99≈228.8136.7137.2886.67≈135.8
饱和光强
Saturation irradiance
(μmol·m-2·s-1)
1 641.71 742.2617.83≈1 6001 740.91 782.5633.1≈1 7501 394.51 326.3395.2≈1 400
叶绿素含量
Chlorophyll content
(mg·m-2)
253.68268.18286.33
本征截面
Eign-absorption cross section
(10-21 m2)
5.41---5.19---4.35---
决定系数
Determined coefficient
1.0000.9980.995-0.9990.9990.997-0.9980.9970.995-
AIC信息准则
Akaike’s information criterion
7.407.349.22-8.347.488.78-7.747.278.46-

Equation (2), Equation (9), Equation (11) see Fig. 1.公式(2)、(9)、(11)同图1

新窗口打开

2.2 三种植物叶片的PSII实际光化学量子效率对光的响应曲线(ΦPSII-I)

图2可知, 3种植物的ΦPSII随光强的增加而非线性下降, 且剑叶金鸡菊和黄荆的下降速度略有差异, 但大狼杷草的ΦPSII随光强的增加而下降的速度最快。由图2还可知, 3个模型都能较好地拟合剑叶金鸡菊、黄荆和大狼杷草的ΦPSII-I曲线, R2超过0.984。此外, 由表2的数据可知, 如果仅从R2值的大小衡量3个模型的优劣, 那么以机理模型的拟合效果最好, 指数模型最差; 如果以AIC值的大小衡量, 那么以负指数模型的拟合效果最好, 拟合效果最差的是指数模型。

图2   剑叶金鸡菊(A)、黄荆(B)和大狼杷草(C)的实际光量子效率对光的响应曲线(ΦPSII-I)(平均值±标准误差)。公式(5)为ΦPSII = ω (1-βI) / (1 + γI), 其中ωΦPSII-I的初始斜率; 公式(8)为ΦPSII = ΦPSIImax e-kwI, 公式(10)为ΦPSII = (Fv/Fm) PARsat / I (1-exp(-I/PARsat))。βγFv/FmΦPSIImaxkwPARsat的定义同图1

Fig. 2   Light-response curves of ΦPSII for Coreopsis lanceolata (A), Vitex negundo (B) and Bideas frondosa (C)(mean ± SE). Equation (5) is ΦPSII = ω (1-βI) / (1 + γI). Here ω is the initial slope of ΦPSII-I; Equation (8) is ΦPSII = ΦPSIImax e-kwI; Equation (10) is ΦPSII = (Fv/Fm) PARsat / I (1-exp(-I/PARsat)). The definition of β, γ, Fv/Fm, ΦPSIImax, kw and PARsat see Fig. 1.

表2   三种模型对三种植物叶片的ΦPSII-I曲线的拟合参数

Table 2   The measured data and fitted values for Coreopsis lanceolata, Vitex negundo and Bidens frondosa using three models

光合参数
Photosynthetic
parameter
剑叶金鸡菊 C. lanceolata黄荆 V. negundo大狼杷草 B. frondosa
公式(5)
Equation (5)
公式(8)
Equation (8)
公式(10)
Equation (10)
测量值
Measured value
公式(5)
Equation (5)
公式(8)
Equation (8)
公式(10)
Equation (10)
测量值
Measured value
公式(5)
Equation (5)
公式(8)
Equation (8)
公式(10)
Equation (10)
测量值
Measured value
ΦPSII-I的初始斜率
Initial slope of ΦPSII-I
(μmol·μmol -1)
0.763ΦPSIImax
(0.771)
Fv/Fm
(0.779)
-0.752ΦPSIImax
(0.759)
Fv/Fm
(0.767)
-0.665ΦPSIImax
(0.678)
Fv/Fm
(0.674)
-
饱和光强
Saturation irradiance
(μmol·m-2·s-1)
1 642.91 744.9824.95≈1 6001 712.71 763.3830.33≈1 7501 406.71 331.2448.50≈1 400
决定系数
Determined coefficient
0.9950.9920.984-0.9960.9940.986-0.9960.9950.991-
AIC信息准则
Akaike’s information criterion
-2.49-2.75-2.03--2.79-3.04-2.20--2.51-3.06-2.46-

Equation (5), Equation (8), Equation (10) see Fig. 2.公式(5)、(8)、(10)同图2

新窗口打开

除了可用公式(9)和(11)拟合ETR-I曲线分别得到这3种植物的ΦPSIImaxFv/Fm的值外, 也可以由公式(8)和(10)拟合ΦPSII-I曲线得到这3种植物的ΦPSIImaxFv/Fm的值(表2)。此外, 由公式(11)和(10)分别拟合ETR-I曲线和ΦPSII-I曲线得到的PARsat有较大的差异。比如, 由公式(11)拟合剑叶金鸡菊的ETR-I曲线得到的PARsat为617.83 μmol·m-2·s-1, 而由公式(10)拟合剑叶金鸡菊的ΦPSII-I曲线得到的PARsat为824.95 μmol·m-2·s-1, 后者是前者的1.33倍。

2.3 3个模型拟合3种植物的ΦPSIImaxFv/Fmα/ω的结果比较

由(8)和(9)式可以估算植物的ΦPSIImax, 由(10)和(11)式可以估算植物的Fv/Fm, 但对比表1表2中的数据可知, 由(8)和(9)式估算同种植物的ΦPSIImax得到的值并不相同, 由(10)和(11)式估算同种植物的Fv/Fm得到的值也不相同。例如, 用(9)式得到剑叶金鸡菊的ΦPSIImax为0.859, 而由(8)式得到剑叶金鸡菊的ΦPSIImax为0.771, 前者是后者的1.114倍。为了确定这些模型估算的ΦPSIImax之间和Fv/Fm之间是否存在差异, 对拟合结果进行了ANOVA检验, 发现用(8)和(9)式分别得到的剑叶金鸡菊和黄荆的ΦPSIImax值之间存在显著差异(表3); 同时也对由(10)和(11)式估算的植物Fv/Fm的值进行了ANOVA检验, 发现用(10)和(11)式分别得到的剑叶金鸡菊(和黄荆)的Fv/Fm值之间也存在显著差异(表3)。此外, 由配对样本t检验得知: 由(2)式和(9)式或(5)式和(8)式得到3种植物的PARsat与实测值之间不存在显著差异, 而与(10)或(11)式得到3种植物的PARsat与实测值之间则存在显著差异(数据没有给出)。

表3   三种模型对三种植物叶片的ΦPSII-I曲线的拟合参数(平均值±标准误差)

Table 3   The fitted values of ΦPSII-I curves for Coreopsis lanceolata, Vitex negundo and Bidens frondosa using three models (mean ± SE)

剑叶金鸡菊 C. lanceolata黄荆 V. negundo大狼杷草 B. frondosa
ΦPSIImaxFv/FmΦPSIImaxFv/FmΦPSIImaxFv/Fm
公式(8) Equation (8)0.859 ± 0.011a-0.834 ± 0.005a-0.678 ± 0.059a-
公式(9) Equation (9)0.771 ± 0.004b-0.759 ± 0.006b-0.675 ± 0.051a-
公式(10) Equation (10)-0.954 ± 0.019a-0.928 ± 0.005a-0.820 ± 0.079a
公式(11) Equation (11)-0.779 ± 0.004b-0.767 ± 0.006b-0.674 ± 0.029a
α/ω0.455 ± 0.007*0.468 ± 0.008*0.474 ± 0.023ns

Equation (8), Equation (9), Equation (10), Equation (11) see Fig. 1 and Fig. 2. Different letters (a and b) indicate statistical significance at p < 0.05 between the same model. *, significant difference (p < 0.05) between α/ω and 0.42 (α'β' = 0.42).公式(8)、(9)、(10)、(11)同图1图2。不同小写字母(a和b)表示同一参数在相同模型间存在显著性差异(p < 0.05)。*, α/ω和0.42 (α'β' = 0.42)之间存在显著性差异(p < 0.05)。

新窗口打开

此外, 对α/ω和0.42 (α'β' = 0.42)之间进行了配对样本t检验, 发现剑叶金鸡菊和黄荆的α/ω和0.42存在显著差异(表3)。由此可知剑叶金鸡菊、黄荆和大狼杷草的α/ω值分别为0.455、0.468和0.474, 这3个值都比LI-6400-40B给出α'β'的推荐值0.42大。

2.4 三种植物天线色素分子的有效光能吸收截面对光的响应曲线(σ'ik-I)

天线色素分子的有效光能吸收截面(σ'ik)表示其在不同光强下天线色素分子捕获光子能力的大小。由图3可知, 剑叶金鸡菊、黄荆和大狼杷草的天线色素分子有效光能吸收截面随光强的增加而下降, 表明这3种植物随光强的增加其捕获光子的能力都在逐渐变小。由图3表1的数据可知, 剑叶金鸡菊和黄荆天线色素分子在高光强(2000 μmol·m-2·s-1)时, 它们的有效光能吸收截面几乎相同, 约为1.85 × 10-21 m2, 表明这两种植物在高光强时依然可以吸收较多的光能。而大狼杷草天线色素分子的有效光能吸收截面在高光强(2000 μmol·m-2·s-1)时约为1.02 × 10-21 m2, 表明在高光强时, 它吸收的光能比剑叶金鸡菊和黄荆少81%左右。

图3   剑叶金鸡菊(A)、黄荆(B)和大狼杷草(C)的有效光能吸收截面对光的响应曲线(σ'ik-I)(平均值±标准误差)。

Fig. 3   Light-response curves of effective light energy absorption cross-section (σ'ik-I) for Coreopsis lanceolata (A), Vitex negundo (B) and Bidens frondosa (C)(mean ± SE).

此外, 对比图3图2可知, 剑叶金鸡菊、黄荆和大狼杷草的ΦPSII-I曲线与它们的σ'ik-I曲线非常相似。

3 讨论和结论

首先, 基于光合电子流对光响应的机理模型构建了一个ΦPSII对光的响应模型。由ΦPSII-I模型可知, ΦPSII不仅与光强有关, 而且还与植物天线色素的内禀特性存在非线性的关系, 即ΦPSII的大小除了与叶绿素的含量有关外, 还与天线色素分子的本征光能吸收截面和激子的传递效率等有关。

其次, 指数模型和负指数模型主要用于藻类(Laws et al., 2002; Ritchie, 2008)、浮游植物(Smyth et al., 2004; Silsbe & Kromkamp, 2012)和少数几种陆生植物(Ritchie, 2008; Ritchie & Bunthawin, 2010), 用于陆地高等植物较少。本研究中, 我们分别用所构建的ΦPSII对光响应的机理模型, 同时用负指数模型和指数模型分别拟合了剑叶金鸡菊、黄荆和大狼杷草的ΦPSII-I响应曲线, 发现这3个模型都可以较好地拟合这3种植物的ΦPSII-I响应曲线, R2最低为0.984, 但由指数模型拟合得到这3种植物的PARsat与实测值相差较大, 约为实测值的一半甚至更少(表2)。由配对样本t检验也可以得知: 由公式(10)或(11)得到3种植物的PARsat与实测值之间存在显著差异, 而由公式(3)和(5)或公式(8)和(9)得到3种植物的PARsat与实测值之间不存在显著差异(数据并没有给出), 这表明由机理模型和指数模型拟合得到的PARsat与实测值相符。此外, 我们用这3个ΦPSII-I响应模型所对应的ETR-I响应模型分别拟合了剑叶金鸡菊、黄荆和大狼杷草的ETR-I响应曲线。如果仅用R2来判断3个模型拟合的优劣, 那么以ETR-IΦPSII-I机理模型为最佳, 指数模型最差; 如果以AIC作为判断标准, 则以负指数模型为最佳, 指数模型为最差。而由图1图2表1表2的拟合结果可知, 只有ETR-I的机理模型和负指数ETR-I的响应模型拟合得到的ETRmaxPARsat与实测值高度符合, 而由指数ETR-I的响应模型得到的ETRmaxPARsat远低于实测值。这些研究结果与前人的结果相似(Laws et al., 2002; Smyth et al., 2004; Silsbe & Kromkamp, 2012)。由此可知, 与机理模型和负指数模型的拟合结果相比, 指数模型不适于研究藻类、浮游植物和陆地植物的ΦPSII-I响应和ETR-I响应问题。

再次, 大量的实验数据表明: ΦPSII随光强的增加而下降是一个普遍现象(Ritchie, 2008; Ritchie & Bunthawin, 2010; Silsbe & Kromkamp, 2012; 冯汉青等, 2014; 安东升等, 2015; 林琭等, 2015; 师生波等, 2015; Ware et al., 2015)。本研究中剑叶金鸡菊、黄荆和大狼杷草的ΦPSII也随光强的增加而下降(图2), 与前人的结果相似。此外, 由图1图2可知, 负指数模型和指数模型可以拟合三种植物的ETR-IΦPSII-I响应曲线, 却不能解释ΦPSII随光强的增加而下降的原因, 对应的电子传递速率对光的响应模型即公式(9)和(11)也不能解释在ETR-I响应曲线中出现的光系统II的动力学下调(White & Critchley, 1999; Ralph & Gademann, 2008)或光抑制(Vijayan & Browse, 2002; Goh et al., 2012; Rochaix, 2014; Murchie et al., 2015)的原因。然而由公式(7)ΦPSIIσik之间的关系可知, 剑叶金鸡菊、黄荆和大狼杷草的ΦPSII随光强的增加而非线性下降的原因来自它们的天线色素分子的有效光能吸收截面随光强的增加而下降(图3), 即天线色素分子吸收光子的能力随光强的增加而减少, 由此3种植物用于光合作用的光能就减少了, 从而导致PSII的实际光化学量子效率降低。

最后, 用公式(8)和(9)分别得到的剑叶金鸡菊和黄荆的ΦPSIImax值之间存在显著差异(表3); 用公式(10)和(11)估算得到剑叶金鸡菊和黄荆的Fv/Fm值之间也存在显著差异(表3)。此外, 剑叶金鸡菊、黄荆和大狼杷草的α/ω值分别为0.455、0.468和0.474, 其中剑叶金鸡菊和黄荆α/ω的值与0.42 (α'β' = 0.42)存在显著差异(表3), 也即这3个值都比LI- 6400-40B给出α'β'的推荐值0.42要大。其原因可能是剑叶金鸡菊、黄荆和大狼杷草叶片的光能吸收系数(β')要大于0.84或光能分配到PSII的系数要大于0.5, 也即由LI-6400-40B给出α'β'的推荐值不适合这3种植物。如果能用球面积分仪测量植物叶片的光吸收系数(β'), 就可以由公式(2)和(5)或由公式(8)和(9)计算叶片吸收的光能是否平均地分配到光系统II和光系统I, 即可以确定α'是否恒等于0.5。这为研究植物光能是否平均分配到光系统II和光系统I提供了一个有效工具。

总之, 本文所构建的ΦPSII-I模型可以很好地拟合3种植物PSII的实际光化学量子效率对光的响应曲线, 并可以很好地说明ΦPSII随光强的增加而非线性下降的原因是捕光色素分子的有效光能吸收截面随光强的增加而下降。同时, 其对应的ETR-I响应模型可以拟合不同植物的电子传递速率对光的响应曲线, 并可以合理地解释植物的PSII动力学下调或光抑制等生理现象。此外, 当ΦPSII-I模型和ETR-I机理模型相结合, 并结合球面积分仪测量植物叶片的光吸收系数时, 就可以研究植物叶片吸收的光能是否平均地分配到PSII和PSI。

The authors have declared that no competing interests exist.

作者声明没有竞争性利益冲突.


参考文献

[1] Akaike H (1973).

A new look at the statistical model identifica- tion

.IEEE Translations Automatic Control, 19, 716-723.

[2] An DS, Cao J, Huang XH, Zhou J, Dou MA (2015).

Application of Lake-model based indices from chlorophyll fluorescence on sugarcane seedling drought resistance study

.Chinese Journal of Plant Ecology, 39, 398-406. (in Chinese with English abstract)
[安东升, 曹娟, 黄小华, 周娟, 窦美安 (2015).

基于Lake模型的叶绿素荧光参数在甘蔗苗期抗旱性研究中的应用

. 植物生态学报, 39, 398-406.]

[本文引用: 3]     

[3] Arnon DJ (1949).

Copper enzymes in isolated chloroplasts. Polyphenoloxidase in Beta vulgaris

. Plant Physiology, 24, 1-15.

[4] Baker NR (2008).

Chlorophyll fluorescence: A probe of photosynthesis in vivo

. Annual of Review Plant Biology, 59, 89-113.

[本文引用: 4]     

[5] Ehleringer J, Pearcy RW (1983).

Variation in quantum yield for CO2 uptake among C3 and C4 plants

.Plant Physiology, 73, 555-559.

[本文引用: 2]     

[6] Feng HQ, Jiao QS, Tian WY, Sun K, Jia LY (2014).

Effects of extracellular ATP on the characteristics of photochemical reaction in bean (Phaseolus vulgaris) leaves under differ- ent light intensities

. Chinese Journal of Plant Ecology, 38, 1117-1123. (in Chinese with English abstract)
[冯汉青, 焦青松, 田武英, 孙坤, 贾凌云 (2014).

不同光强下细胞外三磷酸腺苷对菜豆叶片光化学反应特性的影响

. 植物生态学报, 38, 1117-1123.]

[本文引用: 2]     

[7] Genty B, Briantais JM, Baker NR (1989).

The relationship between the quantum yield of photosynthetic electron transport and quenching of chlorophyll fluorescence

.Biochimica et Biophysica Acta, 990, 87-92.

[本文引用: 2]     

[8] Goh CH, Ko SM, Koh S, Kim YJ, Bae HJ (2012).

Photosyn- thesis and environments: Photoinhibition and repair mech- anisms in plants

.Journal of Plant Biology, 55, 93-101.

[本文引用: 1]     

[9] Govindjee (2002).

A role for a light-harvesting antenna complex of photosystem II in photoprotection

.The Plant Cell, 14, 1663-1668.

[本文引用: 2]     

[10] Krall JP, Edward GE (1992).

Relationship between photo- system II activity and CO2 fixation in leaves

.Physiologia Plantarum, 86,180-187.

[本文引用: 5]     

[11] Laws E, Sakshaug E, Babin M, Dandonneau Y, Falkowski P, Geider R, Legendre L, Morel A, Sondergaard M, Takahashi M, Williams PJ (2002).

Photosynthesis and primary productivity in marine ecosystems: Practical aspects and application of techniques

.Joint Global Ocean Flux Study, 36, 1-77.

[本文引用: 5]     

[12] Lin L, Tang Y, Zhang JT, Yan WL, Xiao JH, Ding C, Dong C, Ji ZS (2015).

Effects of different water potentials on leaf gas exchange and chlorophyll fluorescence parameters of cucumber during post-flowering growth stage

.Chinese Journal of Applied Ecology, 26, 2030-2040. (in Chinese with English abstract)
[林琭, 汤昀, 张纪涛, 闫万丽, 肖建红, 丁超, 董川, 籍增顺 (2015).

不同水势对黄瓜花后叶片气体交换及叶绿素荧光参数的影响

. 应用生态学报, 26, 2030-2040.]

[本文引用: 3]     

[13] Lin MZ, Wang ZW, He LC, Xu K, Cheng DL, Wang GX (2015).

Plant photosynthesis-irradiance curve response to pollution show non-competitive inhibited Michaelis kinet- ics

.PLOS ONE, 10, e0142712. doi: 10.1371/journal. pone. 0142712.

[14] Major KM, Dunton KH (2002).

Variations in light-harvesting characteristics of the seagrass,Thalassia testudinum: Evidence for photoacclimation

. Journal of Experimental Marine Biology and Ecology, 275, 173-189.

[本文引用: 1]     

[15] Murchie EH, Ali A, Herman T (2015).

Photoprotection as a trait for rice yield improvement: Status and prospects

.Rice, 8, 31-40.

[本文引用: 1]     

[16] Nelson N, Yocum CF (2006).

Structure of function of photo- systems I and II

.Annual Review of Plant Biology, 57, 521-565.

[本文引用: 1]     

[17] Niyogi KK, Truong TB (2013).

Evolution of flexible non- photochemical quenching mechanisms that regulate light harvesting in oxygenic photosynthesis

.Current Opinion in Plant Biology, 16, 307-314.

[本文引用: 1]     

[18] Oxborough K (2004).

Imaging of chlorophyll a fluorescence: Theoretical and practical aspects of an emerging technique for the monitoring of photosynthetic performance

.Journal of Experimental Botany, 55, 1195-1205.

[本文引用: 2]     

[19] Panitchayangkoon G, Hayes D, Fransted KA, Caram JR, Harel E, Wen J, Blankenship RE, Engel GS (2010).

Long-lived quantum coherence in photosynthetic complexes at physiological temperature

.Proceedings of the National Academy Sciences of the United States of America, 107, 12766-12770.

[本文引用: 1]     

[20] Ralph PJ, Gademann R (2005).

Rapid light curves: A powerful tool to assess photosynthetic activity

.Aquatic Botany, 82, 222-237.

[本文引用: 2]     

[21] Richter M, Renger T, Knorr A (2008).

A Block equation approach to intensity dependent optical spectra of light harvesting complex II. Excitation dependence of light harvesting complex II pump-probe spectra

.Photosynthesis Research, 95, 119-127.

[本文引用: 1]     

[22] Ritchie RJ (2008).

Fitting light saturation curves measured using modulated fluorometry

.Photosynthesis Research, 96, 201-215.

[本文引用: 9]     

[23] Ritchie RJ, Bunthawin S (2010).

The use of pulse amplitude modulation (PAM) fluorometry to measure photosynthesis in a cam orchid,Dendrobium spp.(D. cv. Viravuth Pink)

. International Journal of Plant Sciences, 171, 575-585.

[本文引用: 8]     

[24] Rochaix JD (2014).

Regulation and dynamics of the light- harvesting system

. Annual of Review Plant Biology, 65, 287-309.

[本文引用: 1]     

[25] Sarovar M, Ishizaki A, Fleming G, Whaley B (2010).

Quantum entanglement in photosynthetic light-harvesting com- plexes

.Nature Physics, 6, 462-467.

[本文引用: 1]     

[26] Shi SB, Li TC, Li M, Liu SZ, Li AD, Ma JP (2015).

Interaction effect analysis of soil drought and strong light on PSII nonphotochemical quenching inKobresia pygmaea leaves

. Plant Physiology Journal, 51, 1678-1686. (in Chinese with English abstract)
[师生波, 李天才, 李妙, 刘世增, 李爱德, 马剑平 (2015).

土壤干旱和强光对高山嵩草叶片PSII反应中心非光化学猝灭的交互影响分析

. 植物生理学报, 51, 1678-1686.]

[本文引用: 1]     

[27] Silsbe GM, Kromkamp JC (2012).

Modeling the irradiance dependency of the quantum efficiency of photosynthesis

.Limnology and Oceanography: Methods, 10, 645-652.

[本文引用: 10]     

[28] Smyth TJ, Pemberton KL, Aiken J, Geider RJ (2004).

A meth- odology to determine primary production and phytoplank- ton photosynthetic parameters from fast repetition rate fluorometry

.Journal of Plankton Research, 26, 1337-1350.

[本文引用: 8]     

[29] Takahashi S, Badger MR (2011).

Photoprotection in plants: A new light on photosystem II damage

.Trends in Plant Sciences, 16, 53-59.

[本文引用: 1]     

[30] Vijayan P, Browse J (2002).

Photoinhibition in mutants of arabidopsis deficient in thylakoid unsaturation

.Plant Physiology, 129, 876-885.

[本文引用: 1]     

[31] Ware MA, Belgio E, Ruban AV (2015).

Photoprotective capacity of non-photochemical quenching in plants acclimated to different light intensities

.Photosynthesis Research, 126, 261-274.

[本文引用: 2]     

[32] Webb WL, Newton M, Starr D (1974).

Carbon dioxide ex- change of Alnus rubra: A mathematical model

. Oecologia, 17, 281-291.

[本文引用: 3]     

[33] White AJ, Critchley C (1999).

Rapid light curves: A new fluo- rescence method to asses the state of the photosynthetic apparatus

.Photosynthesis Research, 59, 63-72.

[本文引用: 1]     

[34] Ye ZP, Hu WH, Xiao YA, Fan DY, Yin JH, Duan SH, Yan XH, He L, Zhang SS (2014).

A mechanistic model of light-response of photosynthetic electron flow and its application

.Chinese Journal of Plant Ecology, 38, 1241-1249. (in Chinese with English abstract)
[叶子飘, 胡文海, 肖宜安, 樊大勇, 尹建华, 段世华, 闫小红, 贺俐, 张斯斯 (2014).

光合电子流对光响应的机理模型及其应用

. 植物生态学报, 38, 1241-1249.]

[本文引用: 1]     

[35] Ye ZP, Robakowski P, Suggett DJ (2013a).

A mechanistic model for the light response of photosynthetic electron transport rate based on light harvesting properties of photosynthetic pigment molecules

.Planta, 237, 837-847.

[本文引用: 4]     

[36] Ye ZP, Suggett DJ, Robakowski P, Kang HJ (2013

b). A mechanistic model for the photosynthesis-light response based on the photosynthetic electron transport of PSII in C3 and C4 species

.New Phytologist, 199, 110-120.

[本文引用: 4]     

/