全株植物经济谱新探索：树皮功能性状纳入揭示亚热带木本植物整体适应策略

本研究由Jinlong Li（第一作者，福建师范大学）、Xiaoping Chen、Karl J. Niklas（康奈尔大学）、Jun Sun、Zhaoying Wang、Quanlin Zhong、Dandan Hu以及通讯作者Dongliang Cheng（福建师范大学）等人合作完成。研究成果于2022年以题为“A whole-plant economics spectrum including bark functional traits for 59 subtropical woody plant species”的研究论文形式，发表在国际知名生态学期刊《Journal of Ecology》（第110卷，第248-261页）上。该项工作得到了中国国家自然科学基金等多个项目的资助。

一、 研究背景与目标

该项研究属于植物生态学领域，核心聚焦于“基于性状的生态学”。该领域旨在通过研究植物的功能性状——即影响其生长、生存、繁殖以及与环境和生态系统互动的可测量形态、生理和物候特征——来预测和解释植物的适应策略、物种分布、群落构建及生态系统功能。

在植物功能性状研究中，一个核心概念是“经济谱”。最广为人知的是“叶片经济学谱”，它揭示了叶片功能性状间普遍存在的协同与权衡关系，例如，叶片氮含量高、比叶面积大的“快速投资-快速收益”策略与叶片寿命长、组织密度高的“保守投资-缓慢收益”策略之间的连续谱。近年来，这一概念已扩展至根系、木质部等其他器官。将植物不同器官的性状整合起来，就构成了“全株植物经济学谱”。PES被认为能更全面地反映植物在资源有限条件下的整体适应策略，即植物如何在生长、防御、繁殖等多方面分配资源。

然而，以往构建PES的研究大多关注叶片、根系和木材（木质部），常常忽略了另一个重要器官——树皮。树皮（位于形成层外的所有组织）承担着多重功能：作为机械支撑和保护层（抵御物理损伤、病虫害及火灾），参与光合作用（尤其在某些物种中），运输光合产物（韧皮部），以及储存水分和养分。因此，树皮同样是植物碳和养分的重要“投资”部位，其性状（如厚度、密度、养分含量）应与其他器官性状存在紧密的关联和权衡。尽管如此，关于“树皮经济学谱”的研究非常匮乏，将其系统性地整合进全株植物经济学谱的研究更是鲜见。

针对这一知识空白，本研究旨在构建一个真正意义上的“全株”经济学谱，将树皮、木材（木质部）和叶片的功能性状进行整合分析。研究以我国武夷山国家级自然保护区内59种代表性的亚热带木本植物（包括常绿和落叶物种）为对象，提出并检验了三个关键科学假说：（1）树皮、木材和叶片的性状相互关联并整合为一个协调的整体；（2）氮和磷养分性状是连接不同器官、形成整合性PES的枢纽性状；（3）常绿和落叶物种在全株水平上表现出不同的资源利用策略，常绿物种偏向保守策略，落叶物种偏向获取策略。

二、 研究方法与流程

本研究是一项系统的野外采样与室内性状测定、数据分析相结合的实证研究，具体流程如下：

步骤一：研究地点与样地设置 研究地点位于中国武夷山国家自然保护区。沿海拔梯度选择了四个森林群落作为研究样地，分别为：常绿阔叶林（海拔1319米）、针阔混交林（海拔1697米）、落叶林（海拔1818米）和高山矮林（海拔1975米）。这种设计旨在捕捉环境梯度下植物性状的变化。

步骤二：物种选择与样品采集 在每个样地内，设立了多个20米×20米的样方。从四个群落中共选取了59种代表性木本植物，涵盖了不同叶习性和生长型（常绿vs.落叶，乔木vs.灌木）。对于每个物种，选择三株具代表性的标准木。从每株标准木的树冠外围，随机采集三根长度约1米、直径约1厘米的向阳小枝。同时，从树冠的上、中、下及四周采集30-50片健康、成熟、无损伤的叶片。所有样品采集后立即放入装有湿纸的密封袋中冷藏保存，以最大程度减少水分损失，并迅速运回实验室进行后续处理。

