未实现的最优设计挑战早期脊椎动物演化的约束模型

《iScience》：Unrealized optimal designs challenge models of constraint in early vertebrate evolution

　　

有颌脊椎动物主导了现代脊椎动物的多样性，而大多数无颌脊椎动物谱系在泥盆纪期间消失。它们的衰落常被部分归因于由坚硬、重装甲身体蓝图所施加的形态功能约束。研究人员通过在翼甲鱼形类（pteraspidomorphs）中整合三维几何形态测量学（three-dimensional geometric morphometrics）、计算流体动力学（computational fluid dynamics, CFD）、帕累托最优性（Pareto optimality）和系统发育比较建模（phylogenetic comparative modeling）来检验这一假说。流体动力学性能在头盾形态空间（morphospace）中变化显著，较晚的谱系倾向于占据更高效的区域。然而，在整个类群的历史中，包含高度最优设计的广阔区域始终未被探索。演化模型拟合几乎不支持阻止进入这些区域的硬表型边界（hard phenotypic bounds）。因此，翼甲鱼形类保留了大量未实现的潜在功能，这表明它们的灭绝发生在完全探索其身体蓝图可用的可行流体动力学解决方案范围之前。更广泛地说，这些结果挑战了认为内在形态功能约束耗尽了无颌脊椎动物演化潜力的解释。

研究背景与科学问题

有颌脊椎动物如今主导着脊椎动物的多样性，而现存的无颌脊椎动物仅包括盲鳗和七鳃鳗。然而，早古生代的化石记录揭示了大量已灭绝的无颌脊椎动物谱系，其中许多被称为甲胄鱼类（ostracoderms），它们不仅缺乏颌部，还拥有坚硬沉重的骨化头盾（headshield）且缺乏成对附肢。长期以来，这种保守的身体蓝图被认为施加了功能约束，阻碍了生态适应性和演化潜力：僵硬的前部骨骼和减少的控制面可能损害了在高要求运动模式下的表现，从而强化了无颌脊椎动物受生态限制且无法适应性多样化的观点。甲胄鱼类经验形态空间中观察到的早期饱和与广泛功能趋同模式也与这一解释一致，暗示内在生物力学限制塑造了其演化轨迹。这种约束常被用来解释大多数无颌脊椎动物在泥盆纪的衰落和最终灭绝。本研究旨在探究无颌脊椎动物身体蓝图的结构和功能特性是否从本质上限制了其可演化性（evolvability），进而增加了它们在生物或非生物因素驱动下被更灵活的有颌类群取代的脆弱性。

研究内容与主要结论

研究人员选择翼甲鱼形类（pteraspidomorphs）作为案例研究对象，这是早期无颌脊椎动物中分类学和形态学多样性最高的类群。他们构建了代表性物种的3D模型，应用几何形态测量学量化形状变异并推导出理论形态空间。通过计算流体动力学模拟水流围绕一系列理论形态的流动，评估了关键流体动力学参数的性能景观。随后采用多目标帕累托排序框架识别在多个性能标准间实现最佳权衡的形态。将最优值投影到形态空间并结合时间和系统发育背景后，发现整个类群历史中存在大片持续未被占据的最优区域。演化模型拟合表明，硬表型边界并非阻止进入这些区域的原因。因此，翼甲鱼形类保留了大量未实现的潜在功能，其灭绝发生在充分探索身体蓝图可用流体动力学解决方案之前。这些结果挑战了内在形态功能约束耗尽无颌脊椎动物演化潜力的传统观点。

研究意义

该研究质疑了经典叙事中的核心论点——即有颌脊椎动物的关键创新解除了早期约束并开启了新的生态机遇。相反，无颌脊椎动物中并未观察到可访问性能景观的饱和，表明所谓的约束可能被夸大甚至从未存在。这提示甲胄鱼类的衰落并非由于身体蓝图强加的不可逾越限制，而是可能涉及偶然性、有限的演化时间或其他外部因素。

关键技术方法

1. 1.
   三维几何形态测量学：对51种翼甲鱼形类头盾进行数字化重建，通过5个固定标志点和4条曲线半标志点定量形态变异，经广义普氏分析叠加后进行主成分分析，构建由PC1和PC2定义的二维理论形态空间。
