一种高血压病肝阳上亢证大鼠模型判别方法及其应用与流程
本发明属于医药技术领域，本发明涉及一种高血压病肝阳上亢证大鼠模型判别方法及其应用，该方法用于评价高血压病肝阳上亢证大鼠模型。
背景技术：
病证结合动物模型，是指在中医药理论指导下，适当结合现代医学理论与实验动物科学知识，分别或同时采用传统中医学病因复制证候动物模型和现代医学病因复制疾病动物模型，令动物模型既能体现疾病的病理特征，又能系统动态地反映此特定阶段所出现的证候特征。高血压病肝阳上亢证是高血压病常见证候类型，具有特征性的发病机制和证候表征组合，因此在研究高血压病肝阳上亢证时，选用既具有高血压病的病理特点又能体现肝阳上亢证阶段性特征的病证结合动物模型，将有助于深入探讨高血压病肝阳上亢证的特点与本质，同时也有助于对治疗高血压病肝阳上亢证相关中药方剂的作用机制进行研究。
为复制高血压病肝阳上亢证动物模型，鄢东红[鄢东红，等.湖南中医学院学报，1999，19(4)：35-38.]等人采用自发性高血压大鼠加灌附子汤的方法，一方面大鼠血压稳定升高，另一方面出现了近似人类肝阳上亢证的表现，如易激惹、饮水多、旋转时间缩短及痛阈降低等。同时应用平肝潜阳的“潜阳方”进行反证，发现上述诸证均有改善。此方法利用附子的“久用必耗伤精血，致肝肾阴虚，阴不潜阳而肝阳偏亢”的药性特点制造的肝阳上亢证，操作简单，重复性好，在动物实验中得到了广泛应用。
目前，对高血压病肝阳上亢证的判定多从以下几个途径进行：模型动物的宏观体征、行为表现，与证候相关的理化指标，以及从中药药物的反证进行推测。目前采集模型动物的四诊表征判定动物证候属性，是较为常用的方法之一。依据高血压病肝阳上亢证临床表现与造模后动物的体征、行为一一对照，进行等效转换，常以检测旋转时间的长短评价是否出现眩晕；检测痛阈值的高低变化评价是否出现头疼；观察易激惹程度的变化评价是否存在烦躁易怒；观察饮水量的多少变化评价是否出现口干的症状；检测面部温度的高低变化评价是否存在面部烘热的表现。在高血压病肝阳上亢证动物模型评价研究时，还常参照临床报道的与高血压病肝阳上亢证表现相关的生化指标进行评价，如儿茶酚胺一类可以反映体内外周交感-肾上腺髓质系统兴奋性的指标。此外，在动物模型建立过程中常选用临床上常用并且疗效确切的治疗该证型的经典复方，通过比较方药对某些具有代表性指标的影响以确定建立的动物模型是否属于该证型。因此，为证实高血压病肝阳上亢证模型的成立，常常选用临床应用上具有代表的方剂，观察其对高血压病肝阳上亢证有关的异常指标的影响，从而说明模型建立的正确性与可靠性。
但当前的评价方法存在诸多弊端。首先，高血压病肝阳上亢证是体内系统病理变化的结果，尽管部分生理生化指标能够一定程度反映证候的特点，但具有相当大的片面性，并且没有特异性检测指标。采用几个理化指标简单叠加以阐释肝阳上亢证，这种对应层面的不一致性决定了难以对证候进行客观化描述。其次，反证高血压病肝阳上亢证动物模型时，应谨慎选取确有其效而又机理明确的方剂才能以此为依据来观察方剂作用模型后的效应，判断所建立的动物模型是否具有高血压病肝阳上亢证的性质，而目前临床上选用的反证方剂多是临床上的经验方，并未深入研究其药理机制。另外，回顾以往应用“以方测证”的方法时，通常只设立了正常组、高血压病肝阳上亢证组即模型组以及测证方剂组，往往缺乏对照组，然而没有对照组就不能简单的将方剂取得的功效判定为该证候的属性。此外，在进行高血压病肝阳上亢证动物模型评价时多关注动物宏观体征与生化指标的改善，对于机体内物质的扰动却关注较少，然而对高血压病肝阳上亢证模型大鼠的内在变化进行整体层面的阐述，探明高血压病肝阳上亢证发生发展过程中所产生的特殊物质群，对模型的评价具有重要意义。
技术实现要素：
为解决现有技术上的不足，本发明提供一种高血压病肝阳上亢证大鼠模型判别方法。
本发明的目的是提供一种判别指标全面，并且包含高血压病肝阳上亢证发生发展过程中所产生的特殊物质群的判别方法。
本发明的技术方案如下：
本发明提供了一种高血压病肝阳上亢证大鼠模型判别方法，包括以下步骤：
（1）复制高血压病肝阳上亢证大鼠模型作为病证组，同时选择自发性高血压大鼠作为模型组,选择正常组作为对照；
（2）实验动物干预，根据人与动物间体表面积折算的等效剂量，给予高血压病肝阳上亢证治疗药物；
（3）检测大鼠指标，包括大鼠血压测量、大鼠痛阈测定、大鼠旋转时间测定、大鼠面温测量、大鼠饮水量测量；
（4）血清样本的收集；
（5）大鼠血清angⅱ、da、ne、e、5-ht的检测；
（6）大鼠血清代谢组学分析；
（7）选取数据，纳入下述高血压病肝阳上亢证大鼠pls回归模型进行检验，pls二维得分图中病证组和模型组分离，则建立的高血压病肝阳上亢证大鼠模型成功，否则不成功；
y=-0.030x28-0.037x23-0.032x2+0.027x12-0.025x27-0.040x24+0.038x11+0.020x56+0.020x46+0.020x48+0.020x57+0.018x58+0.020x55+0.018x52+0.018x50+0.032x15+0.017x49+0.021x14-0.024x26+0.024x59+0.029x13+0.014x47+0.016x54-0.025x29-0.043x6+0.014x51-0.030x1-0.024x36-0.031x7+0.041x17+0.017x53+0.009x45+0.008x44
其中，x1是干预前收缩压数据，x2是干预第1周收缩压数据，x6是干预前舒张压数据，x7是干预第1周舒张压数据，x11是干预前痛阈数据，x12是干预第1周痛阈数据，x13是干预第2周痛阈数据，x14是干预第3周痛阈数据，x15是干预第4周痛阈数据，x17是干预第1周旋转时间数据，x23是干预第2周面温数据，x24是干预第3周面温数据，x26是干预前饮水量数据，x27是干预第1周饮水量数据，x28是干预第2周饮水量数据，x29是干预第3周饮水量数据，x36是phenylethylamine峰强度，x44是tg(18:1(9z)/18:2(9z,12z)，/20:0)[iso6]峰强度，x45是ceramide(d18:1/12:0)峰强度，x46是testosterone峰强度，x47是androsterone峰强度，x48是dihydrotestosterone峰强度，x49是8(r)-hydroperoxylinoleicacid峰强度，x50是tetrahydrodeoxycorticosterone峰强度，x51是pe(16:0/15:0)峰强度，x52是sm(d18:0/16:1(9z))峰强度，x53是chenodeoxycholicacid峰强度，x54是pc(18:4(6z,9z,12z,15z)/，22:6(4z,7z,10z,13z,16z,19z))峰强度，x55是lactosylceramide(d18:1/18:0)峰强度，x56是glucosylceramide(d18:1/26:0)峰强度，x57是3-o-sulfogalactosylceramide，(d18:1/26:1(17z))峰强度，x58是ceramide(d18:1/25:0)峰强度，x59是9,10-dhome峰强度。
优选的，步骤（1）所述的复制高血压病肝阳上亢证大鼠模型方法为，采用每天灌服附子药液方法复制高血压病肝阳上亢证模型。
更优选的，步骤（1）所述的复制高血压病肝阳上亢证大鼠模型方法为，每天灌服附子药液20ml·kg-1的方法复制高血压病肝阳上亢证模型，连续给药6周。
优选的，步骤（2）所述的实验动物干预是，根据人与动物间体表面积折算的等效剂量，计算给予钩藤、天麻或石决明。
更优选的，每日给予钩藤2.29g·kg-1，给予天麻1.15g·kg-1，或每日给予石决明2.29g·kg-1。
优选的，给予钩藤、天麻或石决明每天定时给药，每周给药6d，连续给药4周，根据体重调整给药量。
该方法在高血压病肝阳上亢证大鼠模型判别中的应用。
该方法在高血压病肝阳上亢证中药复方的方-证相关效应物质基础和作用机理研究中的应用。
本发明的有益效果是：
首先，本方法选用生理生化指标全面，并且特异性。
其次，本方法设立高血压病肝阳上亢证组即模型组和病证组，对照观察指标。
再次，本方法不仅关注动物宏观体征与生化指标的改善，还关注高血压病肝阳上亢证发生发展过程中所产生的特殊物质群。
最后，本方法可应用于高血压病肝阳上亢证中药复方的方-证相关效应物质基础和作用机理研究。
附图说明
图1是pca二维得分图(a、b为n-m正负离子模式，c、d为n-ds正负离子模式，e、f为m-ds正负离子模式)。
图2是pls-da二维得分图(a、b为n-m正负离子模式，c、d为n-ds正负离子模式，e、f为m-ds正负离子模式)。
图3是opls-da二维得分图(a、b为n-m正负离子模式，c、d为n-ds正负离子模式，e、f为m-ds正负离子模式)。
图4是n-m、n-ds以及m-ds数据矩阵韦恩图。
图5是各组高血压潜在生物标志物的变化比较（注：与模型组比○p＜0.05，与病证组比●p＜0.05）。
图6是各组肝阳上亢证潜在生物标志物的变化比较（注：与病证组比*p＜0.05。●pc为pc(18:4(6z,9z,12z,15z)/22:6(4z,7z,10z,13z,16z,19z))）。
图7是潜在生物标志物对高血压病肝阳上亢证大鼠roc分析结果。
图8是svm分类评估图。
图9是高血压病肝阳上亢证大鼠vip直方图。
图10是pls二维得分图。
具体实施方式
实施例1判别方法的构建
一、实验材料
（一）实验动物
spf级8周龄雄性shr(自发性高血压大鼠(spontaneoushypertensiverat,shr))56只，体重186.35±8.15g，spf级8周龄雄性wky7只，体重190.15±6.37g，均购自北京维通利华实验动物技术有限公司（scxk（京）：2012-0001）。shr随机分为病证组（dsgroup）、病证钩藤组（dsgtgroup）、病证天麻组（dstmgroup）、病证石决明组（dssjmgroup）、模型组（mgroup）、模型钩藤组（mgtgroup）、模型天麻组（mtmgroup）、模型石决明组（msjmgroup），wky作为正常组（normalgroup，ngroup）。每组7只。
（二）实验药物
附子、钩藤、天麻、石决明药材按照《中华人民共和国药典》（2015年版）规格，一次性购自亳州市沪谯药业有限公司，并鉴定为正品（由山东中医药大学药学院徐凌川教授鉴定）。
①附子灌胃液的标准化制备：附子药材加入8倍量水，浸泡30min，水浴加热提取1h，过滤，残渣加6倍量水再提取1h，先用四层纱布过滤后，再抽滤，合并两次提取液，混合后浓缩为含生药0.1g·ml-1的附子灌胃液，-20℃保存备用。
②钩藤灌胃液的标准化制备：钩藤以75%乙醇8倍量回流提取2次，每次2h，合并提取液，60减压回收，得浓缩为含生药0.1g·ml-1的提取液，置于-20℃保存。
③天麻灌胃液的标准化制备：天麻浸泡0.5h，10倍量水回流提取2次，每次1h，合并提取液，60℃减压回收，得浓缩为含生药0.1g·ml-1的提取液，置于-20℃保存。
④石决明灌胃液的标准化制备：石决明浸泡0.5h，10倍量水回流提取2次，每次1h，60℃减压回收，得浓缩为含生药0.1g·ml-1的提取液，置于-20℃保存。
（三）实验试剂
angⅱ试剂盒（批号：10614r）、da试剂盒（批号：10871r）、ne试剂盒（批号：10953r）、e试剂盒（批号：10105r）、5-ht试剂盒（批号：10801r），购自美国r&dsvstems公司；乙腈、甲醇（色谱级），购自美国merck公司；甲酸（色谱纯），购自fisher公司；超纯水，购于中国屈臣氏有限公司。
（四）实验仪器
大小鼠无创血压仪bp-98a（北京软隆技术有限公司），zxc-a型大鼠压痛仪（山东省医学科学院设备站），手持式红外测温仪（gm550），bb5060uv型co2培养箱（德国heraeus公司），elx50型洗板机（美国biotek公司），multiskango全波长酶标仪（美国thermoscientific公司），ultimate3000超高效液相色谱仪（美国thermo公司），qexactive™hybridquadrupole-orbitrapmassspectrometer（美国thermoscientific公司），微量移液器（芬兰dragon公司），高速台式离心机（美国thermoscientific公司），vortex-genie2涡旋震荡器（美国si公司）。
二、实验方法
（一）高血压病肝阳上亢证大鼠模型的复制方法
大鼠适应性饲养1周后，参照文献方法[鄢东红，等.湖南中医学院学报，1999，19(4)：35-38；熊新贵.长沙：湘雅医院，2009.]，dsgroup、dsgtgroup、dstmgroup、dssjmgroup大鼠按照每天灌服附子药液20ml·kg-1的方法复制高血压病肝阳上亢证模型，连续给药6周，同时正常组与模型组大鼠给予等量的蒸馏水。
（二）实验动物干预方法
高血压病肝阳上亢证大鼠模型复制后，根据人与动物间体表面积折算的等效剂量，计算各组给药量。mgtgroup与dsgtgroup每日给予钩藤2.29g·kg-1，mtmgroup与dstmgroup每日给予天麻1.15g·kg-1，msjmgroup与dssjmgroup每日给予石决明2.29g·kg-1，其余组给予等量蒸馏水。每天定时给药，每周给药6d，连续给药4周，根据体重调整给药量。
（三）大鼠指标检测
①大鼠血压测量：采用无创套尾法检测大鼠安静状态下尾动脉血压。每只大鼠连续测量3次，取平均值作为应测血压值，每周检测一次。
②大鼠痛阈测定：采用大鼠压痛测量仪测定安静状态下大鼠疼痛阈值。每周进行一次测定。
③大鼠旋转时间测定：采用自制旋转测量仪器测定大鼠旋转时间。每周进行一次测定。
④大鼠面温测量：采用红外测温仪，测量安静状态下大鼠面部温度。每周进行一次测量。
⑤大鼠饮水量测量：于早9点每笼供饮用水300ml，24h后测量剩余饮用水量，差值加灌胃量即为当日每笼总饮水量，除以每笼大鼠只数为当日每只大鼠饮水量。每周连续测量三天，取平均值作为当周每日饮水量。
（四）血清样本的收集
末次给药后，所有大鼠禁食不禁水12h，大鼠腹腔注射10%水合氯醛0.3ml·100g-1麻醉，下腔静脉取血，注入不抗凝试管中，静置30min，3500r·min-1离心10min，取上清液。将血清分为两部分，一部分用于elisa试剂盒检测，另一部分用于代谢组学分析，分别置于-80℃冰箱保存。
（五）大鼠血清angⅱ、da、ne、e、5-ht的检测
采用酶联免疫吸附法对大鼠血清angⅱ、da、ne、e、5-ht的含量进行检测，严格按照试剂盒说明书进行操作。
（六）大鼠血清代谢组学分析
a、样本制备
取备用血清室温解冻，摇匀，精密吸取血清100μｌ，加入预冷甲醇300μｌ，涡旋1min后以12000r·min-1离心10min，吸取上清液，经0.22μm微孔滤膜过滤于液相小瓶中，-20℃保存备用
b、色谱条件
液相色谱为ultimate3000超高效液相色谱仪（uplc，thermoscientific），色谱柱规格为thalo-c18柱（2.1mm×100mm，2.7μm）。流速为0.30ml·min-1，柱温保持在45℃，进样器温度为15℃，进样量5μl。流动相a为水（含0.05%甲酸），b为乙腈（含0.05%甲酸），梯度洗脱：0~1min，2%~2%b；1~3min，2%~20%b；3~4min，20%~20%b；4~7min，20%~40%b；7~9min，40%~70%b；9~15min，70%~98%b。，每次进样前以初始流动相平衡3min。
c、质谱条件
质谱为thermoscientific公司的四级杆-静电场轨道阱超高分辨质谱仪（qexactive™hybridquadrupole-orbitrapmassspectrometer），采用正负离子模式进行检测。正离子模式检测条件为：hesi离子源，鞘气45arb，辅助气10arb，裂解电压3.5kv，源内温度300℃，s-lensrflevel为55，质谱采集范围80-1000m/z，分辨率为70000。负离子模式检测条件为：hesi离子源，鞘气40arb，辅助气10arb，裂解电压2.8kv，源内温度320℃，s-lensrflevel为55，质谱采集范围80-1000m/z，分辨率为70000。在正负离子检测模式分析后，对目标代谢物进行分析，获得相关代谢物的ms/ms谱。碰撞能量分别为30ev、45ev、60ev或-30ev、-45ev、-60ev。
d、数据质量控制
同时将所有血清样本各取10μl等量混合，制备成包含所有生物信息的qc。qc包含了所有样品信息，因此最具有代表性，通过分析检测序列中插入qc的主成分分析结果，可以检测从样品前处理到样品检测过程中方法的稳定性和重复性。在大批量检测序列运行之前连续运行5次qc，使系统稳定，之后在样品检测序列中每间隔10针插入1次空白清洗和qc分析。
（七）数据处理
a、统计学处理
采用spss22.0统计软件进行分析，计量资料采用均数±标准差（
b、代谢组学数据处理
将正负离子的原始谱图数据.raw格式文件，利用proteowizard软件转换为.mzxml格式，之后采用r语言进行峰识别、峰对齐、峰校正和rt校正，逐一考察优化每个参数设置，最终获得包括质荷比（m/z）、rt及其峰面积的二维数据矩阵。，之后将对齐后的数据导入metaboanalys3.0软件进行多变量分析。采用pca、pls-da以及opls-da对数据进行分析，采用基于opls-da产生的s-plot结合vip值进行标志物的筛选，并通过vip-plot控制标志物的质量。通过foldchange分析以及t-test进一步筛选变量，保留表达差异2倍以上且p＜0.05的变量，回归原始数据核对后将剩余变量作为显著表达差异的变量，即潜在生物标志物。
由于大部分的内源性小分子代谢物的结构类型、元素组成等信息未知，因而对这些生物标志物的鉴定难度较大。本研究根据生物标志物的精确分子量在hmdb（http://www.hmdb.ca）、metlin（https://metlin.scripps.edu）、kegg（http://www.genome.jp/kegg/）等数据库中进行检索，并利用massfrontier7.0以及相关文献对ms/ms谱进行解析。
三、实验结果
（一）外观性情观察
干预期间，正常组大鼠被毛均匀细密，双目灵活，红亮有神，唇鼻色淡红湿润，笼内较安静，抓取时无明显反抗，反应较灵敏，正常饮食水，大便性状正常。模型组大鼠与病证组大鼠毛色逐渐失去光泽，甚者出现被毛发黄干燥，颈部出现竖毛，尾巴出现中度角化现象，病证组大鼠同笼互相打斗现象明显，抓取时反抗激烈攻击观察者，模型组大鼠在抓取时虽未攻击观察者但有尖叫、惊跳现象。经药物干预后的各组大鼠则随干预时间的延长，性情逐渐温顺，同笼打斗现象逐渐减少。
（二）干预前后血压变化结果
a、收缩压变化结果
各组大鼠干预前后收缩压变化见表1。不同干预时间之间大鼠收缩压存在显著差异（f=296.04，p＜0.001），各药物干预组不同时间点收缩压具有显著差异性（p＜0.001）。各组之间收缩压也存在显著性差异（f=3076.16，p＜0.001）。各组收缩压与干预时间点之间存在交互作用（f=29.57，p＜0.001）。
进步一分析单独效应，结果显示，与模型组相比，给予钩藤、天麻的shr在干预第1周即出现收缩压下降（p＜0.05），随着干预时间的延长各给药组降压效果逐渐加强，以模型钩藤组干预效果最好。与病证组相比，给予钩藤、天麻与石决明的高血压病肝阳上亢证大鼠在第2周开始呈现降压效应（p＜0.05），至干预第4周收缩压降至最低，与病证组比较均具有显著差异（p＜0.001），并以钩藤干预效果最好。与干预前相比，各给药组收缩压均下降明显（p＜0.001）。结果表明钩藤、天麻、石决明均具有降低收缩压的作用，以钩藤最佳，天麻次之。
表1各组大鼠不同时间点收缩压结果（
注：*主效应的f统计量和p值，#交互效应的f统计量和p值；与模型组比，◆p＜0.05，◇p＜0.001；与病证组比，○p＜0.001（单因素方差分析）；与干预前比较，▲p＜0.05，△p＜0.001（单因素重复测量方差分析）。
b、舒张压变化结果
各组大鼠干预前后舒张压变化见表2。不同干预时间之间大鼠舒张压存在显著差异（f=65.07，p＜0.001），各药物干预组不同时间点舒张压具有差异性（p＜0.05）。各组之间舒张压也存在显著性差异（f=380.79，p＜0.001）。各组舒张压与干预时间点之间存在交互作用（f=5.20，p＜0.001）。
进步一分析单独效应，结果显示，与模型组相比，给予钩藤、天麻与石决明的shr在干预第2周时舒张压下降（p＜0.05），随着干预时间的延长降压效果逐渐加强，至干预第4周后与模型组相比具有显著差异性（p＜0.001）。与病证组相比，给予钩藤、天麻的高血压病肝阳上亢证大鼠在干预第1周后舒张压即出现降低（p＜0.05），至干预第4周后各干预组与病证组相比舒张压显著降低（p＜0.001）。与干预前相比，随着干预时间的延长至干预第4周，各药物干预组舒张压持续下降，与干预前比较均具有显著差异性（p＜0.001）。结果表明钩藤、天麻、石决明均具有降低舒张压的作用，以钩藤最佳，石决明次之，但石决明在干预期间对舒张压的调节作用有波动。
表2各组大鼠不同时间点舒张压结果（
注：*主效应的f统计量和p值，#交互效应的f统计量和p值；与模型组比◆p＜0.05，◇p＜0.001；与病证组比，●p＜0.05，○p＜0.001（单因素方差分析）；与干预前比较，▲p＜0.05，△p＜0.001（单因素重复测量方差分析）。
（三）干预前后痛阈变化结果
各组大鼠干预前后痛阈变化见表3。不同干预时间之间大鼠痛阈存在显著差异（f=18.50，p＜0.001），其中病证钩藤组、病证天麻组与病证石决明组不同时间点痛阈具有差异性（p＜0.05）。各组之间痛阈也存在显著性差异（f=156.28，p＜0.001）。各组痛阈与干预时间点之间存在交互作用（f=3.49，p＜0.001）。
进步一分析单独效应，结果显示，在整个干预期间除干预第4周模型钩藤组痛阈值的升高较模型组具有显著差异（p＜0.05），其余时间点给予钩藤、天麻与石决明后shr的痛阈虽出现升高的趋势但无统计学意义（p＞0.05）。与病证组相比，给予钩藤、天麻与石决明3周后痛阈值出现升高（p＜0.05），至干预第4周后与病证组相比具有显著差异性（p＜0.001）。与干预前相比，shr经干预后仅以模型钩藤组在干预第4周痛阈升高较干预前比有统计学意义（p＜0.05）；而病证钩藤组、病证天麻组与病证石决明组痛阈值随着干预时间的延长持续升高，至干预第4周与干预前比较均具有显著差异性（p＜0.001）。结果表明钩藤、天麻、石决明可以提高高血压病肝阳上亢证大鼠的痛阈，钩藤与石决明效果相当。
表3各组大鼠不同时间点痛阈结果（
注：*主效应的f统计量和p值，#交互效应的f统计量和p值；与模型组比，◆p＜0.05；与病证组比，●p＜0.05，○p＜0.001（单因素方差分析）；与干预前比较，▲p＜0.05，△p＜0.001（单因素重复测量方差分析）。
（四）干预前后旋转时间变化结果
各组大鼠干预前后旋转时间变化见表4。不同干预时间之间大鼠旋转时间存在显著差异（f=10.95，p＜0.001），各组之间旋转时间的改善也存在显著性差异（f=48.56，p＜0.001）。各组旋转时间与不同干预时间点之间不存在交互作用（f=1.03，p=0.433）。
进步一分析单独效应，结果显示，与模型组相比，给予天麻、石决明干预4周后shr出现旋转时间延长（p＜0.05），其余干预时间点与模型组相比旋转时间的延长无统计学意义（p＞0.05）。与病证组相比，给予钩藤、天麻、石决明2周后高血压病肝阳上亢证大鼠的旋转时间均出现延长（p＜0.05），至干预第4周后与病证组相比旋转时间的延长具有显著差异性（p＜0.001）。与干预前相比，干预4周后模型天麻组与模型石决明组旋转时间的延长较干预前比较有统计学意义（p＜0.05），但模型钩藤组在整个干预期间的旋转时间变化与干预前均无明显差异（p＞0.05）；病证钩藤组、病证天麻组与病证石决明组经药物干预后旋转时间逐渐延长，至干预第4周，各组旋转时间与干预前比较均具有显著差异（p＜0.001）。结果表明钩藤、天麻、石决明对高血压病肝阳上亢证大鼠的旋转时间的干预作用优于shr，但三者之间效果无显著差异。
表4各组大鼠不同时间点旋转时间结果（
注：*主效应的f统计量和p值，#交互效应的f统计量和p值；与模型组比，◆p＜0.05，◇p＜0.001；与病证组比，●p＜0.05，○p＜0.001（单因素方差分析）；与干预前比较，▲p＜0.05，△p＜0.001（单因素重复测量方差分析）。
（五）干预前后面温变化结果
各组大鼠干预前后面温变化见表5。不同干预时间之间大鼠面温存在显著差异（f=25.29，p＜0.001），各组之间面温也存在显著性差异（f=23.91，p＜0.001）。各组痛阈与干预时间点之间存在交互作用（f=1.58，p=0.029）。
进步一分析单独效应，结果显示，与模型组相比，经钩藤、天麻、石决明干预后，除钩藤干预第3周与第4周以及石决明干预第3周出现明显改善外（p＜0.05），其余时间点药物干预后面温的变化与模型组比较未见差异性（p＞0.05）。与病证组相比，给予钩藤、天麻、石决明1周后高血压病肝阳上亢证大鼠的面温即出现降低（p＜0.05），至干预第4周后降至最低（p＜0.05）。与干预前相比，模型钩藤组、模型天麻组以及模型石决明组在干预4周后面温的下降较干预前比较具有统计学意义（p＜0.05）。病证石决明组在干预1周后即出现面温的降低（p＜0.05），病证钩藤组与病证天麻组在干预2周后出现降低（p＜0.05），至干预第4周与干预前比较均具有显著差异（p＜0.001）。结果表明钩藤、天麻、石决明可降低高血压病肝阳上亢证大鼠的面温，三者效果相当。
表5各组大鼠不同时间面温结果（
注：*主效应的f统计量和p值，#交互效应的f统计量和p值；与模型组比，◆p＜0.05；与病证组比，●p＜0.05，○p＜0.001（单因素方差分析）；与干预前比较，▲p＜0.05，△p＜0.001（单因素重复测量方差分析）。
（六）饮水量变化结果
各组大鼠干预前后饮水量变化见表6。不同干预时间之间大鼠饮水量存在显著差异（f=4.73，p=0.007），各组之间的饮水量也存在显著性差异（f=4.14，p=0.010）。各组饮水量变化与干预时间点之间不存在交互作用（f=1.22，p=0.26）。
进步一分析单独效应，结果显示，与模型组相比，经钩藤、天麻、石决明干预4周后饮水量有所减少（p＜0.05）。与病证组相比，给予天麻、石决明2周后高血压病肝阳上亢证大鼠的饮水量均有不同程度的减少（p＜0.05）。与干预前相比，模型钩藤组、模型天麻组以及模型石决明组在干预各个时间点的饮水量变化未见明显差异（p＞0.05）；而高血压病肝阳上亢证大鼠经钩藤、石决明干预4周后饮水量较干预前具有差异性（p＜0.05）。结果表明钩藤、天麻、石决明对高血压病肝阳上亢证大鼠的饮水量有改善作用，三者作用相当。
表6各组大鼠不同时间饮水量结果（
注：*主效应的f统计量和p值，#交互效应的f统计量和p值；与模型组比，◆p＜0.05；与病证组比，●p＜0.05（单因素方差分析）；与干预前比较，▲p＜0.05，△p＜0.001（单因素重复测量方差分析）。
（七）酶免结果
各组大鼠干预前后酶免变化见表7。与模型组相比，钩藤干预后shr血清中angⅱ含量降低（p＜0.05）；天麻干预后降低了shr血清中da的含量（p＜0.05），但钩藤、天麻、石决明对shr血清中e、ne、5-ht的调整与模型组无显著差异性（p＞0.05）。与病证组比，经钩藤与石决明干预后，除了对高血压病肝阳上亢证大鼠血清中da水平有下调趋势但无显著差异外（p＞0.05），对angⅱ、e、ne以及5-ht均有改善作用（p＜0.05）；天麻则可以调整高血压病肝阳上亢证大鼠血清中angⅱ、e、da、5-ht的含量（p＜0.05）。结果表明，钩藤、天麻与石决明均可以回调高血压病肝阳上亢证大鼠血清中angⅱ、e、5-ht的含量，三种药物不同程度的减少了高血压病肝阳上亢证大鼠体内ca类物质的含量，改善了交感-肾上腺髓质功能。
表7各组大鼠angⅱ、e、ne、da、5-ht结果（
注：与模型组比，◆p＜0.05，◇p＜0.001；与病证组比，●p＜0.05，○p＜0.001。
综上，经钩藤、天麻、石决明干预后，shr与高血压病肝阳上亢证大鼠的血压均出现降低，提示血压指标可为高血压与高血压病肝阳上亢证共性指标。而经钩藤、天麻、石决明对大鼠的痛阈、旋转时间、面温、饮水量以及易激惹程度等宏观指标的改善情况分析，则以高血压病肝阳上亢证大鼠的改善较为明显，故将宏观体征指标作为高血压病肝阳上亢证的指标。因钩藤、天麻、石决明干预后，shr与高血压病肝阳上亢证大鼠在angⅱ、e、ne、da与5-ht生化指标的改善状况上，一方面反映了三者降低血压的药效，另一方面又体现出其改善交感-肾上腺髓质功能的作用，因此将生化指标作为高血压病肝阳上亢证的指标。
（八）高血压病肝阳上亢证代谢标志物数据库构建
由于原始数据分组较多，不同大鼠模型与不同干预药物之间相互影响，难以筛选核心代谢标志物。因此，为分析正常大鼠与shr之间、正常大鼠与高血压病肝阳上亢证大鼠之间以及shr与高血压病肝阳上亢证大鼠之间的代谢产物的差异，区分高血压代谢标志物与高血压病肝阳上亢证代谢标志物，本研究对正负离子模式下原始的九组数据矩阵进行了拆分，分别对正负离子模式下的正常-模型（n-m）数据矩阵、正常-病证（n-ds）数据矩阵、模型-病证（m-ds）数据矩阵进行多变量统计分析以及差异代谢物的筛选，具体研究如下。
a、多变量统计分析
将正负离子模式下n-m、n-ds、m-ds数据矩阵分别导入metaboanalyst3.0中，经数据过滤、归一化（normalizationbysum）和标准化（paretoscaling）处理后，对满足正态性的数据进行分析。
首先应用pca对六个数据矩阵进行分析，结果如图1所示，无论是正离子模式还是负离子模式下，正常组与模型组明显分开，正常组与病证组明显分开，说明两种大鼠模型与正常大鼠代谢模式存在明显差异；同时模型组与病证组分离良好，提示高血压病肝阳上亢证大鼠与shr代谢轮廓不同。进一步采用有监督的pls-da对六个数据矩阵进行分析，结果如图2所示，六个数据矩阵在基于正负离子两种模式下所建立的pls-da二维得分图均实现了完全分离。vip的大小反映了变量对分类的重要程度，vip＞1的变量是潜在差异变量筛选的依据之一，因此本研究同时提取pls-da模型中vip＞1的变量。为获得更大程度的区分内源性代谢产物的差异，我们采用opls-da方法继续对各数据矩阵进行分析，如图3所示，基于opls-da模式清晰的反应了正常大鼠与shr、正常大鼠与高血压病肝阳上亢证大鼠以及shr与高血压病肝阳上亢证大鼠的不同代谢轮廓。所建模式具有较好的解释能力与预测能力，模式识别可靠有效。
b、差异代谢物筛选
通过上述对正常组、模型组、病证组的正确分类，表明大鼠血清中确实存在可以区分高血压和高血压病肝阳上亢证的潜在标志物，然而要确定潜在标志物需要对差异代谢进行仔细筛查提取。本研究首先使用基于opls-da产生的s-plot进行差异变量的筛选，具有较高的p和p(corr)值的变量是对opls-da模型贡献最大的变量。之后对s-plot筛选出的变量的vip值进行筛查，由于vip>1的变量是对结合模型有贡献的变量，因此排除vip<1的变量。同时通过vip-plot排除没有可靠置信区间的变量。对剩余差异变量采用foldchange与t-test相结合的方法，筛选表达相差2倍以上并且p＜0.05的变量为显著差异性变量。最后回归原始数据进行验证，提取筛选出的差异代谢物的正负离子模式下色谱峰来验证结果的准确性，剔除不可靠的代谢物，最终剩余变量为存在明显表达差异的变量，作为潜在标志物进行后续分析。
c、代谢物数据集构建
通过对数据矩阵中的变量进行上述筛选，纳入的最终变量中，n-m数据矩阵差异变量代表了正常大鼠与shr之间的差异代谢物，与高血压病密切相关；n-ds数据矩阵差异变量代表了正常大鼠与高血压病肝阳上亢证大鼠之间的差异代谢物，与高血压病肝阳上亢证密切相关；而m-ds数据矩阵差异变量代表了shr与高血压病肝阳上亢证大鼠之间的差异代谢物，与肝阳上亢证密切相关。但各个差异代谢物之间存在交叉重叠（见图4），因此需要对差异代谢物进行合并同类型。提取三个数据矩阵差异变量之间交叉的46个变量作为高血压病肝阳上亢证代谢物数据集；对n-ds与m-ds数据矩阵中剩余的差异变量进行剔重合并，最终剩余53个差异变量作为肝阳上亢证代谢物数据集。
由前期对药物干预后的大鼠代谢特征进行分析可知，shr与高血压病肝阳上亢证大鼠经钩藤、天麻、石决明干预后，体内的代谢轮廓均发生不同程度的变化，表明药物对shr和高血压病肝阳上亢证大鼠均有干预作用，在降血压的同时调整肝阳上亢证的病理状态。因此，我们通过对药物干预后的高血压病肝阳上亢证代谢数据集和肝阳上亢证代谢物数据集进行比对，筛查出经三种药物干预后变化一致并趋于正常的代谢标志物，将高血压病肝阳上亢证代谢数据集中钩藤、天麻、石决明共同干预调整的标志物作为降血压潜在生物标志物数据集；将肝阳上亢证代谢数据集中钩藤、天麻、石决明共同干预调整的标志物作为平肝潜阳潜在生物标志物数据集。
（九）潜在生物标志物鉴定
本研究中我们对潜在代谢标记物标志物进行二级质谱（ms/ms）实验，一部分代谢标志物是利用massfrontier7.0对获得的二级质谱图进行解析并与代谢组学数据库中的结构进行比对来确定的，一部分代谢标志物过代谢组学数据库信息、文献报道和结构解析来确定。最终我们确定出了14个降血压潜在生物标记物（表8）与23个平肝潜阳潜在生物标记物（表9），钩藤、天麻、石决明对shr以及高血压病肝阳上亢证大鼠降血压潜在生物标志物数据集中的潜在代谢标记物干预情况见图5，钩藤、天麻、石决明对高血压病肝阳上亢证大鼠平肝潜阳潜在生物标志物数据集中的潜在代谢标记物干预情况见图6。
表8降血压潜在标志物
表9平肝潜阳潜在标志物
至此，我们得到了能体现降血压效应的潜在生物标志物与体现平肝潜阳效应的潜在生物标志物，这些标志物即为高血压病肝阳上亢证与肝阳上亢证的潜在生物标志物。由于高血压病肝阳上亢证潜在生物标志物中存在即能够体现高血压病的因子也存在能够体现肝阳上亢证的因子，故将两部分生物标志物合并为高血压病肝阳上亢证潜在生物标志物数据集，共包含37个潜在生物标记物。
为进一步对37个潜在生物标志物判别高血压病肝阳上亢证大鼠模型的能力进行评价，采用roc曲线进行评估,以曲线下面积（auc）的大小反映标志物诊断的准确度。一般而言，auc值越大，标志物临床诊断准确性越高。表10展示了37个潜在生物标志物的auc值、灵敏度和特异度。发现其中有27个潜在生物标志物auc达到0.8以上，并且在最近临界点处有较理想的灵敏度与特异度。
表10高血压病肝阳上亢证潜在标志物roc曲线分析结果
在此基础上，采用支持向量机（supportvectormachine，svm）对拟合的roc曲线进程评估，以进一步验证潜在生物标志物对高血压病肝阳上亢证大鼠模型的评价能力。由于若直接通过auc值选取标志物，则所纳入的标志物是特定数据的最佳生物标志物，从而增加了过拟合的风险，导致评价结果过于乐观，因此对标志物进行了筛选。根据文献知识对27个标志物进行预判，选取pa(16:0/16:0)、tg(18:1(9z)/18:2(9z,12z)/20:0)[iso6]、ceramide(d18:1/12:0)、testosterone、androsterone、dihydrotestosterone、8(r)-hydroperoxylinoleicacid、tetrahydrodeoxycorticosterone、pe(16:0/15:0)、sm(d18:0/16:1(9z))、chenodeoxycholicacid、pc(18:4(6z,9z,12z,15z)/22:6(4z,7z,10z,13z,16z,19z))、lactosylceramide(d18:1/18:0)、glucosylceramide(d18:1/26:0)、3-o-sulfogalactosylceramide(d18:1/26:1(17z))、ceramide(d18:1/25:0)、9,10-dhome、7a-hydroxy-cholestene-3-one、calcidiol、7α-hydroxy-3-oxo-4-cholestenoate、calcitriol共21个生物标志物进行组合。由于本次研究纳入的样本量较少，为了产生平滑的roc曲线，进行100次交叉验证，以结果的平均值产生roc曲线。结果如图7所示，当21个生物标志组合在一起时，其auc为1.00。进一步通过svm用来进行分类评估，由图8可知所选取的21个潜在生物标志物可准确将shr与高血压病肝阳上亢证大鼠区分，两者完全分离未出现交叉。
综上，通过对高血压病肝阳上亢证潜在生物标志物数据集进行roc曲线分析，最终确定27个潜在生物标志物作为评价高血压肝阳上证大鼠模型的核心标志物。
（十）高血压病肝阳上亢证大鼠模型判别系统构建
通过上述实验，我们多层次的采集了高血压病肝阳上亢证大鼠模型的宏观表征、血液生化学以及代谢组学的特征数据，并从中提取了核心判别因子，在此基础上以将所有数据导入simca-p软件构建pls模型，以建立宏观与微观相结合的高血压病肝阳上亢证大鼠模型的判别系统。
首先对变量进行标识，由于易激惹指标为计数资料，其他指标为计量资料，为避免对模型构建造成影响故未纳入，其他所有指标变量标识见表11。
表11高血压病肝阳上亢证模型指标变量标识表
采用vip＞1且具有95%置信区间的变量（见图9）所对应的回归系数建立高血压病肝阳上亢证pls回归模型（系数保留小数点后三位）：y=-0.030x28-0.037x23-0.032x2+0.027x12-0.025x27-0.040x24+0.038x11+0.020x56+0.020x46+0.020x48+0.020x57+0.018x58+0.020x55+0.018x52+0.018x50+0.032x15+0.017x49+0.021x14-0.024x26+0.024x59+0.029x13+0.014x47+0.016x54-0.025x29-0.043x6+0.014x51-0.030x1-0.024x36-0.031x7+0.041x17+0.017x53+0.009x45+0.008x44。
综上，本课题组将中医传统四诊诊疗观念与现代科学技术相结合，多层次的采集大鼠的宏观体征表现、血液生化指标，代谢组学数据，通过对数据的整合分析，得出高血压病肝阳上亢证大鼠模型的判别系统，确定出高血压病肝阳上亢证大鼠模型以血压较高同时伴有痛阈、旋转时间、饮水量与面温指标的明显变化为宏观特征，以体内神经酰胺、睾酮、雄激素、二氢睾酮、四氢脱氧皮质酮、鹅去氧胆酸、磷脂酰乙醇胺、磷脂酰胆碱、鞘磷脂等在内的27个潜在生物标志物异常扰动为微观特征，并借助pls分析方法初步构建了高血压病肝阳上亢证大鼠模型的判别系统。该高血压病肝阳上亢证大鼠模型的判别系统既具有高血压病的病理特点又能体现肝阳上亢证的特征，借助该判别系统将有助于深入探讨治疗高血压病肝阳上亢证中药复方的方-证相关效应物质基础和作用机理。
实施例2判别方法的验证
为了检验所构建的模型对未知样本的预测能力，本研究在建模时引入了7个独立样本作为测试集，分别为来自随机选取的1个病证组样本和6个药物干预组样本，从而保证了所构建模型预测能力的客观性。结果见图10，可见模型组与病证组明显分开，表明shr与高血压病肝阳上亢证大鼠的病理情况不同。pls模型r2y为99.0%，q2为97.9%，表明模型有较好的预测能力且预测准确率较高。同时可以观察到，测试集的1个病证样本准确的落在病证组区域，而6个药物干预组与病证组相互分离，单独分布，反应了所构建模型对未知样本的准确预测。
- 还没有人留言评论。精彩留言会获得点赞！