计算流体力学揭示网纹溞（*Ceriodaphnia dubia*）针对食肉植物 neglecta 狸藻（*Utricularia* × *neglecta*）的可诱导防御的水力机械适应性

《Ecology and Evolution》：Computational Fluid Dynamics Reveal Hydromechanical Adaptations of Ceriodaphnia dubia as an Inducible Defence Against the Carnivorous Plant Utricularia × neglecta

　　

摘要
水蚤科（Daphniidae）水蚤是许多静水生态系统中的重要生态交互者。作为这些系统中数量最丰富的滤食者，它们在初级生产力和更高营养级之间形成了关键联系。由于在淡水食物网中的核心作用以及通过表型可塑性对一系列环境变化作出反应的能力，水蚤科被广泛视为生态研究中的模式生物。它们形成可诱导防御，这代表了一种表型可塑性形式，可降低生物体对许多动物捕食者的脆弱性。就在最近，研究人员在网纹溞（*Ceriodaphnia dubia*）中首次观察到针对食肉水生植物——南方狸藻（*Utricularia* × *neglecta*）的可诱导防御，同时捕食率显著降低。该物种的防御性状尚未被完全理解。研究人员观察到 *C. dubia* 游动速度较慢，从而减少接触机会。此外，其体型更为纤细。这是违反直觉的，因为该植物具有口宽限制（gape-limited），因此超出口宽生长是直观的策略。利用计算流体力学（CFD），研究人员提出了一项关于该可诱导防御功能的新假说。研究人员模拟了典型形态和 *U.* × *neglecta* 暴露形态下水蚤所受的拖曳力。结果表明，暴露于植物环境中的形态所受拖曳力降低，这可能与其第二触角按比例减小从而降低接触植物引发吸捕的触毛的概率协同作用。因此，研究人员提出一种基于形态改变与水动力特性改变相结合的保护性可诱导防御机制。这种新发现的水动力调整防御策略为水蚤作为生态响应生物体的声誉增添了又一例证。

**论文解读：计算流体力学揭示 *Ceriodaphnia dubia* 针对食肉植物 *Utricularia* × *neglecta* 的水动力可诱导防御机制**

**1 研究背景与科学问题**

在许多生态系统中，捕食作用是塑造物种适应性特征的重要选择压力。为应对捕食风险，部分猎物物种进化出了可诱导防御（inducible defence），即仅在存在捕食者信号时才表达的表型改变，包括形态、行为或生活史特征的调整。水蚤科（Daphniidae）是淡水浮游动物的重要类群，在静水生态系统中作为主要滤食者，连接初级生产与更高营养级，被广泛视为生态学模式生物。已有大量研究报道了水蚤针对无脊椎捕食者和鱼类形成头盔、加长尾刺等形态防御，以及昼夜垂直迁移、游泳速度改变等行为防御和繁殖策略调整。然而，针对食肉植物的可诱导防御研究较少。近期研究首次报道了网纹溞（*Ceriodaphnia dubia*）在暴露于食肉植物南方狸藻（*Utricularia* × *neglecta*）后表现出体型减小、游速降低等防御反应，且捕食率显著下降（Kruppert et al. 2022）。但该防御性状的功能机制尚不清楚。由于狸藻捕虫囊的入口尺寸有限（平均高度495 μm，宽度613 μm），传统观点认为针对口径受限捕食者的形态防御应表现为体型增大以超出口径范围。然而，*C. dubia* 在防御形态下反而生长更慢、体型更小，这一反直觉现象亟需解释。研究人员推测，除了降低接触概率之外，形态改变可能通过影响猎物所受的水动力作用力，如拖曳力（drag force），来降低被吸入捕虫囊的风险。为验证这一假说，研究人员利用计算流体力学（Computational Fluid Dynamics, CFD）对两种形态的*C. dubia*三维模型进行吸捕过程的拖曳力模拟。

**2 主要技术方法**

研究使用了已发表的三维形态数据（Kruppert et al. 2022），基于13个对照个体和8个防御个体构建平均形态模型，并利用共聚焦激光扫描显微镜（confocal laser scanning microscope, CLSM）获取触角结构后使用Blender软件虚拟装配到标准静息姿态。CFD模拟采用Ansys ICEM CFD、CFX Pre、CFX Solver和CFD-Post完成（ANSYS 2020 R1）。模拟区域为一个2 mm×2 mm的流体空间，其一侧设有模拟狸藻捕虫囊入口的管道结构（宽度613 μm，高度495 μm），通过移动壁面的数学函数模拟吸捕事件（持续约8 ms，最大流速4 m/s）。水蚤模型分别置于捕虫囊入口前200 μm、400 μm距离处、前上方、前下方、侧方，以及旋转90°（触角朝向入口）等多个位置，比较两类形态所受拖曳力差异。

