基于分形几何的波动性不对称无法追踪潮间带栖息地形成型墨角藻的环境胁迫

《Journal of Phycology》：Fluctuating asymmetry based on fractal geometry does not track environmental stress in an intertidal habitat-forming fucoid

　　

易于测量且成本低廉的生物标志物的可用性，是对生物体环境胁迫监测的一项挑战。尽管有观点认为形态性状中偏离完美对称模式的非定向偏差（称为波动性不对称，FA）可能与环境胁迫直接相关，但这一假设已受到质疑。自相似对称性（分形）是处理FA分析的一种可能途径，但采用该途径的研究极少，且结果尚无定论。本研究评估了栖息地形成型墨角藻Fucus limitaneus的叶状体形态分形描述符是否在多种情境下反映环境胁迫，包括季节变化、潮间带垂直梯度以及入侵物种可能的长期影响。研究结果表明，分形指标在这些对比鲜明的环境条件下并未发生一致性变化，提示在该系统中，基于分形几何的FA不能可靠地追踪环境胁迫。这可能是因为胁迫通过其他生态和生物力学过程表达，而非通过叶状体自相似性的改变。

**基于分形几何的波动性不对称无法追踪潮间带栖息地形成型墨角藻的环境胁迫：论文解读**

环境胁迫的监测是生态保护与管理的重要环节，但缺乏易测且廉价的生物标志物。波动性不对称（fluctuating asymmetry, FA）指形态性状偏离完美对称模式的非定向偏差，理论上可反映个体在发育过程中因环境或遗传胁迫而将资源从对称性维持机制转移的过程，因而被提议为环境胁迫的候选指标。然而，这一假设已受到广泛质疑。自相似对称性（分形）为FA分析提供了新途径，分形维数（fractal dimension, D）可量化结构填充空间的程度，而盒计数回归中的估计标准误（standard error of the estimate, S_yx）被提议作为FA的替代指标。但已有研究结果矛盾，且极少在实际生态系统中得到验证。潮间带大型藻类的形态不仅受发育过程调控，还受波浪、拖曳、断裂等强机械力塑造，机械性修剪和选择性存活可能掩盖发育不稳定性，因此FA指标在这些生物中的适用性需专门评估。栖息地形成型褐藻Fucus limitaneus是直布罗陀海峡潮间带的主要建群种，其扁平分枝叶状体由顶端细胞驱动二叉分支，形成层级自相似结构，是研究该问题的理想模型。该研究发表于《Journal of Phycology》，旨在检验基于分形几何的FA是否能在多个环境情境下追踪Fucus limitaneus的环境胁迫。

研究人员利用两个观测数据集开展研究。第一组数据采集自直布罗陀海峡北侧Tarifa，于2022年2月（寒冷月，生物量最低）和8月（温暖月，繁殖高峰，生物量最高）分别在潮间带中潮区和分布上限采集样本，以比较季节和垂直梯度的效应。第二组为2011–2024年长期监测，覆盖Guadalmesí和Tarifa两个种群，采样均在仲夏进行；该期间直布罗陀海峡藻类群落遭受入侵种Rugulopteryx okamurae的严重干扰，种群覆盖度在2014–2018年间急剧下降且此后未恢复。研究人员采用盒计数法计算叶状体的D和S_yx，并通过单因素方差分析（analysis of variance, ANOVA）和线性模型评估变异。结果显示，季节与高度组合下D和S_yx均无显著差异；长期时间序列中，尽管覆盖度从2011/2014年的38.1%±10.1%降至2018–2024年的11.4%±6.9%，但D和S_yx在各年份间保持稳定。因此，基于分形几何的FA不能可靠追踪该系统中的环境胁迫。这一结论具有重要意义，提示在应用形态学指标时需结合胁迫的实际作用机制，而非简单假设形态与胁迫之间存在线性关系。

关键技术方法方面，样本采集于直布罗陀海峡北侧的Tarifa（季节×高度研究）以及Guadalmesí和Tarifa（长期研究）。在春季大潮低潮期，沿平行于岸线的样带每隔20厘米设置15×15厘米样方，收集样方中心无啃食痕迹的完整叶状体，每个处理15个（季节研究）或21个（长期研究）个体；约15%的叶状体因有草食痕迹被排除。采用盒计数法：将叶状体置于透明醋酸纤维膜间扫描（300 ppi），叠加指数递增的网格（0.125–4 cm2），计数含叶状体的盒子数，线性回归斜率绝对值为D，回归残差均方根为S_yx。为减小量化误差，每个样本随机重定位重复三次取均值；所有扫描图像在2024年由独立观察者以盲法处理，避免确认偏倚。统计使用单因素ANOVA和线性模型，显著性水平设为0.01，并计算效应量η²和η²_p。

研究结果分为两部分。其一，季节和垂直范围内高度对Fucus limitaneus的FA的影响：通过将季节与高度组合为单一因子的单因素ANOVA比较四种组合，D无显著差异（F_3,56=1.36, p=0.265, η²=0.07），总体均值为1.16±0.11；S_yx亦无显著差异（F_3,56=0.13, p=0.939, η²<0.01），总体均值为0.06±0.02。这表明季节和潮间带垂直位置对分形FA无可检测的影响。其二，2011–2024年Fucus limitaneus覆盖度和FA的时间序列：覆盖度在年份间和种群间差异显著，但年份效应远强于种群效应（线性模型，年份：F_4,44=76.65, p<0.001, η²_p=0.87；种群：F_1,44=11.26, p=0.002, η²_p=0.20），平均覆盖度从2011/2014年的38.1%±10.1%降至2018–2024年的11.4%±6.9%。然而，D和S_yx在各年份间均无显著变化（D：Guadalmesí F_4,95=2.23, p=0.071；Tarifa F_4,95=1.98, p=0.103；S_yx：Guadalmesí F_4,95=1.85, p=0.125；Tarifa F_4,95=0.54, p=0.704），总体均值分别为1.37±0.24和0.07±0.02。这些结果表明，尽管种群覆盖度显著下降，存活个体的叶状体自相似结构保持相对稳定。

讨论部分指出，潮间带墨角藻的形态由发育过程与外部机械力共同作用，波浪引起的拖曳、断裂和水动力载荷在塑造藻体形态中占主导地位；形态变异可能源于机械性修剪和个体的选择性存活，而非单纯的发育可塑性。本研究所涵盖的环境情境代表不同的胁迫机制：季节变化主要涉及生理限制（温度和生长动态），垂直梯度涉及脱水和热胁迫，长期衰退则与入侵种Rugulopteryx okamurae的竞争、覆盖、物理掩埋及其对单细胞和少数细胞发育阶段的化感作用有关。分形描述符的稳定性表明，胁迫可能主要通过人口统计过程（如断裂、脱落、竞争排除）选择性移除个体，而非逐步改变存活个体的叶状体结构。因此，研究结果不支持基于分形的FA与Fucus limitaneus环境胁迫之间存在普遍关系。研究结论可翻译为：因此，我们的结果不支持基于分形的FA与Fucus limitaneus环境胁迫之间存在普遍关系。然而，这并不表明不存在胁迫效应，而是表明所选的FA描述符与该系统中胁迫作用的主要途径不一致。与机械性能或种群动态更直接相关的替代性状（如断裂率、组织强度或存活概率）可能为潮间带大型藻类的胁迫提供更有参考价值的指标。此外，分形维数的估计和应用本身存在方法学挑战，不应被视为功能或生态过程的通用指标；但研究人员严格遵循了FA研究的建议（如随机采样、盲法测量），支持了结论的可靠性。