CN107153762A - 利用脂质生物标志物预测糖尿病高血压的方法 - Google Patents

Abstract

本发明公开了一种利用脂质生物标志物预测糖尿病高血压的方法，通过OPLS‑DA模型，筛选出正常人组群和糖尿病高血压组群之间的6个差异性化合物，即与糖尿病高血压相关的脂质生物标志物，并通过构建逻辑回归模型，得到根据这些脂质标志物预测糖尿病高血压的方法，快捷方便，准确度高。通过AIC值初步判断，并进行ROC曲线绘制，准确度达到94.3％。

Description

利用脂质生物标志物预测糖尿病高血压的方法

技术领域

本发明涉及生物技术领域，尤其涉及一种利用脂质生物标志物预测糖尿病高血压的方法，能够准确预测糖尿病高血压。

背景技术

高血压既是糖尿病的风险因素，与其他因素协同促进了糖尿病及其他并发症的发生，也可以是糖尿病的并发症状之一。也可以说糖尿病与高血压经常是相互伴随的，三分之二的糖尿病患者患有高血压，而高血压患者也会出现不同程度的胰岛素抵抗。目前还没有对于高血压作为糖尿证并发症的特异性标志物，而对其进行预测与检测。但是高血压本身的早期标志物主要集中于内皮细胞功能障碍与氧化应激方面，如一氧化氮与前列环素的下降，氧化的低密度脂蛋白与脂质过氧化物的增加，因此，开发快速高效，且可以通过验血就能了解和预测糖尿病并发症高血压的方法和模型，具有重要的临床意义和市场前景。

发明内容

本发明的一个目的是解决至少上述问题，并提供至少后面将说明的优点。

本发明还有一个目的是提供一种利用脂质生物标志物预测糖尿病高血压的方法，能够快捷、准确、高效地预测糖尿病高血压。

本发明还有一个目的是提供一种针对中国人群的利用脂质生物标志物预测糖尿病高血压的方法。

本发明还有一个目的是提供一组针对中国人群的与糖尿病高血压密切相关的化合物。

为了实现根据本发明的这些目的和其它优点，提供了以下技术方案：

一种利用脂质生物标志物预测糖尿病高血压的方法，其中，主要包括以下步骤：

步骤1、筛选出正常人组群和糖尿病高血压组群之间VIP值大于1的排名前6位的差异性化合物，分别为表1所示：

表1

步骤2、利用逻辑回归模型2进行计算，得到TC值，所述逻辑回归模型2的计算公式为：

TC＝15.272+(-0.179)*R751+(-0.751)*R773+(-0.1.315)*R1400+(-1.094)*R1461；

步骤3、根据所得TC值进行判断，TC＝0为否；TC＝1为是。

优选的是，所述的利用脂质生物标志物预测糖尿病高血压的方法中，步骤1中利用OPLS-DA模型对差异性化合物进行筛选，筛选条件VIP>1。

优选的是，所述的利用脂质生物标志物预测糖尿病高血压的方法中，步骤1中筛选的具体方法为：

步骤1.1将样品进行超高效液相色谱和质谱分析，得到脂质组学数据，将正常人组群和糖尿病高血压组群分别计为CK及GXY。

步骤1.2对脂质组学数据进行标准化操作，利用OPLS-DA模型对CK及GXY进行S-plot分布得到S形曲线，并进行强制分组，计算影响CK及GXY分组的变量重要性，即得VIP值。

步骤1.3按照VIP值大于1的标准得到93个化合物，并将该93个化合物作为与糖尿病高血压相关度最高的差异性化合物。

步骤1.4将所得93个化合物按照VIP值大小从高到低排列，取前6位，即得步骤1中所述差异性化合物。

优选的是，所述的利用脂质生物标志物预测糖尿病高血压的方法中，步骤2中TC值或者用逻辑回归模型1进行计算，所述逻辑回归模型1的计算公式为：

TC＝18.661+(-0.484)*R751+(0.475)*R773+(-0.134)*R776+(-0.436)*R904+(-1.445)*R1400+(-1.207)*R1461。

优选的是，所述的利用脂质生物标志物预测糖尿病高血压的方法中，步骤2中TC值或者用逻辑回归模型3进行计算，所述逻辑回归模型3的计算公式为：

