云南松矮生变种(Pinus yunnanensis var. pygmaea)个体在火干扰频发区域更倾向于表现出"迎火(fire-embracing)"适应策略

《Functional Ecology》：Pinus yunnanensis var. pygmaea individuals are more ready to embrace wildfire in more fire-prone areas

　　

摘要：植物功能性状(functional traits)在响应环境异质性中发挥关键作用，并为理解物种适应策略提供理论基础。云南松矮生变种(Pinus yunnanensis var. pygmaea)是中国西南火干扰频发区特有的高韧性物种，其在包括球果迟裂(serotiny)在内的众多功能性状上存在显著种内变异，然而这些功能性状如何演化形成响应异质火干扰制度(fire regime)的适应策略尚未得到充分探讨。研究人员调查了14个P. yunnanensis var. pygmaea种群的25项功能性状，筛选出与球果迟裂显著关联的功能性状并基于其间权衡(trade-offs)构建了迎火指数(fire-embracing index)，随后探究了潜在驱动迎火指数变异的环境因子。结果表明，迟裂程度较高的种群表现出更显著的迎火策略特征——较低叶氮浓度(leaf nitrogen concentration)、较矮植株高度及较高球果粗脂肪含量(cone crude fat content)，此类种群倾向于分布在过火面积(burned area)较大且土壤养分贫瘠的区域。研究结果揭示了P. yunnanensis var. pygmaea迟裂性状与其他功能性状间存在明显权衡关系，为其迎火策略的演化机制提供了新见解，有助于阐明植物如何在火干扰生态系统中发展适应策略。

论文解读：《Functional Ecology》——云南松矮生变种迎火策略的功能性状协同演化与火干扰-环境梯度驱动机制

研究背景与立项依据

植物功能性状(functional traits)是物种与环境互作及适应生物与非生物胁迫的关键指标，不同性状组合可形成促进植物存活与繁殖的生态策略，如叶经济谱(Leaf Economics Spectrum, LES)和Grime的CSR理论。野火是陆地生态系统普遍的选择压力，深刻塑造植物性状与群落结构，许多植物演化出与特定火干扰制度(fire regime)匹配的抗性性状组合。松树常划分为火耐受型(fire-tolerator，具厚树皮隔热、高位树冠)与迎火型(fire-embracer，依靠保护种子如迟裂serotiny或强萌蘖能力在火后再生)。球果迟裂(serotiny)指成熟种子不立即释放而滞留于闭合球果中跨过至少一个生殖周期，是迎火策略的核心特征，但其常与其他功能性状（如叶性状）协同构成适应综合征(syndrome)的机制仍不清楚。云南松(Pinus yunnanensis)为中国西南特有针叶树，该区野火频发，含三个变种，其中灌木型生态型云南松矮生变种(var. pygmaea)特生于火干扰干热河谷，迟裂水平高于其他变种，但种内种群间迟裂存在显著变异，这种变异尤其受火干扰与其他环境因子共同驱动的机制尚未阐明，迟裂如何协同其他功能性状及其权衡机制亦不明确。本研究旨在探明：(1)何种功能性状与球果迟裂共同支撑迎火策略；(2)火干扰及其他环境因子如何塑造种群水平迎火策略的变异。研究人员假设迟裂与结构及资源相关性状协同形成整合的迎火策略以促进火后更新；并假设该策略种内变异主要由火相关环境梯度驱动，火活动(fire activity, FA)更高种群表现更强迟裂及相关性状综合征。

主要技术方法概述

研究人员在中国西南云南松矮生变种自然分布区内选取14个天然种群（云南、四川），每种种群设置20 m × 20 m样方，选4株健康成熟个体测定树高(tree height, HT)，按迟裂比例采集开放与未开放球果，采一年生枝条南北冠层针叶；土壤用五点法取0～20 cm层。测定生长性状（HT）、叶性状[叶长(leaf length, LL)、叶鞘长(leaf sheath length, LSL)、叶束直径(leaf fascicle diameter, LFD)、叶直径(leaf diameter, LD)、叶面积(leaf area, LA)、叶含水量(leaf moisture content, LMC)、单位面积叶质量(leaf mass per area, LMA)、叶热值(leaf calorific value, LCV)、单宁(leaf tannin, LT)、纤维素(leaf cellulose, LC)、粗脂肪(leaf crude fat, LCF)、全碳LTC、全氮LTN、全磷LTP]、球果性状[迟裂水平cone serotiny level, CSL＝留存未开球果数/成熟球果数、球果长CL、宽CW、质量CM、灰分CAC、粗脂肪CCF、热值CCV、全氮CTN、全碳CTC、全磷CTP]及土壤理化性质(pH、容重SBD、重量含水率SGW、全碳STC、有机碳SOC、全氮STN、全磷STP、全钾STK)。环境数据提取19项生物气候变量(WorldClim)、火活动FA与火辐射功率FRP(MODIS/FIRMS, 2001–2024, 0.5°网格)、过火面积burned area, BA(MCD64A1)、火险天气指数Fire Weather Index, FWI(ECMWF)。统计先进行正态性与方差齐性检验，单因素ANOVA与Duncan多重比较检验种群间性状差异，Pearson相关分析性状关联，主成分分析(prcomp包)提取第一主成分PC1作为迎火策略指数(fire-embracing strategy index)，随机森林(randomForest)结合SHAP值(shapviz)量化环境因子相对重要性及方向。

