植物生长过程中有效成分的含量与生长环境密切相关^[1]，生态因素影响中药材质量主要指气候、土壤、地质等生态因子影响药材的结构、性状和药效成分等，即环境饰变导致同种异地药材质量不同^[2]。研究土壤因子对药材药效成分的影响及两者相互关系，筛选药效成分的主要影响因子，能够为评价不同产地药材的质量差异提供依据。土壤因子与药材质量具有一定的相关性，沈晓凤等^[3]研究发现土壤因子有效性影响不同产地丹参药材品质；王利丽等^[4]发现土壤中铁、锰、钾元素显著影响山茱萸药材化学指标，其他土壤元素没有明显影响。

防风Saposhnikovia divaricata (Turcz.) Schischk为伞形科防风属多年生草本植物，以秋季未抽薹植株的根入药，主产于东北及华北地区。防风作为传统中药，在我国有悠久的使用历史，最早收录于《神农本草经》并列为上品，是常用的大宗药材之一^[5]，具有解热、镇痛、抗炎、抗菌、抗肿瘤、抗惊厥等积极的作用。近年来防风野生资源日益匮乏，防风药材需求量的增加使栽培防风成为市场主体，种植区域逐渐扩张^[6]，产地和生产方式发生了较大变化^[7]。2015年版《中国药典》^[8]规定防风以色原酮(包括升麻素苷、升麻素、5−O−甲基维斯阿米醇苷、亥茅酚苷)含量为评价标准。研究表明，不同产地防风的药效成分含量具有一定差异^[9]，这种差异可能与不同产地生态环境有密切关系。本论文旨在研究不同产地栽培防风药材质量差异以及土壤因子与药材质量的关系，通过相关性分析、聚类分析和多重线性回归分析探讨在防风药材质量形成过程中土壤因子是否起到关键作用，为合理栽培种植防风提供理论依据。

1 材料与方法 1.1 材料

2018年10月，在河北、山东、吉林、宁夏、甘肃、辽宁、内蒙古等省（自治区）的16个产地采集防风材料(表1)，16个产地分别为河北承德大屯乡(河北大屯乡)、河北承德金沟屯(河北金沟屯)、山东烟台牟平区(山东烟台)、宁夏固原隆德(宁夏隆德)、吉林长春农大药园(吉林长春)、甘肃临夏和政县(甘肃兰州)、吉林白城洮南(吉林白城)、辽宁朝阳东碾子沟(辽宁东碾子沟)、辽宁朝阳边杖子(辽宁边杖子)、辽宁朝阳黄台子(辽宁黄台子)、内蒙赤峰大明镇(内蒙大明镇)、内蒙赤峰唐家窝铺(内蒙唐家窝铺)、内蒙赤峰汐子镇(内蒙汐子镇)、内蒙赤峰牛营子(内蒙牛营子)、内蒙赤峰杨家营子(内蒙杨家营子)、内蒙赤峰文钟镇(内蒙文钟镇)。每个产地随机采集10株2年生防风及相对应的根际土壤，采集根际土壤是通过抖落的方法，挖取完整的防风根系土壤，将大块、不含根系的土壤抖落，收集附着于防风根系表面的土壤。采集的防风根系经自来水洗净，蒸馏水冲洗后于50 ℃烘干，称取质量，粉碎，备用。土壤样品烘干后过100目筛，储存备用。

1.2 方法 1.2.1 色原酮提取及含量测定

防风根系粉碎后过100目筛，精密称取0.25 g样品加入5 mL的70%(φ)乙醇中，在600 W、40 ℃条件下超声提取45 min，收集滤液。采用Agilent 1100高效液相色谱仪测定含量，色谱柱Diamonsil C₁₈(4.6 mm×200 mm×5 μm)，柱温30 ℃，检测波长254 nm。甲醇−水梯度洗脱：0~15 min，φ(甲醇)20%~45%；15~25 min，φ(甲醇)45%~70%。流速1 mL/min。根据峰面积计算防风4种色原酮(升麻素苷、升麻素、5−O−甲基维斯阿米醇苷、亥茅酚苷)含量^[9]。

1.2.2 土壤因子的测定

