龙胆属(Gentiana)十字组(sect. Cruciata)药用植物在中国气候变化下的潜在适生区预测
《Diversity》:Prediction of Potential Suitable Areas of Medicinal Plants of Gentiana Sect. Cruciata in China Under Climate Change
编辑推荐:
龙胆属(Gentiana)十字组(sect. Cruciata)植物常作为大宗药材使用,具有多种临床治疗功能,但因过度采挖与环境变化野生种群急剧衰退,亟需保护。研究人员收集了该组8个物种的分布点记录,整合42个环境变量,利用“kuenm” R包(v. 3.6.
龙胆属(Gentiana)十字组(sect. Cruciata)植物常作为大宗药材使用,具有多种临床治疗功能,但因过度采挖与环境变化野生种群急剧衰退,亟需保护。研究人员收集了该组8个物种的分布点记录,整合42个环境变量,利用“kuenm” R包(v. 3.6.3)优化最大熵模型(MaxEnt),预测当前及未来气候情景(SSP126与SSP585)下潜在适生区,并分析空间格局变化与分布重心迁移。结果显示所有优化模型平均AUC(受试者工作特征曲线下面积)>0.9,预测精度高。海拔(alt)、年降水量(bio12)、人类足迹指数(hfp)与温度季节性(bio4)主导物种分布,不同物种限制因子组合揭示种间生态位分化。8个龙胆物种呈三类气候响应:气候敏感型G. macrophylla遭受严重不可逆生境丧失且重心远距离迁移;5个稳定适应型(G. straminea、G. crassicaulis、G. tibetica、G. robusta、G. siphonantha)维持稳定核心区并具局部扩张能力;G. dahurica与G. officinalis呈“收缩—恢复”轨迹。青藏高原东部、横断山区、秦岭及甘青川交界为核心保护区。高适应性G. crassicaulis、G. tibetica与G. straminea优先纳入保护与人工繁育。本研究揭示该组药用植物对气候变化的差异化响应与适应潜力,为精准种质保护与圆叶秦艽(Gentianae Macrophyllae Radix)可持续利用提供科学依据。
研究背景方面,山地与高山植物因气候变暖呈现快速向上迁移趋势,但针对高山药用植物(特别是区域特有类群)在不同历史时期气候条件下的动态分布研究仍较匮乏。龙胆属(Gentiana)十字组(sect. Cruciata)含多种秦艽(Gentianae Macrophyllae Radix)基原植物及藏药“解吉(Jieji)”,其野生资源因过度采挖与生态退化急剧下降,G. macrophylla已被列为国家重点保护野生药材。既往物种分布模型(SDM)研究多仅纳入气候变量、忽略人类干扰与土壤属性,且多为单种建模,难以系统揭示该类群对气候变化的综合响应。为此,研究人员以8个具药典或民间药用价值的十字组物种为对象,整合地形、气候、土壤与人类足迹等多源环境数据,采用优化最大熵模型(MaxEnt)开展当前与未来(SSP126、SSP585,2041—2080)适生区模拟与比较分析。
关键技术方法上,研究人员从CVH、GBIF、NSII及野外调查获取8种共2132条分布记录,经去重与5 km空间薄化(thinning)保留独立样点;环境层包含WorldClim v2.1的19个生物气候变量、ArcGIS衍生的海拔(alt)/坡度/坡向、HWSD的19个土壤变量及2009年人类足迹指数(hfp),统一重采样至2.5弧分。利用R包“kuenm”对特征类(FC)与正则化乘数(RM)进行1240种组合交叉验证,按ΔAICc=0选取最优模型;以10次重复平均输出栖息地适宜性指数(HSI),依最大检验敏感度加特异度阈值(MTSPS)二值化,划分不适生/一般/中等/高适生四级,并在ArcGIS与SDMTools中叠加各时期二值栅格量化扩张、稳定与收缩面积,计算分布重心坐标。
研究结果部分,2.1模型优化与精度评估显示8个物种最优MaxEnt模型平均AUC介于0.908±0.008(G. macrophylla)至0.973±0.005(G. robusta),遗漏率<0.05,证实参数优化显著提升转移性与稳定性。2.2主导环境因子表明海拔(alt)、年降水量(bio12)、人类足迹(hfp)、等温性(bio3)、温度季节性(bio4)及最湿月降水(bio13)等为关键变量,但组合因种而异:G. macrophylla受hfp+bio12+alt+bio1驱动;G. straminea、G. crassicaulis、G. siphonantha以alt为核心;G. robusta与G. tibetica依赖bio3与bio4;G. tibetica唯一显著受表土盐基饱和度(t_bs)限制。2.3当前与未来适生区面积显示G. macrophylla当前134.25×104 km2,SSP585-2061—2080大幅收缩;G. straminea、G. crassicaulis、G. tibetica在两种情景下均具扩张趋势;G. dahurica与G. officinalis呈中期收缩、后期回升。2.4空间格局变化指出G. macrophylla稳定区比例降至5.83%、收缩达8.16%;稳定适应型不变区占比多>8.6%,G. crassicaulis达12.16%以上且扩张最多;波动调整型在2061—2080由净丧失转为净扩张。2.5重心迁移显示G. macrophylla重心西北—西南迁移超200 km;其余物种重心摆动<50 km或限于原核心区(如G. robusta、G. siphonantha、G. tibetica)。
讨论部分总结,模型优化克服默认参数过拟合缺陷,AUC>0.9验证可靠性;海拔与降水季节性为高山龙胆分布骨架,hfp兼具轻牧干扰与采样偏倚双重含义;种间生态位分化体现为对地形、水分季相与土壤因子的不同侧重。响应类型归纳为气候敏感型(G. macrophylla)、稳定适应型(5种)与波动调整型(G. dahurica、G. officinalis),差异源于生态位宽度、地形缓冲与排放情景。保护上建议将青藏高原东部、横断山、秦岭、甘青川交界设为核心种质区,对G. macrophylla实施迁地保存与高稳定区再引入,对G. crassicaulis、G. tibetica、G. straminea推进川西、滇西北、藏东人工栽培,并耦合代谢组学评价生境—药效关联。
结论部分翻译:本研究利用R优化的MaxEnt模型系统模拟了龙胆属十字组8种药用植物在当前与未来气候情景下的潜在适生区,通过分析空间格局变化与重心迁移厘清了种间气候适应潜力与未来分布趋势的差异,评估了资源系统的气候脆弱性并识别了关键适生区与高适应物种。海拔、年降水量、人类足迹与温度季节性被确定为塑造8种龙胆分布的首要环境变量,种间存在明显生态位分化。当前气候下各物种共享重叠分布同时保留专属适生区;其对未来气候的差异化响应分为三类:以G. macrophylla为代表的气候敏感型,由G. straminea、G. crassicaulis、G. tibetica、G. robusta与G. siphonantha组成的稳定适应型,以及由G. dahurica与G. officinalis构成的波动调整型。结果为靶向就地保护与原位种质保存提供了坚实理论证据,并为中药资源的可持续利用给予重要指导。论文发表于《Diversity》。