夏枯草Prunella vulgaris L.为唇形科Lamiaceae夏枯草属Prunella植物，又名铁线夏枯、麦夏枯和麦穗夏枯草等，是我国常用中药材，全株可入药，具有清肝明目、消肿散结等功效^[1-2]。夏枯草的主要生物活性成分为三萜、皂苷、酚酸、黄酮及多糖类，具有降压、降糖、抗氧化、抗病毒、抗肿瘤等多种药理活性^[3-7]。夏枯草为药食两用大宗中药材，在中成药生产、临床配方、保健食品和凉茶生产中作原料使用。近年来，由于市场需求量剧增，大量野生夏枯草资源遭到挖掘，致使野生资源保有量下降，市场价格上涨，但人工栽培夏枯草远远不能满足市场需求，野生资源仍是目前夏枯草药材原料的重要来源。夏枯草广泛分布于我国陕西、甘肃、新疆、河南、湖北、湖南、江西、浙江、福建、台湾、广东、广西、贵州、四川及云南等省份，生长于荒坡、草地、溪边及路旁的湿润土地，分布海拔可达3000 m^[2]。目前，对夏枯草的研究主要集中在药用活性物质提取^[8-9]、药理作用^[10-11]、人工栽培种植^[12-13]和种质资源遗传多样性分析^[14]等方面，系统评价不同种源药材质量及影响质量的环境因子构成方面的研究较少。本研究收集了10个野生夏枯草种源，分析其药用成分含量，利用方差分析、相关性分析、聚类分析和隶属函数综合评价法，筛选出主要有效成分影响因子和高品质地理种源，为夏枯草种源筛选、品种选育和人工规模化种植提供参考依据。