通过恒温箱烘烤的方法测定土壤水分，使用酸度计精准测定土壤样品的pH和电导率，使用重铬酸钾容量法(外加热法)测定土壤有机质(Organic matter)含量，通过半微量凯氏法对土壤全氮(Total nitrogen)含量进行消煮测定，通过碱解扩散法测定土壤碱解氮(Alkali-hydrolyzed nitrogen)含量，通过高氯酸−浓硫酸法测定土壤全磷(Total phosphorus)含量，通过碳酸氢钠法^[10]测定土壤有效磷(Available phosphorus)含量，通过火焰光度法^[11]测定土壤有效钾(Available potassium)、有效钙(Available calcium)、有效镁(Available magnesium)、有效铁(Available iron)、有效铜(Available copper)、有效锌(Available zinc)和有效锰(Available manganese)含量。

1.3 数据处理

采用Microsoft Excel 2013和SPSS 25.0统计软件对数据进行整理和分析，相关性分析用Pearson相关系数法,聚类分析和多重线性回归分析分别采用 SPSS 25.0 软件中的聚类分析法和多重线性回归分析法。

2 结果与分析 2.1 防风色原酮含量

防风4种色原酮含量以及色原酮总量的测定结果见表2。不同产地升麻素苷含量、升麻素含量、5−O−甲基维斯阿米醇苷含量、亥茅酚苷含量、色原酮总量、单株产量和单株色原酮总产量均具有显著差异。河北大屯乡、河北金沟屯、山东烟台、吉林白城、内蒙文钟镇色原酮总量较高，宁夏隆德的4种色原酮含量均普遍偏低，未达到药典标准(升麻素苷和5−O−甲基维斯阿米醇苷质量分数大于2.4 mg/g)^[8]。升麻素苷含量最高的产地为河北大屯乡，最低为宁夏隆德；升麻素含量最高的产地为山东烟台，最低为宁夏隆德；5−O−甲基维斯阿米醇苷含量最高的产地是吉林白城，最低为宁夏隆德；亥茅酚苷含量最高的产地为内蒙大明镇，最低为内蒙唐家窝铺。不同产地防风的升麻素苷和5−O−甲基维斯阿米醇苷含量差异较大，可能是决定药材质量的关键物质^[12]。

2.2 防风单株产量及单株色原酮总量

防风单株产量及单株色原酮总量的测定结果见表2。辽宁边杖子和辽宁东碾子沟的单株产量最高，分别为11.025、10.262 g，吉林白城和山东烟台最低，分别为5.777、6.302 g，最高单株产量是最低的1.9倍；河北大屯乡和山东烟台的单株色原酮总量最高，分别为81.329、60.115 mg，宁夏隆德、内蒙牛营子和内蒙大明镇最低，分别为12.394、21.224和21.476 mg，最高是最低的6.6倍。

2.3 防风土壤因子

16个产地的土壤pH在7.46~8.87之间，土壤偏碱性，部分产地的土壤pH较接近；不同产地土壤电导率差异明显，范围为74.30~158.60，电导率最高产地为吉林长春，最低为山东烟台，最高是最低的2.1倍；不同产地土壤水分、有机质和全氮含量差异明显，质量分数分别为3.71%~12.96%、0.23%~4.09%、0.01%~0.32%，含量最高产地均为甘肃兰州，最低分别为内蒙文钟镇、河北金沟屯和内蒙大明镇；土壤中全磷质量分数相对较高，在0.29~1.29 g/kg之间，含量最高产地为河北金沟屯，最低为内蒙牛营子，最高是最低的4.4倍。

土壤中其他9种有效元素分别为有效氮、有效磷、有效钾、有效钙、有效镁、有效铁、有效铜、有效锌和有效锰，不同产地含量差异较为明显，质量分数分别为37.38~282.35 mg/kg、2.77~47.46 mg/kg、45.84~718.85 mg/kg、325.68~1 303.58 mg/kg、83.20~128.46 mg/kg、0.97~31.19 mg/kg、0.55~1.44 mg/kg、0.38~2.49 mg/kg、0.59~5.14 mg/kg。有效氮、有效磷和有效钾含量最高的产地分别为甘肃兰州、辽宁黄台子和吉林长春，含量最低的产地分别为辽宁东碾子沟、山东烟台和河北大屯乡，含量最高产地分别是最低产地的7.6、17.1、15.7倍；有效钙和有效镁含量最高产地均为宁夏隆德，含量最低产地均为山东烟台，含量最高产地分别是最低产地的4.0、1.5倍；有效铁、有效铜、有效锌和有效锰含量最高产地分别为甘肃兰州、辽宁边杖子、山东烟台、河北大屯乡，含量最低产地分别为辽宁黄台子、河北金沟屯、内蒙牛营子、内蒙汐子镇，含量最高产地分别是最低产地的32.1、2.6、6.6、8.7倍(表3)。