TC＝18.309+(-0.490)*R751+(-0.525)*R773+(-0.458)*R904+(-1.410)*R1400+(-1.175)*R1461。

本发明至少包括以下有益效果：

本发明通过筛选，首次发现与糖尿病高血压相关的一组化合物，即与糖尿病高血压相关的脂质生物标志物，并通过构建逻辑回归模型，得到根据这些脂质标志物预测糖尿病高血压的方法，快捷方便，准确度高。通过AIC值初步判断，并进行ROC曲线绘制，准确度高达94.3％。

本发明的其它优点、目标和特征将部分通过下面的说明体现，部分还将通过对本发明的研究和实践而为本领域的技术人员所理解。

附图说明

图1为本发明中所述的OPLS-DA模型的S-plot分布图；

图2为本发明中利用OPLS-DA模型对CK和GXY进行强制分组的结果；

图3为本发明中用来考察筛选出的化合物的火山图；

图4为本发明中ROC曲线图。

具体实施方式

下面结合附图对本发明做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。

应当理解，本文所使用的诸如“具有”、“包含”以及“包括”术语并不配出一个或多个其它元件或其组合的存在或添加。

一种利用脂质生物标志物预测糖尿病高血压的方法，其中，主要包括以下步骤：

步骤1、筛选出正常人组群和糖尿病高血压组群之间VIP值大于1的排名前6位的差异性化合物，分别为表1所示：

表1

步骤2、利用逻辑回归模型2进行计算，得到TC值，所述逻辑回归模型2的计算公式为：

TC＝15.272+(-0.179)*R751+(-0.751)*R773+(-0.1.315)*R1400+(-1.094)*R1461；

步骤3、根据所得TC值进行判断，TC＝0为否；TC＝1为是。

一个优选方案中，步骤1中利用OPLS-DA模型对差异性化合物进行筛选，筛选条件VIP>1。

一个优选方案中，步骤1中筛选的具体方法为：

步骤1.1将样品进行超高效液相色谱和质谱分析，得到脂质组学数据，将正常人组群和糖尿病高血压组群分别计为CK及GXY；

步骤1.2对脂质组学数据进行标准化操作，利用OPLS-DA模型对CK及GXY进行S-plot分布得到S形曲线，并进行强制分组，计算影响CK及GXY分组的变量重要性，即得VIP值；

步骤1.3按照VIP值大于1的标准得到93个化合物，并将该93个化合物作为与糖尿病高血压相关度最高的差异性化合物；

步骤1.4将所得93个化合物按照VIP值大小从高到低排列，取前6位，即得步骤1中所述差异性化合物。

一个优选方案中，步骤2中TC值或者用逻辑回归模型1进行计算，所述逻辑回归模型1的计算公式为：

TC＝18.661+(-0.484)*R751+(0.475)*R773+(-0.134)*R776+(-0.436)*R904+(-1.445)*R1400+(-1.207)*R1461。

一个优选方案中，步骤2中TC值或者用逻辑回归模型3进行计算，所述逻辑回归模型3的计算公式为：

TC＝18.309+(-0.490)*R751+(-0.525)*R773+(-0.458)*R904+(-1.410)*R1400+(-1.175)*R1461。

本发明通过筛选，首次发现与糖尿病高血压相关的一组化合物，即与糖尿病高血压相关的脂质生物标志物，并通过构建逻辑回归模型，得到根据这些脂质标志物预测糖尿病高血压的方法，快捷方便，准确度高。通过AIC值初步判断，并进行ROC曲线绘制，准确度高达94.3％。

实施例1

材料和方法

1.实验对象(均选自中国人)：59例糖尿病患者(无其他并发症)，111例糖尿病高血压患者，抽取静脉血5ml。准确量取100μL的血液，加入0.9mL的提取液(100％异丙醇)，转入2mL离心管中，漩涡振荡10s以上，超声10min，然后在-20度冰箱中冷冻1小时，取出后在室温下漩涡振荡，用冷冻离心机10000rpm离心10min，然后取上清液1mL到玻璃进样瓶中，保存在冰箱中待测。

2.主要仪器

