中药材天麻又名赤箭、定风草，为兰科(Orchidaceae)植物天麻(Gastrodia elata)的干燥块茎，具有息风止痉、平抑肝阳、祛风通络之功效，常用于治疗小儿惊风，癫痫抽搐，破伤风，头痛眩晕，手足不遂，肢体麻木和风湿痹痛等病症^[1]，在中国已有2 000多年的药用历史。现代药理研究表明，天麻还具有镇静催眠、镇痛、抗惊厥、抗氧化、改善记忆力、增强免疫和降血压等作用^[2]。天麻的主要活性成分有天麻素、对羟基苯甲醇和巴利森苷类等^[3-5]，而浸出物也是评价其药材质量的重要指标。1987年天麻被列入《中国珍稀濒危保护植物名录》^[6]，2021年又被收录于《国家重点保护野生植物名录》^[7]，为国家二级重点保护野生植物。长期以来，因受生态环境和人为因素的影响，野生天麻的生长发育和繁殖机能逐渐减退，人工栽培呈必然趋势。已有研究表明，气温在20-25 ℃时天麻生长较快，而在30 ℃以上时，天麻的生长便会受到抑制^[8]；乌天麻(G.elata f.glauca)适宜生长在海拔较高的云南和贵州，而红天麻(G.elata f.elata)则主要分布在海拔相对较低的陕西、河南等地^[9]。此外，土质疏松，利水透气的微酸性或中性土壤较适宜天麻生长^[10]。目前有关生态环境因子对天麻内在品质相关成分含量影响的研究还很有限，探明生态环境因子与天麻品质形成之间的关系，对天麻资源的可持续利用及道地性形成机制的解析具有重要意义。本研究对在陕西省汉中市宁强县采集到的48份天麻样品进行浸出物及主要有效成分含量测定并建立指纹图谱，利用ArcGIS 10.2和SPSS 20.0软件分析气候、土壤、地形地貌等77个生态环境因子与成分含量之间的相关性，总结天麻品质相关成分在小尺度(县域)空间范围内的变异规律，为天麻资源的可持续开发利用提供参考。

1 材料与方法 1.1 仪器

LC-20AT型高效液相色谱仪(日本岛津公司)，Diamonsil C₁₈型色谱柱(250 mm×4.6 mm，5 μm，北京迪马科技有限公司)，Quintix35-1CN型十万分之一电子天平[赛多利斯科学仪器(北京)有限公司]，BSA124S型万分之一电子分析天平[赛多利斯科学仪器(北京)有限公司]，HH-4A型数显单控单列恒温电热水浴锅[国华(常州)仪器制造有限公司]，KQ-500DE型超声波清洗器(昆山舒美超声仪器有限公司)，SHB-IIIS型循环水真空泵(北京世纪予华仪器有限公司)，350T型德蔚高速多功能粉碎机(广州市渲龙厨具有限公司)，DHG-9123A型鼓风干燥箱(上海一恒科学仪器有限公司)。

1.2 试剂

天麻素(批号：T10M9F55562)、对羟基苯甲醇(批号：H21D6Q7813)、巴利森苷A(批号：P24S11F124204)、巴利森苷B(批号：A13GB157725)、巴利森苷C(批号：D06N11S130122)、巴利森苷E(批号：D08D11S133702)均购自上海源叶生物科技有限公司，纯度≥98%。其他试剂还有乙腈[色谱纯，赛默飞世尔科技(中国)有限公司]、磷酸[色谱纯，福晨(天津)化学试剂有限公司]、水(娃哈哈纯净水)、无水乙醇(分析纯，天津市致远化学试剂有限公司)、甲醇(分析纯，天津市致远化学试剂有限公司)。

1.3 样品

于2021年10-12月在陕西省汉中市宁强县18个镇采集林下栽培的新鲜天麻，记录经纬度和海拔，随机从中抽取分析样品。具体样品信息见表 1。样品经北京中医药大学孙志蓉教授鉴定为天麻(Gastrodia elata Bl.)的新鲜块茎，将块茎放入蒸锅中，保持水沸腾蒸制30 min，待天麻无白心时取出，室内自然冷却，随后将天麻样品放入鼓风干燥箱中于60 ℃烘干至恒重，粉碎，过4号筛后干燥保存备用。