通过文献检索和实地调查，本研究于2018—2019年的5—6月间收集野生夏枯草植株，10个地理种源采样点的基本情况见表1，表1中的气候资源数据来源于中国自然资源数据库( http://www.data.ac.cn/zy/show/shi3.asp)，经纬度和海拔数据在调查时实地记录。收集时要求植株长势正常，处于盛花期，无病虫害，无明显缺陷，每个地理种源收集15株以上。每个地理种源设置3个采集点，各采集点间隔10 m，每个采集点采集5株。采集的夏枯草样品在60 ℃条件下烘干粉碎，过60目孔筛，然后各采集点称取10 g样品，将相同地理种源的样品均匀混合，干燥保存，供药用有效成分——多糖和黄酮含量测试。

表 1 Table 1

表 1 夏枯草10个野生种源采样点地理位置及主要气候因子 Table 1 Geographical location and main climatic factors of ten wild population sampling points of Prunella vulgaris 种源 Provenance 经度 Longitude 纬度 Latitude 海拔/m Altitude 年平均气温/℃ Annual mean air temperature 年日照时长/h Annual sunshine hour 年降水量/mm Annual precipitation 江苏宜兴 Yixing, Jiangsu E119°48′37″ N31°20′48″ 18 16.7 1 807.5 1 805.4 浙江丽水 Lishui, Zhejiang E119°04′21″ N28°38′57″ 450 17.8 1 676.6 1 568.4 安徽黄山 Huangshan, Anhui E118°18′12″ N29°42′09″ 759 15.5 1 750.3 1 670.1 江西分宜 Fenyi, Jiangxi E114°36′45″ N27°39′50″ 250 17.2 1 535.3 1 643.6 湖南衡阳 Hengyang, Hunan E112°36′20″ N26°56′12″ 613 18.2 1 688.9 1 510.8 湖北宜昌 Yichang, Hubei E111°33′42″ N30°44′18″ 533 16.9 1 710.5 1 215.6 广东清远 Qingyuan, Guangdong E112°20′11″ N25°10′10″ 648 18.0 2 290.0 1 329.0 广西桂林 Guilin, Guangxi E110°35′28″ N25°26′32″ 309 18.9 1 670.0 1 949.5 福建三明 Sanming, Fujian E117°29′18″ N26°33′24″ 548 18.2 1 727.1 1 688.0 福建宁化 Ninghua, Fujian E116°31′49″ N25°59′47″ 527 16.5 1 757.0 1 750.0

表 1 夏枯草10个野生种源采样点地理位置及主要气候因子 Table 1 Geographical location and main climatic factors of ten wild population sampling points of Prunella vulgaris

UV-2450紫外−可见光分光光度计(日本岛津公司)，M17035数显电热恒温水浴锅(北京中西远大科技有限公司)，JA10003N数显电子分析天平(上海精其仪器有限公司)，FQ-1042变频超声波清洗器(杭州法兰特超声波科技有限公司)，LC-100PHP型高效液相色谱仪(上海伍丰科学仪器有限公司)，1810-B型石英自动双重纯水整流器(常州万合仪器制造有限公司)。葡萄糖、芦丁、苯酚、乙醇、浓硫酸、亚硝酸钠、甲醇、氯化铝、氢氧化钠等均为分析纯，试验用水为超纯水。

夏枯草多糖含量的测定参照席与斌等^[15]的硫酸−苯酚比色法，黄酮含量的测定参照邓斌等^[16]的亚硝酸钠−硝酸铝比色法。在测定夏枯草多糖和黄酮的含量时，每个地理种源称取1 g样品测定多糖和黄酮的含量，每种成分含量的每次测定值取3个样品的平均值，重复5次，取各次测得的平均值作为多糖和黄酮的含量。

模糊数学隶属函数分析方法能够避免单一指标评价的片面性问题，更加全面客观地评定结果，为此本文采用隶属函数综合评价法对夏枯草主要药用活性成分含量进行综合评价。隶属函数值[ $U({x_i}) $ ]按公式 $U({x_i}) = ({x_i} - {x_{\min }})/({x_{\max }} - {x_{\min }})$ 计算，其中， ${x_i}$ 为指标测定值， ${x_{\max }}$ 和 ${x_{\min }}$ 分别为该指标的最大值和最小值。使用Microsoft Office Excel 2016整理数据和制作相关图表，使用SPSS19.0社会统计分析软件中的One-way ANOVA模块进行方差分析，比较不同种源夏枯草多糖和黄酮含量的差异，并进行相关性分析和聚类分析。

不同地理种源夏枯草多糖和黄酮含量的方差分析结果见表2。结果显示，多糖和黄酮含量在不同地理种源之间存在显著差异(P<0.05)，不同地理种源夏枯草多糖的质量分数为70.45~120.39 mg·g⁻¹，黄酮的质量分数为34.40~59.04 mg·g⁻¹。在夏枯草多糖含量方面，湖北宜昌种源的含量最高，质量分数为120.39 mg·g⁻¹，显著高于其他地理种源；江西分宜种源的多糖含量次之，江苏宜兴种源的多糖含量位列第3，质量分数分别为116.33和113.00 mg·g⁻¹；浙江丽水种源的多糖含量最少，质量分数仅为70.45 mg·g⁻¹，显著低于其他地理种源。在夏枯草黄酮含量方面，广西桂林种源的含量最高，质量分数为59.04 mg·g⁻¹，显著高于其他地理种源；浙江丽水种源的黄酮含量次之，湖南衡阳种源的黄酮含量位列第3，质量分数分别为56.03和55.90 mg·g⁻¹；广东清远种源的黄酮含量最低，质量分数仅为34.40 mg·g⁻¹，显著低于其他地理种源。

表 2 Table 2

表 2 不同地理种源夏枯草多糖和黄酮含量1) Table 2 Polysaccharide and flavonoid contents in Prunella vulgaris from different geographical provenances w/(mg·g−1) 种源 Provenance 多糖 Polysaccharide 黄酮 Flavonoid 种源 Provenance 多糖 Polysaccharide 黄酮 Flavonoid 江苏宜兴 Yixing, Jiangsu 113.00±12.79e 53.72±2.08f 湖北宜昌 Yichang, Hubei 120.39±6.36f 36.84±1.92b 浙江丽水 Lishui, Zhejiang 70.45±4.85a 56.03±1.62g 广东清远 Qingyuan, Guangdong 98.75±3.35c 34.40±2.12a 安徽黄山 Huangshan, Anhui 97.63±6.85c 50.44±1.18e 广西桂林 Guilin, Guangxi 105.47±6.43d 59.04±1.83h 江西分宜 Fenyi, Jiangxi 116.33±6.44ef 47.66±1.30d 福建三明 Sanming, Fujian 96.52±3.54bc 46.79±1.88d 湖南衡阳 Hengyang, Hunan 106.44±3.83d 55.90±1.39g 福建宁化 Ninghua, Fujian 91.95±3.82b 42.92±2.25c 　1) 相同指标数据后的不同小写字母表示不同种源间差异显著(P<0.05，Duncan’s法) 　1) Different lowercase letters after the data of the same index represent significant differences among different geographical provenances (P<0.05, Duncan’s method)

表 2 不同地理种源夏枯草多糖和黄酮含量1) Table 2 Polysaccharide and flavonoid contents in Prunella vulgaris from different geographical provenances

相关性分析结果表明，夏枯草多糖含量和黄酮含量呈现极显著负相关(r=−0.212，P=0.009)，说明不同地理种源夏枯草多糖和黄酮这2种次生代谢物的积累不可同时兼得，多糖和黄酮的代谢机制存在相互抑制的可能。因此，在筛选夏枯草高品质种源时应该综合考虑这2个指标，而在临床使用或药物栽培生产方面可根据需要选择目标地理种源。

