环境塑造灰毡毛忍冬适宜生境与品质：基于最大熵模型、高效液相色谱、化学计量学与基因表达分析的见解

《Plants》：Environmental Shaping Suitable Habitats and Quality of Lonicera macranthoides Hand.?Mazz.: Insights from MaxEnt, HPLC, Chemometrics, and Gene Expression Analysis Nan Xu, Canfeng Long, Meixin Zhou, Weijia Wang, Qiang Zeng, Yingying Shen, Pan Wu, Liqun Rao, Guoping Peng and Qiming Wang

　　

灰毡毛忍冬（Lonicera macranthoides Hand.?Mazz.）是中国珍稀药用植物，在全球范围内均有应用。本研究旨在采用最大熵（MaxEnt）模型预测其中国境内适宜生境，并通过高效液相色谱（HPLC）与化学计量学评估环境变量对指标成分（绿原酸、灰毡毛忍冬苷B、川续断皂苷B）的影响。此外，为探究环境与品质关联的分子机制，研究人员开展了室内外胁迫预实验，采用实时荧光定量PCR（qRT?PCR）分析绿原酸生物合成基因的表达特征。结果显示，最干月降水量是影响分布的最关键变量。当前适宜区主要位于北纬21°至33°之间。末次盛冰期（LGM）生境范围更广，而中全新世（MH）发生收缩。未来情景预测表明，SSP585情景下将出现生境丧失，而SSP126情景下该趋势在2090年可得到部分缓解。适宜生境中绿原酸与皂苷含量更高，且与土壤、海拔及降水显著相关。值得注意的是，室外复合胁迫（低温与低光照）显著调控LmPAL、LmCHS、LmCHI、LmC4H、LmCCoAOMT及LmANS的表达。本研究为中国灰毡毛忍冬的保护、可持续栽培及胁迫导向育种提供了科学依据。

论文解读

灰毡毛忍冬是中国南方重要的药用植物，其干燥花具有较高的经济价值，现代药理研究证实其具有免疫调节、抗炎抗菌、降糖降脂等功效。目前野生资源较丰富但质量不均，药材来源以栽培为主，主产区包括湖南、贵州、重庆、四川、云南等地。不同地理单元会影响药用植物的化学成分，进而决定其品质，《中国药典》以绿原酸、灰毡毛忍冬苷B、川续断皂苷B作为品质评价指标。现有研究多聚焦于灰毡毛忍冬的化学成分分析、药理学及组学层面关键基因与代谢物挖掘，针对其物种分布的研究较少，且气候变化影响其次生代谢产物的分子机制尚不明确。同时，植物在自然生境中常面临低温、干旱或光照波动等复合胁迫，这些环境信号如何在分子层面调控次生代谢，进而影响最终药材品质，仍需系统解析。为此，研究人员整合生态位建模、植物化学分析与分子生物学技术，开展灰毡毛忍冬适宜生境预测、品质与环境关联分析及分子机制探索，研究成果发表于《Plants》。

为开展该研究，研究人员采用了四个核心技术体系：一是基于最大熵（MaxEnt）模型的物种分布模拟，整合34份实地采样样本与329条公开分布记录，筛选21个关键环境变量，以AUC（受试者工作特征曲线下面积）与TSS（真实技巧统计）评估模型精度，划分适宜生境等级并分析不同气候情景下的分布动态；二是高效液相色谱（HPLC）结合蒸发光散射检测器（ELSD）测定不同生境样品中三类指标成分的含量；三是采用层次聚类分析（HCA）、主成分分析（PCA）及相关性分析等化学计量学方法，解析生境适宜性与成分积累的关联；四是开展室内外复合胁迫实验，通过实时荧光定量PCR（qRT?PCR）检测苯丙烷通路关键基因的相对表达水平。

2.1 模型精度评估

研究人员共纳入34份代表性产区样本与329条公开分布记录构建MaxEnt模型，选取21个环境变量以降低过拟合风险。结果显示，各重复运行的平均AUC值均大于0.95，TSS值介于0.842至0.856之间，符合AUC>0.9、TSS>0.8的高精度预测标准，表明模型可准确模拟灰毡毛忍冬的适宜生境分布。

2.2 环境变量分析

影响灰毡毛忍冬分布的关键变量共8个，分别为最干月降水量（Bio_14）、全球干旱指数（Ai_v3_yr）、海拔（Elev）、温度季节性（Bio_4）、表土碳酸钙含量（T_caco3）、坡度（Slope）、最低月平均UV?B辐射（UVB4）及土地覆盖类型（Gm_lc_v3），其贡献率与置换重要性合计分别达93.3%和57.9%，其中Bio_14的贡献率最高。响应曲线显示，除T_caco3与Gm_lc_v3外，其余变量的适宜范围均呈单峰正态分布，明确了各环境因子的适宜区间与最优值，例如最干月降水量适宜范围为0~190 mm，最优值为50.31 mm。

2.3 当前气候情景下的分布特征

当前气候条件下，灰毡毛忍冬主要分布于中国南方北纬21°~33°区域，集中分布在北纬23°~31°之间。总适宜生境面积为659.307×10³km²，其中高适宜区面积207.78×10³km²，集中于黔东、湘西、渝东南、鄂西南、桂北、粤北、闽东北、浙南及四川盆地东部；中适宜区面积451.52×10³km²，分布范围更广；低适宜区面积329.44×10³km²，呈扩散状环绕高中适宜区。