1.4 生态环境因子

月度气象数据(包括平均气温、降水量、水气压、风速、太阳辐射，时间尺度为1970-2000年)以及海拔数据均源自全球气候数据库(http://www.worldclim.org/)。土壤数据(土壤沙含量、黏土含量)来源于联合国粮食和农业组织土壤门户网站的世界土壤数据库(HWSD)土壤数据集(v1.2，https://www.fao.org/soils-portal/en/)，中国区域的土壤数据由第二次全国土地调查提供的《1∶100万中华人民共和国土壤图》(1995年编制)得到。坡度数据使用ArcGIS 10.2软件的投影栅格、按掩膜提取、坡度分析工具处理海拔数据得到。所有数据空间分辨率为30 s(约1 km×1 km)。

1.5 浸出物含量测定

参照《中华人民共和国药典(一部)》(2020年版)^[1](以下简称《中国药典》)的方法测定。取供试品2-4 g，精密称定，置于100-250 mL锥形瓶中，精密加水50-100 mL，密塞，称定质量，静置1 h后，连接回流冷凝管，加热至沸腾，并保持微沸1 h。自然冷却后，取下锥形瓶，密塞，再称定质量，用水补足减失的质量，摇匀，用干燥滤器滤过，精密量取滤液25 mL，置于已干燥至恒重的蒸发皿中，水浴蒸干后，于105 ℃干燥3 h，置干燥器中冷却30 min，迅速精密称定质量。

1.6 有效成分含量测定 1.6.1 色谱条件

采用Diamonsil-C18 (250 mm×4.6 mm，5 μm)色谱柱；流动相为乙腈(A)-0.1%磷酸水(B)，梯度洗脱(0-10 min，3%-10% A；10-15 min，10%-12% A；15-25 min，12%-18% A；25-40 min，18% A；40-42 min，18%-95% A)；流速为0.8 mL/min；检测波长为220 nm；柱温为30 ℃；进样量为10 μL。

1.6.2 对照品溶液的制备

取对照品适量，精密称定，置于10 mL容量瓶中，加乙腈-水(3∶97)混合溶液溶解定容，配制成质量浓度分别为天麻素(100 μg/mL)、对羟基苯甲醇(30 μg/mL)、巴利森苷A (362.6 μg/mL)、巴利森苷B (101.8 μg/mL)、巴利森苷C (102 μg/mL)、巴利森苷E (102.2 μg/mL)的混合对照品溶液。

1.6.3 供试品溶液的制备

精密称定天麻样品粉末0.5 g，置于具塞锥形瓶中，加入50%甲醇25 mL，称定质量，超声处理(500 W，40 kHz) 30 min，自然冷却后，用水补足减失的质量，4 000 r/min离心5 min，取上清液过0.22 μm微孔滤膜，即得供试品溶液。将48份天麻供试品的指纹图谱数据导入“中国色谱指纹图谱相似度评价系统(2004版)”，以S10(大安镇新明村)样品为参照峰，采用中位数法和自动匹配法，获取共有峰并建立宁强县天麻共有的对照指纹图谱，计算供试品相对于对照指纹图谱相似度值。

1.6.4 标准曲线的制备

精密量取1.6.2节制备的混合对照品溶液，用乙腈-水(3∶97)混合溶液逐级稀释至0.00、1.25、2.50、5.00、10.00、20.00倍，制成一系列不同浓度的混合对照品溶液，各取10 μL注入高效液相色谱仪。以1.6.1节的色谱条件测定各对照品峰面积，以对照品浓度为横坐标(X)，以峰面积为纵坐标(Y)进行线性回归(表 2)。

1.6.5 精密度试验

精密量取1.6.2节制备的混合对照品溶液，以1.6.1节的色谱条件测定各对照品峰面积，重复6次。其中天麻素、对羟基苯甲醇、巴利森苷E、巴利森苷B、巴利森苷C与巴利森苷A的峰面积相对标准偏差(Relative Standard Deviation，RSD)分别为0.26%、0.20%、0.06%、0.05%、0.06%、0.11%。

1.6.6 稳定性试验

精密量取1.6.3节制备的供试品溶液，分别于配制后0、6、12、24、36 h各进样1次，以1.6.1节的色谱条件测定峰面积。其中天麻素、对羟基苯甲醇、巴利森苷E、巴利森苷B、巴利森苷C与巴利森苷A的峰面积RSD分别为1.32%、2.18%、1.83%、1.98%、2.21%与2.25%，表明供试品溶液36 h内稳定性良好。

1.6.7 重复性试验

精密称取同一天麻样品粉末2.0 g，共6份。按1.6.3节下制备成供试品溶液，以1.6.1节的色谱条件测定各供试品峰面积并计算各成分质量分数。其中天麻素、对羟基苯甲醇、巴利森苷E、巴利森苷B、巴利森苷C与巴利森苷A的质量分数RSD分别为2.60%、1.17%、1.39%、1.42%、1.54%、1.43%。

1.6.8 加样回收率试验

精密称定已知成分含量的天麻供试品0.5 g，共6份，按照要求制备供试品溶液，分别加入各对照品溶液，进样测定，记录峰面积，计算平均加样回收率及RSD，结果见表 3。

1.7 数据分析

使用ArcGIS 10.2软件的采样工具对48个采样点的环境数据进行提取，其中海拔匹配情况见图 1。使用SPSS 20.0软件对数据进行方差分析及皮尔逊相关性分析，运用Microsoft Office Excel 2020、Adobe Illustrator CS6、OriginPro 2015完成数据和图表分析制作。

2 结果与分析 2.1 指纹图谱建立及相似度分析

通过比较对照品色谱图(图 2)，判定共有成分分别为A：天麻素，B：对羟基苯甲醇，C：巴利森苷E，D：巴利森苷B，E：巴利森苷C，F：巴利森苷A，与《中国药典》规定的特征图谱保留峰一致，说明建立的检测方法稳定、可靠。48批天麻样品的共有峰为6个，共有峰保留时间RSD较小(0.09%-0.17%)，共有峰的峰面积RSD较大(18.50%-52.59%)，说明48份天麻样品中同一成分的含量差异较大。由表 4可看出，不同采样点之间天麻的特征图谱相似度较高，除S37(舒家坝镇文家河村)样品相似度低于0.9(为0.843)，其他样品相似度均为0.948-0.998，表明宁强县内天麻有效成分基本一致。

2.2 不同采样点天麻浸出物及有效成分含量聚类分析

对48份天麻样品的浸出物及6种有效成分的含量作聚类分析，结果表明：汉源街道办黄坝驿村天麻的天麻素、巴利森苷B、巴利森苷C含量最高，分别为1.05%、0.79%、0.30%；阳平关镇张家河村天麻的对羟基苯甲醇含量最高(0.18%)；高寨子街道办古城村天麻的巴利森苷A含量最高(1.95%)；高寨子街道办何家院子村天麻的巴利森苷E含量最高(0.93%)，巴利森苷类成分总含量以高寨子街道办古城村最高(3.57%)，约为大安镇龙泉村(1.52%)的2.35倍，其次是汉源街道办黄坝驿村(3.54%)。数据标准化后的聚类分析结果表明，总体来看不同采样点聚为3类。6种成分聚为4类，其中天麻素与巴利森苷C、巴利森苷A、巴利森苷B聚为一类，浸出物、对羟基苯甲醇、巴利森苷E各为一类(图 3)。

2.3 天麻浸出物及有效成分含量

48份天麻样品的浸出物及有效成分含量统计分析结果见表 5，其中天麻素的含量为0.252 6%-1.050 7%，对羟基苯甲醇为0.027 4%-0.182 7%，巴利森苷E为0.378 9%-0.928 3%，巴利森苷B为0.312 2%-0.787 8%，巴利森苷C为0.087 3%-0.301 2%，巴利森苷A为0.379 2%-1.949 0%，浸出物为16.05%-30.59%。含量变化幅度对羟基苯甲醇＞天麻素＞巴利森苷A＞巴利森苷C＞巴利森苷B＞巴利森苷E＞浸出物，表明浸出物及巴利森苷E含量稳定性高，离散程度小，对羟基苯甲醇及天麻素稳定性低，离散程度大。县内所有样点的浸出物均符合《中国药典》要求。不同采样点天麻浸出物及6种有效成分含量的差异显著(P＜0.001, 表 6)，其中6种有效成分总量以汉源街道办黄坝驿村天麻最高，达4.65%，约为大安镇龙泉村天麻(1.81%)的2.57倍；汉源街道办黄坝驿村红天麻浸出物含量最高，达30.59%，禅家岩镇张家坝村乌天麻浸出物含量最低，为16.05%。

2.4 天麻浸出物及有效成分含量与生态环境因子的相关性分析

将天麻浸出物及有效成分含量与采样点相应的77个生态环境因子进行Spearman相关性分析和显著性检验，结果表明：浸出物含量与年平均气温、年平均水气压及土壤沙含量均呈显著或极显著正相关，与年平均降水量、年平均风速、海拔、土壤淤泥含量均呈显著负相关。6种有效成分总量与年平均气温、年平均降水量、年平均风速均呈显著正相关，与海拔、坡度呈显著负相关。巴利森苷类成分含量与环境因素的关系和6种有效成分总量表现一致，而天麻素+对羟基苯甲醇总量仅与坡度呈显著负相关关系(图 4)。随着土壤沙含量和淤泥含量的升高，6种有效成分总量整体呈先升后降的趋势(图 5)。由图 6可以看出，对羟基苯甲醇与2月降水，8、9月太阳辐射，交换性盐基，土壤黏土含量及土壤的阳离子交换能力呈显著正相关；巴利森苷A与1、12月降水，3、7、8月风速，土壤淤泥含量呈显著负相关，与黏性层土壤的阳离子交换能力呈显著正相关；巴利森苷E与坡度呈显著负相关；巴利森苷B与黏性层土壤的阳离子交换能力呈显著正相关。根据相关性大小发现，生态环境因子与天麻浸出物及有效成分含量的关系为气候＞海拔＞土壤，气候因子中影响较大的为12月降水量、4月风速和12月平均气温。与坡度相比，海拔影响较大。土壤因子中影响最大的为土壤淤泥含量，其次为土壤沙含量。

3 讨论

天麻的人工培育会受到多种因素的影响，如促进种子萌发的萌发菌、为其提供直接营养的蜜环菌、栽培的环境和技术等。本研究选择栽培条件相对一致的样点(萌发菌、蜜环菌和栽培技术相同)采集天麻样品，并结合采样点生态环境因子数据分析发现，宁强县天麻浸出物及有效成分含量与生态环境因子的关系为气候＞海拔＞土壤。其中，气候因子中影响较大的为年平均气温、年平均降水量和年平均风速；海拔的影响大于坡度；土壤因子中影响最大的为土壤淤泥含量，其次为土壤沙含量。

气候是地球上某一地区多年时段大气的一般状态，是多个生态环境因子的综合表现^[11]。其中，光照、温度、水分对中药材的分布与品质的形成最为重要，如防己(Stephania tetrandra)中粉防己碱和防己诺林碱含量与年平均气温呈正相关^[12]；欧洲莨菪(Scopolia carniolica)在干旱时阿托品含量达1%，而在湿润环境下只有0.3%-0.5%^[13]。本研究发现，宁强县属山地暖温带湿润季风气候，县域内所产天麻的6种有效成分总量与年平均气温呈显著正相关，天麻浸出物含量与1-5月、10-12月降水量，1-12月风速呈显著或极显著负相关，原因可能是一定气候区域内气温影响了蜜环菌和天麻的生长及二者的伴生关系，有利于生物量积累和品质形成，而春季和秋冬季的降水直接影响林地土壤水分和空气湿度，风速会增加林地通透性并调节气温，同时也会在少雨年份或季节加大蒸发量而加剧干旱程度，影响萌发菌、蜜环菌和天麻的生长。中药资源具有明显的空间和地域分布规律^[13]，药材品质的形成与产地的地理位置、地形地貌、海拔等因子密切相关，如纹党(Codonopsis pilosula)中浸出物含量与海拔呈正相关^[14]。本研究中天麻浸出物含量与海拔却呈显著负相关，可能是由于高海拔主要种植乌天麻，其麻形偏小，水分含量低，浸出物含量相对较高；但是低海拔红天麻个大，浸出物绝对含量较高。已有研究证实中药材在不同采收时间里的有效成分含量存在差异^[15]。本研究选择天麻成熟期进行采收，不同采样点的天麻成熟期具有差异，导致天麻成熟期时间不一的原因可能包括环境因素和人为因素(栽培技术)，本研究在保证栽培条件(所用菌种和栽培技术)一致的前提下进行采收，以此来探究环境因素对天麻质量的影响。土壤是药用植物生长的基础，土壤的质地、养分、酸碱性等均会对中药材的品质产生影响^[16]，如最适合金银花(Lonicera japonica)的土壤类型是中性或稍碱性的沙壤土^[17]。本研究发现，天麻6种有效成分总量与土壤淤泥含量呈显著负相关，与土壤沙含量呈显著正相关。因天麻适宜生长在较疏松的沙土或沙壤土中^[18]，土壤通透性不佳不利于蜜环菌生长，且易染杂菌，从而影响天麻生长。此外，本研究发现天麻在生长过程中受到环境胁迫(如气温、水气压较高等)时，在外观形态上会变细长、生物量降低、浸出物含量增加，当处于顺境(如气温较低、风速较快等)时，天麻会变短粗、生物量增加、浸出物及巴利森苷A含量降低，同样体现了中药材的“逆境效应”^[19]。宁强县天麻的6种有效成分总量和海拔、坡度呈负相关，高海拔下气温低，可能不利于有效成分的积累，这与6种有效成分含量与年平均气温呈正相关的结果吻合。

位于宁强县巴山北坡高中山区的毛坝河镇、禅家岩镇、二郎坝镇及巴山镇的气候特点为气温低、风速大、促生产月降水量大、光照稍不足，浸出物含量相对较低，6种有效成分含量相对居中，在天麻栽培过程中应注意调节林分郁闭度，采取覆盖保温和防积水等措施；宁强县北部秦岭南坡中山的代家坝镇、巨亭镇、大安镇、太阳岭镇以及县南部的二郎坝镇、巴山镇的气候特点为7-8月降水量相对较少，光照充足，气温较低，2-4月、9-11月降水量稍不足，浸出物和6种有效成分含量相对较低，种植天麻时应注意在天麻生长前期和中后期适当补水，适当调节林分郁闭度；宁强县中部代家坝镇、胡家坝镇、汉源街道办、高寨子街道办、铁锁关镇、阳平关镇的气候特点为气温较高、光照较为充足、风速较低、降水量较为充沛，浸出物和6种有效成分含量相对较高，栽培时应当注意林分郁闭度、温度和湿度的调控，保持通风以保证产量；宁强县西部低山区青木川镇、广坪镇、安乐河镇、燕子砭镇的气候特点为气温高，促生产月降水量少，风速低，气候较不利于天麻增产，浸出物和6种有效成分含量相对居中，栽培时应当注意调控温度和湿度。

4 结论

陕西省汉中市宁强县天麻浸出物含量与年平均气温、年平均水气压、土壤沙含量均呈显著或极显著正相关，与年平均降水量、年平均风速、海拔、土壤淤泥含量均呈显著负相关；6种有效成分总量与年平均降水量、年平均气温、年平均风速均呈显著正相关，与海拔、坡度呈显著负相关。生态环境因子与天麻浸出物及有效成分含量的关系为气候＞海拔＞土壤，气候因子中影响较大的为12月降水量、4月风速和12月平均气温。与坡度相比，海拔影响较大。土壤因子中影响最大的为土壤淤泥含量，其次为土壤沙含量。

本研究从小尺度层面上探究生态环境因子对天麻品质相关成分含量的影响，研究结果可为宁强县优质天麻培育及天麻资源合理开发利用提供参考，同时为其他行政区域面积大、地形地貌复杂县域的中药材产业规划提供参考。在具体研究过程中应考虑到所用生态环境因子数据库的实时性和分析方法的准确性。