海洋鱼类最大寿命的体型、功能形态与环境相关性

《Journal of Fish Biology》：Body size, functional form and environmental correlates of maximum longevity in marine-occurring fishes

　　

海洋鱼类最大寿命在不同物种间差异显著，但其与体型的关系形式以及广泛环境梯度的作用仍不明确。研究人员分析了757个物种，以最大长度为主要体型指标，完整体重仅在敏感性分析中使用。二次长度模型较线性比较模型改善了拟合，并捕捉到体型–寿命关联中广泛的衰减现象。在控制该曲率后，最大深度与log10最大寿命呈正相关（每1000 m估计值=0.0692，95% CI 0.0312至0.1072，p<0.001）。纬度估计在基于分布范围指标和模型设定间存在差异，且在剔除整个科和目时纬度项表现较不稳定。所有分组或连续交互对比均未达到定义的FDR阈值q<0.05。FishTree匹配605个物种但排除典型数据集中的所有软骨鱼；因此分支长度系统发育广义最小二乘（PGLS）提供的是匹配子集敏感性而非完整系统发育覆盖。5000次物种级bootstrap重采样的线性对比与非线性分位数敏感性均未指示广泛的上尾陡峭化。仅真骨鱼、宽纲调整及科聚类分析均保留正的最大深度系数，而纬度在聚类不确定性下较不稳定。因此体型与最大深度提供了与最大寿命最一致的宏观生态关联。纬度仍具探索性，数据未识别因果机制。

研究背景与立项依据

海洋鱼类最大寿命整合了生长、成熟与死亡调度，是比较生活史研究的核心性状。体型常被视为最强比较相关变量，但跨谱系衰老与存活差异显著，单一指数无法描述全尺寸范围；同时最大年龄记录依赖 ageing method 与验证质量，且深度、纬度等环境梯度并不等同温度、氧可得性、生态死亡率或采样历史。已有研究指出深海鱼类常具慢生长、晚成熟与长寿命，但机制未明。鉴于此，研究人员在《Journal of Fish Biology》报道了该项工作，旨在明确体型–寿命关系的功能形式，并评估最大深度与纬度代理在控制曲率后的宏观生态贡献。

主要技术方法概述

研究人员自 FishBase 25.04 经 rfishbase 组装队列，从36,132条记录经咸水、正最大长度、正最大寿命、最大深度与纬度端点完整案例过滤得757种规范数据集；体型以最大长度（log₁₀转换）为主，体重经长度–重量式重建仅作敏感性。环境预测包括连续最大深度（每1000 m）与自北南界重建的绝对分布中点纬度（每10°）；系统发育以 fishtree R 包匹配605种并数值修正枝长满足超度量，用 ape/caper 跑 PGLS（Pagel's λ 最大似然估计）。功能形式比较采用线性、中心化二次与低复杂度 GAM（mgcv），AIC/嵌套F/残差平滑 RMS 判定二次保留；分布尾部用 quantreg 做 τ=0.10–0.95 线性与非线性分位数回归，5000次物种 bootstrap 取百分位区间；分类敏感性含仅真骨鱼、宽 Class 固定效应、科聚类 CR2 标准误（clubSandwich）；内部五折验证分物种随机、科阻断与目阻断。

研究结果

3.1 数据集覆盖与分析集

规范集757种中605种匹配 FishTree，152种未匹配；TL 记录595种；体重来源分直接424、长重推导261与未明72。深度组373/273/111，纬度中点组314/438/5（高纬组仅描述）。FishTree 匹配集含595真骨鱼与10软骨硬鳞类，排除全部75软骨鱼、3七鳃鳗、1腔棘鱼，界定 PGLS 为匹配子集敏感性。

3.2 非线性体型–寿命关联

二次模型 AIC=497.486，较低线性降43.345，调整 R² 由0.317升至0.356，残差平滑 RMS 自0.0971降至0.0378；中心化线性项0.5812（p<0.001），二次项?0.3233（p<0.001），中央共同支持区二次方向不变。导数零点近487 cm，仅15种超出且位于共同支持外，极端下降尾不解释为“大体型寿命降低”。GAM 灵活形状敏感一致但不作独立生物模型。

3.3 控制曲率后环境关联

主 OLS 中最大深度每1000 m 系数0.0692（p<0.001，调整 R²=0.379），折合原尺度约+17% 期望寿命；中央共同支持与仅中纬子集均为正。绝对中点纬度每10°为0.0200（p=0.001），但旧版端点均值/最大端点/严格共同支持区间含零，故纬度仅作探索。PGLS（n=605，λ=0.887）深度0.1062、纬度0.0397均正；体重敏感性中深度仍正，印证非长度专属。

3.4 诊断分组交互与分类敏感性

所有深度组间（epi-meso、epi-bathy、meso-bathy）与低减中纬斜率对比 FDR q 均≥0.216，连续交互（长度×深度、长度×纬度）未达 q<0.05。仅真骨鱼、宽纲调整、科聚类 CR2 下深度系数皆正（0.0832/0.0720/0.0692），纬度在科聚类中区间含零（0.0200，p=0.316）。75种软骨鱼子集深度亦正但支持弱。

3.5 系统发育、分位数与验证诊断

PGLS 匹配子集各项目符号同 OLS 方向，λ=0.887，调整 R²=0.352，分组对比不显著。5000次线性分位数 bootstrap 中 β_0.90?β_0.50=?0.1406（p=0.128），β_0.95?β_0.50=?0.2173（p=0.288），上尾未陡于中位；二次分位数局部斜率区间含零。预测 R² 随加深度由0.352→0.368、加纬度→0.376（物种随机折）；科阻断下加纬度反降至0.273，目阻断降至0.277，示纬度在分类阻断时不稳。

讨论浓缩

体型是最稳相关但关系为广域衰减二次而非单指数，极大体型尾端因样本稀而不支持“转折下降”；最大深度正向关联跨 OLS/PGLS/分类敏感/共同支持均稳，但深度是广谱生态代理（低温、低光、资源、捕食、采样）而非已辨明因果通路；分组斜率弱源于分类切段与曲线局部斜率混淆，共同曲率移除后无 FDR 显著交互；纬度依赖分布汇总法、在共同支持/科聚类/分类阻断下不稳，且中点代理不等于出现加权热暴露；PGLS 仅校正605匹配种协方差，未补软骨鱼缺失亦未愈环境测量限；记录极值受采样努力与 ageing 偏倚放大。

结论部分翻译

海洋鱼类最大寿命与体型呈正但非线性关系，简约二次模型捕捉观测范围内的广泛衰减。最大深度在体型之外表现出最一致的正向环境关联。纬度在基于分布范围指标、共同支持限制、聚类不确定性与分类阻断诊断间变动，故仍具探索性。分组斜率对比与连续交互在 FDR 校正后未显示中心化长度系数的稳健环境变异，分位数分析未指示广泛上尾陡峭化。这些结果描述有界的宏观生态关联，而非已检验的因果机制。