**3 研究结果**

**3.1 体型、体宽和触角长度的形态测量结果**
在孵化后第3天和第6天，防御个体的体长和体宽均显著小于对照个体（第6天体长：对照728 ± 105 μm vs 防御529 ± 53 μm，p＜0.0001；体宽：对照500 ± 90 μm vs 防御300 ± 40 μm，p＜0.0001）。相对体宽（体宽/体长）在第6天显著减小（对照0.683 ± 0.04 vs 防御0.566 ± 0.03），表明水蚤暴露于狸藻后身体形态变得更为细长。防御个体第6天绝对触角长度（背支、腹支、基段）和触角刚毛长度总和均显著短于对照个体，但相对触角长度（触角长度/体长）在多数比较中无显著差异，提示触角的减小与体长基本成比例。根据三维模型估算，对照模型体积为0.000131 cm3，防御模型为0.000026 cm3，前者约为后者的5倍。防御个体模型在体长为对照73%的情况下，体积减少至对照的约20%，表明除线性尺寸缩小外，体形也发生了显著变化。

**3.2 不同位置下拖曳力的 CFD 模拟结果**
模拟结果显示，两种形态所受拖曳力在捕虫囊吸捕的初始阶段（约前2.5–3 ms）均接近零，随后急剧上升并在峰值流速时刻达到最大值。在所有模拟位置中，位于捕虫囊入口正前方时两种模型承受的拖曳力最大，侧方位置最小。具体而言：
（1）入口正前方位置（距离约200 μm）：对照模型达到最大拖曳力1.76×10?? N，防御模型为9.88×10?? N，约为对照的56%。
（2）有一定距离的正前方位置（约400 μm）：对照模型最大拖曳力3.28×10?? N，防御模型为1.87×10?? N，约为对照的57%。
（3）旋转90°（触角朝向入口）位置：该位置产生两种形态间最大差异，防御模型拖曳力（9.64×10?? N）仅为对照模型（5.15×10?? N）的约20%。
（4）入口上方：对照模型最大拖曳力1.73×10?? N，防御模型为8.71×10?? N，约为对照的一半。
（5）入口下方：对照模型最大拖曳力1.73×10?? N，防御模型为1.1×10?? N，约为对照的60%。
（6）入口侧方：对照模型最大拖曳力为1.81×10?? N，防御模型为1.35×10?? N，约为对照的75%。
在吸捕主阶段的所有时间点，防御模型所受拖曳力均低于对照模型，多数情况下仅为对照的一半或更少。

**4 讨论与结论**

本研究通过计算流体力学模拟首次证明，*C. dubia*在*U.* × *neglecta*存在下形成的形态改变——包括整体体型缩小、触角按比例缩短以及体积显著减小——能够有效降低猎物在狸藻捕虫囊吸捕水流中所承受的拖曳力。水蚤的触角伞状结构（antenna parachute）在低雷诺数条件（10–100）下水动力上可视为连续表面，尽管触角刚毛间存在微小间隙。吸捕水流作用于表面积较大的触角结构时，触角起到类似“帆”或“降落伞”的作用，因而触角的缩短直接减少了被捕水流“捕获”的面积，从而降低拖曳力。研究还表明，防御形态的此种水动力优势在不同空间位置（入口前方、上方、下方、侧方）均存在，特别是在触角朝向吸捕入口的姿态下，拖曳力降幅最高达约80%。结合此前观察到的防御个体游速更慢、捕食率显著下降的行为和生态学证据，形态缩小与水动力特性调整协同降低了*C. dubia*被吸入捕虫囊的风险。
研究还指出了几个限制因素。其一，模拟中的猎物模型保持固定位置，未模拟实际吸捕过程中猎物被吸入并向入口加速的运动，因此模拟结果是对真实受力情况的近似估计。其二，真实的捕虫囊吸捕压力差（最高约0.176 bar）取决于距上次吸捕事件的时间等因素，本研究未将这一变量纳入模拟。其三，研究使用单一克隆（clone S04）的*C. dubia*，其他克隆或水蚤物种的普适性尚需进一步验证。
结论：本研究的CFD模拟证实，*C. dubia*在*U.* × *neglecta*存在下表现出的形态防御（整体体型缩小、触角按比例缩短）显著降低了猎物在吸捕水流中所受的拖曳力，不同位置下降幅可达20%至80%。该水动力调整机制与先前报道的游速降低、行为回避等防御策略协同作用，共同降低了捕食风险。这一发现首次揭示了形态改变可作为一种水动力防御策略对抗食肉植物捕食者，拓展了水蚤可诱导防御研究的理论范畴，也提示其他小型猎物可能通过类似形态调整来应对吸捕型捕食者。论文发表在《Ecology and Evolution》。