厘清阿尔卑斯岩羚羊角生长的个体发育与环境驱动因素

《Mammalian Biology》：Disentangling ontogenetic and environmental drivers of horn growth in Alpine chamois

　　

环境条件与遗传因素相互作用，共同塑造性选择特征。在温带有蹄类动物中，角以离散的年度增量生长，这为将特定年份的环境条件与发育轨迹联系起来提供了机会。研究人员将252只阿尔卑斯岩羚羊（Rupicapra r. rupicapra）的节段水平角数据与气候及遥感植被生产力（增强型植被指数，EVI）数据相结合，并应用路径建模框架，以量化性别、环境变异和补偿性生长如何共同决定总角长。正如预期，角生长主要发生在生命早期。第一个可测量的生长节段（L₂）解释了大部分个体间差异，并受性别和当年植被生产力的强烈影响。后续节段对总变异的贡献逐渐减小，但并非相互独立。特别是，研究人员检测到L₂与L₃之间存在条件依赖性补偿性生长：早期生长减弱的个体在植被生产力高的条件下倾向于表现出更强的恢复，而在环境条件较差时补偿性生长较弱。这表明补偿早期生长不足的能力取决于环境条件。因此，角大小既反映了生命早期建立的发育约束，也反映了后期年份中环境介导补偿的机会。由于角节段代表带时间戳的年度生长增量，它们为重建环境变异性如何塑造表型发育提供了框架。解读总角长需要理解哪些生命阶段主导其变异，这对理解性选择特征的进化和野生动物管理具有直接影响。

**解构阿尔卑斯岩羚羊角生长的个体发育与环境驱动因素：基于路径模型的整合分析**

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

在生态学与进化生物学中，理解内在因素与外在因素如何协同塑造形态性状的发育，始终是一项核心议题。对有蹄类动物而言，角作为典型的性选择特征，其形态与大小不仅受雄性间交配竞争等进化力量驱动，更受到资源获取与分配策略等近端机制的深刻影响。温带有蹄类动物的角每年形成清晰的生长环，这一特性使得研究者得以将特定年份的环境条件与个体的发育轨迹逐一对应，为精细解析角生长的调控机制提供了理想的天然实验体系。

然而，角生长节段之间并非相互独立。年增长量通常遵循强烈的年龄结构化轨迹——早期生命阶段投入高，角增量大，随后逐年递减。在此基础上，个体间还存在显著差异：某些物种中，早期生长较快的个体在后续年份仍保持相对较高的生长速率，形成生长轨迹的正相关；另一些物种则表现出补偿性生长，即早期生长相对不足的个体在后续年份出现加快生长的现象，形成负相关。北山羊、马鹿等物种的研究揭示了不同的相关模式，而岩羚羊的相关研究则呈现出不一致的结果：部分种群存在部分恢复生长，即在生命第三年补偿前两年的生长不足，但也有研究发现补偿生长缺失，提示补偿能力可能高度依赖当地环境条件。

在此背景下，本研究聚焦于阿尔卑斯岩羚羊（Rupicapra r. rupicapra），旨在回答一个关键问题：性别、气候、栖息地生产力与补偿性生长如何共同塑造角生长的个体发育轨迹？具体而言，研究人员假设早期生长对成体总角长具有不成比例的巨大贡献，并进一步假设补偿性生长仅发生在生命早期，且其强度受栖息地生产力调节——有利环境允许部分补偿，而低生产力环境制约恢复。该研究以德国贝希特斯加登国家公园为研究区域，结合角节段数据与环境变量，运用路径建模框架系统验证上述假说。研究论文发表于《Mammalian Biology》。

**主要技术与方法概要**

