牛膝药材的红外指纹图谱建立及多元统计分析
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
牛膝药材的红外指纹图谱建立及多元统计分析
作者：贾豪雷益铭张维方雷敬卫杨春静李莹莹谢彩侠
来源：《中国药房》2022年第02期
中圖分类号 R284 文献标志码 A 文章编号 1001-0408（2022）02-0153-07
DOI 10.6039/j.issn.1001-0408.2022.02.05
摘要目的建立不同产地牛膝药材的红外指纹图谱，并进行多元统计分析。
方法采用Spectrum for Window 3.02和OMNIC 9.2软件建立61批牛膝药材样品的红外指纹图谱;以红外指纹图谱共有峰的相对峰高为变量，采用Excel 2016软件进行正态分布分析，采用SPSS 22.0软件进行聚类分析和主成分分析并计算综合得分，采用SIMCA 14.1软件进行正交偏最小二乘法-判别分析，以变量重要性投影（VIP）>1为标准，筛选影响牛膝药材成分质量的标志性波数。
结果 61批牛膝药材样品红外光谱图的相关系数为0.967 2～0.997 7;共有13个共有峰。
正态分布分析结果显示，河南产与河北产牛膝药材共有峰相对峰高的正态分布曲线未有交叉，河
南产与内蒙古产牛膝药材的正态分布曲线存在交叉。
聚类分析结果显示，当组间距离为15
时，61批牛膝药材样品可聚为3类，其中N1～N12聚为一类，N13～N45聚为一类，N46～
N61聚为一类。
主成分分析结果显示，前3个主成分的累计方差贡献率为91.121%;河南省焦作市驾步村产牛膝药材（编号N40）的综合评分最高（2.39），河北省安国市新安村产牛膝（编号N4）的综合评分最低（-2.89）。
正交偏最小二乘法-判别分析结果显示，61批牛膝药材样品可分为3类，其中N1～N12为一类，N13～N28为一类，N29～N61为一类;共筛选出7个影响药材样品质量的标志性波数，其VIP值从大到小对应的波数依次为1 059、927、2 933、813、1 732、1 128、3 367 cm-1，其中1 732 cm-1处为皂苷类成分的特征吸收峰，1 059、1 128、927 cm-1处为糖苷类成分的特征吸收峰。
结论红外指纹图谱结合正态分布分析、聚类分析、主成分分析和正交偏最小二乘法-判别分析可用于鉴别不同产地牛膝药材。
关键词牛膝;红外指纹图谱;正态分布分析;聚类分析;主成分分析;正交偏最小二乘法-判别分析;产地
Establishment of infrared fingerprints and multivariate statistical analysis of Achyranthes bidentata
JIA Hao1，2，LEI Yiming3，ZHANG Weifang1，2，LEI Jingwei1，2，YANG
Chunjing2，4，LI Yingying1，2，XIE Caixia1，2 （1. School of Pharmacy， Henan University of Chinese Medicine， Zhengzhou 450046， China; 2. Henan Provincial Engineering Technology Research Center for TCM Quality Control and Evaluation， Zhengzhou 450046， China; 3. College of Orthopedics and Traumatology， Henan University of Chinese Medicine， Zhengzhou 450046，China; 4. Section of Clinical Pharmacy， Pharmacy Department， the Third Affiliated Hospital of Henan University of Chinese Medicine， Zhengzhou 450046， China）
ABSTRACT OBJECTIVE To establish the infrared fingerprints of Achyranthes bidentata from different producing areas， and to conduct multivariate statistical analysis. METHODS The infrared fingerprints of 61 batches of A. bidentata samples were established by Spectrum for Window 3.02 and OMNIC 9.2 software. Taking the relative peak height of common peaks of infrared fingerprint as the variable， the normal distribution analysis was carried out by Excel 2016 software; SPSS 22.0 software was used for cluster analysis and principal component analysis， and the comprehensive score was calculated; the orthogonal partial least squares-discriminant analysis was carried out by SIMCA 14.1 software， and the marker wave numbers affecting the quality of A. bidentata were screened by taking the variable importance in projection （VIP） >1 as the standard. RESULTS The correlation coefficients of infrared spectra of 61 batches of A. bidentata samples were 0.967 2-0.997 7; there were 13 common peaks. The results of normal distribution analysis showed that the normal distribution curve of relative peak height of common peaks for A. bidentata from Henan and Hebei did not cross， and the normal distribution curve of A. bidentata from Henan and Inner Mongolia crossed. The results of cluster analysis showed that when the distance between groups was 15， 61 batches of A. bidentata samples could be clustered into 3 categories， including N1-N12 were clustered into one
category， N13-N45 were clustered into one category， and N46-N61 were clustered into one category. The results of principal component analysis showed that the cumulative variance contribution rate of the first three principal components was 91.121%; comprehensive score of A. bidentata （number N40） in Jiabu village， Jiaozuo City， Henan Province was the highest
（2.39）， and that of A. bidentata （number N4）in Xin’an village， Anguo City， Hebei Province was the lowest （-2.89）. The results of orthogonal partial least squares-discriminant analysis showed that 61 batches of A. bidentata samples were divided into three categories，including N1-N12 were clustered into one category， N13-N28 were clustered into one category and N29-N61 were clustered into one category. Seven marker wave numbers affecting the quality were selected. The corresponding wave numbers of VIP from large to small were 1 059， 927， 2 933，813， 1 732， 1 128 and 3 367 cm-1， 1 732 cm-1 was the characteristic obsorption peak of saponins，1 059， 1 128， 927 cm-1 were the characteristic obsorption peaks of glycosides. CONCLUSIONS Infrared fingerprint combined with normal distribution analysis， cluster analysis，principal component analysis and orthogonal partial least squares-discriminant analysis can be used to identify A. bidentata from different producing areas.
KEYWORDS Achyranthes bidentata; infrared fingerprint; normal distribution analysis; cluster analysis; principal component analysis; orthogonal partial least squares-discriminant analysis; producing areas
牛膝始载于《神农本草经》，为苋科植物牛膝Achyranthes bidentata Bl.的干燥根[1]，具有补肝肾、强筋骨、活血化瘀的功效，常用于临床治疗高血压、冠心病、心绞痛、哮喘等[2]。
牛膝主要含有皂苷类、甾酮类、多糖类等化合物[3]，具有抗生育、抗肿瘤、抗衰老、抗炎及抗骨质疏松等药理作用[4-5]。
因牛膝产于古怀庆府（今河南焦作一带），为河南四大怀药之一，故又称为怀牛膝[6]。
目前，我国已形成三大牛膝产区，即内蒙古赤峰、河北安国和河南焦作[7]。
关于牛膝的研究主要集中在β-蜕皮甾酮、人参皂苷Ro等成分的含量测定，高效液相色谱（HPLC）指纹图谱的建立以及炮制工艺的优化等方面[8-12]。
2020年版《中国药典》（一部）规定采用薄层色谱法对牛膝中β-蜕皮甾酮和人参皂苷Ro进行鉴别，采用HPLC法对牛膝中β-蜕皮甾酮进行含量测定[13]。
但由于牛膝药材的化学成分复杂，仅对单一成分进行定量分析，具有一定的局限性;加之HPLC法多成分检测的成本较高，方法建立及测试费时、操作繁琐，所需流动相（甲醇、乙腈等）具有一定的毒性，应用受限[14]。
红外光谱法具有灵敏度高、特征性强、检测快速、无损等特点[15]，同时该法具有一定的整体性、客观性和科学性，既能反映药材内部所有复杂成分峰的叠加，又能反映药材所含不同成分及其比例，故可全面体现药材的质量[16]。
随着红外光谱技术的发展，其分辨率及图谱识别能力均有了新的突破，已成为中药快速鉴别及质量评价的一种有效手段[17-19]。
美日韩等国药典均已将红外光谱技术作为药物鉴定的重要方法[20]。
化学模式识别分析可以数字化地表达光谱信息，其结合红外光谱
能更加客观地评价中药材的质量[21]。
基于此，本研究建立了上述三大产地牛膝药材的红外指纹图谱，同时结合化学模式识别进行分析，旨在为快速鉴别牛膝提供参考。
1 材料
1.1 主要仪器
本研究所用主要仪器有Spectrum 100型傅里叶变换红外分光光度计、Spectrum for Window 3.02软件（美国Pekin Elmer公司），FW-4A型粉末压片机（天津市拓普仪器有限公司），FW-100型高速万能粉碎机、101-3AB型点热恒温鼓风干燥箱（北京中兴伟业仪器有限公司），ME204E/OL型万分之一天平[梅特勒-托利多仪器（上海）有限公司]等。
1.2 主要药品与试剂
25S-牛膝甾酮对照品（批号3737，纯度≥98.0%）、25R-牛膝甾酮对照品（批号3736，纯度≥98.0%）、竹节参皂苷Ⅳa对照品（批号2665，纯度≥98.0%）、β-蜕皮甾酮对照品（批号1889，纯度≥98.0%）均购自上海诗丹德标准技术服务有限公司;溴化钾（光谱纯，批号
C12045707）购自天津市科密欧化学试剂有限公司;无水乙醇（分析纯，批号20210720）购自天津市致远化学试剂有限公司;水为纯净水。
牛膝药材共61批（编号N1～N61），其中N1～N12购自河北省安国市药材市场，N13～N28购自内蒙古自治区赤峰市药材市场，N29～N61采自河南省焦作市牛膝产地。
所有药材经河南中医药大学药学院陈随清教授鉴定，均为苋科植物牛膝A. bidentata Bl.的干燥根。
取牛膝药材，除去须根、泥沙，捆成小把，晒至干皱，将顶端切齐，晒干后，用水润透，除去残留芦头，于55 ℃烘干，粉碎，过200目筛，用自封袋密封，保存至干燥器，备用。
61批牛膝药材样品的来源信息见表1。
2 方法与结果
2.1 红外光谱的建立
2.1.1 样品片的制备取牛膝药材样品粉末2 mg，加入干燥溴化钾200 mg，研磨，混匀;取混合均匀的样品适量，置于专用压片模具中，以8 MPa压制30 s，得均匀半透明的薄片，即供试样品片。
同法制备溴化钾空白片和各对照品样品片。
2.1.2 测定条件光谱扫描范围为400～4 000 cm-1，扫描前扣除CO2和H2O的干扰，扫描次数为16次/s，扫描速度为0.2 cm/s，光谱分辨率为4 cm-1，实验室温度为20～25 ℃，相对湿度为25%～35%。
2.1.3 精密度试验取牛膝药材样品（编号N1）粉末，按“2.1.1”项下方法制备供试样品片，按“2.1.2”项下条件连续测定6次，将所得红外光谱图导入OMNIC 9.2软件。
结果显示，红外光谱图的相关系数为0.997 2～0.997 4（n=6），表明方法精密度良好。
2.1.4 重复性试验取牛膝药材样品（编号N1）粉末，共6份，按“2.1.1”项下方法制备供试样品片，按“2.1.2”项下条件测定，将所得红外光谱图导入OMNIC 9.2软件。
结果显示，红外光谱图的相关系数为0.995 5～0.998 1（n=6），表明方法重复性良好。
2.1.5 稳定性试验取牛膝药材样品（编号N1）粉末，按“2.1.1”项下方法制备供试样品片，分别于室温下放置0、30、60、90、120、150 min时按“2.1.2”项下条件测定，将所得红外光谱图导入OMNIC 9.2软件。
结果显示，红外光谱图的相关系数为0.991 3～0.997 4（n=6），表明样品在室温下放置150 min内稳定性良好。
2.1.6 红外光谱数据分析取61批牛膝药材样品粉末，按“2.1.1”项下方法制备供试样品片，按“2.1.2”项下条件测定，采用Spectrum for Window
3.02软件采集上述药材样品的红外光谱图，得到61批牛膝药材的红外指纹图谱，采用OMNIC 9.2软件计算红外光谱图的相关系数。
结果显示，61批牛膝药材样品红外光谱图的相关系数为0.967 2～0.997 7;所得紅外光谱图的峰形、峰位、峰高基本相似，表明不同产地牛膝药材中含有相似的化学成分，所含特征峰及其形状基本一致，但指纹区（400～1 350 cm-1）存在一定差异。
结果见图1。
KEYWORDS Achyranthes bidentata; infrared fingerprint; normal distribution analysis; cluster analysis; principal component analysis; orthogonal partial least squares-discriminant analysis; producing areas
牛膝始載于《神农本草经》，为苋科植物牛膝Achyranthes bidentata Bl.的干燥根[1]，具有补肝肾、强筋骨、活血化瘀的功效，常用于临床治疗高血压、冠心病、心绞痛、哮喘等[2]。
牛膝主要含有皂苷类、甾酮类、多糖类等化合物[3]，具有抗生育、抗肿瘤、抗衰老、抗炎及抗骨质疏松等药理作用[4-5]。
因牛膝产于古怀庆府（今河南焦作一带），为河南四大怀药之一，故又称为怀牛膝[6]。
目前，我国已形成三大牛膝产区，即内蒙古赤峰、河北安国和河南焦作[7]。
关于牛膝的研究主要集中在β-蜕皮甾酮、人参皂苷Ro等成分的含量测定，高效液相色谱（HPLC）指纹图谱的建立以及炮制工艺的优化等方面[8-12]。
2020年版《中国药典》（一部）规定采用薄层色谱法对牛膝中β-蜕皮甾酮和人参皂苷Ro进行鉴别，采用HPLC法对牛膝中β-蜕皮甾酮进行含量测定[13]。
但由于牛膝药材的化学成分复杂，仅对单一成分进行定量分析，具有一定的局限性;加之HPLC法多成分检测的成本较高，方法建立及测试费时、操作繁琐，所需流动相（甲醇、乙腈等）具有一定的毒性，应用受限[14]。
红外光谱法具有灵敏度高、特征性强、检测快速、无损等特点[15]，同时该法具有一定的整体性、客观性和科学性，既能反映药材内部所有复杂成分峰的叠加，又能反映药材所含不同成分及其比例，故可全面体
现药材的质量[16]。
随着红外光谱技术的发展，其分辨率及图谱识别能力均有了新的突破，已成为中药快速鉴别及质量评价的一种有效手段[17-19]。
美日韩等国药典均已将红外光谱技术作为药物鉴定的重要方法[20]。
化学模式识别分析可以数字化地表达光谱信息，其结合红外光谱能更加客观地评价中药材的质量[21]。
基于此，本研究建立了上述三大产地牛膝药材的红外指纹图谱，同时结合化学模式识别进行分析，旨在为快速鉴别牛膝提供参考。
1 材料
1.1 主要仪器
本研究所用主要仪器有Spectrum 100型傅里叶变换红外分光光度计、Spectrum for Window 3.02软件（美国Pekin Elmer公司），FW-4A型粉末压片机（天津市拓普仪器有限公司），FW-100型高速万能粉碎机、101-3AB型点热恒温鼓风干燥箱（北京中兴伟业仪器有限公司），ME204E/OL型万分之一天平[梅特勒-托利多仪器（上海）有限公司]等。
1.2 主要药品与试剂
25S-牛膝甾酮对照品（批号3737，纯度≥98.0%）、25R-牛膝甾酮对照品（批号3736，纯度≥98.0%）、竹节参皂苷Ⅳa对照品（批号2665，纯度≥98.0%）、β-蜕皮甾酮对照品（批号1889，纯度≥98.0%）均购自上海诗丹德标准技术服务有限公司;溴化钾（光谱纯，批号
C12045707）购自天津市科密欧化学试剂有限公司;无水乙醇（分析纯，批号20210720）购自天津市致远化学试剂有限公司;水为纯净水。
牛膝药材共61批（编号N1～N61），其中N1～N12购自河北省安国市药材市场，N13～N28购自内蒙古自治区赤峰市药材市场，N29～N61采自河南省焦作市牛膝产地。
所有药材经河南中医药大学药学院陈随清教授鉴定，均为苋科植物牛膝A. bidentata Bl.的干燥根。
取牛膝药材，除去须根、泥沙，捆成小把，晒至干皱，将顶端切齐，晒干后，用水润透，除去残留芦头，于55 ℃烘干，粉碎，过200目筛，用自封袋密封，保存至干燥器，备用。
61批牛膝药材样品的来源信息见表1。
2 方法与结果
2.1 红外光谱的建立
2.1.1 样品片的制备取牛膝药材样品粉末2 mg，加入干燥溴化钾200 mg，研磨，混匀;取混合均匀的样品适量，置于专用压片模具中，以8 MPa压制30 s，得均匀半透明的薄片，即供试样品片。
同法制备溴化钾空白片和各对照品样品片。
2.1.2 测定条件光谱扫描范围为400～4 000 cm-1，扫描前扣除CO2和H2O的干扰，扫描次数为16次/s，扫描速度为0.2 cm/s，光谱分辨率为4 cm-1，实验室温度为20～25 ℃，相对湿度为25%～35%。
2.1.3 精密度试验取牛膝药材样品（编号N1）粉末，按“2.1.1”项下方法制备供试样品片，按“2.1.2”项下条件连续测定6次，将所得红外光谱图导入OMNIC 9.2软件。
结果显示，红外光谱图的相关系数为0.997 2～0.997 4（n=6），表明方法精密度良好。
2.1.4 重复性试验取牛膝药材样品（编号N1）粉末，共6份，按“2.1.1”项下方法制备供试样品片，按“2.1.2”项下条件测定，将所得红外光谱图导入OMNIC 9.2软件。
结果显示，红外光谱图的相关系数为0.995 5～0.998 1（n=6），表明方法重复性良好。
2.1.5 稳定性试验取牛膝药材样品（编号N1）粉末，按“2.1.1”项下方法制备供试样品片，分别于室温下放置0、30、60、90、120、150 min时按“2.1.2”项下条件测定，将所得红外光谱图导入OMNIC 9.2软件。
结果显示，红外光谱图的相关系数为0.991 3～0.997 4（n=6），表明样品在室温下放置150 min内稳定性良好。
2.1.6 红外光谱数据分析取61批牛膝药材样品粉末，按“2.1.1”项下方法制备供试样品片，按“2.1.2”项下条件测定，采用Spectrum for Window
3.02软件采集上述药材样品的红外光谱图，得到61批牛膝药材的红外指纹图谱，采用OMNIC 9.2软件计算红外光谱图的相关系数。
结果显示，61批牛膝药材样品红外光谱图的相关系数为0.967 2～0.997 7;所得红外光谱图的峰形、峰位、峰高基本相似，表明不同产地牛膝药材中含有相似的化学成分，所含特征峰及其形状基本一致，但指纹区（400～1 350 cm-1）存在一定差异。
结果见图1。
KEYWORDS Achyranthes bidentata; infrared fingerprint; normal distribution analysis; cluster analysis; principal component analysis; orthogonal partial least squares-discriminant analysis; producing areas
牛膝始载于《神农本草经》，为苋科植物牛膝Achyranthes bidentata Bl.的干燥根[1]，具有补肝肾、强筋骨、活血化瘀的功效，常用于临床治疗高血压、冠心病、心绞痛、哮喘等[2]。
牛膝主要含有皂苷类、甾酮类、多糖类等化合物[3]，具有抗生育、抗肿瘤、抗衰老、抗炎及抗骨质疏松等药理作用[4-5]。
因牛膝产于古怀庆府（今河南焦作一带），为河南四大怀药之一，故又称为怀牛膝[6]。
目前，我国已形成三大牛膝产区，即内蒙古赤峰、河北安国和河南焦作[7]。
关于牛膝的研究主要集中在β-蜕皮甾酮、人参皂苷Ro等成分的含量测定，高效液相色谱（HPLC）指纹图谱的建立以及炮制工艺的优化等方面[8-12]。
2020年版《中国药典》（一部）规定采用薄层色谱法对牛膝中β-蜕皮甾酮和人参皂苷Ro进行鉴别，采用HPLC法对牛膝中β-蜕皮甾酮进行含量测定[13]。
但由于牛膝药材的化学成分复杂，仅对单一成分进行定量分
析，具有一定的局限性;加之HPLC法多成分检测的成本较高，方法建立及测试费时、操作繁琐，所需流动相（甲醇、乙腈等）具有一定的毒性，应用受限[14]。
红外光谱法具有灵敏度高、特征性强、检测快速、无损等特点[15]，同时该法具有一定的整体性、客观性和科学性，既能反映药材内部所有复杂成分峰的叠加，又能反映药材所含不同成分及其比例，故可全面体现药材的质量[16]。
随着红外光谱技术的发展，其分辨率及图谱识别能力均有了新的突破，已成为中药快速鉴别及质量评价的一种有效手段[17-19]。
美日韩等国药典均已将红外光谱技术作为药物鉴定的重要方法[20]。
化学模式识别分析可以数字化地表达光谱信息，其结合红外光谱能更加客观地评价中药材的质量[21]。
基于此，本研究建立了上述三大产地牛膝药材的红外指纹图谱，同时结合化学模式识别进行分析，旨在为快速鉴别牛膝提供参考。
1 材料
1.1 主要仪器
本研究所用主要仪器有Spectrum 100型傅里叶变换红外分光光度计、Spectrum for Window 3.02软件（美国Pekin Elmer公司），FW-4A型粉末压片机（天津市拓普仪器有限公司），FW-100型高速万能粉碎机、101-3AB型点热恒温鼓风干燥箱（北京中兴伟业仪器有限公司），ME204E/OL型万分之一天平[梅特勒-托利多仪器（上海）有限公司]等。
1.2 主要药品与试剂
25S-牛膝甾酮对照品（批号3737，纯度≥98.0%）、25R-牛膝甾酮对照品（批号3736，纯度≥98.0%）、竹节参皂苷Ⅳa对照品（批号2665，纯度≥98.0%）、β-蜕皮甾酮对照品（批号1889，纯度≥98.0%）均购自上海诗丹德标准技术服务有限公司;溴化钾（光谱纯，批号
C12045707）购自天津市科密欧化学试剂有限公司;无水乙醇（分析纯，批号20210720）购自天津市致远化学试剂有限公司;水为纯净水。
牛膝药材共61批（编号N1～N61），其中N1～N12购自河北省安国市药材市场，N13～N28购自内蒙古自治区赤峰市药材市场，N29～N61采自河南省焦作市牛膝产地。
所有药材经河南中医药大学药学院陈随清教授鉴定，均为苋科植物牛膝A. bidentata Bl.的干燥根。
取牛膝药材，除去须根、泥沙，捆成小把，晒至干皱，将顶端切齐，晒干后，用水润透，除去残留芦头，于55 ℃烘干，粉碎，过200目筛，用自封袋密封，保存至干燥器，备用。
61批牛膝药材样品的来源信息见表1。
2 方法与结果
2.1 红外光谱的建立
2.1.1 样品片的制备取牛膝药材样品粉末2 mg，加入干燥溴化钾200 mg，研磨，混匀;取混合均匀的样品适量，置于专用压片模具中，以8 MPa压制30 s，得均匀半透明的薄片，即供试样品片。
同法制备溴化钾空白片和各对照品样品片。
2.1.2 测定条件光谱扫描范围为400～4 000 cm-1，扫描前扣除CO2和H2O的干扰，扫描次数为16次/s，扫描速度为0.2 cm/s，光谱分辨率为4 cm-1，实验室温度为20～25 ℃，相对湿度为25%～35%。
2.1.3 精密度试验取牛膝药材样品（编号N1）粉末，按“2.1.1”项下方法制备供试样品片，按“2.1.2”项下条件连续测定6次，将所得红外光谱图导入OMNIC 9.2软件。
结果显示，红外光谱图的相关系数为0.997 2～0.997 4（n=6），表明方法精密度良好。
2.1.4 重复性试验取牛膝药材样品（编号N1）粉末，共6份，按“2.1.1”项下方法制备供试样品片，按“2.1.2”项下条件测定，将所得红外光谱图导入OMNIC 9.2软件。
结果显示，红外光谱图的相关系数为0.995 5～0.998 1（n=6），表明方法重復性良好。
2.1.5 稳定性试验取牛膝药材样品（编号N1）粉末，按“2.1.1”项下方法制备供试样品片，分别于室温下放置0、30、60、90、120、150 min时按“2.1.2”项下条件测定，将所得红外光谱图导入OMNIC 9.2软件。
结果显示，红外光谱图的相关系数为0.991 3～0.997 4（n=6），表明样品在室温下放置150 min内稳定性良好。
2.1.6 红外光谱数据分析取61批牛膝药材样品粉末，按“2.1.1”项下方法制备供试样品片，按“2.1.2”项下条件测定，采用Spectrum for Window
3.02软件采集上述药材样品的红外光谱图，得到61批牛膝药材的红外指纹图谱，采用OMNIC 9.2软件计算红外光谱图的相关系数。
结果显示，61批牛膝药材样品红外光谱图的相关系数为0.967 2～0.997 7;所得红外光谱图的峰形、峰位、峰高基本相似，表明不同产地牛膝药材中含有相似的化学成分，所含特征峰及其形状基本一致，但指纹区（400～1 350 cm-1）存在一定差异。
结果见图1。
KEYWORDS Achyranthes bidentata; infrared fingerprint; normal distribution analysis; cluster analysis; principal component analysis; orthogonal partial least squares-discriminant analysis; producing areas
牛膝始载于《神农本草经》，为苋科植物牛膝Achyranthes bidentata Bl.的干燥根[1]，具有补肝肾、强筋骨、活血化瘀的功效，常用于临床治疗高血压、冠心病、心绞痛、哮喘等[2]。
牛膝主要含有皂苷类、甾酮类、多糖类等化合物[3]，具有抗生育、抗肿瘤、抗衰老、抗炎及抗骨质疏松等药理作用[4-5]。
因牛膝产于古怀庆府（今河南焦作一带），为河南四大怀药之一，故又称为怀牛膝[6]。
目前，我国已形成三大牛膝产区，即内蒙古赤峰、河北安国和河南焦作[7]。
关于牛膝的研究主要集中在β-蜕皮甾酮、人参皂苷Ro等成分的含量测定，高效液相色谱（HPLC）指纹图谱的建立以及炮制工艺的优化等方面[8-12]。
2020年版《中国药典》（一部）规定采用薄层色谱法对牛膝中β-蜕皮甾酮和人参皂苷Ro进行鉴别，采用HPLC法對牛膝中β-蜕皮甾酮进行含量测定[13]。
但由于牛膝药材的化学成分复杂，仅对单一成分进行定量分析，具有一定的局限性;加之HPLC法多成分检测的成本较高，方法建立及测试费时、操作繁琐，所需流动相（甲醇、乙腈等）具有一定的毒性，应用受限[14]。
红外光谱法具有灵敏度高、特征性强、检测快速、无损等特点[15]，同时该法具有一定的整体性、客观性和科学性，既能反映药材内部所有复杂成分峰的叠加，又能反映药材所含不同成分及其比例，故可全面体现药材的质量[16]。
随着红外光谱技术的发展，其分辨率及图谱识别能力均有了新的突破，已成为中药快速鉴别及质量评价的一种有效手段[17-19]。
美日韩等国药典均已将红外光谱技术作为药物鉴定的重要方法[20]。
化学模式识别分析可以数字化地表达光谱信息，其结合红外光谱能更加客观地评价中药材的质量[21]。
基于此，本研究建立了上述三大产地牛膝药材的红外指纹图谱，同时结合化学模式识别进行分析，旨在为快速鉴别牛膝提供参考。
1 材料
1.1 主要仪器
本研究所用主要仪器有Spectrum 100型傅里叶变换红外分光光度计、Spectrum for Window 3.02软件（美国Pekin Elmer公司），FW-4A型粉末压片机（天津市拓普仪器有限公司），FW-100型高速万能粉碎机、101-3AB型点热恒温鼓风干燥箱（北京中兴伟业仪器有限公司），ME204E/OL型万分之一天平[梅特勒-托利多仪器（上海）有限公司]等。
1.2 主要药品与试剂
25S-牛膝甾酮对照品（批号3737，纯度≥98.0%）、25R-牛膝甾酮对照品（批号3736，纯度≥98.0%）、竹节参皂苷Ⅳa对照品（批号2665，纯度≥98.0%）、β-蜕皮甾酮对照品（批号1889，纯度≥98.0%）均购自上海诗丹德标准技术服务有限公司;溴化钾（光谱纯，批号
C12045707）购自天津市科密欧化学试剂有限公司;无水乙醇（分析纯，批号20210720）购自天津市致远化学试剂有限公司;水为纯净水。
牛膝药材共61批（编号N1～N61），其中N1～N12购自河北省安国市药材市场，N13～N28购自内蒙古自治区赤峰市药材市场，N29～N61采自河南省焦作市牛膝产地。
所有药材经河南中医药大学药学院陈随清教授鉴定，均为苋科植物牛膝A. bidentata Bl.的干燥根。
取牛膝药材，除去须根、泥沙，捆成小把，晒至干皱，将顶端切齐，晒干后，用水润透，除去残留芦头，于55 ℃烘干，粉碎，过200目筛，用自封袋密封，保存至干燥器，备用。
61批牛膝药材样品的来源信息见表1。
2 方法与结果
2.1 红外光谱的建立
2.1.1 样品片的制备取牛膝药材样品粉末2 mg，加入干燥溴化钾200 mg，研磨，混匀;取混合均匀的样品适量，置于专用压片模具中，以8 MPa压制30 s，得均匀半透明的薄片，即供试样品片。
同法制备溴化钾空白片和各对照品样品片。
2.1.2 测定条件光谱扫描范围为400～4 000 cm-1，扫描前扣除CO2和H2O的干扰，扫描次数为16次/s，扫描速度为0.2 cm/s，光谱分辨率为4 cm-1，实验室温度为20～25 ℃，相对湿度为25%～35%。
2.1.3 精密度试验取牛膝药材样品（编号N1）粉末，按“2.1.1”项下方法制备供试样品片，按“2.1.2”项下条件连续测定6次，将所得红外光谱图导入OMNIC 9.2软件。
结果显示，红外光谱图的相关系数为0.997 2～0.997 4（n=6），表明方法精密度良好。
2.1.4 重复性试验取牛膝药材样品（编号N1）粉末，共6份，按“2.1.1”项下方法制备供试样品片，按“2.1.2”项下条件测定，将所得红外光谱图导入OMNIC 9.2软件。
结果显示，红外光谱图的相关系数为0.995 5～0.998 1（n=6），表明方法重复性良好。
2.1.5 稳定性试验取牛膝药材样品（编号N1）粉末，按“2.1.1”项下方法制备供试样品片，分别于室温下放置0、30、60、90、120、150 min时按“2.1.2”项下条件测定，将所得红外光谱图导入OMNIC 9.2软件。
结果显示，红外光谱图的相关系数为0.991 3～0.997 4（n=6），表明样品在室温下放置150 min内稳定性良好。
2.1.6 红外光谱数据分析取61批牛膝药材样品粉末，按“2.1.1”项下方法制备供试样品片，按“2.1.2”项下条件测定，采用Spectrum for Window
3.02软件采集上述药材样品的红外光谱图，得到61批牛膝药材的红外指纹图谱，采用OMNIC 9.2软件计算红外光谱图的相关系数。
结果显示，61批牛膝药材样品红外光谱图的相关系数为0.967 2～0.997 7;所得红外光谱图的峰形、峰位、峰高基本相似，表明不同产地牛膝药材中含有相似的化学成分，所含特征峰及其形状基本一致，但指纹区（400～1 350 cm-1）存在一定差异。
结果见图1。