基于多相荧光动力学曲线特征...质生物毒性综合表征参数构建_胡丽.pdf

第 卷，第期 光谱学与光谱分析 ，年月 ，基于多相荧光动力学曲线特征的水质生物毒性综合表征参数构建胡丽，殷高方，赵南京，傅强合肥师范学院物理与材料工程学院，安徽 合肥 中国科学院安徽光学精密机械研究所，中国科学院环境光学与技术重点实验室，安徽 合肥 摘要基于藻类光合抑制效应的水质生物毒性检测方法，具有响应快速、测量简捷等优点。然而现有十多种光合荧光参数都源于 光合反应中心，对典型的 光合作用抑制剂，均表现出响应灵敏度差或无响应。利用毒性胁迫会引起藻类荧光动力学曲线发生变形、且变化程度与毒性强度成正比的特性，直接以多相叶绿素荧光动力学曲线形态为分析对象，在准确定位特征位点的基础上，构建基于曲线分段抑制的综合表征参数 ，并与常用的光合荧光参数 以及性能指标 进行毒性响应时间、响应灵敏度和稳定性的对比分析。结果表明，综合参数 对于 抑制剂阿特拉津，即在全样本表现出显著性抑制，最低检测限和 分别为 和 ，对 比 、分 别 降 低了 、和 、。对于典型 抑制剂百草枯，综合参数 在 表现出显著性抑制，对应的检测限和 分别为 和 ，而参数 无响应，具有时间响应优势；在长时抑制下，对应的最低检测限和 分别为 和 ，对比 、分别降低了 、和 、。上述结论表明了综合参数 作为响应指标可用于检测 及 抑制剂的生物毒性，且表现出良好的响应灵敏度和稳定性，解决了藻类光合抑制法对 抑制剂响应灵敏度低的问题，且对于统一不同种类毒性物质的毒性响应指标提供了关键参数。关键词生物光学；生物毒性；光合抑制效应；综合参数；毒性响应中图分类号：文献标识码：（）收稿日期：，修订日期：基金项目：国 家 重 点 研 发 计 划 项 目（），安 徽 省 科 技 重 大 专 项（），安 徽 省 杰 出 青 年 科 学 基 金 项 目（），国家自然科学基金项目（），安徽省教育厅重点项目（），中国科学院环境光学与技术重点实验室开放基金项目（），光电探测科学与技术安徽高校联合重点实验室项目（），年高校优秀青年骨干教师国内访学研修项目（皖教工委函 号），高校优秀拔尖人才培育项目（）资助作者简介：胡丽，女，年生，合肥师范学院物理与材料工程学院副教授 ：通讯作者 ：引言藻类作为生态系统初级单细胞生物，对毒物敏感、世代周期短、易于培养、可直接观察细胞水平上的中毒症状，是理想的生物毒性受试生物。区别于生长抑制方法采用细胞密度、细胞产量和生长速率为反应终点，光合荧光参数能够实时反映活体细胞光合作用状态，是更为快速敏感的藻类生物毒性测试反应终点。叶绿素荧光动力学过程根据其变化趋势和发生时间分为快相荧光和慢相荧光两个阶段，目前快相荧光机理研究已形成较为完整的体系，现有十多种光合荧光参数都源于快相荧光，这些参数能够反映捕光色素吸收、反应中心、等电子受体氧化还原状态，但无法分析 及其后面的光合能流过程。慢相荧光已被证明与 反应中心状态、非光化学猝灭，合成，卡尔文循环等过程密切相关，但由于慢相荧光形成机理至今尚不明晰、缺乏相关理论模型及表征参数，至今其所蕴含的光合作用信息仍无法有效使用。因此，以现有光合荧光参数为表征的藻类光合抑制法对百草枯、敌草快等典型的 光合作用抑制剂，均表现出响应灵敏度差或无响应。鉴于此，本文另辟蹊径，利用毒性胁迫会引起藻类荧光动力学曲线发生变形、且变化程度与毒性强度成正比这一特点，直接以多相叶绿素荧光动力学曲线形态为分析对象，构建一种基于曲线分段抑制的水体生物毒性综合表征参数 （），并将 与常用的光合荧光参数 以及 构建的性能指标 ，在毒性响应时间、响应灵敏度和稳定性三个方面进行对比分析，检验综合参数 作为响应指标在检测 及 抑制剂生物毒性中的有效性和稳定性。基于曲线分段抑制的综合参数构建多相叶绿素荧光动力学曲线包含、七个相点，如图（）所示。毒性物质作用不同特征位点，导致被作用的特征位点处电子传递链被打断、电子受体被阻塞，正常光合作用过程被阻断，对应的荧光动力学曲线特征相点的位置及相点之间的荧光曲线线型发生变化，如图（）所示。基于此，本文从以下两个方面构建藻类光合作用状态的综合表征参数 ：（）准确跟踪毒性物质胁迫下的多相叶绿素荧光动力学曲线、相位特征时间，这是构建基于曲线分段抑制综合表征参数的前提。（）在准确定位多相叶绿素荧光动力学曲线各相位的基础上，最大程度提取、整合各相之间曲线形变信息，这是构建综合表征参数的关键。图（）藻类叶绿素荧光诱导动力学曲线示意图（仿 ）、（）阿特拉津和百草枯胁迫下的荧光动力学曲线 （），（）多相叶绿素荧光动力学曲线相位跟踪通过曲线起点、终点和最大值分别确定、和相位；、是曲线上的瞬态相位，当光合作用位点的供体侧氧化与受体侧还原速率平衡时，曲线上升或下降趋势变缓，出现相位点，即曲线拐点，因此通过分析曲线上的局部曲率最大值，便能够快速定位各相位点出现的时间。曲线曲率具体计算过程如下：首先，曲线采用（）式所示的阶多项式进行拟合，随着多项式阶数的增加，拟合误差减小并趋于平稳时的值选为最终的阶数期望值。（）式（）中，为荧光强度、横坐标是时间的对数、是多项式系数、是多项式次数、为拟合多项式的阶数。然后，利用式（）和式（）求得拟合曲线的一阶导数和二阶导数，再通过式（）计算出曲线上各点的曲率，最后对曲率取对数让曲率变化更加明显，从而实现、相位准确定位。（）（）（）（）（）综合表征参数 构建毒性物质作用藻类某些特征位点或光合作用过程时，最终都会引起与位点或过程相关的荧光上升或者下降，进而导致、曲线下方包围面积发生变化，且相面积变化程度直接反映了毒性物质对藻类的抑制程度。鉴于此，通过对比空白样本，得到、相面积的抑制率，并对所有相面积抑制率数据项加和，获得基于曲线分段抑制的综合表征参数 ，该参数理论上可全面反映毒性物质对藻类所有光合作用位点或过程的抑制程度。图藻类光合抑制效应的综合表征参数构建示意图 （）曲线分段相面积为了消除实验条件（如受试藻浓度、荧光探测灵敏度）对荧光绝对强度的影响，进而影响曲线分段相面积计算，首先采用最大值归一化方法对荧光曲线进行归一化处理；其次，由于曲线各相间时长存在量级上的差异，为了平衡各相间曲线线型表征光合抑制效应的贡献率，将时间取对数；最后通过式（）计算曲线分段相面积，式中（）是归一化后的荧光强度，为对数时间，、是曲线起止相位，例如 相面积 是 相位到相位曲线段下方与坐标轴围成的面积。（）（）光谱学与光谱分析第 卷（）曲线分段抑制率及综合参数 曲线分段抑制率利用式（）计算，为曲线在 相之间的抑制率，为空白样本 相之间的曲线包围面积；为实验样本 相之间的曲线包围面积。（）为了防止数据项加和过程中抑制率的正负抵消一部分曲线形变信息，先对各分段抑制率取绝对值再进行数据项加和。同时，为解决曲线分段面积大小不同带来的权重差异，利用曲线上部分包围面积抑制率和下部分包围面积抑制率的均值，构建综合表征参数 。（）（）（）式（）中，是综合表征参数，表示 相之间上部分包围面积抑制率，表示 相之间下部分包围面积抑制率，其中 相包含 、段。实验测量与结果 实验测量实验以常见淡水藻蛋白核小球藻为受试生物，选取两种典型的作用位点抑制剂：阿特拉津（，抑制位点在 中心的处）、百草枯（，抑制位点在 反应中心处）。多相叶绿素荧光动力学曲线通过可变光脉冲诱导叶绿素荧光分析仪（，中科院安徽光机所研制）测量（该技术采用光电直流检测方式，点和点噪音干扰较大，本文仅分析 点，并加入 特征点）。蛋白核小球藻选种及扩大培养方法同参考文献 ，毒性物质用去离子水配制，并根据实验需要进一步稀释成不同浓度工作液，阿特拉津浓度梯度为 、和 ，百 草 枯 浓 度 梯 度 为 、和 ，每个样品设置个平行样。实验结果（）时间响应灵敏度对两种毒性物质分别进行 、和 胁迫实验，根据上文相位跟踪分析方法，得到蛋白核小球藻在阿特拉津和百草枯抑制下的相位位置。结果表明，两种毒性物质对应的相位、位置基本相同，分别位于、，其中点位置随浓度变化差异较大，与已有文献结论一致。在准确定位特征位点的基础上，利用单因素方差分析（）方法计算不同抑制时间下实验组与空白组之间的显著性差异，获得参数 、时间响应特性。结果如表所示，在阿特拉津抑制 后，参数 、在全浓度范围内均表现出显著性差异（），因此，选择 作为阿特拉津抑制实验的时间节点。百草枯抑制 后，、在 及以上浓度上表现出显著性差异（），但 在全浓度范围内均未表现出显著性差异（这与已有研究结果一致，对 抑制剂响应灵敏度差或无响应）；抑制后，、在全浓度范围表现出显著性差异（）。因此，满足百草枯中高浓度样液快速测量需求，可以选择 作为时间节点，采用综合参数 、性能指标 作为响应指标；对于中低浓度样液，可以选择作为时间节点。综上，参数 、均可以作为阿特拉津和百草枯的毒性响应参数，但 、在 抑制剂百草枯的毒性检测中，具有时间响应优势。图可变光脉冲诱导叶绿素荧光分析仪结构示意图 （）浓度响应灵敏度分析参数 、在阿特拉津 抑制、百草枯长时抑制和 短时抑制的剂量效应关系，符合生物毒性实验反“”形，以阻滞增长模型（）为基础进行非线性拟合，均大于 （如 图 图 所示）。在此基础上，建立低浓度毒性物质的线性剂量效应关系（如图图所示），根据国际纯粹与应用化学联合会（）对最低检测限的定义，利用式（）计算以参数 、为响应指标的毒性物质最低检测限，且以最低检测限表示不同参数对毒性效应的响应灵敏度。（）（）式（）中，为空白样品多次测量的标准偏差（本研究中空白样品重复测量 次）；为线性剂量效应曲线的斜率。以参数 、为响应指标的阿特拉津、百草枯最低第期 胡丽等：基于多相荧光动力学曲线特征的水质生物毒性综合表征参数构建检测限如表所示。结果表明，对于典型 和 抑制剂，以综合参数 为响应指标所获得的检测限最低。其中，阿特拉津 抑制，与常用参数 和性能指标 相比，检测限分别降低了 、；对于百草枯 抑制，的检测限比 降低了 ，参数 无响应；对于百草枯 抑制，对比 和 ，检测限分别降低了 、。因此，以综合参数 为毒性响应指标，更有利于提高藻类光合抑制法水体生物毒性的检测灵敏度。表 、实验组与空白组之间的差异性分析结果 ，（）（）（）注：“”、“”分别代表试验组与空白组在 和 水平上有差异图阿特拉津 抑制下（）、（）、（）对应的剂量效应关系 ，对应空白样的 分别为 ，（），（），（）图百草枯 抑制下（）、（）、（）对应的剂量效应关系 ，对应的空白样 分别为 ，（），（），（）光谱学与光谱分析第 卷图百草枯 抑制下（）、（）对应的剂量效应关系 ，对应的空白样 分别为 ，（），（）表 、为响应指标的毒性物质最低检测限 ，（）（）（）（）（）响应稳定性为了对比参数 、作为响应指标的稳定性，计算平行样相对标准偏差 ，用以表示毒性检测的稳定性。结果如表所示，阿特拉津 毒性胁迫下，的 平均值是 ，相对 和 ，分别降低了 、。百草枯 短时抑制下，对于 及以上浓度样本，的 平均值是 ，相对 降低了 ；百草枯长时胁迫下，的 平均值是 ，相对于 和 ，分别降低了 、。上述结果表明，稳定性最好的参数是 ，其次是表 、对两种毒性物质的 ，（）（）（）（），稳定性相对较差的是 。这与参数的构建方式密切相关，采用多个数据项加和的方式，而 采用三个数据项连乘的方式，连乘将放大单个数据波动对整个环节的影响，数据稳定性降低；参数 由初始荧光和最大荧光决定，受单点荧光强度影响较大，相比综合参数 ，更容易受测量环境干扰。结论直接以多相叶绿素荧光动力学曲线形态为分析对象，在准确定位特征相位的基础上，构建基于曲线分段抑制的综合表征参数 ，并与常用的光合荧光参数 以及性能指标 进行毒性响应时间、响应灵敏度和稳定性的对比分析。结果表明，对于 抑制剂阿特拉津，综合参数 表现优异，即 在 全 样 本 表 现 出 显 著 性 抑 制，相 对 于 、，最 低 检 测 限 和 分 别 降 低 了 、和 、。同时，对于典型 抑制剂百草枯，综合参数 在 即表现出显著性抑制，而参数 无响应，具有时间响应优势；在长时抑制下，相对于 、，最低检测限和 分别降低了 、和 、，表现出良