2.4 防风色原酮含量与土壤因子的相关性分析

防风色原酮含量与土壤因子的相关性分析见表4。升麻素苷含量与土壤pH显著负相关(P<0.05)，与电导率极显著负相关(P<0.01)，与有效磷含量显著负相关(P<0.05)，与有效锰含量显著正相关(P<0.05)；升麻素含量与土壤pH、有效磷含量、有效钙含量显著负相关(P<0.05)，与有效锌含量显著正相关(P<0.05)；5−O−甲基维斯阿米醇苷含量与土壤有效磷含量显著负相关(P<0.05)；亥茅酚苷含量与土壤pH显著负相关(P<0.05)，与有效铁含量显著正相关(P<0.05)；色原酮总量与土壤电导率、有效磷含量、有效钙含量显著负相关(P<0.05)，与有效锰含量显著正相关(P<0.05)。结果说明：防风升麻素苷含量可能受土壤pH、电导率、有效磷含量和有效锰含量影响，升麻素含量可能受土壤pH、有效磷含量、有效钙含量和有效锌含量影响，5−O−甲基维斯阿米醇苷含量可能受土壤有效磷含量影响，亥茅酚苷含量可能受土壤pH和有效铁含量影响，色原酮总量可能受电导率、有效磷含量、有效钙含量和有效锰含量影响。防风色原酮含量与土壤因子的相关性分析充分表明根际土壤因子参与并影响药材品质的形成，色原酮含量与土壤因子具有一定的相关性^[13-15]。

2.5 防风单株产量及单株色原酮总量与土壤因子的相关性分析

防风单株产量及单株色原酮总量与土壤因子的相关性分析见表4。单株产量与土壤因子不具有显著相关性；单株色原酮总量与土壤有效钙含量显著负相关(P<0.05)，与有效锰含量显著正相关(P<0.05)。结果说明：防风单株产量可能不受土壤因子影响，单株色原酮总量可能受有效钙和有效锰含量影响。

2.6 防风色原酮总量与土壤因子的聚类分析

以防风色原酮总量和根际土壤因子为评价指标进行聚类分析，结果见图1。在欧氏距离15~20之间被分为I、II、III大类，宁夏隆德单独聚为II大类，吉林长春单独聚为III大类，其余产地聚为I大类；缩小欧氏距离至5~10之间又被分为A、B、C、D 4类，I大类分为A、B 2类，山东烟台单独聚为B类；继续缩小欧氏距离至0~5，A类又分为①和② 2小类，①包括9个产地，分别是辽宁边杖子、辽宁黄台子、内蒙汐子镇、辽宁东碾子沟、内蒙唐家窝铺、河北大屯乡、河北金沟屯、内蒙杨家营子和内蒙大明镇，②包括4个产地，分别是甘肃兰州、内蒙文钟镇、吉林白城和内蒙牛营子。聚类分析结果显示，以欧氏距离5~10为标准，16个产地分为4个产区，吉林长春产区、宁夏隆德产区、山东烟台产区和其他产区(包括内蒙、河北、辽宁等)。原因可能是药材质量差异和地理位置导致的土壤元素含量差异，山东烟台和吉林长春药材质量较高且距离内蒙药材产区较远单独聚为两类；宁夏隆德药材质量最低，同样距离内蒙药材产区较远聚为一类；其他产地药材质量相似、地理位置相近聚为一类。聚类分析结果充分说明土壤因子影响药材质量^[16-18]。

2.7 防风色原酮总量与根际土壤因子多重线性回归分析