2. 2.
   计算流体动力学模拟：对49个理论形态模型在两种缩放方式（等雷诺数和等体积）下，使用ANSYS Fluent进行稳态不可压缩三维流场模拟，提取阻力、升力、升阻比及俯仰和偏航力矩五个参数，覆盖0°至40°的俯仰和偏航角度偏差。
3. 3.
   多目标帕累托最优性分析：基于五个生物力学标准——低阻力、高升力、高升阻比、俯仰和偏航轴的稳定力矩，对理论形态进行帕累托排序，识别非支配的最优解。
4. 4.
   系统发育比较建模：利用100棵时间校准的系统发育树，拟合布朗运动、有界布朗运动、奥恩斯坦-乌伦贝克等连续性状演化模型，通过赤池信息准则比较模型支持度，检验未占据最优区域的演化可及性。

研究结果

头盾形状的系统发育形态空间

理论形态空间主要由两个主成分定义，PC1（37.4%）捕捉头盾内外侧扩张的变化，PC2（21.2%）描述背腹厚度变异。实际异甲鱼类物种占据相对宽广的区域，集中在中间PC1和PC2得分附近，避开最极端的形态配置。主要类群形成不同的聚类，沿PC1方向依次排列，显示非随机的形状空间占据模式。

跨理论形态空间的流体动力学性能

性能景观揭示了一致趋势：阻力沿PC1轴降低；升力随PC1降低而增加，在体积归一化条件下随PC2增加而增加；升阻比沿对角线梯度变化，随PC2增加而增加、随PC1增加而降低。俯仰和偏航旋转力矩在形态空间大部分区域为负值，表明整体稳定倾向。实际形态聚集在与有利流体动力学状态相关的区域——低阻力、低至中等升力、中等升阻比和稳定力矩。与独立经验CFD数据的比较证实，PC1-PC2简化形态空间保留了有意义的流体动力学信号。

流体动力学最优性景观与演化轨迹

多标准帕累托框架一致识别出一个连续的、位于PC2正端的最优区域。帕累托等级数量为6至9级，但实际形态大多局限于前四或五级区域，仅有少数晚期分化类群延伸至最高最优性区域边缘。结合系统发育和时间信息后发现，多个谱系中较年轻的物种相对于其祖先持续向更优区域移动。

检验未占据最优区域的演化可及性

演化模型拟合几乎不支持限制进入未占据最优区域的硬边界。沿PC1，无界布朗运动在所有100棵系统发育树中均优于有界模型；沿PC2，无界模型同样获得更强支持。奥恩斯坦-乌伦贝克模型在PC1上是最佳支持模型，其推断的最优值位于占据区域内部而非极端表型附近，表明吸引力朝向内部最优而非受硬表型极限约束。PC2上模型不确定性较大，但总体而言有界模型未获得积极证据支持。

讨论与结论

研究结果表明，不同形态空间区域表现出显著对比的流体动力学性能，实际形态趋向于占据提供有利性能特征的中心区域。然而，大片最优区域在整个类群历史中始终未被占据。尽管多个谱系显示出向最优性演化的趋势，但已实现与潜在流体动力学性能之间存在持续差距。这些发现挑战了无颌脊椎动物身体蓝图施加不可逾越功能约束的传统观点。存在一片连续且相邻的最优区域未被占据，且演化模型不支持硬表型边界的限制，表明这些区域在演化上是可及的。虽然可能存在发育、解剖或其他功能约束，但它们未留下可检测的大尺度演化信号。因此，研究结果更符合以下情景：甲胄鱼类的衰落反映了未实现的适应性潜力，而非身体蓝图限制导致的形态空间饱和。这质疑了有颌脊椎动物创新解除早期约束的核心论点，提示无颌脊椎动物中可访问性能景观的饱和缺失意味着那些假设的约束可能被夸大。总之，翼甲鱼形类头盾流体动力学测试的方面表明，可行的功能设计可能在演化上仍可触及，甲胄鱼类的历史不仅受约束和替代塑造，也受未实现的适应性潜力影响。