气候驱动的中国铁皮石斛生境适宜性与植物化学组成变异

《Agronomy》：Climate-Driven Variation in Habitat Suitability and Phytochemical Composition of Dendrobium officinale in China

　　

铁皮石斛（Dendrobium officinale）是一种具有高治疗价值的濒危药用兰花。了解其适宜生长区域及药用化合物积累最多的地点，对保护与栽培均具有重要意义。本研究中，研究人员利用物种分布模型（MaxEnt）预测了中国铁皮石斛的适宜生境，并结合野外采样和紫外-可见分光光度法（UV-Vis spectrophotometry），分析了贵州省四个产区植物中生物活性化合物的含量。研究结果表明，铁皮石斛的潜在分布主要受水分相关变量约束，尤其是最干季度降水量、年降水量、水汽压以及植被活动（归一化植被指数（NDVI））。总适宜面积估计约为21.4×10⁴ km²。对贵州省四个地点样品的植物化学分析提示，多糖、生物碱和总酚的积累可能对环境因素产生不同响应。生物碱和总酚与降水相关变量（bio12, bio17）呈强负相关，而多糖与这些变量呈正相关。然而，这些相关性基于小样本（n=4），应仅作为针对贵州省的探索性假设，而非确立的关系。它们指出，干旱季水分可利用性可能对不同类别次生代谢产物产生相反效应，但该推测需扩大采样加以验证。总体而言，本研究提供了一个空间框架，可为栽培选址和保护提供参考，而非一个完全验证的精准农业工具。

**论文解读：气候驱动的铁皮石斛生境适宜性与植物化学组成变异**

**研究背景与问题**

铁皮石斛（*Dendrobium officinale*）是兰科多年生附生草本植物，兼具观赏、药用及食用价值，在传统中医中被视为珍贵药材，其多糖、生物碱、总酚等活性成分具有抗炎、抗氧化、免疫调节等药理作用。然而，由于野生资源濒危（IUCN列为濒危物种），人工栽培成为主要来源，但不同产区品质差异显著，普遍存在生境适宜性下降、产量低、活性成分含量不足等问题。根本原因在于栽培基地缺乏科学选址和标准化管理，环境变量如何同时调控大尺度生境适宜性和局地植物化学品质的关键知识空白尚未填补。因此，本研究旨在通过整合物种分布模型与化学分析，探索环境因子对铁皮石斛分布及活性成分积累的联合调控机制，为保护与可持续利用提供科学依据。论文发表在《Agronomy》。

**主要技术方法**

研究人员主要采用以下关键技术：
1）**物种分布模型（MaxEnt）**：基于128个野生种群分布记录（经空间稀化至1 km分辨率，数据来源包括GBIF、CVH、NSII等数据库，1980–2020年），结合26个环境变量（19个生物气候变量、海拔、坡度、坡向、太阳辐射、水汽压、NDVI、人类足迹指数），通过Pearson相关分析（|r|≥0.85）筛选11个生物气候变量，并利用方差膨胀因子（VIF）评估多重共线性，最终保留18个预测因子。
2）**模型优化与评价**：使用“kuenm”包进行参数优化（测试正则化乘子0.5–4，29种特征组合），以部分ROC、5%遗漏率及AICc为标准选择最优模型（RM=4, FC=LH），采用空间分块交叉验证（10次重复），最终模型AUC=0.899，TSS=0.78。
3）**植物化学分析**：于2020年11月采集贵州省安龙、松桃、平塘、德江四个栽培基地的样品（每个基地3个10 m×10 m样方，每样方随机采5株），采用80%乙醇提取后，利用紫外-可见分光光度计（UV-Vis）按2020年版《中国药典》测定多糖、生物碱、总酚含量（三次技术重复），并通过单因素方差分析（ANOVA）和LSD事后检验比较区域差异，Pearson相关分析评估化合物含量与环境变量的关系。

**研究结果**

**3.1 模型优化与性能验证**
通过128个分布记录和筛选变量，确定最优参数组合（RM=4, FC=LH），最终模型AUC=0.899（良好），TSS=0.78，表明模型预测能力优良。

**3.2 当前铁皮石斛适宜生境分布**
模型预测适宜生境主要集中在中国南方和东南部，包括台湾省。最高适宜区位于湖南、广东、广西、福建、浙江。总适宜面积约21.4×10⁴ km²，其中高适宜区34,100 km²（15.94%），中适宜区132,100 km²（61.76%），低适宜区47,700 km²（22.30%）。