2.1.1冷冻离心机：型号D3024R，Scilogex公司，美国

2.1.2漩涡振荡器：型号MX-S，Scilogex公司，美国

2.1.3高分辨质谱仪：ESI-QTOF/MS；型号：Xevo G2-S Q-TOF；厂家：Waters

2.1.4超高效液相色谱：UPLC；型号：ACQUITYUPLC I-Class系统；厂家：Waters

2.1.4数据采集软件：MassLynx4.1；厂家：Waters

2.1.5分析鉴定软件：Progenesis QI；厂家：Water

2.1.6作图软件：EZinfo；HemI；Simca-P

3.主要试剂

甲醇、乙腈、甲酸、甲酸铵、亮氨酸脑啡肽、甲酸钠。厂家均为Fisher。

4.实验设置

以提取液(100％异丙醇)作为空白对照(Blank)样品；从每个样品的进样品种取出100μL混合到新的进样品种作为质控(QC)样品；正式样品按照每组样品间隔进样，例如先是Blank1，然后是QC1，然后是W1，然后是W2这个顺序，然后再进行下一轮进样。

液相方法

色谱柱：ACQUITY UPLC CSH C18 Column,1.7μm,1mm X50mm,1/pkg[186005292]；

柱温：55度

流速：0.4mL/min

流动相：

A：ACN/H2O(60％/40％)，含有10mM甲酸铵和0.1％甲酸

B：IPA/ACN(90％/10％)，含有10mM甲酸铵和0.1％甲酸

(注：ACN为乙腈，IPA为异丙醇)

进样体积：0.2μL

洗脱程序：

质谱方法

数据采集方式：MSe；分子量扫描范围：50-1500m/z；分辨率模式(轮廓图)。

正负离子模式各采集一次。

离子源：电喷雾电离源(ESI)

毛细管电压：3KV

锥孔电压为：25V

碰撞能：15-60V

源温度：120度

脱溶剂温度：500度

锥孔气体速度：50L/h

脱溶剂气体速度：500L/h

扫描时间：0.2s

使用亮氨酸脑啡肽(m/z 556.2771，正离子；554.2615，负离子)进行实时校正。使用甲酸钠进行校正。

脂质组学数据分析

Progenesis QI软件(Waters,Massachusetts,USA)用于结果分析，抽提非靶向脂质分子的特征峰，进行比对和筛选。同时，用QC(quantifying control)和Blank(空白)来筛选背景数据。最终的数据，导入EZinfo 3.0，并进行Principal Component Analysis(PCA)分析,Orthogonal signal correction Partial Least Square DiscriminationAnalysis(OPLS-DA)建模,Variable Importance in Projection(VIP)的计算，同时得到S-Plot和火山图(coefficients vs.VIP spots)。其中，逻辑回归模型以及ROC曲线通过R语言进行建设和绘制。

结果描述

糖尿病(无其他并发症)与糖尿病高血压人群脂质差异物质的鉴定。

我们首先建立OPLS-DA模型，对糖尿病(无其他并发症)(CK)和糖尿病高血压(GXY)组群进行分类，并研究造成他们出现差异的原因。我们看到在OPLS-DA模型中，利用相关性(correlation)和协方差(covarience)的p值作出的S-plot形成了S形曲线，OPLS-DA模型的S-plot分布如图1所示，更多的差异脂质显示有极大的升高。

利用OPLS-DA将两组数据进行强制分组，结果如图2所示。

计算影响CK和GXY分组的变量重要性，即VIP(Variable Importance inProjection)值。一共筛选到93个化合物，他们的VIP值大于1；我们挑选VIP贡献率大的6个化合物。

我们通过VIP>1筛选出来的化合物在S-plot上用红色方框标出，发现他们均匀的分布在两侧。同时，我们利用火山图，来考察筛选出的化合物的分布，如图3所示，发现筛选出来的化合物都分布在火山图的外围。这些结果都说明，借由OPLS-DA模型，我们成功筛选出造成CK和GXY差异的化合物。

建立逻辑回归模型以及ROC曲线。

我们对脂质组学的数据进行标准化操作(Rproject：scale)。然后将VIP>1的排名前6的数据提取出来，建立逻辑回归模型和ROC曲线。