2.4 过去气候情景下的分布变化

末次盛冰期（LGM）的适宜生境面积达849.23×10³km²，较当前增加189.93×10³km²，分布范围向北、向南扩张，向西、向东收缩。中全新世（MH）适宜生境面积降至613.28×10³km²，较LGM减少235.42×10³km²，整体向南、向东、向西扩张但总面积收缩，呈带状分布。

2.5 未来气候情景下的分布预测

SSP126情景（低排放）下，2050S、2070S适宜生境面积先降后升，2090S面积达766.14×10³km²，恢复至接近LGM水平且高于当前。SSP585情景（高排放）下，适宜生境面积持续下降，2090S仅为351.06×10³km²，高适宜区缩减最为显著，未出现缓解趋势。

2.6 不同时期分布格局变化

当前分布质心位于湖南省祁东县（111.84° E，26.89° N）。LGM时期质心较当前向东北移动50.25 km，MH时期向东南移动70.84 km，整体随温度变化呈现“先南移后北移”的规律。未来情景下，分布质心总体向更寒冷的北方地区迁移，且CO₂排放越高，北移距离越远。

2.7 指标成分含量分析

所有样品的指标成分均符合《中国药典》要求，高、中适宜区的绿原酸、灰毡毛忍冬苷B、川续断皂苷B含量均显著高于低适宜区。遵义样品绿原酸（53.61 mg·g^?1）与川续断皂苷B（6.50 mg·g^?1）含量最高，怀化样品灰毡毛忍冬苷B含量最高（90.43 mg·g^?1）；长沙样品绿原酸（29.23 mg·g^?1）与灰毡毛忍冬苷B（53.77 mg·g^?1）含量最低，昭通样品川续断皂苷B含量最低（3.67 mg·g^?1）。

2.8 层次聚类、主成分与相关分析

层次聚类（HCA）将34份样品分为两类，低适宜区样品聚为一类，适宜区样品聚为另一类，与环境差异高度吻合。主成分分析（PCA）前两主成分累计解释方差99.75%，可清晰区分三类生境样品，低适宜区样品差异尤为显著。相关性分析显示，绿原酸积累主要受土壤（表土黏土含量、表土碳酸钙含量、底土pH、表土电导率）、海拔与最干月降水量调控；两类皂苷积累主要与土壤因子、土地覆盖类型及土壤有效水含量相关。

2.9 室内外环境下的基因表达模式

以室内环境（25 ℃，20000 lux）为对照，室外复合胁迫（8 ℃，低日照）下，灰毡毛忍冬苯丙烷通路关键基因LmPAL、LmCHS、LmCHI、LmC4H、LmCCoAOMT的表达量显著下调，LmANS表达量显著上调，表明低温与低光照复合胁迫可显著调控苯丙烷通路的转录水平，在分子层面影响次生代谢。

讨论与结论

讨论部分指出，最干月降水量是影响灰毡毛忍冬分布的首要因子，其贡献率达64.2%，体现了该物种对干旱胁迫的适应性，但充足的水分供应仍是必要条件，这与药用植物分布受降水驱动的规律一致。地形因子（海拔、坡度）通过影响水热分配与土壤结构间接调控分布，UVB4的贡献则印证了其喜光的生物学特性。品质分析证实适宜生境可促进活性成分积累，低适宜区成分含量较低，而适度环境胁迫可能进一步诱导次生代谢产物合成，该现象在其他药用植物中也有报道。分子机制层面，复合胁迫下苯丙烷通路上游基因下调、LmANS上调，可能与MYB?bHLH?WD40复合物特异性调控花青素分支有关，揭示了环境信号调控品质形成的分子基础。过去与未来分布模拟表明，温度是驱动生境变迁的核心因素，LGM低温促进分布扩张，MH高温导致收缩；未来高排放情景下生境将持续丧失，低排放情景下可在2090年逐步恢复，且分布质心总体北移。基于研究结果，研究人员提出保护策略：优先保护南方山地野生种质资源，建立种质保存库；选择海拔约500 m的向阳坡地种植，改良沙质土壤并完善排水设施；未来可基于胁迫响应分子机制开展抗逆优质品种选育。

研究结论可归纳为：本研究整合MaxEnt模型、HPLC、化学计量学与qRT?PCR技术，系统阐明了灰毡毛忍冬适宜生境的关键环境变量（最干月降水量、全球干旱指数、海拔、温度季节性、表土碳酸钙含量、坡度、最低月平均UV?B、土地覆盖类型）及其适宜区间。当前适宜区位于中国南方北纬21°~33°，LGM生境扩张、MH收缩，未来高排放情景下生境持续丧失，低排放情景下2090年可部分恢复。适宜生境中绿原酸、灰毡毛忍冬苷B、川续断皂苷B含量更高，积累过程受土壤、海拔、降水调控。室外低温低光照复合胁迫显著调控LmPAL、LmCHS、LmCHI、LmC4H、LmCCoAOMT、LmANS的表达，从转录水平揭示了环境胁迫对苯丙烷通路及次生代谢的调控作用。该研究构建了多学科融合的研究框架，为灰毡毛忍冬的资源保护、可持续栽培及抗逆育种提供了科学依据。