夏枯草2个主要药用成分含量与各地理种源环境因子的相关性分析见表3，相关性顺序为年降水量(r = 0.793)>经度(r = 0.673)>年日照时长(r = 0.666)>海拔(r = 0.660)>年平均气温(r = 0.327)>纬度(r = 0.321)。其中，经度与夏枯草多糖含量呈极显著负相关(r = −0.427，P<0.01)，与黄酮含量呈极显著正相关(r = 0.246，P<0.01)；年日照时长与夏枯草黄酮含量呈极显著负相关(r = −0.557，P<0.01)，与多糖含量相关性不显著；年降水量和年平均气温与夏枯草黄酮含量呈现极显著正相关(r =0.674、0.239，P<0.01)，与多糖含量相关性不显著；纬度与夏枯草多糖含量呈现极显著正相关(r = 0.247，P<0.01)，与黄酮含量相关性不显著；海拔与夏枯草多糖含量和黄酮含量均呈现极显著负相关(r = −0.295、−0.365，P<0.01)。

表 3 Table 3

表 3 夏枯草多糖和黄酮含量与环境因子的相关性1) Table 3 Correlation of the contents of Prunella vulgaris polysaccharide and flavonoid with environmental factors 成分 Ingredient 项目 Project 经度 Longitude 纬度 Latitude 海拔 Altitude 年平均气温 Annual mean air temperature 年日照时长 Annual sunshine hour 年降水量 Annual precipitation 多糖 Polysaccharide r −0.427** 0.247** −0.295** −0.088 −0.109 −0.121 P 0.000 0.002 0.000 0.285 0.186 0.141 黄酮 Flavonoid r 0.246** 0.074 −0.365** 0.239** −0.557** 0.674** P 0.002 0.366 0.000 0.003 0.000 0.000 总和 Total 0.673 0.321 0.660 0.327 0.666 0.793 　1) “**”表示在0.01水平显著相关(Pearson相关) 　1) “**” indicates significant correlation at 0.01 level (Pearson correlation)

表 3 夏枯草多糖和黄酮含量与环境因子的相关性1) Table 3 Correlation of the contents of Prunella vulgaris polysaccharide and flavonoid with environmental factors

依据多糖和黄酮含量测定结果，运用模糊数学隶属函数法对各地理种源夏枯草品质进行综合评价。多糖和黄酮均赋予一样的权重，根据隶属函数计算公式确定各指标的隶属函数值，然后依据平均隶属函数值大小对各地理种源进行排序，进而综合评价各地理种源夏枯草品质优劣。各种源综合评价从优到劣的排序为广西桂林、江苏宜兴、湖南衡阳、江西分宜、安徽黄山、湖北宜昌、福建三明、浙江丽水、福建宁化和广东清远(表4)。广西桂林种源的综合评价排名最高(平均隶属函数值为0.850 6)，其黄酮含量为所有种源中最优，多糖含量亦在平均值(101.69 mg·g⁻¹)以上；广东清远种源综合评价结果最差(平均隶属函数值为0.283 3)，其黄酮含量为所有种源中最低，多糖含量亦在平均值(101.69 mg·g⁻¹)以下。

表 4 Table 4

表 4 不同种源夏枯草品质的隶属函数评价 Table 4 Evaluation of membership functions of Prunella vulgaris from different provenances 种源 Provenance 多糖含量隶属函数值 Polysaccharide content membership function value 黄酮含量隶属函数值 Flavonoid content membership function value 平均隶属函数值 Average membership function value 平均隶属函数值排序 Average membership function value ranking 广西桂林 Guilin, Guangxi 0.701 2 1.000 0 0.850 6 1 江苏宜兴 Yixing, Jiangsu 0.852 0 0.784 1 0.818 1 2 湖南衡阳 Hengyang, Hunan 0.720 7 0.872 6 0.796 7 3 江西分宜 Fenyi, Jiangxi 0.918 7 0.538 1 0.728 4 4 安徽黄山 Huangshan, Anhui 0.544 3 0.651 0 0.597 7 5 湖北宜昌 Yichang, Hubei 1.000 0 0.099 0 0.549 5 6 福建三明 Sanming, Fujian 0.522 0 0.502 8 0.512 4 7 浙江丽水 Lishui, Zhejiang 0 0.877 8 0.438 9 8 福建宁化 Ninghua, Fujian 0.430 5 0.345 8 0.388 2 9 广东清远 Qingyuan, Guangdong 0.566 7 0 0.283 3 10

表 4 不同种源夏枯草品质的隶属函数评价 Table 4 Evaluation of membership functions of Prunella vulgaris from different provenances