**3.3 影响铁皮石斛资源分布的主要环境因子**
Jackknife分析显示，最干季度降水量（bio17）是主导因子，贡献率83.3%；其次为年降水量（bio12, 10.6%）、NDVI（1.6%）和水汽压（VAPR, 1.4%）。响应曲线表明适宜阈值：bio17≥38.86 mm，bio12≥885.13 mm，NDVI≥129，VAPR≥1.13 kPa。总体而言，降水、植被覆盖和大气湿度是分布的主要决定因素。

**3.4 贵州省铁皮石斛生物活性化合物谱的区域变异**
四个地区化合物含量差异显著（p<0.05）。安龙样品的水分、醇提物、黄酮、总酚、生物碱含量最高；松桃醇提物和生物碱最低，但多糖含量最高（305.06 mg/g）；平塘各项中等；德江多糖最低（68.25 mg/g）。表明地理起源显著影响植物化学谱。

**3.5 生物活性化合物对季节性降水模式的差异化响应**
Pearson相关分析显示：生物碱与年降水量（bio12, r=-0.89）、最干季度降水量（bio17, r=-0.99, p=0.001）及人类足迹指数（HF, r=-0.99, p=0.001）呈强负相关；总酚也呈类似负相关（bio17: r=-0.97, p=0.017）；而多糖则与这些变量呈正相关（bio17: r=0.90, p=0.10；HF: r=0.99, p=0.001）。NDVI与所有化合物相关性较弱（|r|<0.45）。基于Cook's距离分析，无异常值驱动。这些模式提示干旱季水分可利用性可能对初级代谢与次级代谢产生相反调控效应，但因样本量小（n=4），应视为探索性假设。

**总结与讨论**

讨论部分强调：模型预测性能良好，但受限于数据质量与空间尺度；核心适宜区与传统道地产区（福建、浙江）部分重叠，差异可能源于发生数据更全面、新增NDVI与VAPR变量以及参数配置差异。化学组分与环境因子的相关分析揭示，生物碱和总酚与降水变量负相关，多糖正相关，暗示干旱胁迫可能促进防御性次生代谢产物积累，而水分充足利于多糖合成。但需注意，化学数据来自栽培种群，MaxEnt模型基于野生分布记录，两者环境信号可能因管理干预而不同。人类足迹指数与降水变量相似的相关模式可能反映人类活动区多位于水分充足的河谷低地，而非直接因果。未来需扩大采样、整合土壤与微气候测量，并开展生理与分子机制研究。

**研究结论翻译**
本研究中，研究人员探索了环境因素如何影响铁皮石斛的生境适宜性和药用品质。结果表明，该物种在中国的潜在分布主要由少数关键生物气候变量决定，尤其是最干季度降水量和年降水量。研究人员还发现，生物碱、总酚和总黄酮含量与降水相关变量呈显著负相关，而多糖与这些变量呈正相关。这种相反模式表明，干旱季水分可利用性可能对铁皮石斛的初级与次级代谢产生差异调控效应。然而，由于植物化学采样仅局限于贵州省四个地点，这些气候-化合物关系应视为该区域的初步假设，而非全国性结论。一个有趣的发现是冷干季的干旱胁迫可能促进生物碱和酚类积累，这指向环境胁迫与次生代谢之间的复杂相互作用。通过结合经过验证的国家尺度生境适宜性模型（基于128个分布记录和空间交叉验证）与探索性植物化学数据（来自贵州四个地点），研究人员识别出高生境适宜性与高生物活性化合物积累可能共存的检验假设。这尚未成为精准栽培的验证框架，而是概念验证，展示了如何利用物种分布模型（SDM）输出为靶向植物化学采样生成假设。未来需扩大空间覆盖范围以验证这些初步相关性。然而，这一“双标准”构想目前仍属推测，因为植物化学部分缺乏预测制图所需的空间覆盖和统计效力。该框架为未来精准栽培研究提供了起点，但需通过扩大采样进一步验证。此外，由于化学数据来自栽培种群而MaxEnt模型基于野生分布记录，整合制图（图10）应视为概念性图示而非直接预测工具。展望未来，值得探索的方向包括：扩大植物化学采样的地理范围与环境梯度；研究水胁迫等环境信号与次生代谢产物生产的生理与分子机制；开发包含土壤、生物相互作用等更全面的模型；通过田间栽培试验验证空间预测。总之，这一濒危药用植物的可持续保护与利用，有赖于从传统经验栽培向数据驱动、气候信息导向、品质导向的方法转变，以在变化的世界中保护其珍视的“道地”品质。