研究结果

3.1 Variation in adopting fire-embracing strategy within P. yunnanensis var. pygmaea

几乎所有被测性状在14个种群间存在显著种内变异；CSL地理分布沿纬度、经度、海拔无单调趋势但区域间有中等变异。CSL与球果粗脂肪(CCF)、全磷(CTP)显著正相关，与球果宽(CW)、长(CL)、全氮(CTN)显著负相关；与叶热值(LCV)、单宁(LT)、含水量(LMC)显著正相关，与叶全氮(LTN)、叶束直径(LFD)、全磷(LTP)、LMA、叶长(LL)显著负相关；树高(HT)与CSL显著负相关。对CSL及其显著关联功能性状做PCA，前两轴解释83.6%总变异(PC1=76.9%, PC2=6.7%)，CSL与PC1显著正相关，故PC1被解释为迎火策略指数，高分值代表强迎火策略（高迟裂、矮株、低叶氮低LMA、高球果粗脂肪等），低分值为弱迎火策略。

3.2 Drivers of the fire-embracing strategy index within P. yunnanensis var. pygmaea

环境因子在种群间具异质性。迎火策略指数与BA、土壤pH、全钾(STK)、容重(SBD)显著负相关，与SOC、全碳(STC)、重量含水率(SGW)显著正相关。对环境显著关联变量做PCA，PC1解释75.6%变异——PC1高值对应高STK、SBD、pH、BA及低SOC、STC、SGW，可解释为火干扰偏强、干燥、贫养、高容重高pH土壤的综合环境指数(火频发/贫养端)。迎火策略指数与该环境PC1呈显著正相关(R²=0.897, p<0.001)：营养较丰富、火干扰较少区域(如ZY、FSC种群)迎火指数低(高叶氮、低迟裂)；火干扰更频发区域(如QF、YB、XJK种群)迎火指数高(矮株、薄叶、低叶养分、高迟裂)。随机森林SHAP分析显示土壤全钾STK( Mean |SHAP|=0.56)影响最大，依次为pH(0.48)、SOC(0.47)、STC(0.43)、SBD(0.38)、SGW(0.24)、BA(0.05)；STK与pH正向贡献(碱性贫钾贫养土迎火指数高)，SOC、STC、SGW负向贡献，BA亦正向贡献。表明土壤理化性质、水分有效性及火干扰共同塑造种内迎火策略变异，土壤属性解释力超过单一火干扰指标。

讨论与结论翻译总结

讨论指出P. yunnanensis var. pygmaea迟裂与矮株、高球果粗脂肪等性状协同变异(co-variation)，说明其适应策略源于多性状联合演化而非单性状变化，构成整合的迎火适应综合征，支持将宏观策略框架延伸至种内水平检测种群适应综合征。火不仅直接选择且通过改变土壤(降低有机质、改变pH与水分)间接重塑选择压，该物种迎火性状由火干扰与火频发生境非生物条件(贫养、偏碱、低湿)协同作用驱动，土壤特性解释力优于气候与火指标本身，印证植物适应是干扰与环境背景交互选择的结果。研究局限含未测光合、呼吸及地下性状，未区分遗传本地适应与表型可塑性，遥感火产品分辨率限制微尺度火动态捕捉。

结论(Conclusion)翻译：

P. yunnanensis var. pygmaea各种群间迟裂性状与其他功能性状存在明显权衡。表现更显著迎火策略（较高迟裂、较低叶氮、较低LMA、较矮株高、较高球果粗脂肪CCF及全碳CTC含量）的种群倾向于占据历史过火面积较大且土壤更贫养的区域。结果表明云南松矮生变种的迎火策略演化是火干扰制度及相关环境因子塑造多功能性状协同变异的结果。未来整合基因组工具研究此种内变异将深化对火干扰生态系统中该策略演化过程的理解。