评估吡虫啉、温度和食物可得性对水生昆虫联合效应的DEB-TKTD模型构建与应用

《Aquatic Toxicology》：A DEB–TKTD model for assessing the combined effects of imidacloprid, temperature and food availability on aquatic insects

为解决如何评估多重环境应激因子对水生生物毒性的问题，研究人员结合动态能量预算(DEB)理论，开发了针对摇蚊(C. riparius)和四节蜉蝣(C. dipterum)的毒代-毒效(TKTD)模型。研究揭示了暴露时长是比暴露起始时间更关键的效应驱动因子，在短期暴露中温度是主要敏化因素，而在长期暴露中食物可得性更为重要，为理解野外复杂环境下的农药风险评估提供了有力工具。

　　

在洁净的实验室中，研究者们通常小心翼翼地将水生昆虫暴露于单一浓度的农药之下，以评估其对环境的毒性风险。然而，在真实的河流、湖泊之中，这些生物面临的挑战远不止于此。温度的波动、食物的短缺，以及其他化学物质的混杂，共同构成了一个复杂多变的“压力场”。尤其对于新烟碱类农药，如吡虫啉，其在水生昆虫中的毒性评估存在巨大变异性，不同研究报道的半数效应浓度(EC₅₀)可相差数十倍。这背后，究竟是哪些环境“幕后推手”在起作用？它们是如何联手放大或改变农药的毒性效应？将这些“幕后”因素纳入考量，是连接“理想”实验室数据与“复杂”现实环境，进行精准生态风险评估的关键所在。

为了回答这些问题，来自希腊克里特大学的一个研究团队在《Aquatic Toxicology》上发表了一项研究。他们以动态能量预算(Dynamic Energy Budget, DEB)理论为基础，开发了整合毒代动力学与毒效动力学(Toxicokinetic–Toxicodynamic, TKTD)的数学模型，并选取了两种对农药敏感且具有生态指示作用的水生昆虫——摇蚊(Chironomus riparius)和四节蜉蝣(Cloeon dipterum)作为模式生物，以新烟碱类杀虫剂吡虫啉为目标化合物。通过模拟不同环境变量（温度、食物可得性、暴露时机与时长）与毒物浓度的组合情景，该研究成功揭示了驱动毒性效应及数据变异性的核心因素，为理解多重应激下的毒性机制和提升风险评估的可靠性提供了新颖而有力的工具。

本研究所采用的关键方法包括：1. DEB-TKTD模型开发与参数化：基于物种特有的DEB生理学模型，耦合毒代动力学(TK)与毒效动力学(TD)模块，用以描述生物能量分配、化学物质体内动力学及其对摄食、活动和生存的影响。2. 实验数据拟合：利用针对C. riparius和C. dipterum幼虫的吡虫啉暴露实验数据（包括摄食抑制、活动性和存活率），对模型的TKTD参数（如无效应浓度c₀、损伤速率常数k_d、致死速率b_M、耐受浓度c_T等）进行估计和验证。3. 多情景模拟分析：在已校准的模型基础上，设定并运行涵盖不同温度、食物供应水平、暴露起始时间、暴露持续时间和毒物降解速率的数百种组合情景，计算并比较不同条件下的半数致死浓度(LC₅₀)等关键毒性指标。4. 蒙特卡洛模拟：为解释文献中报道的C. dipterum的效应浓度存在高度变异的现象，研究引入了五种变异来源（包括暴露前温度、个体初始长度、最大同化率、致死速率和耐受浓度的个体间差异），通过500次蒙特卡洛模拟，重现了观测到的效应浓度变异范围。

研究结果

3.1. 观察到的模式与模型预测

模型成功拟合了吡虫啉对两种昆虫摄食和活动性的影响数据。对于C. riparius，更长的预暴露时间会降低导致摄食减少20%所需的浓度。相比之下，其活动性在所有预暴露时间下受影响较小。对于C. dipterum，由于其损伤动力学更快(k_d= 1.42 d^-1)，摄食和活动性在相似的浓度下都表现出明显影响。而在C. riparius中，摄食抑制效应比活动性抑制效应更强。慢性暴露数据对确定C. dipterum的无效应浓度c₀和损伤速率k_d至关重要。