步骤三：功能性状测定 本研究共计测定了25个功能性状，分为三类： 1. 树皮性状（6个）：从采集的小枝上测量。首先测量带皮和去皮（仅木质部）后小枝的直径与体积，以此计算树皮厚度和树皮体积。树皮样品烘干后称重，计算总树皮密度。使用元素分析仪和连续流动分析仪测定树皮干样中的碳、氮、磷含量，并计算碳氮比。 2. 木材性状（7个）：从小枝的木质部样品中获得。木材密度的计算方法与树皮密度类似。采用滑动切片机将小枝横切成厚度为20微米的薄片，经染色脱水后，在光学显微镜下观察并拍照。利用Photoshop软件分析图像，测量并计算导管直径和导管密度。同样测定木材的碳、氮、磷含量和碳氮比。 3. 叶片性状（12个）：测量鲜重、饱和鲜重、叶面积（使用爱普生V39扫描仪）和干重。基于这些基础数据计算一系列性状：比叶面积（干重/面积）、叶片干物质含量（干重/饱和鲜重）、叶片组织密度（干重/(面积×厚度)）和叶片相对含水量。使用便携式光合作用测定系统（Li-6800, Licor）在自然光强和控温条件下，测量叶片的光合速率和暗呼吸速率。使用同位素质谱仪测定叶片的稳定碳同位素比率，该值可间接反映植物的水分利用效率。同样测定叶片的碳、氮、磷含量和碳氮比。

步骤四：土壤性质测定 为了解研究地点的环境背景，在各样地采用网格法采集了0-10厘米深度的土壤样品，测定其总碳、总氮和总磷含量。

步骤五：数据分析流程 研究采用了一套严谨、多维度的统计分析方法，以检验假说并揭示数据背后的模式： 1. 相关性网络分析：首先，采用Pearson相关分析检验所有25个性状之间的两两关系（排除同一组织内的性状），并将这些显著性关系（p < 0.05）构建成“性状相关网络”。使用Gephi软件进行可视化，节点大小表示某一性状与其他性状连接的数目（连接度），用以识别枢纽性状。 2. 主成分分析：分别对树皮性状、木材性状、叶片性状以及全部25个全株性状进行主成分分析。PCA是降维和识别主要变异轴的强大工具。在本研究中，PC1轴（解释变异最大的轴）通常被解释为反映植物从“获取策略”到“保守策略”的连续谱。通过计算各性状与PC1、PC2的相关性，可以判断每个性状在该经济学谱上的位置（例如，高氮含量通常与“获取”端关联）。通过比较常绿与落叶物种在PC1轴上的得分（t检验），可以检验它们在整体策略上的差异。 3. 标准化主轴回归分析：这是一种适用于研究两个变量间功能关系的回归方法，特别考虑了测量误差。本研究使用SMA分析来探究BES、WES、LES以及PES各自的PC1得分之间的数量关系（即它们的PC1轴如何相互关联），并检验常绿和落叶物种在这些关系上是否存在斜率或截距的差异。 4. 多元因子分析：用于同时评估树皮、木材、叶片三组性状集之间的整体关联强度。

三、 主要研究结果

1. 性状间的广泛关联与枢纽性状：相关网络分析清晰地显示，树皮、木材和叶片的大多数功能性状之间存在显著的相关性。更关键的是，氮和磷养分性状（如树皮磷含量、叶片磷含量、叶片氮含量、木材氮含量等）在网络中具有最高的“连接度”。这意味着它们是连接不同器官性状网络的核心枢纽。此外，叶片组织密度也与众多其他性状高度相关。这一结果为假说1和假说2提供了强有力的支持，表明不同器官的性状确实协同变化，而N、P养分代谢可能是驱动这种协同的关键生理过程。

2. 各器官及全株水平经济学谱的构建：

   • 树皮经济学谱：BES的PC1轴解释了48.1%的变异。树皮氮含量、磷含量（反映活跃的代谢和运输功能）位于“获取”端，而树皮厚度、碳含量和碳氮比（反映结构投资和防御）位于“保守”端。然而，常绿和落叶物种在BES的PC1轴上没有表现出显著差异。
   • 木材经济学谱：WES的PC1轴解释了44.8%的变异。导管密度高、木材碳氮比高（可能反映更安全但效率较低的水分运输和更耐用的结构）位于“保守”端，而导管直径大、木材氮磷含量高（反映高效的水分和养分运输）位于“获取”端。常绿物种显著地聚集在保守一端，而落叶物种则聚集在获取一端。
   • 叶片经济学谱：LES的PC1轴解释了36.2%的变异，结果符合经典的叶片经济学谱理论。比叶面积高、叶片组织密度高、相对含水量高、碳含量高的性状组合（通常与叶片寿命长、光合能力中等相关）位于“保守”端；而光合速率、暗呼吸速率、氮磷含量高的性状组合（与快速生长相关）位于“获取”端。同样，常绿物种显著偏向保守策略，落叶物种显著偏向获取策略。
   • 全株植物经济学谱：整合所有25个性状的PES，其PC1轴解释了28.1%的变异。PCA图清晰地显示，常绿物种和落叶物种分别聚集在谱系的两端，支持了假说3。