研究基于2020至2024年间在贝希特斯加登国家公园干预区内猎获的118只≥5岁个体（70雄、48雌）的角生长数据。测量分成L₂至L₅共4个节段，其中L₂代表生命前两年合并生长量，L₃至L₅为各年度增量；样本涵盖2006至2023年。环境数据包括区内气象站春季与初夏温度、MODIS卫星250米分辨率增强型植被指数（EVI）生长季面积下曲线值（AUC）、种群密度及采集山谷。统计上采用分段结构方程模型（piecewiseSEM），对总角长与各节段均取对数变换，以弹性系数解析各节段的相对贡献，并检验L₂与EVI对L₃的交互效应，量化条件依赖性补偿生长。

**研究结果**

**总角长与年度生长节段的弹性关系**
对数尺度下，L₂的弹性系数为0.79，即L₂每增加1%，总角长约增加0.79%；而L₃、L₄、L₅的弹性分别仅为0.13、0.05和0.03。标准化系数进一步印证：L₂（β=0.87）远大于L₃（β=0.31）、L₄（β=0.22）和L₅（β=0.14）。原始数据表明，L₂平均占总角长的约74%（±6%标准差）。由此得出明确结论：个体间总角长的差异主要源自生命早期的生长表现，后续节段仅贡献有限附加变异。

**个体与环境变量的间接效应**
性别是总角长变异最强的间接驱动因素，雄性各节段生长率均高于雌性，对总角长的间接效应β=0.498。采集山谷之间亦存在显著空间异质性（β=0.290），最大与最小角长的山谷间差异可能反映栖息地结构、饲草组成或微气候的差别。种群密度未检测到显著效应。当年春季至夏季温度对总角长显示出负间接效应（β=-0.110），主要通过抑制L₂实现；而生命第一年EVI对L₂和总角长呈负效应（β=-0.094），当年EVI则对L₂呈正效应，对总角长的正间接效应为β=0.094。这些结果说明角生长对营养资源可利用性的反应具有明确的年龄依赖性。

**补偿性生长的环境依赖性**
L₂与EVI的交互项在L₃生长模型中呈显著负效应（β=-0.073）。条件效应图揭示：在高植被生产力年份，早期生长（L₂）较差的个体在L₃期间表现出更强的恢复性生长；而在植被生产力较低的年份，这种补偿显著减弱，甚至完全消失。这一发现证实补偿性生长并非岩羚羊固有的生长轨迹特性，而是受环境条件动态调控的可塑响应。

**讨论与结论整合**

研究结果清晰表明，L₂节段对总角长的主导地位使角大小差异主要定型于生命早期，后期节段的增量贡献相对有限，这与岩羚羊及近缘物种的既有认知一致。性别的强效应可能反映幼年个体在建立优势等级过程中对角长度的需求，尽管目前缺乏直接证据将早期角长与交配成功联系。山谷间差异提示栖息地质量的空间斑块性通过营养景观影响角发育，而密度效应的缺失可能与研究期间种群密度波动较小、未达密度制约阈值有关。温度负效应可能的机制包括高温降低觅食效率以避免热应激，以及暖春导致物候提前、冷季适应型高营养饲草资源减少。EVI的正向作用则从营养机制层面得到支持：角生长涉及骨核心矿化（需要钙、磷）和角蛋白鞘合成（需要含硫氨基酸及锌、铜等微量元素），植被生产力综合反映了这些营养素的供给潜力。

讨论中尤为值得关注的是补偿性生长的环境依赖性。本研究中高生产力条件下L₂与L₃呈负相关，意味着早期生长劣势可在有利环境中被部分弥补；而低生产力条件下二者无显著相关，说明资源限制既阻碍了优势延续，也阻断了补偿通道。这种条件依赖的发育缓冲机制表明，岩羚羊的角生长轨迹是个体发育动态与外部生态机会相互作用的产物。

**研究结论**：早期角生长（L₂）对成体总角长变异的贡献最大，角大小本质上反映生命早期建立起来的发育约束。然而，后期补偿性生长的机会在环境条件有利时显著增强，在资源受限条件下则受到抑制。这种环境依赖性补偿表明，岩羚羊角生长受个体发育轨迹与外部生态条件共同调节，理解这一动态过程对于解读性选择特征的进化机制以及制定野生动物管理策略均具有重要的理论和实践意义。