本发明属于检测分析技术领域，尤其涉及一种基于液质联用代谢组学分析玛咖标志性代谢产物的方法。
背景技术：
玛咖(lepidiummeyeniiwalp.),十字花科(cruciferae)独行菜属(lepidium)植物，原产于海拔3500～4500米的南美安第斯山区,已经有2000多年的栽培历史，具有食用和药用价值。自20世纪90年代以来，美国、日本、德国等许多国家尝试引进和培育玛咖。2002年，玛咖被引入中国。在过去的十多年中，这种植物已经在中国的许多地区，如云南、新疆、吉林和四川等，成功引进和种植。
玛咖的营养成分合理，种类繁多，主要包括：蛋白质、多种氨基酸、脂肪酸、纤维素、糖类、钙、锌、铁等矿物质及多种维生素等。玛咖含有多种生物活性成分，已经通过高效液相色谱—紫外分光光度计、lc-ms、gc-ms、核磁共振波谱法等方法来分析和鉴定，包括玛咖烯、玛咖酰胺、生物碱、甾醇和芥子油苷及其分解产物异硫氰酸苄酯。玛咖含有的生物活性物质使之具有特殊的药理功效，例如增强性功能、提高生育能力、抗疲劳、抗癌、抗氧化、促进生长、缓解良性前列腺增生、缓解女性更年期症状的发生等。同时，玛咖的下胚轴具有红色、黑色和黄色等不同表型，不同表型的玛咖表现出生物活性多样性。
代谢组学及其相关应用：代谢组学是继基因组学和蛋白质组学之后新近发展起来的一门学科，是对细胞、组织、以及其他生物在一特定生理时期内所有低分子量代谢产物同时进行研究的一门新学科。代谢组学发展迅速并且涉及多个领域，比如疾病诊断、毒理学、营养科学、环境学、医药研制开发和植物学等与人类健康护理密切相关的领域。代谢组学研究方法有两种：一种是代谢物指纹分析(metabolomicfingerprinting)，即通过不同方法得到代谢产物的核磁峰或质谱峰，了解代谢产物的化学结构，建立识别不同化合物特征的分析方法；另一种方法是代谢轮廓分析(metabolomicprofiling),研究人员假定一条特定的代谢途径，并对此进行更深入的研究。代谢组学分析包括样品收集与前处理、样品的分析检测、原始数据的处理、化学统计学分析和标志性代谢产物分析。其中样品分析检测技术包括核磁共振技术、质谱技术、气质联用技术和液质联用技术。化学统计学分析主要包括聚类分析、主成分分析、偏最小二乘分析和正交-偏最小二乘分析。
zhangl等人为了评价玛咖样品的质量，采用高效液相色谱(hplc)和液相色谱串联质谱(lc-ms)分析了15批不同地理来源的玛咖下胚轴。在玛咖下胚轴中鉴别出三种芳香族芥子油苷：glucosinalbin(gsb)、glucotropaeolin(gtl)和glucolimnanthin(glh)。hplc图谱显示15批玛咖都含有这三种芥子油苷，但含量有差异。此外，根据玛咖的地理来源和芥子油苷含量，使用hca(等级聚类分析)来评估玛咖质量差异性，结果表明，三种芳香族芥子油苷(gsb，gtl和glh)含量可以作为评价玛咖材料质量的化学分类学标记物。
zhaoj等人为了深入了解玛咖的栽培历史(之前是否种植过玛咖)和生长地点对玛咖(lepidiummeyeniiwalpers)下胚轴组成变化的影响，应用nmr分析结合化学计量学分析研究不同玛咖的代谢物变异性。为了检测代谢物的差异，对在两个不同地点(距离较远)栽培的玛咖、相同地区相近的两块区域但栽培历史不同的栽培玛咖和相同地点栽培的具有不同颜色(黄色，粉红色，紫色和铅色)的玛咖下胚轴进行了分析。1hnmr谱获得的数据用于主成分分析(pca)，进而分析不同栽培条件下种植的玛咖的分离度。通过nmr分析共鉴定了16种代谢物，并且使用单变量统计分析评估与玛咖胚轴的颜色类型和生长条件相关的代谢物水平的变化。pca结果显示不同地点栽培的玛咖在得分图上可以较好分开，pca载荷图显示两个地点玛咖的代谢物不同。单变量统计分析的结果表明，两个种植地点的玛咖中的16种代谢物含量有显著性差异。栽培历史也会对玛咖代谢物产生差异，但差异小于种植地点。由颜色类型引起的代谢物差异小于由栽培因素引起的代谢物差异。
技术实现要素：
本发明针对上述现有技术存在的不足，提供一种基于液质联用代谢组学分析玛咖标志性代谢产物的方法。
本发明解决上述技术问题的技术方案如下：一种基于液质联用代谢组学分析玛咖标志性代谢产物的方法，步骤如下：
(1)玛咖提取液的制备
将玛咖与甲醇混合，经超声提取后，离心，所得上清液过孔径0.22μm滤膜，即得玛咖提取液，置于-80℃备用；
(2)玛咖提取液预处理
将步骤(1)的玛咖提取液用甲醇稀释15-20倍，经过稀释后的玛咖提取液过孔径0.22μm滤膜，得滤液；
(3)分析检测
将步骤(2)的滤液采用uhplc-qexactiveorbitraplc-ms上样分析检测，检测参数为：
色谱柱为acquityuplcbehc18柱2.1×50mm，1.7μm；
正谱条件：洗脱液a相为含0.1％(v:v)甲酸的乙腈溶液，b相为0.1％(v:v)甲酸水溶液；梯度洗脱流程为：0-2min1％a，2-3.25min1％-5t％a，3.25-4.25min5％a，4.25-7.75min5％-55％a，7.75-9.75min55％-90％a，9.75-14.75min90％a，14.75-15min90％-1％a，15-18min1％a；
负谱条件：洗脱液a相为含0.1％(v:v)甲酸的乙腈溶液，b相为10mm醋酸铵；梯度洗脱流程为：0-2min1％a，2-3.25min1％-5％a，3.25-4.25min5％a，4.25-7.75min5％-55％a，7.75-9.75min55％-90％a，9.75-14.75min90％a，14.75-15min90％-1％a，15-18min1％a；
流速：0.3ml/min；进样体积：2μl；1min之前和16min之后不进质谱；
质谱条件：一级质谱分辨率70000(fwhm)，鞘气40arb，辅助气10arb，反吹气0arb，正谱喷雾电压3.5kv,负谱喷雾电压3.1kv，毛细管温度320℃，辅助气温度350℃，扫描范围75-1125，扫描模式fullms；
(4)compounddiscover数据处理
使用compounddiscover软件对原始数据进行提取、分析、整理，通过峰对齐、去卷积、降噪、归一化处理，得到原始数据矩阵；在compounddiscover数据处理过程中设置参数保留时间偏差为0.2min，质量偏差为5ppm，信噪比最大窗口为3，信号强度最大偏差30％，基础离子[m+h]+1、[m-h]-1，未知元素组成设定c90、h190、k、n5、na、o15、p3、s5；
(5)simca软件数据处理
将步骤(4)的原始数据矩阵导入simca14.0进行分析，使用主成分分析pca、偏最小二乘判别pls-da和正交-偏最小二乘判别opls-da分析数据，运用置换排列实验检验；在正交-偏最小二乘判别opls-da中导出的第一主成分的变量权重重要性排序值vip，结合原始数据矩阵中的物质在两种玛咖中峰面积比值ratio和p值，vip值大于1、ratio值大于2或小于0.5且p值小于0.01的代谢物即认定为玛咖标志性代谢产物；
(6)玛咖标志性代谢产物的二级质谱信息
将步骤(5)的玛咖标志性代谢产物进行二级质谱分析，得到二级质谱信息；
质谱条件：二级质谱分辨率17500(fwhm)，扫描模式ms/ms，hcd高能碰撞池碰撞能量nce：30、50、100、150；
二级质谱保留时间与一级质谱保留时间的允许偏差为0.2min，保留时间对应的物质即为不同表型玛咖的标志性代谢产物。
其中，步骤(1)中玛咖与甲醇的用量比为1g：20ml混合，超声提取时间为20min。
本发明的特点和有益效果在于：
本发明的分析方法通过甲醇提取玛咖代谢物，前处理后进行uhplc-qexactiveorbitraplc-ms分析，uhplc-qexactiveorbitraplc-ms数据使用compounddiscover处理，再将数据导入simca14.0软件，使用主成分分析(pca)、偏最小二乘判别分析(pls-da)和正交-偏最小二乘判别分析(opls-da)分析数据，寻找差异性代谢产物；使用置换排列检验方法检验pls-da模型建立有效性，结果显示模型建立有效。
本发明的分析方法操作简单，寻找不同表型玛咖代谢产物，从而用于不同表型玛咖的区分鉴别。
本发明中，uhplc-qexactiveorbitraplc-ms是thermofisher公司的超高效液相与质谱联用仪器。
uhplc-ms：超高效液相色谱与质谱两用技术。
compounddiscover是thermofisher公司开发的与lc-ms配套的数据处理软件。
simca软件为多元统计分析软件。
pca：主成分分析(principalcomponentsanalysis)，是一种无监督的模式统计方法。通过正交变换将一组可能存在相关性的变量转换为一组线性不相关的变量，转换后的这组变量叫主成分。
pls-da：偏最小二乘判别(partialleastsquaresprojectiontolatentstructure-discriminantanalysis)，是一种有监督的模式统计方法，利用偏最小二乘法对数据结构进行投影分析。
opls-da：正交-偏最小二乘判别分析(orthogonal-pls-da)，是将正交信号校正方法(orthogonalsignalcorrection,osc)与pls-da进行结合、从而对pls-da进行修正的分析方法。
附图说明
图1为compounddiscover数据处理工作流程图；
图2是负谱的pca得分图；
图3是负谱的pls-da得分图；
图4是负谱的opls-da得分图；
图5是负谱的排列实验模型图；
图6是正谱的pca得分图；
图7是正谱的pls-da得分图；
图8是正谱的opls-da得分图；
图9是正谱的排列实验模型图；
图10是449号物质的二级质谱图；
图11是655号物质的二级质谱图。
具体实施方式
以下结合实例对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
本发明实施例中使用的仪器、试剂和玛咖样品如下：
1、仪器与软件：
2、试剂：
3、玛咖样品由烟台新时代健康产业有限公司提供，以下为20种玛咖样品表型和产地信息：
实施例1
一种基于液质联用代谢组学分析玛咖标志性代谢产物的方法，步骤如下：
(1)玛咖提取液的制备
称取20种不同表型玛咖的粉末样品各1g，分别加入20ml甲醇混合，分别经超声提取20min后，3900rpm转速下离心20min，所得上清液过孔径0.22μm滤膜，即得玛咖提取液；
(2)玛咖提取液预处理
将步骤(1)的玛咖提取液用甲醇稀释20倍，经过稀释后的玛咖提取液过孔径0.22μm滤膜，得滤液；
(3)分析检测
将步骤(2)的滤液各1ml注入进样瓶，采用uhplc-qexactiveorbitraplc-ms上样分析检测，检测参数为：
色谱柱为acquityuplcbehc18柱2.1×50mm，1.7μm；
正谱条件：洗脱液a相为含0.1％(v:v)甲酸的乙腈溶液，b相为0.1％(v:v)甲酸水溶液；梯度洗脱流程为：0-2min1％a，2-3.25min1％-5t％a，3.25-4.25min5％a，4.25-7.75min5％-55％a，7.75-9.75min55％-90％a，9.75-14.75min90％a，14.75-15min90％-1％a，15-18min1％a；
负谱条件：洗脱液a相为含0.1％(v:v)甲酸的乙腈溶液，b相为10mm醋酸铵；梯度洗脱流程为：0-2min1％a，2-3.25min1％-5％a，3.25-4.25min5％a，4.25-7.75min5％-55％a，7.75-9.75min55％-90％a，9.75-14.75min90％a，14.75-15min90％-1％a，15-18min1％a；
流速：0.3ml/min；进样体积：2μl；1min之前和16min之后不进质谱；
质谱条件：一级质谱分辨率70000(fwhm)，鞘气40arb，辅助气10arb，反吹气0arb，正谱喷雾电压3.5kv,负谱喷雾电压3.1kv，毛细管温度320℃，辅助气温度350℃，扫描范围75-1125，扫描模式fullms；
(4)compounddiscover数据处理
使用compounddiscover软件对原始数据进行提取、分析、整理，通过峰对齐、去卷积、降噪、归一化处理，得到原始数据矩阵；在compounddiscover数据处理过程中设置参数保留时间偏差为0.2min，质量偏差为5ppm，信噪比最大窗口为3，信号强度最大偏差30％，基础离子[m+h]+1、[m-h]-1，未知元素组成设定c90、h190、k、n5、na、o15、p3、s5；
(5)simca软件数据处理
将步骤(4)的原始数据矩阵在第一列前插入“序号”一列，对该列填充序列，导入simca14.0软件，经过行列转换，将“序号”一行及第一列设置为primaryid，ratio值、p值、保留时间和分子式行设置为secondaryid，保存；使用主成分分析pca、偏最小二乘判别pls-da和正交-偏最小二乘判别opls-da分析数据，找出可能的不同表型玛咖标志性代谢产物，并采用置换排列检验方法来检验模型是否有效；
(6)玛咖标志性代谢产物的二级质谱信息
将步骤(5)的玛咖标志性代谢产物进行二级质谱分析，得到二级质谱信息；
质谱条件：二级质谱分辨率17500(fwhm)，扫描模式ms/ms，hcd高能碰撞池碰撞能量nce：30、50、100、150；
二级质谱保留时间与一级质谱保留时间的允许偏差为0.2min，保留时间对应的物质即为不同表型玛咖的标志性代谢产物。
采用compounddiscover软件对原始数据进行处理，图1是compounddiscover数据处理工作流程图，导入文件后进行峰提取、峰对齐、未知物的检测、未知物组成预测等步骤，但由于在compounddiscover处理之前没有进行二级质谱分析，所以未知物组成预测并不准确。
compounddiscover处理后得到的原始数据矩阵进行序号填充处理后导入simca14.0，然后进行行列转换、设置primaryid和secondaryid后保存为excel格式，然后进行pca、pls-da和opls-da处理，并对pls-da进行排列置换检验，结果如图2-图9所示，其中图2-图5是负谱数据，图6-图9是正谱数据。
图2和图6是pca得分图，在pca得分图上我们可以看出黄色、黑色和红色三种表型玛咖的分组趋势不明显，主成分分析属于无监督分析，分析结果可以呈现组内差异和组间差异，由于本实验选择的同种表型玛咖的产地不同，而产地不同也会引起玛咖中的物质发生改变，使得同种表型玛咖产生组内差异，影响三种表型玛咖分组结果。
图3和图7是pls-da得分图，图4和图8是opls-da得分图，pls-da和opls-da属于有监督分析方法，对样品进行指定分组，消除组内差异，有利于发现组间差异，从而寻找标志性代谢产物。opls-da与pls-da相比，可以进一步强化组间差异。图3和图7的pls-da得分图可以看出三种表型玛咖基本可以分开，而图4和图8的opls-da得分图则将三种表型玛咖完全分开。负谱数据的opls-da得分图可以看出，紫色玛咖和黑色玛咖在第一主成分上可以分开，紫色玛咖位于第一主成分的负轴，黑色玛咖位于第一主成分正轴；黑色、紫色玛咖与黄色玛咖在第二主成分可以完全分开，黑色、紫色玛咖位于第二主成分负轴，黄色玛咖位于第二主成分正轴。正谱数据的opls-da得分图可以看出，黄色玛咖和紫色玛咖在第一主成分上可以分开，黄色玛咖位于第一主成分的负轴，紫色玛咖位于第一主成分正轴；紫色、黄色玛咖与黑色玛咖在第二主成分可以完全分开，紫色、黄色玛咖位于第二主成分负轴，黑色玛咖位于第二主成分正轴。
图5和图9是pls-da置换检验结果，当q2在y轴的截距小于0、r2在y轴的截距大于0、且q2和r2的最终值接近1，说明模型建立比较好。由图5和图9可看出正负谱原始模型建立有效。
在正交-偏最小二乘判别(opls-da)中导出的第一主成分的变量权重重要性排序值(vip)，结合原始数据矩阵中的该物质在两种玛咖中峰面积比值(ratio)、p值和正负谱数据量大小，在正谱中选出vip大于1、ratio值大于20或小于0.05并且p值小于0.01的代谢物，在负谱中选出vip大于2、ratio值大于2或小于0.5并且p值小于0.01的代谢物。将这些代谢产物进行二级质谱分析，在二级质谱中保留时间与一级质谱保留时间对应的物质确定为不同表型玛咖标志性代谢产物，保留时间允许偏差为0.2min，其二级质谱信息如以下表1-3所示。
表1、表2和表3分别为黄色、紫色和黑色玛咖的标志性代谢物质二级质谱信息。经过对比文献中玛咖代谢物质的质谱数据，确定黑色玛咖中序号为449的标志性代谢物为n-苄基-15z-烯二十四酰胺，二级质谱图如图10所示。黄色和紫色玛咖中序号为655的标志性代谢物为3-甲氧基苯基乙酸，二级质谱图如图11所示。由于目前对于玛咖中代谢物质定性的数据库和文献较少，449和655以外的物质没有找到与之相对应的二级质谱信息，所以没有定性结果。
表1黄色玛咖标志性代谢产物二级质谱信息
表2紫色玛咖标志性代谢产物二级质谱信息
表3黑色玛咖标志性代谢产物二级质谱信息
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。