本发明涉及一种肝损伤诊断
技术领域：
：，涉及包含生物小分子及基因的肝损伤生物标志物、方法及应用。
背景技术：
：：国内的流行病学调查显示，中草药制剂是引起药物性肝损伤排名前三位的常见药物之一。雷公藤多苷是目前在临床上用于治疗类风湿性关节炎、银屑病等免疫性疾病的重要中药制剂，其临床疗效显著，但也是发生毒性作用较多的中药制剂之一，以急性肝损伤最为严重。目前临床诊断指标与缺陷：目前临床尚无统一、公认的药物性肝损伤的诊断标准，仅根据丙氨酸氨基转移酶(alt)、天门冬氨酸氨基转移酶(ast)、碱性磷酸酶(alp)等传统指标进行判断，特异性较差，很难完全反映肝损伤的严重程度，并且无法鉴别肝损伤类型与诱导因素。另外，病理组织学检查具有侵袭性，也无法在肝损伤发生之前进行及时准确的预警。传统指标的滞后性使其不适合作为早期预警与疗效评价的标准，不能做到有效的预防与预后判断。技术实现要素：本发明的一个目的是解决至少上述问题和/或缺陷，并提供至少后面将说明的优点。本发明还有一个目的是提供用于早期发现与预警药物性肝损伤的生物标志物。本发明再有一个目的是提供用于早期发现与预警药物性肝损伤的方法。本发明另有一个目的是提供生物标志物在科学研究、制备药物性肝损伤诊断试剂盒或诊断设备中的应用。为此，本发明提供的技术方案为：用于早期发现与预警药物性肝损伤的生物标志物，所述生物标志物为生物小分子或者差异基因，所述生物小分子包括：乙烯乙酰甘氨酸、2-甲基-3戊酮酸、3-吲哚丁酸、lpc(20:2)和lpc(22:6)。优选的是，所述的用于早期发现与预警药物性肝损伤的生物标志物中，所述差异基因包括：乙醇脱氢酶adhs基因和多巴脱羧酶ddc基因nm_001270852.1/nm_001270853.1/nm_012545.4，所述乙醇脱氢酶adhs基因包括adh4nm_017270.2，adh6nm_001012084.1和adh7nm_134329.1。优选的是，所述的用于早期发现与预警药物性肝损伤的生物标志物中，所述药物性肝损伤为雷公藤多苷导致的肝损伤。用于早期发现与预警药物性肝损伤的方法，检测样品中生物标志物的浓度，所述生物标志物为生物小分子或差异基因，所述生物小分子包括乙烯乙酰甘氨酸、2-甲基-3戊酮酸、3-吲哚丁酸、lpc(20:2)和lpc(22:6)，当检测生物小分子的浓度时，利用液相色谱串联质谱法进行定量检测，若乙烯乙酰甘氨酸比正常浓度升高4倍以上、2-甲基-3戊酮酸比正常浓度升高1.5倍以上、3-吲哚丁酸比正常浓度下降0.65倍以下、lpc(20:2)比正常浓度下降0.7倍以下、lpc(22:6)比正常浓度下降0.7倍以下，则判定该样品发生药物性肝损伤。优选的是，所述的方法中，当样品中的乙烯乙酰甘氨酸比正常浓度升高4.07倍以上、2-甲基-3戊酮酸比正常浓度升高1.55倍以上、3-吲哚丁酸比正常浓度下降0.67倍以下、lpc(20:2)比正常浓度下降0.71倍以下、lpc(22:6)比正常浓度下降0.73倍以下，则判定该样品发生药物性肝损伤。优选的是，所述的方法中，所述差异基因包括乙醇脱氢酶adhs基因和多巴脱羧酶ddc基因nm_001270852.1/nm_001270853.1/nm_012545.4，所述乙醇脱氢酶adhs基因包括adh4nm_017270.2，adh6nm_001012084.1和adh7nm_134329.1，当检测样品中差异基因的表达量时，若adh4基因的表达量为正常值的4.2倍以上，ddc基因的表达量为正常值的1.2倍以上，adh6基因的表达量为正常值的47倍以上，adh7基因的表达量为正常值的0.8倍以上，则判定该样品发生药物性肝损伤。优选的是，所述的方法中，所述药物性肝损伤为雷公藤多苷导致的肝损伤。优选的是，所述的方法中，当检测生物小分子的浓度时，所述样品采用动物血清。优选的是，所述的方法中，当检测样品差异基因的表达量时，所述样品为动物肝脏组织。所述的生物标志物在科学研究、制备药物性肝损伤诊断试剂盒或诊断设备中的应用。优选的是，所述的应用中，药物性肝损伤诊断试剂盒或诊断设备应用于雷公藤多苷导致的肝损伤的诊断中。本发明至少包括以下有益效果：本发明基于代谢组学与lncrna基因芯片技术的学科交叉特色，获得的代谢物类生物标志物与基因类生物标志物具有如下优势：1、通过lncrna→mrna→代谢终产物的代谢通路对生物标志物进行了有效验证。2、生物标志物的变化发生在药物性肝损伤之前，相比于传统指标的滞后性，具有预警的作用。3、有利于降低药物性肝损伤的发病率和病死率，为临床药物性肝损伤高危人群筛查、预防干预和个体化防治及精准用药提供重要依据，将有助于改善我国药物性肝损伤防控现状，有效缓解患者机体及经济的双重负担。本发明的其它优点、目标和特征将部分通过下面的说明体现，部分还将通过对本发明的研究和实践而为本领域的技术人员所理解。附图说明图1a、1b、1c、1d、1e和1f依次分别为雷公藤多苷大鼠肝损伤模型造模六周过程中，每周对照组大鼠与相应低剂、高剂组模型大鼠肝酶指标谷草转氨酶(ast)、谷丙转氨酶(alt)、总胆红素(tbil)、γ-谷胺酰转肽酶(ggt)、碱性磷酸酶(alp)和甘油三酯(tg)的变化；图2为本发明其中一个实施例中的雷公藤多苷大鼠肝损伤模型肝组织he染色病理切片各自处理结果图，2a高剂量一周组、2b高剂量两周组、2c高剂量三周组、2d高剂量四周组、2e高剂量五周组、2f高剂量六周组、2g低剂量一周组、2h低剂量两周组、2i低剂量三周组、2j低剂量四周组、2k低剂量五周组、2l低剂量六周组，2m高剂量对照组，2n低剂量对照组；图3a和3b分别为本发明其中一个实施例中正和负离子模式下色谱总离子流图；图4a和4b分别为给药前对照组与给药第一至第五周高剂量组在正离子检测模式(4a)和负离子检测模式下(4b)的opls-da得分图(●：给药前对照组；■：给药一周组；▲：给药两周组；给药三周组；◆：给药四周组；给药五周组)；图5a为对照组与高剂量组在正离子模式下的opls-da得分图(●代表对照组，▲代表高剂量组)，5b为对照组与高剂量组在正离子模式下的pls-da模型经过100次置换验证的结果；图6a为对照组与高剂量组在负离子模式下的opls-da得分图((●代表对照组，▲代表高剂量组)，6b为对照组与高剂量组在负离子模式下的pls-da模型经过100次置换验正的结果；图7a和7b为13个药物性肝损伤潜在生物标志物的roc曲线，7a为药物性肝损伤模型中4个上调代谢物的roc曲线，7b为药物性肝损伤模型中9个下调代谢物的roc曲线；图8a和8b为药物性肝损伤模型大鼠血清中发生上调的2中代谢物的含量变化，图8a为2-甲基-3-戊酮酸的含量变化，图8b为乙烯乙酰甘氨酸的含量变化；图9a、9b、9c为药物性肝损伤模型大鼠血清中发生下调的3种代谢物的含量变化，9a为lpc(22:6)的含量变化、9b为lpc(20:2)的含量变化、9c为3-吲哚丁酸的含量变化；图10为药物性肝损伤相关代谢通路分析的可视化结果；图11为本发明另一个实施例中差异表达基因的火山图(灰色点：p>0.05的基因；绿色点：fc<2，p≤0.05的基因；红色点：fc≥2，p≤0.05的显著上调的差异基因；蓝色点：fc≤-2，p≤0.05的显著下调的差异基因)；图12为本发明另一个实施例中给药前对照组(bk)与高剂量第四周组(h4)的差异基因分层聚类图。图13为本发明另一个实施例中上调差异表达基因相关的20条生物学通路；图14为本发明另一个实施例中下调差异表达基因相关的20条生物学通路；图15为本发明中生物小分子与差异基因的示意关联图。具体实施方式下面结合附图对本发明做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。应当理解，本文所使用的诸如“具有”、“包含”以及“包括”术语并不配出一个或多个其它元件或其组合的存在或添加。基于代谢组学与lncrna基因芯片技术的学科交叉特色，获得的代谢物类生物标志物与基因类生物标志物具有如下优势：1、通过lncrna→mrna→代谢终产物的代谢通路对生物标志物进行有效验证。2、生物标志物的变化发生在药物性肝损伤之前，相比于传统指标的滞后性，具有预警的作用。3、有利于降低药物性肝损伤的发病率和病死率，为临床药物性肝损伤高危人群筛查、预防干预和个体化防治及精准用药提供重要依据，将有助于改善我国药物性肝损伤防控现状，有效缓解患者机体及经济的双重负担。本发明提供用于早期发现与预警药物性肝损伤的生物标志物，所述生物标志物为生物小分子或者差异基因，所述生物小分子包括：乙烯乙酰甘氨酸、2-甲基-3戊酮酸、3-吲哚丁酸、lpc(20:2)和lpc(22:6)。在本发明的其中一个实施例中，作为优选，所述差异基因包括：乙醇脱氢酶adhs基因和多巴脱羧酶ddc基因nm_001270852.1/nm_001270853.1/nm_012545.4，所述乙醇脱氢酶adhs基因包括adh4nm_017270.2，adh6nm_001012084.1和adh7nm_134329.1。在本发明的其中一个实施例中，作为优选，所述药物性肝损伤为雷公藤多苷导致的肝损伤。本发明提供用于早期发现与预警药物性肝损伤的方法，检测样品生物标志物的浓度，所述生物标志物为生物小分子或差异基因，所述生物小分子包括乙烯乙酰甘氨酸、2-甲基-3戊酮酸、3-吲哚丁酸、lpc(20:2)和lpc(22:6)，当检测生物小分子的浓度时，利用液相色谱串联质谱法进行定量检测，若乙烯乙酰甘氨酸比正常浓度升高4倍以上、2-甲基-3戊酮酸比正常浓度升高1.5倍以上、3-吲哚丁酸比正常浓度下降0.65倍以下、lpc(20:2)比正常浓度下降0.7倍以下、lpc(22:6)比正常浓度下降0.7倍以下，则判定该样品发生药物性肝损伤。在上述方案中，作为优选，当样品中的乙烯乙酰甘氨酸比正常浓度升高4.07倍以上、2-甲基-3戊酮酸比正常浓度升高1.55倍以上、3-吲哚丁酸比正常浓度下降0.67倍以下、lpc(20:2)比正常浓度下降0.71倍以下、lpc(22:6)比正常浓度下降0.73倍以下，则判定该样品发生药物性肝损伤。在本发明的其中一个实施例中，作为优选，所述差异基因包括乙醇脱氢酶adhs基因和多巴脱羧酶ddc基因nm_001270852.1/nm_001270853.1/nm_012545.4，所述乙醇脱氢酶adhs基因包括adh4nm_017270.2，adh6nm_001012084.1和adh7nm_134329.1，当检测样品中差异基因的表达量时，若adh4基因的表达量为正常值的4.2倍以上，ddc基因的表达量为正常值的1.2倍以上，adh6基因的表达量为正常值的47倍以上，adh7基因的表达量为正常值的0.8倍以上，则判定该样品发生药物性肝损伤。在本发明的其中一个实施例中，作为优选，所述药物性肝损伤为雷公藤多苷导致的肝损伤。在上述方案中，作为优选，当检测生物小分子的浓度时，所述样品采用动物血清。在上述方案中，作为优选，当检测样品差异基因的表达量时，所述样品为动物肝脏组织。所述的生物标志物在科学研究(比如药物性肝损伤动物模型的建立与评价；药物或化合物肝脏毒性的筛选与评价等)、制备药物性肝损伤诊断试剂盒或诊断设备中的应用。在本发明的其中一个实施例中，作为优选，所述的应用中，药物性肝损伤诊断试剂盒或诊断设备应用于雷公藤多苷导致的肝损伤的诊断中。为使本领域技术人员更好地理解本发明，现提供如下的实施例进行说明：实施例1雷公藤多苷致肝损伤大鼠模型的血清代谢组学研究2材料与方法2.1仪器与材料戴安3000系列超高效液相色谱仪(ultrahighperformanceliquidchromatography，uhplc，california，usa)，q-exactive四级杆-轨道阱高分辨串联质谱仪(thermoscientific，waltham，ma，usa)，以及其他本领域常用设备及试剂。雷公藤多苷片(批号：150302；规格：10mg*50片/瓶，福建汇天生物药业有限公司)；实验动物：6～8周龄wistar大鼠114只，体重160～180g，购自维通利华实验动物技术有限公司，许可证编号：scxk(京)2012-0001。2.2样本的收集与制备2.2.1大鼠肝损伤模型建立与样本采集雷公藤多苷大鼠肝损伤模型在18.9mg/kg的剂量下连续灌胃42天可成功建立肝损伤模型。大鼠给药分为三组：正常对照组、低剂量组(9.5mg/kg，根据成人临床最大使用剂量1.5mg/kg换算)和高剂量组(18.9mg/kg)。并根据给药时长具体分为：低剂量第一周组至第六周组(l1～l6)；高剂量第一周组至第六周组(h1～h6)；给药前空白组(bk0)及各周相对应的对照组(bk1～bk6)，每小组各6只大鼠，共114只大鼠。分别于给药第7日、第14日、第21日、第28日、第35日和第42日处死每周相应的对照及高、低剂量组大鼠。每只大鼠下腔静脉取血4ml，4℃下3500rpm离心10min。将分离的血清分为三份，第一份4℃保存进行血清生化指标检测，第二份-80℃冻存进行血清代谢组检测，剩余-80℃冻存备用。完整分离肝脏并称重，将肝组织分为三份，第一份用10％中性甲醛溶液固定以制作肝组织病理切片，第二份液氮冻存进行差异基因表达分析，剩余-80°冻存备用。2.2.2样本的制备将-80℃冻存的血清样本放置于4℃下融化，2500rpm涡旋混合3min。精密吸取血清样本40μl，并加入300μl冰乙腈，2500rpm涡旋5min，使其充分混合后，于4℃条件下10000rpm离心10min，沉淀蛋白。吸取上清液200μl，36℃离心浓缩1h至干燥。上样检测前，向经离心浓缩的样本中加入含2％乙腈的水溶液200μl进行复溶，2500rpm涡旋混合3min使其全部溶解，在4℃下10000rpm离心10min，取上清液并经96孔板过滤后进行uplc-ms/ms分析。2.3样本分析检测条件2.3.1色谱条件色谱柱：watersacquityuplchsst3(1.8μm，2.1×100mm)，与保护柱acquityuplchsst3vanguard(1.8μm)。柱温：30℃。流速：0.25ml/min。进样体积：5μl。流动相：a为含0.1％甲酸的水溶液；b为乙腈；线性梯度洗脱，进样前用起始流动相平衡色谱柱8min，具体梯度洗脱条件如下：表1-1.梯度洗脱条件流动相a：水，含0.1％甲酸；mobilephasea：0.1％formicacidinwater；流动相b：乙腈；mobilephaseb：acetonitrile。2.3.2质谱条件离子源为加热esi源，采用正、负离子检测模式检测。所使用的各种气路均为氮气。正离子检测模式质谱参数如下：负子检测模式质谱参数如下：2.4酶学指标统计学处理采用spss16.0软件对每周对照组、低剂量及高剂量组大鼠的血清生化结果进行统计分析，采用单因素方差分析(one-wayanova)进行组间比较，以p<0.05表示差异具有统计学意义。2.5多变量统计分析数据处理2.5.1数据预处理分别采用正、负离子检测模式下的lc-ms/ms方法对正常对照组大鼠及肝损伤组大鼠的血清样本进行分析。针对获得的lc-ms/ms原始数据，通过数据格式转换软件massmatrixmsdatafileconversion，将原始.raw格式的数据文件转换为.mzxml格式。再将其导入开源数据处理软件xcms进行峰识别、峰过滤和峰对齐，获得包括质荷比(m/z)、保留时间及峰面积在内的二维数据阵[21-22]。2.5.2多变量统计分析将上述二维数据阵输入多变量统计分析软件simca-p13.0(versionab，sweden)，在数据分析前进行mean-centering中值化及pareto-scaling标度化处理，消除噪音及假峰对数据分析产生的影响。采用正交偏最小二乘判别分析(orthogonalprojectionstolatentstructuresdiscriminantanalysis，opls-da)建立模式识别模型，并采用偏最小二乘判别分析(partialleastsquares-discriminantanalysis，pls-da)对所建立的模型进行交叉验证分析。2.5.3差异代谢物的筛选通过opls-da模型中变量的vip(variableimportanceonprojection)值和s-plot载荷图，挑选所有对分组有贡献的差异变量(vip>1.0)。经过带jack-knifed置信区间进行可靠性验证，并使用原始变量轮廓图删除两组数据严重重叠的变量。通过对差异变量进行组间t-test(p<0.05),保证筛选出的变量在组间的差异具有统计学意义。进一步使用pearson协相关分析计算差异变量间的二级偏相关系数；对于相关系数r>0.8的变量，结合camera分析结果以及提取离子流色谱图或离子的精确质量数和保留时间，剔除同位素离子、加和离子或碎片离子，获得准确分子离子。2.6潜在生物标志物的结构鉴定通过高分辨的ms及ms/ms谱分析、同位素丰度比和质谱裂解特征，结合网络数据库检索：hmdb(http://www.hmdb.ca/)，massbank(http://www.massbank.jp/)，metlin(http://metlin.scripps.edu/)andkegg(http://www.kegg.jp/)对发现的差异代谢物进行结构鉴定。并通过与标准品的保留时间和ms/ms谱对照分析，最终确定差异代谢物的结构。3结果与讨论3.1酶学指标分析通过测定常用临床肝功能指标谷草转氨酶(ast)、谷丙转氨酶(alt)、总胆红素(tbil)、γ-谷胺酰转肽酶(ggt)、碱性磷酸酶(alp)和甘油三酯(tg)来衡量大鼠肝损伤的进展情况。图1a-1f中为造模六周过程中，每周对照组大鼠与相应低剂、高剂组模型大鼠肝酶指标的变化。如图所示，与对照组相比，ast指标在第一周至第四周的低、高剂量组均有升高趋势；高剂量组的第一周和第四周中ast指标升高具有统计学意义；低剂量组的第二周至第四周ast指标具有明显升高；两个剂量组在第五周和第六周均无明显变化。alt指标在第一周和第二周的低、高剂量组同样显示出升高趋势，且低、高剂量组在第一周alt的升高均具有统计学意义；第三周至第六周alt变化趋于平缓，且与对照组相比有降低的趋势，这可能与动物的自身耐受有关。tbil指标在第三周的高剂量组中发生了明显的升高，alp指标在第一周的高剂量组中具有显著的增高。ggt和tg指标在第一周至第六周中未出现明显变化。血清中肝酶指标ast和alt的变化与文献报道基本相符。3.2病理组织学结果分析雷公藤多苷致大鼠肝损伤的he染色病理切片可见图2。根据文献报道，雷公藤多苷引起的肝脏损伤具有时间及剂量依赖性。由图2a至2n中可以观察到，造模六周后对照组大鼠的肝组织无明显病理改变，肝小叶及细胞结构正常、肝索排列整齐；高剂量组在第一周至第四周出现肝细胞肿胀，大部分肝细胞出现小泡性的脂肪变。第五周和第六周高剂量组出现肝索结构紊乱，肝小叶和汇管区可见炎性细胞浸润，肝细胞出现脂肪变性。低剂量组从第一周至第四周的肝细胞结构基本正常，视野中有个别细胞内出现脂肪变；第五周和第六周出现脂肪变的肝细胞较前四周稍多。综上所述，在六周造模时间内高剂量组对大鼠肝脏造成了比较明显的损害，而低剂量组的大鼠肝脏基本正常，随用药时间增加而在第五和第六周出现轻度改变。3.3数据可靠性评价在正、负离子检测模式下，对高、低剂量组肝损伤大鼠及对照组大鼠的血清样本进行lc-ms检测，总离子流色谱图(totalionchromatography，tic)如图3a和3b所示。血清样本中的诸多代谢物能够得到良好的分离，保留时间在30min以内，适用于大样本检测分析要求。3.4血清中代谢物的轮廓分析在正离子检测模式下共获得4019个变量，在负离子检测模式下获得2064个变量。将正、负离子模式下获得的二维数据阵分别导入simca-p软件进行多变量统计分析。综合病理切片及血清酶学结果，提示高剂量组大鼠在给药的六周中发生了明显的肝脏损害，选取前五周的数据进行多变量分析。为了更加清晰的观察雷公藤多苷致肝损伤大鼠在给药前及给药的五周中血清代谢物整体轮廓的变化轨迹，选择给药前对照组及高剂量组的样本同时进行opls-da模型分析。图4a和4b所示lc-(+)esi/ms谱所获得的opls-da分析结果，表明给药前对照组及第一至五周高剂量组共6个样本组可达到比较明显的聚类及分组，其中给药第四、五周组的样本比给药第一至三周组的样本更加偏离对照组，表明随着肝损伤发生时间的增加，其病理损害逐渐加重。对lc-(-)esi/ms(图4b)谱所获得的数据阵进行同样的统计分析，结果显示六组间同样能够得到比较明显的聚类及分组。通过对给药前对照组及高剂量组的大鼠血清样本进行整体代谢轮廓分析后，发现高剂量组在给药第一周后即在代谢水平上出现明显的扰动，与病理组织的观察结果相一致。因此，后期将选择高剂量组所有样本与对照组所有样本的lc-ms谱数据进行opls-da多变量统计分析，有利于寻找可用于药物性肝损伤早期发现的潜在生物标志物。3.5模式识别模型的建立及验证选用pls-da和opls-da模型进行有监督的数据分析，并采用置换验证(permutationtest)避免pls-da模型拟合过度以证明模型的可靠性。进一步结合opls-da模型中的变量vip值，带jack-knife置信区间的载荷图及原始变量轮廓图对差异代谢物进行筛选。为寻找对照组与高剂量组大鼠血清中的差异代谢物，本研究选用opls-da模型对血清样本的lc-(±)esi/ms数据进行有监督的数据分析。如图5a所示，lc-(+)esi/ms谱获得的数据得到的opls-da模型含有1个预测成分和2个正交成分，46.2％的变量[r2(x)]可用于解释74.5％的组间差异[r2(y)]，经过交叉验证显示预测能力为52.2％[q2y]。r2(y)与q2y之间的差值为0.223(一般不超过0.2～0.3)，q2y值大于50％，表示模型预测能力良好。为了避免模型过度拟合，对与opls-da模型主成分数相同的pls-da模型进行置换验证，图5b为其经过100次建模的置换验证结果，获得的r2y(绿色圆形)与真实模型的r2y值共同构成的回归线截距为0.375(理论上应小于0.3～0.4)，q2y(蓝色方块)与真实模型的q2y值共同构成的回归线截距为-0.192(理论上应小于0.05)。由lc-(-)esi/ms谱获得的数据得到的opls-da模型(图6a)含有1个预测成分和2个正交成分，61.5％的变量[r2(x)]可用于解释64.5％的组间差异[r2(y)]，经过交叉验证显示预测能力为49.4％[q2y]。r2(y)与q2y之间的差值为0.151。对主成分数为3的pls-da模型进行置换验证，图6b为其经过100次建模的置换验证结果，获得的r2y(绿色圆形)与真实模型的r2y值共同构成的回归线截距为0.388，q2y(蓝色方块)与真实模型的q2y值共同构成的回归线截距为-0.227。以上验证结果表明，血清样本在正、负离子检测模式下获得的数据结果，经过多变量统计分析后，数据模型符合相关参数标准，能够有效的进行潜在生物标志物的筛选。3.6潜在生物标志物的筛选(1)本研究采用viplist与s-plot载荷图相结合来挑选所有对分组有贡献的差异变量。首先通过s-plot载荷图初步筛选差异变量；再结合viplist，挑选vip值大于1的变量，当vip值大于1时则认为该变量对于模型分组的贡献高于平均水平[31]，进一步通过jack-knifed置信区间验证这些变量的可靠性，删除跨越零值的变量。经过上述筛选步骤，在正离子检测模式下获得252个差异变量，在负离子检测模式下获得个170差异变量。(2)对上述差异变量在对照组和高剂量组中的检测强度(峰面积)进行独立样本t检验，删除两组间差异无统计学意义(p>0.05)的变量；结合原始代谢轮廓图，删除组间数据严重交叉的变量。经过上述筛选，在正离子检测模式下获得152个差异变量，在负离子检测模式下获得129个差异变量。(3)通过手动提取原始数据图谱中的差异变量色谱峰，以验证结果的准确性，将不可靠的差异代谢物剔除，在正离子检测模式下获得84个差异变量，在负离子检测模式下获得38个差异变量。(4)进一步对差异变量进行pearson协相关分析，计算各差异变量的二级偏相关系数。对于相关系数r>0.8的变量，结合提取离子流色谱图或离子的精确质量数和保留时间，以及r语言中camera方法给出的信息，判断相同保留时间的离子是否来源于相同代谢物，剔除同一代谢物的加合离子、同位素离子及碎片离子，获得准确分子离子。经过以上筛选，最终在正离子检测模式下获得个30差异变量，在负离子检测模式下获得25个差异变量，详细信息可见于表1-2。表1-2.在正、负离子模式下对照组与肝损伤高剂量组的潜在生物标志物a.foldchangewascalculatedbytheratioofmeanvalueofcontrolgrouptohighdosegroup.3.7潜在生物标志物的结构鉴定对可能的生物标志进行结构鉴定时，首先需确定潜在生物标志物的[m+h]+或[m-h]-分子离子，通过其精确质量数结合同位素丰度比法初步推测可能生物标志物的分子组成；结合代谢物的高分辨ms、ms/ms谱，代谢物结构的质谱裂解特征和网上代谢物数据库(hmbd、metlin、kegg等)检索，推断可能生物标志物的可能结构。采用上述方法及思路对本研究发现的其他可能生物标志物进行分析，共鉴定出13个潜在生物标志物的结构，包括：吡啶甲酸、2-甲基-3-戊酮酸、l组氨醇、乙烯乙酰甘氨酸、谷氨酸、7-甲基鸟嘌呤、5-羟色胺、3-吲哚丁酸、龙胆酸及四种溶血磷脂酰胆碱lpc(18:2)、lpc(20:3)、lpc(20:2)和lpc(22:6)，具体信息可见表1-3。其余一些未确定结构的生物标志物，主要是由于其在血清中含量较低，未能获得理想的ms/ms谱图，或者是在代谢物数据库中未能获得足够的信息确认其结构。3.8潜在生物标志物的诊断敏感性和特异性评价本研究采用接收者工作特征曲线(receiveroperatingcharacteristiccurve，roc)对已鉴定出的13个潜在生物标志物的诊断性能进行评价。roc曲线下面积(areaunderthecurve，auc)可反映诊断指标区分阳性和阴性诊断能力的大小。auc值越大，则表示诊断准确性越高，通常认为当0.5<auc<0.7时，诊断价值较低；在0.7<auc<0.9时，诊断价值中等；当auc≥0.9时，诊断价值较高。采用spss16.0软件进行roc曲线分析，图7a为在肝损伤大鼠血清中发生上调的4个潜在生物标志物的roc曲线，其中乙烯乙酰甘氨酸和2-甲基-3-戊酮酸的auc分别为0.840和0.727，表示其具有一定的诊断价值；图7b为肝损伤大鼠血清中发生下调的9个潜在生物标志物的roc曲线，其中3-吲哚丁酸、lpc(20:2)和lpc(22:6)的auc分别为0.828、0.735和0.715，表示其作为药物性肝损伤潜在标志物的诊断具有一定的准确性。进一步分析上述5个潜在标志物在正常对照组及低剂量组中大鼠血清中的变化。由病理结果可知，低剂量组在第一至第四周均未发生明显病理改变，而第五周出现少量细胞脂肪变，因此将结合与肝损伤关系更敏感的alt指标和5个潜在标志物在低剂量组第一至第五周中的含量进行比较。图8a和8b所示为低剂量组大鼠血清中发生上调的两个代谢物的含量变化图。与对照组相比，乙烯乙酰甘氨酸和2-甲基-3-戊酮酸在低剂量组第一至第五周中的含量均明显升高，而alt指标未发生显著改变；表明乙烯乙酰甘氨酸和2-甲基-3-戊酮酸在肝损伤发生时比alt指标更为灵敏，能够在肝实质损伤发生前进行预警，可能为药物性肝损伤早期发现的潜在生物标志物。图9a、9b、9c所示为低剂量组大鼠模型血清中发生下调的三个代谢物的含量变化图。与对照组相比，3-吲哚丁酸、lpc(20:2)和lpc(22:6)在低剂量第四、五周中的含量均显示明显下降，alt指标未发生显著改变，提示三种代谢物相联合可更加有效的预警肝损伤，可能成为药物性肝损伤早期发现的潜在生物标志物。3.9药物性肝损伤相关的代谢通路分析为进一步分析药物性肝损伤发生发展机制可能涉及的代谢通路，利用所发现的13个潜在生物标志物进行了相关的代谢通路分析。metaboanalyst3.0(www.metaboanalyst.ca)是一款专门为代谢组学提供综合数据分析的在线网站，其中的通路分析模块(pathwayanalysismodule)以京都基因与基因组百科全书(kyotoencyclopediaofgenesandgenomes，kegg)代谢途径数据库为依托，可将富集分析(enrichimentanalysis)与通路拓扑结构分析(pathwaytopologyanalysis)的结果相整合，确定与研究结果最相关的代谢通路，并将其可视化以直观分析。将已鉴定出的13个潜在药物性肝损伤生物标志物的英文名称输入至代谢通路分析模块，并与kegg等数据库中标准化合物信息进行匹配，选择所要分析的物种rattusnorvegicus(rat)通路数据库进行检索，最终共发现13条代谢通路可能与药物性肝损伤的发生机制相关，涉及谷氨酸盐代谢，氮素代谢，类脂化合物代谢，组氨酸代谢，丁酸甲酯代谢，丙氨酸、天门冬氨酸和谷氨酸代谢，谷胱甘肽代谢，甘油磷脂代谢，色氨酸代谢，酪氨酸代谢，精氨酸和脯氨酸代谢，氨酰基-trna生物合成与卟啉和叶绿素代谢，结果可见表1-4。代谢通路的可视化结果如图10所示，当p值越小时，节点的颜色越深；代谢通路影响值(pathwayimpact，pi)越大时，节点的半径越大，其所代表的通路与药物性肝损伤的相关性越大。发现的相关13条代谢通路中，影响值pi>0.1的共有4条，即谷氨酸盐代谢，丙氨酸、天门冬氨酸和谷氨酸代谢，类脂化合物代谢，色氨酸代谢，表明这4条通路与药物性肝损伤的发生密切相关。表1-4.药物性肝损伤潜在生物标志物的代谢通路分析a.therawpwasoriginalpvaluederivedfromenrichimentanalysis.4小结本研究采用基于lc-ms/ms的代谢组学方法，对雷公藤多苷致肝损伤模型大鼠与正常对照大鼠的血清进行分析。通过多元统计分析，建立opls-da模式识别模型；通过viplist、jack-knifed置信区间等一系列筛选，从高剂量肝损伤组大鼠与正常对照组大鼠的血清中获得55个差异代谢物。通过高分辨ms谱和ms/ms谱分析，对潜在生物标志物的结构进行分析，鉴定出13个差异代谢物的结构，其中4个代谢物：2-甲基-3戊酮酸、l组氨醇、乙烯乙酰甘氨酸及谷氨酸在肝损伤大鼠血清中发生明显上调；9个代谢物：吡啶甲酸、7-甲基鸟嘌呤、5-羟色胺、3-吲哚丁酸、龙胆酸及四种溶血磷脂酰胆碱lpc(18:2)、lpc(20:3)、lpc(20:2)和lpc(22:6)发生下调。采用roc曲线对差异代谢物诊断敏感性及特异性进行分析，发现乙烯乙酰甘氨酸、2-甲基-3戊酮酸、3-吲哚丁酸、lpc(20:2)和lpc(22:6)这5个代谢物作为药物性肝损伤早期发现的潜在生物标志物具有一定的诊断价值。进一步通过已鉴定的潜在生物标志物进行了相关代谢通路分析，发现药物性肝损伤大鼠体内谷氨酸盐代谢，丙氨酸、天门冬氨酸和谷氨酸代谢，类脂化合物代谢，色氨酸代谢4条通路发生明显紊乱，表明上述代谢通路与药物性肝损伤的发生显著相关。实施例2雷公藤多苷致肝损伤大鼠模型的肝脏差异基因表达谱研究2材料与方法2.1仪器与材料均为本领域常用仪器及材料2.2样本的收集本研究所使用的样本为雷公藤多苷致肝损伤模型的大鼠肝脏组织，本实施例中采用实施例1中的大鼠肝脏组织。结合肝酶指标与病理结果，选取给药前对照组大鼠和药物性肝损伤高剂量组大鼠给药第一周、第二周和第四周的样本进行差异基因表达分析。从每周的6个样本中各随机抽取3个样本进行基因芯片检测。抽取的样本编号为：对照组bk1、bk5、bk6；高剂量第一周组h1-3、h1-5、h1-6；高剂量第二周组h2-4、h2-5、h2-6；高剂量第一周组h4-2、h4-5、h4-6。所有样本在提取及检测前均保存于液氮。2.3肝组织rna的提取与质量检测肝组织总rna利用nanodropnd2000进行定量，并经agilent2100型生物分析仪及1％琼脂糖电泳检测rna完整性和纯度。rna完整指数(rnaintegritynumber，rin)大于7且28s/18s>0.7时，表示rna完整性好；测定rna在260nm和280nm的吸光度值，od260/280值处于2.0左右表示提取的rna纯度较高，可用于基因芯片分析。rna质检合格后，样本的标记、芯片的杂交以及洗脱参照芯片标准流程。肝脏总rna的提取及基因芯片表达谱检测均由上海欧易公司完成。2.4基因芯片检测分析将肝脏总rna反转录成双链cdna，再进一步合成用生物素cyanine-3-ctp(cy3)标记的crna。本研究使用的芯片为agilentratlncrnaarray。将标记好的crna和芯片杂交，洗脱后利用agilentscannerg2505c扫描芯片获得原始探针信号。采用featureextraction软件(version10.7.1.1,agilenttechnologies)处理原始图像并提取原始数据，并利用genespring软件(version13.1,agilenttechnologies)对原始数据进行quantile标准化处理。2.5统计学处理与生物信息学分析将标准化后的数据进行统计学处理，组间比较采用独立样本t检验，以p≤0.05为差异有统计学意义。利用t检验的p值和倍数变化(foldchange，fc)值进行差异表达基因筛选，筛选的标准为上调或者下调倍数变化fc≥2.0且p<0.05。进一步对差异表达基因进行非监督层次聚类，利用热图的形式展示差异基因在不同样本间的表达模式。对差异基因进行kegg富集分析以判定差异基因主要影响的生物学通路。2.6差异表达基因mrna荧光定量pcr检测分析针对主要影响的生物学通路上的差异基因进行pcr定量检测分析，获得准确定量结果。具体步骤包括：rna提取：提取总rna，利用nanodrop2000分光光度计(thermoscientific，usa)测定浓度及od260/od280，琼脂糖凝胶电泳检测rna完整性。逆转录：利用hiscriptiiqrtsupermixforqpcr(+gdnawiper)(vazyme，r223-01)将待测rna逆转录成cdna。逆转录体系：第一步：总rna，0.5μg；4×gdnawipermix，2ul；；nuclease-freeh2o加至8μl，反应程序：42℃2min。第二步：加入5×hiscriptiiqrtsupermixiia，2μl，反应程序：25℃10min，50℃30min，85℃5min。总反应体积10ul。逆转录完毕后加入90μlnuclease-freeh2o储存在-20℃冰箱备用。引物设计：引物采用rochelcpds2软件设计并由北京擎科新业生物技术有限公司合成。荧光定量pcr：利用greenpcrkit试剂盒(qiagen，germany)在480ⅱ型荧光定量pcr仪(roche，swiss)上进行反应。体系：greenpcrmastermix，5μl；10μmforwardprimer，0.2μl；10μmreverseprimer，0.2μl；cdna，1μl；nuclease-freeh2o，3.6μl。pcr程序：95℃5min；95℃10s，60℃30s，40个循环。循环结束后利用熔解曲线检测产物特异性：从60℃缓慢升温至97℃，每℃采集5次荧光信号。计算表达量：获得扩增曲线和熔解曲线，采用2-δδct法(δct＝ct目的基因-ct内参基因，δδct＝δct实验组样品-δct对照组样品)，计算表达量与相对表达量。3结果与讨论3.1肝组织提取总rna的浓度和质量检测表2-1为样本总rna的浓度及质量检测结果。12个样本的总rna提取结果良好，总rna的吸光度od260/280值均在(2.0±0.2)之间，28s/18s>1.5且rin≥7，表示rna的完整性和纯度均符合芯片检测要求。表2-1.总rna的浓度及质量检测结果3.3差异表达基因的筛选基因芯片的分析结果利用t检验得到的差异显著性p值和标准化信号值的差异倍数值进行筛选，标准为fc≥2.0且p≤0.05。筛选后h4组与bk组中共有1803条mrna存在表达差异，其中表达上调759个，表达下调1044个。图11为差异表达基因的火山图，x轴是差异倍数以2为底取对数的值，y轴是p-value值以10为底取负对数的值。图12为给药前对照组(bk)与高剂量第四周组(h4)的差异表达基因的分层聚类图。3.4差异表达基因相关通路分析利用kegg数据库对差异基因进行相关通路(pathway)分析，并且用统计检验的方法计算每个pathway条目中差异基因富集的显著性，以p<0.05表示差异基因在该pathway中出现了显著性的富集。通过kegg数据库的富集分析，对差异表达的1803个基因进行pathway功能分析，图13和图14为上调及下调表达基因富集分析排名前20的相关生物学通路。富集结果显示，上调差异表达基因显著参与生物学通路17条，其中涉及氨基酸代谢通路3条，包括精氨酸生物合成、氨基酸生物合成及丙氨酸、天门冬氨酸和谷氨酸盐代谢；下调差异表达基因显著参与生物学通路8条，其中涉及1条氨基酸代谢通路，即酪氨酸代谢。具体通路信息见于表2-2和表2-3。其中酪氨酸代谢通路，丙氨酸、天门冬氨酸和谷氨酸盐代谢通路。表2-2.上调差异基因显著富集的17条生物学通路table2-2.20biologicalpathwaysofrelatedup-regulatedgenesa.thedescriptionofbiologicalpathwayb.thenumberoftargetgenesinthispathwayc.significantpvaluesofenrichmentinpathwayd.thepvaluesaftercalibrationbybenjamini-hochbergmethod表2-3.下调差异基因显著富集的8条生物学通路a.thedescriptionofbiologicalpathwayb.thenumberoftargetgenesinthispathwayc.significantpvaluesofenrichmentinpathwayd.thepvaluesaftercalibrationbybenjamini-hochbergmethod通过前期代谢组学的研究结果，确定了tyrosinemetabolism代谢通路，可以通过对酪氨酸代谢通路的基因进行检测以作为佐证。3.5差异表达基因mrna荧光定量pcr检测分析结合代谢组学获得的研究结果，针对酪氨酸代谢通路上差异表达的基因：乙醇脱氢酶adhs(adh6、adh4、adh7)与多巴脱羧酶ddc基因进行pcr定量检测分析，获得健康组与药物性肝损伤模型组基因的准确定量结果。表2-5为基因的具体名称与引物信息。表2-5基因的具体名称与引物信息具体表达量与相对表达量及相关倍数关系如表2所示。表2.差异mrna荧光定量pcr结果4小结本研究采用基因芯片技术结合生物信息学分析，对雷公藤多苷致肝损伤大鼠模型的肝脏组织进行差异基因表达谱初步研究。针对空白对照组大鼠和药物性肝损伤组大鼠间的差异表达基因，以fc≥2.0且p≤0.05为筛选标准，共检出1803个差异表达基因，其中759个表达上调，1044个表达下调。对差异表达基因进行基于kegg数据库的通路(pathway)功能分析。富集结果显示：上调差异表达基因主要参与生物学通路17条，其中涉及氨基酸代谢通路3条，包括精氨酸生物合成、氨基酸生物合成及丙氨酸、天门冬氨酸和谷氨酸盐代谢；下调差异表达基因主要参与生物学通路8条，其中涉及1条氨基酸代谢通路，即酪氨酸代谢。进一步针对酪氨酸代谢通路关注的基因，获得对应分组中的表达量，并进行mrna荧光定量pcr验证，验证生物标志物的可靠性提供依据。综合实施例1和2，由图15可见，实施例1中发现的生物标志物共鉴定出13个潜在生物标志物的结构，包括：吡啶甲酸、2-甲基-3-戊酮酸、l组氨醇、乙烯乙酰甘氨酸、谷氨酸、7-甲基鸟嘌呤、5-羟色胺、3-吲哚丁酸、龙胆酸及四种溶血磷脂酰胆碱lpc(18:2)、lpc(20:3)、lpc(20:2)和lpc(22:6)，涉及谷氨酸盐代谢，氮素代谢，类脂化合物代谢，组氨酸代谢，丁酸甲酯代谢，丙氨酸、天门冬氨酸和谷氨酸代谢，谷胱甘肽代谢，甘油磷脂代谢，色氨酸代谢，酪氨酸代谢，精氨酸和脯氨酸代谢，这其中大多涉及氨基酸代谢；并且，基因lncrna代谢通路分析结果显示确有氨基酸代谢通路相关，同时对酪氨酸代谢通路相关基因(adh4等)进行定量验证，确有差异；这表明上述13个潜在生物标志物与相关差异基因(adh4等)确有关联。实施例3一种高灵敏度、高通量液相色谱串联质谱法定量动物血清的生物小分子试剂盒及其制备以及基因芯片技术检测动物基因标志物试剂盒及其制备，由特定浓度的生物小分子标准曲线、质控样本、样本沉淀试剂、内标液、流动相调节液组成；以及rna的提取及基因芯片表达谱检测。包括如下步骤：(1)生物小分子(乙烯乙酰甘氨酸、2-甲基-3戊酮酸、3-吲哚丁酸、lpc(20:2)和lpc(22:6))储备液、标志物母液、标准曲线样本、质控样本的配制；(2)样本沉淀试剂的配制；(3)流动相调节液甲酸的配制。血清的生物小分子定量检测方法，包括：(1)血清样本预处理：血清于4℃条件下融化，涡旋混匀，通过乙腈沉淀蛋白，加入内标液，涡旋混匀，高速离心，取上清液。(2)液相分离：a.采用c18、c8或氰基键和硅胶作固定相；b.流动相：乙腈；甲醇；乙酸铵；甲酸铵；等度或梯度洗脱；控制流速。其所使用的流动相为：乙腈；甲醇；乙酸铵；甲酸铵；反相色谱等度或梯度分离。(3)质谱测定：电喷雾离子源；大气压化学电离源；正离子方式。利用液相色谱串联质谱法定量检测样品肝脏组织血清中生物标志物的浓度，当检测血清中生物小分子的浓度时，若血清中的乙烯乙酰甘氨酸比正常浓度升高4倍以上、2-甲基-3戊酮酸比正常浓度升高1.5倍以上、3-吲哚丁酸比正常浓度下降0.65倍以下、lpc(20:2)比正常浓度下降0.7倍以下、lpc(22:6)比正常浓度下降0.7倍以下，则判定该样品发生药物性肝损伤实施例4基因芯片技术检测动物基因标志物试剂盒及其制备，包括如下步骤：(1)rna的提取：总rna利用nanodropnd2000进行定量，并经agilent2100型生物分析仪及1％琼脂糖电泳检测rna完整性和纯度。(2)差异基因(adh6、adh4、adh7与ddc)及内参基因引物序列、pcr曲线，溶解曲线的配制；(3)测定rna在260nm和280nm的吸光度值；(4)基因芯片分析血清的差异基因检测方法，包括：(1)基因芯片检测分析：使用的芯片agilentratlncrnaarray。将标记好的crna和芯片杂交，洗脱后利用agilentscannerg2505c扫描芯片获得原始探针信号。采用featureextraction软件(version10.7.1.1,agilenttechnologies)处理原始图像并提取原始数据，并利用genespring软件(version13.1,agilenttechnologies)对原始数据进行quantile标准化处理。(2)mrna检测与数据处理：rt-pcr将以rna为模板的cdna合成同pcr结合在一起，提供了一种分析基因表达的快速灵敏的方法。当检测样品肝脏组织血清中差异基因的表达量时，若adh4基因的表达量为正常值的4.2倍以上，ddc基因的表达量为正常值的1.2倍以上，adh6基因的表达量为正常值的47倍以上，adh7基因的表达量为正常值的0.8倍以上，则判定该样品发生药物性肝损伤。实施例5所述的生物标志物在科学研究(比如药物性肝损伤动物模型的建立与评价；药物或化合物肝脏毒性的筛选与评价等)、制备药物性肝损伤诊断试剂盒或诊断设备中的应用，药物性肝损伤诊断试剂盒或诊断设备应用于雷公藤多苷导致的肝损伤的诊断中这里说明的模块数量和处理规模是用来简化本发明的说明的。对本发明的用于早期发现与预警肝损伤的生物标志物、方法及应用的应用、修改和变化对本领域的技术人员来说是显而易见的。尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的领域，对于熟悉本领域的人员而言，可容易地实现另外的修改，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的图例。当前第1页12当前第1页12