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

1.   材料与方法

1.1   样品采集与数据收集

通过文献检索和实地调查，本研究于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  夏枯草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

下载: 导出CSV 

| 显示表格

1.2   仪器与试剂

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

1.3   测定方法

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

1.4   数据分析

模糊数学隶属函数分析方法能够避免单一指标评价的片面性问题，更加全面客观地评定结果，为此本文采用隶属函数综合评价法对夏枯草主要药用活性成分含量进行综合评价。隶属函数值[ $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.   结果与分析

2.1   不同地理种源夏枯草主要药用成分含量差异

不同地理种源夏枯草多糖和黄酮含量的方差分析结果见表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.  Polysaccharide and flavonoid contents in Prunella vulgaris from different geographical provenances w/(mg·g⁻¹)

种源 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)

下载: 导出CSV 

| 显示表格

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

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.  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)

下载: 导出CSV 

| 显示表格

2.3   不同种源夏枯草主要药用成分含量的综合评价

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

表  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

下载: 导出CSV 

| 显示表格

2.4   不同地理种源夏枯草聚类分析

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

图  1  10个地理种源聚类结果

Figure  1.  Clustering results of ten geographical provenances

下载: 全尺寸图片 幻灯片

3.   讨论与结论

3.1   不同地理种源夏枯草的品质差异

夏枯草广泛分布在我国各地，集中分布于淮河流域及长江中下游地区^[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的地理种源。

3.2   不同地理种源夏枯草药用成分含量与环境因子的相关性

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

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

3.3   夏枯草地理种源选择及利用

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