3.2. 多重应激的效应

通过情景模拟，研究揭示了不同环境因子的相对重要性：

• •
  暴露时长：在恒定浓度暴露下，LC₅₀值随暴露时间延长而降低，并在最初几天下降最剧烈，随后趋于稳定。高温加速了这种平台期的到来。
• •
  温度与食物的相对作用：在短期暴露中，温度是影响敏感性的主要因素，温度从15°C升至30°C可使LC₅₀降低约70%，表明高温加剧了急性毒性。在长期暴露中，食物可得性的影响变得更为关键，食物匮乏会使LC₅₀降低约20%，而温度的影响减弱。这是因为温度通过影响生物的所有生理速率（生长、呼吸等）以及毒物的摄取、消除和解毒过程，全方位地调控毒性反应。
• •
  暴露起始时间：生物在孵化后越早暴露，对毒物越敏感（LC₅₀值越低），显示出体型依赖的敏感性。在长期食物匮乏条件下生长的个体，由于体内储备能量（Scaled reserve density）低，会影响毒物的分配系数，从而改变其内部浓度动态。在短期暴露前经历长期食物匮乏的个体，可能因体型较小而表现出不同的毒性反应模式，这凸显了暴露前的生理状态对毒性结果的重要性。
• •
  模拟观测到的变异性：针对文献中C. dipterum的96小时-EC₅₀值从1 μg/L到40.9 μg/L的巨大差异，研究通过蒙特卡洛模拟引入了环境（暴露前温度）、生理（初始长度、最大同化率）和化学敏感性（致死速率b_M、耐受浓度c_T）三个方面的个体间变异。模拟结果表明，仅引入这五种变异源，模型就能重现文献中报道的宽范围效应浓度。其中，致死速率b_M的个体间差异是造成观测变异性的主要来源，其次是无效应浓度c₀的差异。
• •
  时变浓度暴露：模拟毒物浓度随时间衰减（如半衰期9天）的更真实场景。结果显示，在低降解率（即毒物持续存在）时，初始浓度是决定生存率的关键因素。C. riparius对浓度的微小变化反应更为剧烈，生存率呈陡降趋势，而C. dipterum的响应则更为平缓。在较高温度下，导致生存率降低所需的浓度阈值也更低。

结论与讨论

本研究通过构建和运用DEB-TKTD模型，系统评估了环境因子（温度、食物）如何与吡虫啉暴露相互作用，影响两种关键水生昆虫的毒性响应。核心结论是：在恒定暴露浓度下，暴露时长（急性或慢性）是比暴露起始时间更关键的效应驱动因子。在短期暴露中，温度是主要的敏化因素；而在长期暴露中，食物可得性变得更为重要。两种物种虽然都对吡虫啉高度敏感，但响应模式不同：C. riparius对浓度变化表现出剧烈而突然的响应，而C. dipterum的响应更为渐进，且对环境变异表现出更高的敏感性。

这项研究的重要意义在于：首先，它揭示了实验室毒性数据产生高度变异性的潜在机制，指出个体生理状态、暴露前环境历史以及遗传/适应差异是主要贡献者。这解释了为何使用野外采集个体进行实验时，结果往往难以重复。其次，所开发的DEB-TKTD模型为预测多重环境应激下的毒性效应提供了强大的定量工具。模型能够整合生物的能量状态、生长动态和毒物动力学，从而更真实地模拟野外复杂条件。最后，研究结果对改进生态风险评估和实验设计具有直接指导价值。它强调在实验室毒性测试中，需要更好地控制和报告关键条件，如食物供应、生物来源和暴露前驯化历史，以减少数据变异性，提高不同研究间的可比性和对野外风险预测的可靠性。通过弥合可控实验室条件与复杂自然环境之间的鸿沟，这项工作朝着更精准、更稳健的环境风险评估迈出了关键一步。