对10个夏枯草地理种源进行聚类分析(图1)，基于多糖和黄酮含量将10个种源划分为3大类：第1大类仅包括浙江丽水种源，其多糖含量最低；第2大类包括安徽黄山、福建三明、福建宁化和广东清远种源，其多糖含量和黄酮含量较接近平均值；第3大类包括湖南衡阳、广西桂林、江苏宜兴、湖北宜昌和江西分宜种源，其多糖和黄酮含量较高。进一步分析发现，地理位置较为接近的聚为一类，说明环境因子对夏枯草多糖和黄酮的代谢积累有较大影响。同时，聚类的结果与综合评价结果基本一致，说明种源和产地对夏枯草品质的影响较为明显。

图 1 10个地理种源聚类结果 Fig. 1 Clustering results of ten geographical provenances

夏枯草广泛分布在我国各地，集中分布于淮河流域及长江中下游地区^[2]。通过对分布区域内10个夏枯草野生地理种源的主要药用成分含量进行对比研究，发现夏枯草多糖和黄酮含量在不同地理种源之间存在显著差异，这与刘光敏等^[17]对不同产地夏枯草属、皮胜玲等^[18]和郭巧英等^[19]对不同产地夏枯草野生品种有效成分含量的研究结果相吻合。药用成分含量是药材品质的客观反映，本研究发现，夏枯草多糖含量最高的种源是湖北宜昌(质量分数为120.39 mg·g⁻¹)，比所有地理种源含量的平均值(101.69 mg·g⁻¹)高18.39%；夏枯草黄酮含量最高的种源是广西桂林(质量分数为59.04 mg·g⁻¹)，比所有地理种源含量的平均值(48.37 mg·g⁻¹)高22.06%。隶属函数评价中，综合品质最好的种源是广西桂林，其平均隶属函数值为0.850 6，综合品质排名2、3、4位的种源依次为江苏宜兴、湖南衡阳和江西分宜，其平均隶属函数值均在0.7以上，明显好于其他地理种源。因此，依据本文研究结果，在夏枯草产业化种植中，若是以夏枯草多糖含量为品质衡定指标，建议选用湖北宜昌种源；若是以夏枯草黄酮含量为品质衡定指标，建议选用广西桂林种源；若是兼顾考虑多糖和黄酮含量，建议选用研究结果中推荐的排名前4的地理种源。

10个研究种源均分布在我国热带和亚热带地区，气候湿润，光照充足，经度范围在E110°35′28″~E119°48′37″，纬度范围在N25°10′10″~N31°20′48″，海拔范围为18~759 m，这些环境因子的差异是导致夏枯草多糖和黄酮含量差异的外在原因。本研究发现，经度、纬度和海拔对夏枯草多糖含量有极显著影响，经度、海拔、年平均气温、年日照时长和年降水量对夏枯草黄酮含量有极显著影响。经度和海拔与夏枯草多糖含量呈极显著负相关，与黄酮含量分别呈极显著正相关和极显著负相关。随着经度的增加，夏枯草中多糖含量降低，黄酮含量急剧增加；随着海拔的不断升高，夏枯草中多糖和黄酮含量均下降；随着纬度的增加，夏枯草多糖含量不断增加，黄酮含量变化不明显；随着年平均气温和年降水量不断增加，夏枯草黄酮含量显著增加，多糖含量变化不明显；随着年日照时长增加，夏枯草黄酮含量显著降低，多糖含量变化不明显。

通过聚类分析，10个地理种源可分为3大类：第1大类为浙江丽水种源，第2大类为安徽黄山、福建三明、福建宁化和广东清远种源，第3大类包括湖南衡阳、广西桂林、江苏宜兴、湖北宜昌和江西分宜种源。聚类结果基本依据地理距离的远近，说明环境因子对夏枯草黄酮和多糖的次生代谢合成影响较大。

本研究发现，夏枯草多糖和黄酮含量呈现显著的负相关关系(r=−0.212，P<0.01)。在初生代谢比较旺盛的条件下，多糖大量积累，皂苷、黄酮类、生物碱等次生代谢合成受到抑制的现象在多花黄精、黄芪和金线莲等多种中药材中也有出现^[20-24]。经纬度和海拔是影响夏枯草多糖含量的主要环境因子，经度、海拔、年平均气温、年日照时长和年降水量是影响夏枯草黄酮含量的主要环境因子，这些环境因子在影响多糖和黄酮合成方面的趋势表现不一。因此，在选择夏枯草产业开发区域时，应综合考虑各环境因子对次生代谢的影响，才能确保生产高品质药材。就高品质药材种源选择而言，本文在综合考虑2个主要药用成分含量的基础上，使用隶属函数综合评价法，对10个地理种源的综合品质进行排名，筛选出品质最优的地理种源为广西桂林种源，本研究结果可为夏枯草引种栽培和遗传改良提供选择依据和参考。