以防风色原酮总量为因变量、土壤因子为自变量进行多重线性回归分析，表5为拟合模型的拟合优度情况简报，重要指标R²(决定系数)反映了回归贡献的相对程度^[19]，模型1含有变量有效磷含量，R²=0.365,说明有效磷含量可以解释36.5%的防风色原酮总量信息；模型2含有变量有效磷和有效锰含量，R²=0.568，说明这2个变量可以解释56.8%的防风色原酮总量信息；模型3含有变量有效磷、有效锰和全磷含量，R²=0.718，说明有效磷、有效锰和全磷含量3个变量可以解释71.8%的防风色原酮总量信息。模型3为最终的拟合结果，有效磷、有效锰和全磷含量3个变量可以有效解释防风色原酮总量信息，其他变量对防风色原酮信息的解释较小，无统计学意义，未引入模型中。模型偏回归系数(表5)是对线性回归分析是否具有统计学意义的检测，有效磷含量t=−0.579，P<0.01；有效锰含量t=−0.671，P<0.01；全磷含量t=−0.445，P<0.05；说明变量的各项偏回归系数均具有统计学意义,回归方程为Y=7.252−0.132X₈−1.206X₁₅−3.631X₇。图2为因变量累计概率和模型预测值累计概率间的正态P-P图，其散点基本围绕参考直线均匀分布，说明试验数据符合线性回归模型的适用条件，具有一定独立性、正态性和方差齐性。

3 讨论与结论

生态因子是影响药材药效成分含量的重要因素，包括气候因子、土壤因子及地理因素等。适宜的生态环境能促进药材的生长发育,有利于有效成分的合成和积累,从而提高药材质量。周国富等^[20]研究生态因子对黄芩5种主要指标成分的影响时发现黄芩有效成分含量与气候因子密切相关，尤其是温度和降水量。卫昊等^[21]在黄芩7种黄酮类有效成分含量的研究中发现，随海拔升高，黄酮类成分含量呈显著上升趋势。这些结果充分说明气候因子对药用植物质量有重要影响。

土壤因子是植物生长发育的基础，土壤的物理、化学性质以及所含的各种元素对药材生长发育及次生代谢物积累均具有重要影响^[22]。土壤中的磷元素对植物吸收营养具有积极作用，能够促进植物碳水化合物的合成和运输，提高土壤中磷含量，增强植物抗旱、抗寒能力。土壤中的锰元素是植物生长发育不可缺少的微量元素之一，锰可以促进植物光合作用，调节植物体内氧化还原反应，含量充足时可以增强植物病害抗性。防风本身就具有抗旱、抗寒能力，土壤中磷元素含量过高可能会导致磷元素不能被防风根系有效吸收，从而影响防风药效成分的形成；土壤中锰元素含量少会导致防风光合作用受到影响，增加患病几率，从而导致药材药效成分含量降低。孙晶波等^[23]研究了根际土壤化学性质对防风色原酮和香豆素含量的影响，发现速效磷是影响防风有效成分的关键因素，锰含量增加虽然可以促进防风色原酮和香豆素类成分的形成，但却会抑制防风对锰元素的摄取与累积，促进防风对钙元素的吸收。本研究结果与其结论部分相似，防风色原酮含量与部分土壤因子具有一定相关性，包括pH、电导率、有效磷含量、有效锰含量等。以防风色原酮总量为因变量、土壤因子为自变量进行多重线性回归分析，发现有效磷、有效锰、和全磷含量3个变量可以解释71.8%的防风色原酮总量信息，进一步验证了有效磷和有效锰在防风药材质量形成过程中起到的重要作用。辽宁黄台子、内蒙文钟镇、吉林长春和内蒙大明镇等产地有效磷含量较高，内蒙文钟镇、内蒙汐子镇和辽宁东碾子沟等产地有效锰含量较低，建议这些产地在施肥时适当减少磷肥的使用，增加锰肥的比例。

不同地区不同土壤类型可能导致防风有效成分含量及产量的显著差异，本文通过不同产地防风单株产量及单株色原酮总量与土壤因子的相关性分析，发现单株产量与土壤因子不具相关性，单株色原酮总量与有效钙含量显著负相关，与有效锰含量显著正相关，说明土壤中有效钙含量高可能会降低防风单株色原酮总量，有效锰含量高可能会提高防风单株色原酮总量。本研究测定了不同产地防风药材产量及色原酮含量，总结出5个种植防风具有经济效益的地区，分别是河北大屯乡、河北金沟屯、山东烟台、吉林白城和内蒙文钟镇。这5个产地的防风色原酮含量、单株产量、单株色原酮总量都普遍偏高。另外，本研究通过对不同产地防风色原酮总量以及土壤因子进行聚类分析，将16个产地分为四大产区：吉林长春产区、山东烟台产区、宁夏隆德产区、其他产区(包括内蒙、河北、辽宁等产地)，合理有效地使防风药材质量相似或相近的产地聚为一类。

因此，不同产地防风有效成分的差异与根际土壤环境密切相关，本研究为选择防风栽培环境以及适当施肥提供了必要的理论依据。