公式中，TC：是否患有高血压，0为否，1为是。

模型1：TC＝18.661+(-0.484)*R751+(0.475)*R773+(-0.134)*R776+(-0.436)*R904+(-1.445)*R1400+(-1.207)*R1461。经过ROC曲线绘制，模型1的AUC值为0.942，准确度为94.2％。

模型3：TC＝18.309+(-0.490)*R751+(-0.525)*R773+(-0.458)*R904+(-1.410)*R1400+(-1.175)*R1461。经过ROC曲线绘制，模型1的AUC值为0.943，准确度为94.3％。

经过ROC曲线绘制，模型2的AUC值为0.923，准确度为92.3％。

模型2：TC＝15.272+(-0.179)*R751+(-0.751)*R773+(-0.1.315)*R1400+(-1.094)*R1461。

经过ROC曲线绘制，模型3的AUC值为0.943，准确度为94.3％。

针对以上三个逻辑回归模型，进行ROC曲线绘制，如图4所示，模型1为M1，AUC值达到0.942；模型2为M2，AUC值为0.943；模型3为M3，AUC值达到0.923。我们发现模型2最靠近左上角定点，同时AUC值最高，最终，确定模型2为基于脂质指标来预测糖尿病血脂异常较好的预测模型。也可以根据所测样品数据情况，选择模型1或者模型3进行计算预测，或者优先选择模型2进行计算，同时使用1和3模型进行辅助验证，协同分析预测。

本发明利用OPLS-DA模型，首次发现与糖尿病高血压相关的6个biomarker。利用逻辑回归模型以及ROC曲线，首次提出了基于脂质生物标志物的糖尿病高血压预测模型。

尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的领域，对于熟悉本领域的人员而言，可容易地实现另外的修改，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的图例。

Claims (5)

1.一种利用脂质生物标志物预测糖尿病高血压的方法，其特征在于，主要包括以下步骤：

步骤1、筛选出正常人组群和糖尿病高血压组群之间VIP值大于1的排名前6位的差异性化合物，分别为表1所示：

表1

步骤2、利用逻辑回归模型2进行计算，得到TC值，所述逻辑回归模型2的计算公式为：

TC＝15.272+(-0.179)*R751+(-0.751)*R773+(-0.1.315)*R1400+(-1.094)*R1461；

步骤3、根据所得TC值进行判断，TC＝0为否；TC＝1为是。

2.如权利要求1所述的利用脂质生物标志物预测糖尿病高血压的方法，其特征在于，步骤1中利用OPLS-DA模型对差异性化合物进行筛选，筛选条件VIP>1。

3.如权利要求2所述的利用脂质生物标志物预测糖尿病高血压的方法，其特征在于，步骤1中筛选的具体方法为：

步骤1.1将样品进行超高效液相色谱和质谱分析，得到脂质组学数据，将正常人组群和糖尿病高血压组群分别计为CK及GXY；

步骤1.2对脂质组学数据进行标准化操作，利用OPLS-DA模型对CK及GXY进行S-plot分布得到S形曲线，并进行强制分组，计算影响CK及GXY分组的变量重要性，即得VIP值；

步骤1.3按照VIP值大于1的标准得到93个化合物，并将该93个化合物作为与糖尿病高血压相关度最高的差异性化合物；

步骤1.4将所得93个化合物按照VIP值大小从高到低排列，取前6位，即得步骤1中所述差异性化合物。

4.如权利要求1所述的利用脂质生物标志物预测糖尿病高血压的方法，其特征在于，步骤2中TC值或者用逻辑回归模型1进行计算，所述逻辑回归模型1的计算公式为：

TC＝18.661+(-0.484)*R751+(0.475)*R773+(-0.134)*R776+(-0.436)*R904+(-1.445)*R1400+(-1.207)*R1461。

5.如权利要求1所述的利用脂质生物标志物预测糖尿病高血压的方法，其特征在于，步骤2中TC值或者用逻辑回归模型3进行计算，所述逻辑回归模型3的计算公式为：

TC＝18.309+(-0.490)*R751+(-0.525)*R773+(-0.458)*R904+(-1.410)*R1400+(-1.175)*R1461。