3. 不同经济学谱间的数量关系及叶习性影响：SMA回归分析揭示了BES、WES、LES和PES的PC1轴之间存在显著的相关关系，表明它们是相互整合的。然而，这种整合关系在常绿和落叶物种中存在有趣的异同：

   • 共性：在大多数关系中（如叶片PC1 vs. 树皮PC1， 全株PC1 vs. 木材PC1等），常绿和落叶物种的回归斜率没有显著差异。这表明，对于这两个类群，不同器官策略之间的比例变化关系是相似的，即存在一个普遍的整合框架。
   • 差异：在某些关系中，两类物种表现出不同的关联模式。例如，仅在落叶物种中，木材PC1与树皮PC1呈显著负相关；而仅在常绿物种中，叶片PC1与木材PC1呈显著正相关。此外，在叶片PC1与树皮PC1、全株PC1与树皮PC1的关系中，两类物种具有相同的斜率但截距不同。这些截距差异意味着，在相同的比例关系下，常绿和落叶物种达到了不同“绝对值”的策略组合。例如，在给定的树皮策略下，常绿物种总是比落叶物种表现出更保守的叶片策略（更高的截距）。

四、 研究结论与意义

本研究的核心结论是：成功构建了一个包含树皮功能性状的、更为完整的亚热带木本植物全株经济学谱。研究证实： 1. 整合性：树皮、木材和叶片的性状并非独立，而是通过紧密的关联与权衡，形成了一个协调的全株适应策略框架。其中，氮和磷的代谢与分配扮演了连接各器官的“枢纽”角色。 2. 策略分化：常绿物种和落叶物种在该PES上表现出明显的生态策略分化。常绿物种整体上采用以高结构投资、低营养代谢速率、长寿命器官为特征的“保守型”策略，以适应持续的环境压力（如低温、贫瘠）。而落叶物种则采用以高资源获取速率、快速生长为特征的“获取型”策略，以便在生长季内最大化资源利用。 3. 树皮功能：研究将树皮这一长期被忽视的多功能器官正式纳入了经济学谱框架，揭示了其性状（如养分含量）在整体策略中的重要作用。

科学价值：这项研究极大地推进了我们对植物整体适应策略的理解。它超越了传统的叶片中心视角，强调了从全株水平考察性状间协调与权衡的重要性。研究结果支持了“植物作为一个整体对其环境做出协调响应”的观点，并为此提供了来自亚热带森林生态系统的实证数据。提出的“氮磷作为枢纽性状”的观点，为理解植物内部资源分配与整合的生理机制提供了新的视角。

应用与理论意义：该PES框架有助于预测不同策略的植物如何响应全球变化（如气候变化、氮沉降），以及这些响应如何进一步影响生态系统的结构和功能（如碳循环、养分循环）。研究结果对于理解物种共存机制（不同策略的物种如何分享生态位）、森林动态模拟以及生物多样性保护都具有重要参考价值。例如，了解常绿与落叶物种的策略差异，有助于预测在环境变化下森林群落组成的可能演变方向。

五、 研究亮点

1. 视角新颖，填补空白：这是为数不多系统性地将树皮功能性状整合进全株植物经济学谱的研究之一，弥补了该领域的一个重要短板。
2. 假说驱动，验证严谨：研究基于清晰的科学假说展开，并运用了多种先进的统计方法（相关网络、PCA、SMA等）从不同角度进行交叉验证，结论坚实可靠。
3. 枢纽性状的识别：通过相关网络分析，明确指出了氮、磷养分性状在整合不同器官功能中的核心作用，为后续机制研究指明了方向。
4. 揭示叶习性的整体策略分化：不仅证实了常绿与落叶物种在叶片水平上的策略差异，更将这种差异扩展到了木材和全株水平，深化了对不同生活史类型植物适应策略的认识。
5. 研究设计周全：样地设置覆盖海拔梯度，物种选择具有代表性（59种），性状测定全面（25个），数据质量高，为研究结论的普适性提供了良好基础。

六、 其他有价值的讨论

研究在讨论部分还对一些有趣的现象进行了解释。例如，为何常绿和落叶物种在树皮经济学谱上没有显著差异？作者推测，可能的原因包括：（1）本研究的树皮样品来自小枝，其性状可能无法完全代表主干树皮的特性；（2）树皮作为保护组织，其性状（尤其是外层树皮厚度）可能更多地受局地干扰（如火灾频率）等特定环境因素驱动，而在本研究相对温和的亚热带山地环境中，这种选择压力可能不强，导致两类物种在树皮策略上趋同。这一思考为未来研究提出了新的问题，例如需要比较不同器官（主干 vs. 枝条）的树皮性状，以及在更强环境胁迫下（如干旱区、火干扰频繁区）的树皮经济学谱。