2. 宁强天麻研究所有限责任公司, 陕西汉中 724400;
3. 宁强县发展和改革局科技开发中心, 陕西汉中 724400;
4. 宁强县中药材产业发展中心, 陕西汉中 724400
2. Ningqiang Tianma Research Institution Limited Liability Company, Hanzhong, Shaanxi, 724400, China;
3. Science and Technology Development Center, Ningqiang County Development and Reform Bureau, Hanzhong, Shaanxi, 724400, China;
4. Ningqiang County Traditional Chinese Medicine Industry Development Center, Hanzhong, Shaanxi, 724400, China
天麻为兰科(Orchidaceae)植物天麻(Gastrodia elata Bl.)的干燥块茎,具有息风止痉、平抑肝阳、祛风通络之功效,对头晕头痛、镇静催眠、癫痫抽搐等有显著疗效。天麻的主要产地分布在云南、湖北、贵州、陕西、安徽、四川等地[1, 2],其中陕西省汉中市为西天麻的道地产区,所产天麻大多为红天麻(G.elata Bl.f.elata)[3]。汉中天麻主要产自宁强县、勉县和略阳县,不同产地天麻的产量和质量有所不同,故本研究采集了以上3个地区的红天麻作为研究材料进行测定与分析。
药材的外观性状,如形状、颜色、气味等,往往与其内在质量密切相关,但主观评判的标准难以统一。如今,借助电子感官仪器可以将药材外观特征量化,使评判更加客观与直观,为药材质量与视觉、嗅觉等信息的联系提供了可行性[4, 5]。在市场中,人们普遍依据天麻外观性状对其质量进行评判,但是人为评判的主观性过强,缺乏准确性[6-8]。电子鼻与分光测色仪可以准确、客观地量化药材气味与颜色,常用于中药材的鉴别与质量评价[9-12]。《中华人民共和国药典(一部)》(2020年版)(以下简称2020年版《中国药典》)对天麻药材中水分、浸出物、天麻素(Gastrodin)和对羟基苯甲醇(P-hydroxybenzylalcohol)的含量做了明确要求[2],另外巴利森苷(Parishins)在天麻中具有显著药理活性且含量丰富[13],所以本研究选择以上指标进行化学成分含量测定。
主成分分析(PCA)法是一种降维数据的多元统计分析方法,在中药研究中应用广泛[14-16]。赵梦利等[17]利用PCA法对鸡血藤(Spatholobi Caulis)饮片数据降维以观察差异。巫晓霞等[18]利用PCA法对不同产地茯苓进行质量评价。熵权TOPSIS(Technique for Order Preperence by Similarity to an Ideal Solution)法可充分利用原始数据,且可以排除主观因素的干扰,更客观地筛选出综合功能特性优良的样本,适用于中药质量综合评价[19-21]。高伟城等[22]利用熵权TOPSIS法对不同产地枇杷叶质量进行综合评价,李雨昕等[23]、史丛晶等[24]利用熵权法的灰色关联法-TOPSIS法分别对不同产地生品、蒸品三七,知母及盐知母进行质量评价。本研究将PCA法和熵权TOPSIS法相结合,先利用PCA法提取主成分,再通过熵权TOPSIS法计算权重并排名,能有效减少冗余的评价指标,提高多指标评价的直观性与科学性[25, 26]。
综上,本研究对在陕西省汉中市宁强县、勉县和略阳县采集到的天麻样品进行外观性状、气味、颜色、含水量、浸出物和主要有效成分含量测定,在传统形状指标与生物量的基础上,结合感官量化指标与有效成分含量,利用PCA法结合熵权TOPSIS法对天麻质量进行综合评价,探讨各指标在天麻质量评价中的贡献,为汉中天麻质量标准建立及天麻资源的开发利用提供参考。
1 材料与方法 1.1 材料于2020年10-11月在陕西省汉中市宁强县、勉县、略阳县采集27份红天麻,经北京中医药大学孙志蓉教授鉴定为兰科植物天麻Gastrodia elata Bl.的块茎。清洗后统一加工、粉碎,过65目筛,备用。样品信息见表 1。
市县 Cities and Counties |
乡镇 Towns |
编号 Number |
村 Village |
海拔/m Altitude/m |
栽培模式 Cultivation pattern |
采收时间 Collect time |
Ningqiang County | Yanzibian Town | 1 | Caojiagou Village | 962.1 | Cultivation under forest | 2020-10-04 |
3 | Panjiaba Village | 1 007.2 | Cultivation under forest | 2020-10-04 | ||
4 | Xinchangjie Village | 1 300.0 | Cultivation under forest | 2020-11-10 | ||
5 | Shenjiaping Village | 1 068.7 | Cultivation under forest | 2020-10-08 | ||
6 | Jindongzi Village | 630.7 | Cultivation under forest | 2020-10-08 | ||
Da′an Town | 7 | Mandongzi Village | 734.7 | Cultivation under forest | 2020-10-07 | |
8 | Tielugou Village | 677.3 | Greenhouse cultivation | 2020-10-07 | ||
9 | Longquan Village | 895.0 | Cultivation under forest | 2020-10-07 | ||
Yangpingguan Town | 13 | Wujiaba Village | 850.0 | Cultivation under forest | 2020-10-08 | |
15 | Dabahe Village | 674.0 | Cultivation under forest | 2020-11-10 | ||
Shujiaba Town | 16 | Songshugou Village | 1 050.0 | Cultivation under forest | 2020-12-11 | |
Taiyangling Town | 18 | Zhaojia Village | 1 063.7 | Cultivation under forest | 2020-10-09 | |
19 | Qinglinju Village | 1 083.3 | Cultivation under forest | 2020-10-09 | ||
20 | Tangjiaying Village | 1 058.1 | Cultivation under forest | 2020-10-09 | ||
Juting Town | 21 | Longgangba Village | 775.5 | Flat ground cultivation | 2020-10-07 | |
22 | Juting Village | 960.0 | Cultivation under forest | 2020-10-07 | ||
25 | 867.2 | Cultivation under forest | 2020-10-05 | |||
Tiesuoguan Town | 30 | Xiaogou Village | 940.0 | Cultivation under forest | 2020-11-30 | |
Mian County | Xinpu Town | mx-31 | Yujiadong Village | 910.0 | Cultivation under forest | 2020-10-08 |
mx-32 | Yujiahe Village | 1 100.0 | Cultivation under forest | 2020-11-30 | ||
Changgouhe Town | mx-33 | Changgouhe Village | 860.0 | Cultivation under forest | 2020-11-30 | |
Tonggousi Town | mx-34 | Liuba Village | 780.0 | Cultivation under forest | 2020-11-30 | |
Zhangjiahe Town | mx-35 | Bamiao Village | 1 300.0 | Cultivation under forest | 2020-11-30 | |
Lueyang County | Guo Town | ly-36 | 890.0 | Cultivation under forest | 2020-10-08 | |
Xihuaiba Town | ly-37 | Dagou Village | 1 050.0 | Cultivation under forest | 2020-11-11 | |
ly-38 | Liangjiahe Village | 1 089.0 | Cultivation under forest | 2020-11-11 | ||
Hengxianhe Town | ly-39 | Shizhuanggou Village | 900.0 | Cultivation under forest | 2020-11-11 |
1.2 仪器与试剂
小型粉碎机(CK-20,北京金洋利科技发展有限公司), 电磁炉(21FS37-20110632,九阳股份有限公司), 数显游标卡尺(DL91150,得力集团有限公司), 鼓风干燥箱(DHG-9123A,上海一恒科技有限公司), 分光测色仪(CM-5,日本柯尼达美能达公司), 超快速气相电子鼻(HeraclesⅡ,法国Alpha Mos公司), 自动进样器(HS-100,法国Alpha Mos公司), 高效液相色谱仪(2498,美国Waters科技有限公司), C18色谱柱(X Bridge,美国Waters科技有限公司), 十万分之一电子天平(XS205, 梅特勒-托利多仪器上海有限公司), 电热恒温水浴锅(HH-4A,国华仪器制造有限公司), 数控超声波清洗器(KQ-500DE,昆山市超声仪器有限公司)。
天麻素对照品(B21243,纯度≥98%),对羟基苯甲醇对照品(B20326,纯度≥98%),巴利森苷A、B、C、E对照品(B20923、B20912、B2091、B20914,纯度≥98%),以上6种对照品均为上海源叶生物科技有限公司生产;乙腈、无水乙醇、磷酸(色谱级,Thermo Fisher Scientific公司);纯水(娃哈哈饮用纯净水,杭州娃哈哈集团有限公司)。
1.3 方法 1.3.1 外观形状与生物量测定采集的鲜天麻除去腐坏、虫蛀、不完整的块茎,洗去外部泥土,待表面水分晾干后,观察记录天麻的螺环纹数,测量天麻块茎的长度、宽度,称量鲜重,再将其蒸至内无白心,冷却至室温后放入烘箱,55 ℃烘干后称量干重。
1.3.2 颜色色度值测定用分光测色仪进行颜色测定。测定波长为350-750 nm;测定口径为8 mm;光源为脉冲氙灯;标准偏差ΔE*ab≤0.4;数据处理软件为Spectra Magic NX。称取天麻粉末适量,置于分光测定仪的培养皿中,平行测定3次,查看并记录色度空间的明亮度(L*)、红绿色(a*)、黄蓝色(b*),由公式(1)求得总色度值E*:
$ E^*=\sqrt{\left(a^*\right)^2+\left(b^*\right)^2+\left(L^*\right)^2}。$ | (1) |
天麻样品用高速粉碎机粉碎后,过24目筛,取上部较粗颗粒,按表 2的条件进行检测。
条件 Condition |
参数 Parameter |
Sample bottle | 20 mL |
Sample weight | 1.5 g |
Time of shaking | 20 min |
Temperature of shaking | 40 ℃ |
Injection volume | 1 000 μL |
Sample injection rate | 125 μL/s |
Sample injection time | 13 s |
FID gain | 12 |
Initial temperature (Collecting trap) | 40 ℃ |
Shunt | 10 mL/min |
Final temperature | 240 ℃ |
Initial temperature (Column temperature) | 50 ℃ |
Temperature programming | 3-250 ℃(21 s) |
Acquisition time | 110 s |
Detector temperature | 260 ℃ |
1.3.4 水分和浸出物含量测定
水分含量参照2020年版《中国药典》通则中的水分测定法(通则8302第二法烘干法)测定。浸出物含量参照2020年版《中国药典》通则中醇溶性浸出物测定下的热浸法(通则2201)测定。
1.3.5 天麻素和对羟基苯甲醇含量测定① 色谱条件。色谱柱:Waters X Bridge C18(250 mm×4.6 mm,5 μm);以乙腈为流动相A,0.05%磷酸溶液为流动相B,A: B(3:97)等度洗脱;检测波长:220 nm;理论塔板数不低于5 000。
② 对照品溶液的制备。取天麻素和对羟基苯甲醇对照品适量,精密称定,加乙腈-水(3∶97)混合,制成天麻素53 μg/mL、对羟基苯甲醇28 μg/mL的混合溶液,超声使其充分溶解,即得。
③ 供试品溶液的制备。精密称定过65目筛的天麻粉末2 g,置具塞锥形瓶中,精密加入50%稀乙醇50 mL,称重,超声处理(功率120 W,频率40 kHz)30 min,放冷,称重,用50%稀乙醇补足减失的重量,滤过,精密量取续滤液10 mL,放置水浴锅中浓缩近无醇味,加入乙腈-水(3∶97)混合溶液溶解,转移至25 mL容量瓶中,用乙腈-水(3∶97)混合溶液稀释至刻度,摇匀滤过,取续滤液,备用。
④ 标准曲线的建立。分别取5、10、15、20、25 μL的标准品溶液,按照天麻素和对羟基苯甲醇含量测定的色谱条件测定分析,以峰面积为纵坐标(y),浓度为横坐标(x),得到线性回归方程,标准曲线见表 3。
化合物 Components |
回归方程 Regression equation |
R2 | 线性范围/(μg/mL) Linear range/(μg/mL) |
Gastrodin | y=620248x-472.2 | 1 | 0.026 5-0.132 5 |
P-hydroxybenzylalcohol | y=577347x-17186 | 1 | 0.014 0-0.070 0 |
⑤ 样品含量测定。精密吸取对照品和供试品各10 μL,按照1.3.5节①的色谱条件测定,得到天麻素和对羟基苯甲醇的峰面积并计算二者含量。
1.3.6 巴利森苷类成分含量测定① 色谱条件。色谱柱:Waters X Bridge C18 (250 mm×4.6 mm,5 μm);以乙腈为流动相A,0.1%磷酸溶液为流动相B,梯度洗脱,具体见表 4;流速:0.8 mL/min,柱温:30 ℃;检测波长:220 nm;理论板数不低于5 000。
时间/min Time/min |
流动相A/% Mobile phase A/% |
流动相B/% Mobile phase B/% |
0-10 | 3-10 | 97-90 |
10-15 | 10-12 | 90-88 |
15-25 | 12-18 | 88-82 |
25-40 | 18 | 82 |
40-42 | 18-95 | 82-5 |
② 对照品溶液制备。取巴利森苷A、B、C、E对照品适量,精密称定,加50%甲醇混合溶液制成巴利森苷A 97 μg/mL、巴利森苷B 65 μg/mL、巴利森苷C 30 μg/mL、巴利森苷E 50 μg/mL的混合溶液,适当超声溶解后,得到对照品溶液。
③ 供试品溶液的制备。取天麻药材粉末约0.5 g,精密称定,置具塞锥形瓶中,精密加入50%甲醇25 mL,超声处理(功率500 W,频率40 kHz) 30 min,放冷,摇匀,滤过,量取续滤液,备用。
④ 标准曲线的建立。分别取5、10、15、20、25 μL的标准品溶液,按照色谱条件测定分析,以峰面积为纵坐标(y),浓度为横坐标(x),得到线性回归方程,标准曲线见表 5。
化合物 Components |
回归方程 Regression equation |
R2 | 线性范围/(μg/mL) Linear range/(μg/mL) |
Parishin A | y=993498x+10188 | 0.999 | 0.048 5-0.242 5 |
Parishin B | y=729616x-18983 | 1.000 | 0.032 5-0.162 5 |
Parishin C | y=262075x+2174.6 | 0.997 | 0.015 0-0.075 0 |
Parishin E | y=479631x-2687.7 | 1.000 | 0.025 0-0.125 0 |
⑤ 样品含量测定。吸取对照品和供试品各10 μL,按照1.3.6节①的色谱条件测定,得到巴利森苷A、B、C、E的峰面积并计算各巴利森苷的含量。
2 结果与分析 2.1 不同产地天麻的外观形状与生物量27份汉中天麻的块茎长、块茎粗、螺环纹数、鲜重与干重分别为77.6-160.0 mm、26.3-60.8 mm、11-19条、50.0-235.0 g、8.2-37.1 g,略阳县天麻的各项指标均值最高,分别为135.5 mm、51.4 mm、16.6条、176.0 g、27.3 g,运用SPSS 20.0以及Excel 2020进行数据分析以及图表制作,其中各地区天麻块茎长均值无显著差异,略阳县天麻在块茎粗、鲜重、干重均值上与宁强县天麻有显著差异,与勉县天麻差异不显著,勉县天麻与宁强县天麻差异不显著,见图 1。由图 2,勉县天麻的螺纹环数与宁强县天麻、略阳县天麻差异显著,宁强县天麻的螺纹环数与略阳县天麻差异不显著。
2.2 不同产地天麻的颜色色度值
27份汉中天麻的L*和E*值变化明显,如表 6所示,二者均为勉县>略阳县>宁强县,说明勉县天麻在外观上比宁强县、略阳县天麻明亮,总色度值较高;勉县天麻的L*值较大,平均值为86.68,外观明亮度较高。a*值以宁强县较高,b*值以略阳县较高,说明通过色度值大小可以区分汉中市不同地区的天麻。结合图 3,可以发现不同地区天麻色度值变化趋势不同,其中L*和E*值变化最为明显,说明不同地区红天麻的外观颜色明亮度和总色度值存在差异。
地区 Regions |
L*(D65) | a*(D65) | b*(D65) | E*(D65) |
宁强县 Ningqiang County |
79.56±0.63 | 2.68±0.15 | 16.44±0.28 | 81.31±0.60 |
勉县 Mian County |
86.68±0.71 | 2.31±0.32 | 15.63±1.48 | 85.57±2.89 |
略阳县 Lueyang County |
81.54±3.90 | 2.09±0.78 | 16.72±0.72 | 84.76±2.20 |
2.3 不同产地天麻的气味响应值
27份汉中天麻的电子鼻气味分析如图 4所示,不同地区天麻的气味响应值根据产区的不同,在峰形一致的基础上,峰强度有一定差异。
可根据气味响应值将汉中市不同地区天麻大致区分为4类(图 5),宁强县大部分天麻样品主要聚为一类,分布在黄色和绿色区域;勉县、略阳县天麻则在各个分类中均有分布,大部分聚集在蓝色区域。3个地区天麻样品的分布未见明显规律,聚类分析中代表点相对分散,气味响应值特征点较为相似。
2.4 不同产地天麻的成分含量 2.4.1 水分和浸出物含量
汉中天麻的水分、浸出物含量均达到了2020年版《中国药典》标准(图 6)。各地区天麻的水分含量为5.6%-9.4%,浸出物含量为15.4%-28.3%,其中,略阳县天麻水分、浸出物含量均值最高,分别为8.0%和22.0%。
2.4.2 天麻素和对羟基苯甲醇含量
由图 7可知,不同地区天麻的天麻素含量为0.13%-0.65%,其中宁强县天麻的天麻素平均含量最高,为0.35%。不同地区天麻的对羟基苯甲醇含量为0.02%-0.38%,其中勉县天麻的对羟基苯甲醇平均含量较高,为0.20%。不同地区天麻的天麻素和对羟基苯甲醇的总含量为0.26%-0.68%,均达到了2020年版《中国药典》规定标准,宁强县天麻的2种主要有效成分平均总含量最高,达到了0.44%。此外,勉县天麻的对羟基苯甲醇和天麻素的含量接近,2种成分含量比例和其他地区差异显著。
2.4.3 巴利森苷类成分含量
不同地区间天麻中巴利森苷E、B、C、A含量测定结果如图 8所示,不同地区天麻巴利森苷E、B、C、A均值含量均为宁强县天麻最高,分别为7.65、5.87、2.10和10.67 mg/g,总含量为26.30 mg/g。
2.5 PCA-熵权TOPSIS综合评价体系建立
利用PCA选取质量评价主成分,结合熵权法算得权重后,通过DPS数据统计软件采用TOPSIS法对天麻样品整体质量进行综合评价排名。
2.5.1 PCA提取主要指标选择外观性状指标:块茎长、块茎宽、鲜重、干重、螺环纹数、色度值(L*、a*、b*、E*);气味响应值指标:15、17、18、24 s峰值;化学指标:水分、浸出物、天麻素、对羟基苯甲醇,巴利森苷E、B、C、A含量共21个指标。根据各指标差异系数以及相关性最终选出16个指标建立评价体系,数据标准化后进行主成分分析,如表 7、表 8所示。
样品 Sample |
C1:块茎长 C1:Length of tuber |
C2:块茎粗 C2:Thick of tuber |
C3:鲜重 C3:Fresh weight |
C4:干重 C4:Dry weight |
C5:a* | C6:15 s峰值 C6:15 s intensity |
C7:17 s峰值 C7:17 s intensity |
C8:18 s峰值 C8:18 s intensity |
C9:24 s峰值 C9:24 s intensity |
1 | -0.457 | 0.216 | -0.485 | -0.332 | 0.778 | -0.961 | -1.196 | -0.645 | -1.168 |
3 | -1.224 | -1.297 | -1.611 | -1.477 | 1.293 | -0.605 | -0.858 | -0.509 | -0.514 |
4 | -0.430 | 1.689 | 0.918 | 0.803 | -0.902 | 1.728 | 0.544 | -0.145 | 2.117 |
5 | -1.161 | -1.332 | -1.680 | -1.371 | 1.087 | 0.413 | 0.811 | -0.334 | -1.005 |
6 | -1.424 | -0.493 | -0.679 | -0.129 | -0.959 | -1.067 | -1.110 | -0.681 | -1.004 |
7 | 0.563 | -0.410 | -0.456 | -0.430 | 0.755 | -0.168 | -0.485 | -0.413 | 1.218 |
8 | 0.216 | -0.848 | -0.785 | -0.973 | 0.047 | 0.719 | -0.938 | -0.556 | 0.013 |
9 | -0.315 | -2.053 | -1.506 | -1.185 | 0.607 | 0.554 | 0.085 | -0.422 | 1.165 |
13 | -2.187 | 0.748 | 0.729 | 0.571 | -0.285 | 0.392 | 0.163 | -0.429 | 1.213 |
15 | -0.057 | -0.500 | -0.269 | -0.605 | -0.868 | 1.622 | 0.202 | -0.192 | -0.066 |
16 | 0.616 | -0.840 | -0.694 | -0.754 | 0.047 | 0.612 | 0.123 | -0.292 | -0.307 |
18 | -0.294 | -0.209 | -0.362 | -0.372 | -0.079 | 0.582 | 3.681 | -0.193 | 1.165 |
19 | -0.273 | 0.494 | 0.410 | 0.336 | 1.487 | -0.481 | -0.643 | -0.500 | -0.812 |
20 | 0.794 | 0.725 | 0.743 | 1.717 | 0.767 | -0.507 | 0.191 | -0.555 | 0.410 |
21 | 2.146 | 0.914 | 0.953 | 1.154 | 0.778 | -0.671 | -0.971 | -0.690 | 0.099 |
22 | 1.257 | -0.209 | -0.083 | -0.494 | 0.161 | -1.390 | -0.634 | -0.661 | -1.075 |
25 | -1.587 | -1.557 | -1.052 | -0.927 | -0.719 | 0.802 | 0.285 | 0.729 | -1.043 |
30 | -0.316 | -0.418 | -1.004 | -1.049 | 0.127 | 2.630 | -0.025 | -0.621 | -0.791 |
mx-31 | 0.348 | 0.417 | 0.953 | 0.768 | -0.186 | -0.763 | -0.882 | -0.341 | -0.409 |
mx-32 | 0.024 | 0.046 | 0.515 | 0.232 | -0.110 | -1.127 | -0.578 | -0.677 | -0.938 |
mx-33 | 0.169 | 0.630 | 0.602 | 0.389 | -0.281 | -0.534 | 0.365 | -0.397 | -0.403 |
mx-34 | 0.005 | 0.375 | 0.712 | 0.167 | -2.004 | 0.638 | 1.195 | 0.161 | 0.497 |
mx-35 | 0.164 | 0.276 | 0.412 | 0.066 | -1.169 | -1.361 | 0.752 | 1.359 | 1.417 |
ly-36 | -0.194 | -0.375 | -0.209 | -0.602 | 2.459 | 0.411 | 0.825 | 1.903 | 1.914 |
ly-37 | 1.696 | 1.851 | 2.012 | 2.100 | -0.216 | -0.008 | 0.336 | 3.223 | -0.956 |
ly-38 | 1.496 | 2.029 | 2.151 | 2.239 | -1.154 | -0.761 | -0.290 | -0.275 | -0.572 |
ly-39 | 0.427 | 0.130 | -0.235 | 0.158 | -1.459 | -0.699 | -0.949 | 2.156 | -0.164 |
样品 Sample |
C10:水分 C10:Moisture |
C11:天麻素 C11:Gastrodin |
C12:对羟基苯甲醇 C12:P-hydroxybenzylalcohol |
C13:巴利森苷E C13:Parishin E |
C14:巴利森苷B C14:Parishin B |
C15:巴利森苷C C15:Parishin C |
C16:巴利森苷A C16:Parishin A |
1 | -0.263 | -0.238 | -0.226 | 1.224 | 0.096 | 0.613 | -0.249 |
3 | -0.927 | -0.073 | -0.921 | -0.771 | -0.279 | -0.269 | -0.800 |
4 | 0.135 | 0.503 | -0.921 | -0.396 | -1.855 | -1.351 | -1.401 |
5 | 0.907 | 1.902 | 0.468 | 1.367 | 2.314 | 2.959 | 2.878 |
6 | -0.642 | -0.814 | -0.782 | 1.176 | -0.985 | -0.657 | -0.822 |
7 | 1.517 | 0.338 | -0.087 | -0.031 | 0.323 | 0.364 | 0.025 |
8 | -0.459 | -0.485 | -0.643 | 0.469 | -0.559 | -0.561 | 0.138 |
9 | -1.798 | 0.421 | 0.329 | 0.486 | 1.569 | 1.191 | 1.617 |
13 | -1.235 | 0.091 | -0.365 | 1.062 | 0.930 | 0.747 | 0.654 |
15 | -0.890 | 2.807 | -1.198 | -0.229 | -0.391 | -0.552 | -0.196 |
16 | -1.310 | -0.402 | -0.782 | -0.841 | 0.551 | 0.258 | 0.624 |
18 | -0.840 | -0.238 | 0.746 | -0.977 | 0.610 | 0.475 | 0.373 |
19 | -0.505 | 0.174 | 0.468 | -0.170 | 0.645 | 1.747 | -0.796 |
20 | 0.490 | -0.402 | 0.190 | 0.414 | 1.239 | 0.998 | 1.070 |
21 | 1.929 | 0.009 | 0.051 | 2.363 | -0.255 | -0.181 | -0.216 |
22 | 1.140 | 0.585 | -0.226 | -1.323 | 0.748 | 0.610 | 0.181 |
25 | -0.866 | -0.402 | -0.782 | -0.688 | 0.335 | 0.086 | 0.087 |
30 | -0.127 | 2.231 | -0.087 | -0.952 | -0.056 | -0.056 | -0.252 |
mx-31 | -1.299 | -0.814 | 0.468 | 0.248 | 0.236 | 0.019 | 1.358 |
mx-32 | -0.131 | -0.978 | 0.885 | 0.982 | -0.769 | -1.033 | -0.811 |
mx-33 | 0.902 | -0.814 | 0.746 | 0.231 | -0.652 | -0.860 | -0.583 |
mx-34 | 1.095 | -0.731 | 0.329 | -1.364 | -0.598 | -0.689 | 0.346 |
mx-35 | 0.912 | -1.472 | 3.801 | -1.681 | -0.948 | -1.183 | -0.499 |
ly-36 | -0.266 | -1.143 | 1.023 | -0.416 | -1.836 | -1.269 | -1.725 |
ly-37 | 0.400 | 0.503 | -0.921 | 0.938 | 0.470 | 0.087 | 0.687 |
ly-38 | 1.056 | -0.731 | -0.782 | 0.089 | 0.692 | -0.184 | -0.237 |
ly-39 | 1.076 | 0.174 | -0.782 | -1.208 | -1.576 | -1.311 | -1.451 |
对各项指标数据进行KMO和Bartlett检验,KMO值为0.586>0.5,P<0.05,说明数据适宜进行主成分分析。如表 9所示,提取出6个主成分,累计方差贡献率达到82.415%,代表了绝大部分原始指标信息。由表 10可知,第1个主成分中,块茎长、块茎粗、鲜重、干重载荷值较大,因此第1主成分可作为形状与生物量的综合体现,定性为形状与生物量因子;同理,第2个主成分中巴利森苷E、B、C、A载荷值较大,定性为巴利森苷类因子;第3个主成分中气味响应值载荷值较大,定性为气味因子;第4个主成分中天麻素载荷值较大,定性为天麻素因子;第5个主成分中水分载荷值较大,定性为水分因子;第6个主成分中颜色值载荷值较大,定性为颜色因子[27]。
主成分 Principal component |
特征根 Characteristic root |
方差解释率/% Variance explained rate/% |
累积方差贡献率/% Cumulative variance contribution rates/% |
1 | 4.442 | 27.763 | 27.763 |
2 | 2.921 | 18.256 | 46.018 |
3 | 1.936 | 12.102 | 58.120 |
4 | 1.726 | 10.786 | 68.906 |
5 | 1.128 | 7.052 | 75.958 |
6 | 1.033 | 6.457 | 82.415 |
评价指标 Evaluation index |
主成分1 Principal component 1 |
主成分2 Principal component 2 |
主成分3 Principal component 3 |
主成分4 Principal component 4 |
主成分5 Principal component 5 |
主成分6 Principal component 6 |
C1 | 0.564 | 0.413 | 0.074 | 0.047 | 0.437 | 0.222 |
C2 | 0.824 | 0.340 | 0.163 | 0.227 | -0.195 | 0.067 |
C3 | 0.841 | 0.371 | 0.212 | 0.139 | -0.209 | -0.046 |
C4 | 0.781 | 0.499 | 0.168 | 0.142 | -0.236 | -0.039 |
C5 | -0.407 | 0.196 | -0.100 | -0.393 | -0.052 | 0.698 |
C6 | -0.380 | -0.368 | 0.359 | 0.652 | -0.106 | 0.226 |
C7 | -0.064 | -0.288 | 0.856 | -0.023 | -0.050 | -0.056 |
C8 | 0.421 | -0.251 | 0.250 | 0.066 | 0.263 | 0.012 |
C9 | 0.186 | -0.459 | 0.531 | -0.145 | -0.305 | 0.393 |
C10 | 0.531 | 0.268 | 0.061 | -0.057 | 0.581 | 0.217 |
C11 | -0.493 | 0.127 | 0.114 | 0.655 | 0.295 | 0.261 |
C12 | 0.171 | -0.149 | 0.39 | -0.760 | 0.115 | 0.005 |
C13 | -0.070 | 0.643 | -0.257 | -0.014 | -0.424 | 0.284 |
C14 | -0.572 | 0.687 | 0.342 | -0.063 | 0.036 | -0.163 |
C15 | -0.657 | 0.632 | 0.256 | -0.107 | 0.034 | 0.017 |
C16 | -0.533 | 0.567 | 0.445 | -0.037 | 0.025 | -0.239 |
如表 11所示,将载荷系数除以对应特征根的平方根,得到主成分特征向量,构建主成分函数[以形状与生物量因子为例,式(2)],并求得主成分矩阵[28]。
$ \begin{array}{r} F_{\text {shape } \& \text { biomass }}=0.268 \times \mathrm{C} 1+0.391 \times \mathrm{C} 2+ \\ 0.399 \times \mathrm{C} 3+0.370 \times \mathrm{C} 4-0.193 \times \mathrm{C} 5-0.180 \times \end{array} \\ \begin{aligned} & \mathrm{C} 6-0.030 \times \mathrm{C} 7+0.200 \times \mathrm{C} 8+0.088 \times \mathrm{C} 9+ \\ & 0.252 \times \mathrm{C} 10-0.234 \times \mathrm{C} 11+0.081 \times \mathrm{C} 12-0.033 \times \\ & \mathrm{C} 13-0.272 \times \mathrm{C} 14-0.312 \times \mathrm{C} 15-0.253 \times \mathrm{C} 16 。\end{aligned} $ | (2) |
评价指标 Evaluation index |
形状与生物量因子 Shape and biomass factors |
巴利森苷类因子 Parishin factors |
气味因子 Odor factor |
天麻素因子 Gastrodin factor |
水分因子 Water content factor |
颜色因子 Color Factor |
C1 | 0.268 | 0.242 | 0.053 | 0.036 | 0.411 | 0.218 |
C2 | 0.391 | 0.199 | 0.117 | 0.173 | -0.184 | 0.066 |
C3 | 0.399 | 0.217 | 0.152 | 0.106 | -0.197 | -0.045 |
C4 | 0.370 | 0.292 | 0.121 | 0.108 | -0.222 | -0.039 |
C5 | -0.193 | 0.115 | -0.072 | -0.299 | -0.049 | 0.686 |
C6 | -0.180 | -0.215 | 0.258 | 0.496 | -0.100 | 0.223 |
C7 | -0.030 | -0.169 | 0.615 | -0.017 | -0.047 | -0.055 |
C8 | 0.200 | -0.147 | 0.180 | 0.051 | 0.248 | 0.012 |
C9 | 0.088 | -0.269 | 0.381 | -0.110 | -0.287 | 0.387 |
C10 | 0.252 | 0.157 | 0.044 | -0.043 | 0.547 | 0.214 |
C11 | -0.234 | 0.075 | 0.082 | 0.499 | 0.278 | 0.257 |
C12 | 0.081 | -0.087 | 0.280 | -0.579 | 0.108 | 0.005 |
C13 | -0.033 | 0.376 | -0.185 | -0.011 | -0.399 | 0.280 |
C14 | -0.272 | 0.402 | 0.246 | -0.048 | 0.034 | -0.161 |
C15 | -0.312 | 0.370 | 0.184 | -0.082 | 0.032 | 0.016 |
C16 | -0.253 | 0.332 | 0.320 | -0.028 | 0.024 | -0.235 |
2.5.2 熵权TOPSIS法
运用熵权法[29]计算主成分矩阵权重(表 12),结合DPS统计软件进行TOPSIS分析。其中,颜色因子、巴利森苷类因子和气味因子所占权重较大,说明上述指标对评价结果影响较大。表 13为评价结果,从中可以看出宁强县、略阳县天麻综合质量较优。
主成分指标 Principal component index |
熵值 Entropy |
差异系数 Diversity factor |
权重/% Weight/% |
Shape and biomass factors | 0.967 5 | 0.032 5 | 13.85 |
Parishin factors | 0.954 0 | 0.046 0 | 19.64 |
Odor factor | 0.957 8 | 0.042 2 | 18.02 |
Gastrodin factor | 0.970 4 | 0.029 6 | 12.61 |
Water content factor | 0.974 7 | 0.025 3 | 10.81 |
Color factor | 0.941 2 | 0.058 8 | 25.07 |
样品 Sample |
D+ | D- | 相对接近程度CI | 名次 Ranking |
1 | 0.096 5 | 0.072 6 | 0.429 4 | 17 |
3 | 0.107 3 | 0.057 6 | 0.349 4 | 25 |
4 | 0.074 8 | 0.097 6 | 0.566 1 | 6 |
5 | 0.080 8 | 0.097 5 | 0.546 6 | 7 |
6 | 0.129 6 | 0.045 4 | 0.259 5 | 27 |
7 | 0.063 0 | 0.100 0 | 0.613 3 | 2 |
8 | 0.094 2 | 0.067 2 | 0.416 2 | 18 |
9 | 0.089 6 | 0.073 8 | 0.451 7 | 13 |
13 | 0.095 1 | 0.074 2 | 0.438 5 | 15 |
15 | 0.085 5 | 0.081 0 | 0.486 4 | 10 |
16 | 0.103 4 | 0.061 3 | 0.372 2 | 23 |
18 | 0.096 7 | 0.083 7 | 0.464 1 | 12 |
19 | 0.079 5 | 0.084 2 | 0.514 4 | 9 |
20 | 0.072 4 | 0.097 3 | 0.573 5 | 5 |
21 | 0.057 6 | 0.122 5 | 0.680 2 | 1 |
22 | 0.096 4 | 0.076 4 | 0.442 0 | 14 |
25 | 0.129 2 | 0.048 2 | 0.271 9 | 26 |
30 | 0.079 4 | 0.085 6 | 0.518 5 | 8 |
mx-31 | 0.111 6 | 0.067 0 | 0.375 3 | 22 |
mx-32 | 0.106 0 | 0.061 2 | 0.365 8 | 24 |
mx-33 | 0.092 9 | 0.070 9 | 0.433 0 | 16 |
mx-34 | 0.115 3 | 0.071 7 | 0.383 4 | 21 |
mx-35 | 0.114 8 | 0.072 5 | 0.387 1 | 19 |
ly-36 | 0.080 9 | 0.119 4 | 0.596 1 | 3 |
ly-37 | 0.071 5 | 0.103 7 | 0.592 0 | 4 |
ly-38 | 0.099 1 | 0.088 9 | 0.473 0 | 11 |
ly-39 | 0.107 8 | 0.067 1 | 0.383 7 | 20 |
Note:“D+” defines the distance between the i-th evaluation object and the maximum value, “D-”defines the distance between the i-th evaluation object and the minimum value. |
3 讨论
本研究在天麻传统性状指标的基础上,结合有效成分含量以及颜色、气味的定量数据对汉中不同地区天麻进行了比较分析与质量评价。在综合评价中为了有效消除主观影响,本研究运用PCA-熵权TOPSIS法构建了天麻质量综合评价体系[30],且筛选出了3个重要评价指标,但由于在略阳县、勉县采集的样品较少,天麻样品范围存在局限性,还需进一步完善与验证。
在外观性状上,3个地区所产天麻的块茎粗、螺环纹数、鲜重和干重均有一定的差异;颜色明亮度和总色度值有一定差别,通过色度值大小可以区分汉中市不同地区的天麻;不同地区天麻样品的气味响应值特征点较为相似。与李巧玲等[30]研究的不同产区天麻相比,汉中3县红天麻与湖北宜昌、安徽岳西、重庆南川红天麻块茎长度相近。略阳县红天麻鲜重及干重、勉县红天麻鲜重较安徽岳西红天麻低,较湖北宜昌、重庆南川高。宁强县红天麻鲜重较安徽岳西、重庆南川低,较湖北宜昌高。勉县红天麻干重较安徽岳西低,较湖北宜昌、重庆南川高。宁强县红天麻干重较安徽岳西低,较湖北宜昌高,与重庆南川相近[31]。
在化学成分上,汉中天麻中水分、浸出物含量均符合2020年版《中国药典》标准,宁强县天麻的天麻素和对羟基苯甲醇的总含量在3县天麻中最高,平均总含量达到了0.44%,勉县、略阳县总含量为0.40%左右,与其他文献研究结果基本一致[32, 33]。勉县天麻中天麻素和对羟基苯甲醇含量接近,2种成分含量比例和其他地区差异显著。3县天麻中宁强县天麻的巴利森苷E、B、C、A及总含量均为最高。与张卫等[33]研究的贵州天麻相比,汉中3县红天麻水分含量与贵州红天麻相近,浸出物含量更低。与薛华丽等[34]研究的雪峰山区域不同产地天麻相比,汉中3县红天麻浸出物含量与湖南洪江、溆浦以及湖北罗田等地红天麻相近,天麻素及对羟基苯甲醇总含量与湖南红天麻相似,较湖北红天麻含量高。此外,与周媛等[35]研究的安徽金寨红天麻相比,汉中3县红天麻的天麻素含量较低,对羟基苯甲醇、巴利森苷E含量较高;宁强县红天麻的巴利森苷B含量较高,勉县、略阳县红天麻较低;巴利森苷C含量方面,宁强县红天麻与安徽金寨天麻相近,勉县、略阳县较低。与王庆等[36]研究的西南不同产区红天麻相比,宁强县红天麻的天麻素含量较云南昭通、四川广元红天麻低,较贵州大方红天麻高,与贵州德江红天麻相近,勉县、略阳县较低;对羟基苯甲醇含量方面,勉县红天麻较高,宁强县、略阳县红天麻较四川广元红天麻低,较贵州德江红天麻高,与云南昭通、贵州大方红天麻相近。
田孟华等[37]表示可将巴利森苷A、巴利森苷E作为乌天麻(G.eleta Bl.f.glauca S.Chow)与其他天麻品种差异性分析的标志物。巴利森苷类化合物有生理活性和多种药理作用,且与天麻素之间存在相互转化的关系[4, 13],间接提高天麻素含量。张琦[5]在分析引种前后天麻6种主要活性成分含量变化时,发现引种后天麻中天麻素、巴利森苷A、B、C含量有明显的相关性。由此可见,在对天麻开展质量评价时,可将巴利森苷类成分含量纳入评价指标。在本研究建立的评价体系中,颜色因子、气味因子、巴利森苷类因子权重较大,说明在质量评价中应重视上述指标。
4 结论本研究采用PCA结合熵权TOPSIS法对汉中天麻整体质量进行了综合评价,探究了天麻外观性状与内在成分指标所占权重,发现颜色、气味、巴利森苷类因子权重较大,应注重上述指标在质量评价中的重要性,在此评价体系下,宁强县、略阳县天麻排名靠前。本研究为天麻质量评价提供了参考,也为进一步的分析研究提供了理论基础,从而使评估天麻质量的因素可以被更加科学地把握。
[1] |
张进强, 肖承鸿, 周涛, 等. 神奇的天麻[J]. 中国食品药品监管, 2021(2): 114-119. |
[2] |
国家药典委员会. 中华人民共和国药典: 一部[M]. 2020年版. 北京: 中国医药科技出版社, 2020.
|
[3] |
黄璐琦, 詹志来, 郭兰萍. 中药材商品规格等级标准汇编[M]. 北京: 中国中医药出版社, 2019.
|
[4] |
田紫平, 肖慧, 冯舒涵, 等. 天麻有效成分巴利森苷的降解规律分析[J]. 中国实验方剂学杂志, 2017, 23(23): 18-21. |
[5] |
张琦. 天麻共生菌生物学特性及天麻引种前后成分变化的研究[D]. 北京: 北京协和医学院, 2020.
|
[6] |
高原, 方妍, 单梦瑶, 等. 基于色差原理分析不同产地桑白皮有效成分含量与颜色的相关性[J]. 中国药房, 2021, 32(2): 213-219. |
[7] |
王杨, 甄臻, 隆毅, 等. 基于颜色变化的酒白芍质量标准及炮制工艺研究[J]. 中药材, 2020, 43(9): 2141-2145. |
[8] |
宿莹, 李翟, 侯晓琳, 等. 基于色差原理分析龙胆有效成分含量与颜色的相关性[J]. 中国实验方剂学杂志, 2019, 25(13): 151-156. |
[9] |
宿莹, 侯晓琳, 刘战, 等. 基于色差原理分析黄柏有效成分含量与颜色的相关性[J]. 中药材, 2019, 42(8): 1766-1770. |
[10] |
张冬月, 吴浩善, 李思雨, 等. 电子鼻技术鉴别通关藤真伪及其道地性[J]. 中成药, 2018, 40(10): 2228-2233. |
[11] |
AMKOR A, EL BARBRI N. A measurement prototype based on gas sensors for detection of pesticide residues in edible mint[J]. Journal of Food Measurement and Characterization, 2021, 15(1): 170-180. DOI:10.1007/s11694-020-00617-8 |
[12] |
PENG L, ZOU H Q, BAUER R, et al. Identification of Chinese herbal medicines from zingiberaceae family using feature extraction and cascade classifier based on response signals from e-nose[J]. Evidence-Based Complementray and Alternative Medicine: ECAM, 2014, 2014: 963035. |
[13] |
张琦, 王兆丰, 王秋颖. 天麻中巴利森苷类化合物的研究概述[J]. 中国现代中药, 2020, 22(1): 148-153. |
[14] |
覃桂, 葛锦蓉, 周锐, 等. 基于PCA-Logistic回归分析的艾叶药材等级评价方法研究[J]. 医药导报, 2023, 42(3): 317-321. |
[15] |
孔浩, 郭庆梅, 王慧慧, 等. 主成分分析法在中药质量评价中的应用[J]. 辽宁中医杂志, 2014, 41(5): 890-892. |
[16] |
石岩, 魏锋, 马双成. 关于主成分分析在中药和天然药物分析研究中应用的探讨[J]. 中国中药杂志, 2018, 43(14): 3031-3035. |
[17] |
赵梦利, 刘妍如, 宋忠兴, 等. 基于成分-抗氧化活性相关的鸡血藤饮片等级评价研究[J]. 中草药, 2020, 51(4): 943-949. |
[18] |
巫晓霞, 黄群, 李羿. 基于多元统计学的茯苓质量综合评价[J]. 化学研究与应用, 2022, 34(10): 2472-2476. |
[19] |
郝佳旭, 李元增, 范晓, 等. 基于熵权TOPSIS法综合评价草果质量[J]. 中国药房, 2022, 33(17): 2087-2092. |
[20] |
张梦奇, 李彬, 王琼芬, 等. 基于熵权逼近理想解排序法的西红花质量综合评价[J]. 医药导报, 2023, 42(12): 1825-1830. |
[21] |
吴红伟, 李东辉, 张育贵, 等. 基于熵权TOPSIS模型对不同加工方式黄芪药材质量的综合评价[J]. 中国药学杂志, 2021, 56(16): 1325-1331. |
[22] |
高伟城, 王小平, 陈荣珠, 等. 基于化学模式识别和熵权TOPSIS法的不同品种枇杷叶品质分析研究[J]. 甘肃中医药大学学报, 2023, 40(5): 33-39. |
[23] |
李雨昕, 邢娜, 白浩东, 等. 基于熵权法的灰色关联法-TOPSIS法对不同产地三七及其炮制品质量的评价研究[J]. 中草药, 2023, 54(4): 1252-1259. |
[24] |
史丛晶, 白浩东, 李雨昕, 等. 基于熵权法结合灰色关联法和TOPSIS法对不同产地知母及盐知母质量的评价研究[J]. 中国现代应用药学, 2023, 40(7): 873-880. |
[25] |
姚芳. 基于PCA-EWM-TOPSIS模型的绿色食品企业财务风险评估[J]. 广东石油化工学院学报, 2022, 32(6): 81-86. |
[26] |
李俊杰, 李钢铁, 麻云霞, 等. 基于PCA- TOPSIS的沙地不同林龄榆树土壤养分评价[J]. 中国水土保持科学(中英文), 2022, 20(5): 17-23. |
[27] |
李伟, 郜海燕, 陈杭君, 等. 基于主成分分析的不同品种杨梅果实综合品质评价[J]. 中国食品学报, 2017, 17(6): 161-171. |
[28] |
邵明灿, 胡花丽, 王毓宁, 等. 基于主成分分析法分析一氧化氮对绿芦笋贮藏品质的影响[J]. 食品科学, 2012, 33(20): 318-322. |
[29] |
肖璞, 吴潇威. 基于熵值法的评价模型构造与应用[J]. 现代计算机, 2021, 27(30): 50-54. |
[30] |
李巧玲, 肖波, 邓才富, 等. 不同产区天麻种质资源收集保存、评价及优选[J]. 西南大学学报(自然科学版), 2020, 42(12): 54-59. |
[31] |
仲瑞雪, 陈韵蓓, 段小云, 等. 天麻商品等级与天麻素含量相关性研究[J]. 亚太传统医药, 2016, 12(17): 27-30. |
[32] |
冯永辉, 马朝英. 陕滇两地天麻中天麻素的含量比较[J]. 光谱实验室, 2013, 30(4): 1998-2000. |
[33] |
张卫, 刘春艳, 谢宇, 等. 贵州天麻的品质评价研究[J]. 亚太传统医药, 2019, 15(8): 66-69. |
[34] |
薛华丽, 田玉桥, 邱小燕, 等. 雪峰山区域不同产地天麻化学成分含量比较研究[J]. 中国现代医药杂志, 2022, 24(12): 23-26. |
[35] |
周媛, 周娜, 王泽旭, 等. 安徽金寨红天麻、红乌天麻高效液相色谱法指纹图谱与化学成分比较[J]. 安徽医药, 2023, 27(4): 683-690. |
[36] |
王庆, 李丹丹, 陈艾萌, 等. 西南不同产区3种天麻变型主要化学成分含量比较[J]. 中草药, 2018, 49(11): 2646-2652. |
[37] |
田孟华, 袁天军, 周瑞, 等. 不同产地及变型天麻有效成分差异性分析[J]. 中成药, 2020, 42(7): 1824-1829. |