基于GC-MS获取转基因和非转基因玉米代谢差异的方法与流程

本发明涉及分析化学领域，具体涉及基于gc-ms获取转基因和非转基因玉米代谢差异的方法。
背景技术：
：随着转基因技术的飞速发展以及转基因玉米的大面积推广应用，在关注转基因玉米所带来的巨大的社会、经济和生态效益的同时，转基因玉米及其产品的安全性问题也引起了世界范围内的广泛关注，在发展我国农业生物技术的同时，对转基因玉米及其产品的安全性加强监督非常重要。对转基因玉米的安全性评价主要从两个方面来考虑：一个是环境安全性，一个是食用安全性。随着转基因玉米大面积商业化种植，人们对其潜在的生态环境危害尤其是转基因食品对消费者的食用安全性表示担忧。1993年经济合作与发展组织提出了实质等同性原则，即如果转基因食品在组成、营养价值、拟议用途上没有变化，那么可以认为它与非转基因食品是等同的。此后实质等同性原则被广泛应用于转基因产品的安全性评价、审批、监管。目前，代谢组学方法已成为转基因产品安全评估的重要工具。在代谢组学的研究中，已经建立的方法包括气-质联用(gc-ms)、液相层析串联质谱(lc-ms)、核磁共振和傅里叶变换质谱(ft-ms)等分析技术。其中，gc-ms分析方法长期以来被用于植物和微生物的代谢物谱图分析，可被视为对代谢物分析的黄金标准。这项技术的主要优点是具有高灵敏度和分析覆盖面宽广，是代谢组学研究的核心方法之一，适合分析挥发性有机物和亲水性代谢物。但目前由于gc-ms预处理过程分多步进行，样品准备过程中容易引入误差；气质分析平台精确度低，且操作复杂，费时费力，成本较高，使得一次性获得的化合物信息较少，不能有效地完成转基因非预期效应的高通量检测。技术实现要素：本发明的目的在于，提供基于gc-ms获取转基因和非转基因玉米代谢差异的方法；即，建立一种基于气相色谱-高分辨率质谱联用系统进行玉米叶片代谢轮廓分析，并对转基因和非转基因玉米叶片代谢差异进行研究的方法。所述方法采用agilent7890b气相色谱与lecopegasusht质谱联用分析平台对玉米叶片提取液进行分析，得到玉米叶片代谢轮廓，然后用多变量统计方法区分转基因和非转基因玉米数据，获得对代谢表型差异有重要贡献的化合物。与传统方法相比，该分析平台更适合快速色谱分析，可以快速可靠地识别、定量和确认更多化合物；更重要的是，其检测分辨率高，保证谱峰扫描点的完整性，不遗漏任何细节，实现高效可靠的代谢物信息识别，得到的结果更真实可靠。所述方法包括如下步骤：1)取转基因玉米和非转基因玉米叶片，粉碎至碎末状，然后进行提取，并加入内标物，获得含有内标物的待测样品；2)将所述待测样品进行gc-ms检测，获得所述转基因玉米和非转基因玉米的代谢物原始谱图信息；3)采用lc-ms数据处理软件对步骤2)所获得的代谢物原始谱图信息进行提取，然后再采用正交偏最小二乘法-判别分析的统计方法进行分析，获得代谢指纹数据；4)根据步骤3)获得代理指纹数据进行建模，获得opls-da模型；5)根据所述opls-da模型，获得转基因和非转基因玉米的差异代谢物。其中，所述转基因玉米为转双斑萤叶甲mhsnf7基因玉米；所述非转基因玉米为tj806非转基因玉米。本发明所述方法具有检测分辨率高、操作简单的优点，在提高灵敏度的同时又节省了时间，降低了成本。所述方法可以为深入研究转基因玉米的安全性提供根据。为进一步提高该方法灵敏度和效率等优势，本发明提供下述优选方案。本发明进一步提出的，步骤1)具体为：取5～6叶期的玉米嫩叶，采用液氮，研磨至碎末状，获得碎末样品；将所述碎末样品加入提取液中，并加入内标物；在40～50hz的条件下，研磨3～5min；然后在冰水浴中超声4～6min；重复研磨、冰水浴的步骤3次，再低温静置20～40min，高速离心后，弃掉上清液体，保留下层沉淀，将沉淀物干燥处理，然后再进行孵育处理，获得含有内标物的待测样品；更优选的，所述提取液为体积比为3.5～4.5:1的甲醇和水。优选的，所述碎末样品于所述提取液的固液比为100～150:1(g/l)。优选的，所述内标物为质量比为1:1.5～2.5核糖醇和l-2-氯苯丙氨酸。所述核糖醇的浓度为0.4～0.6mg/ml，所述l-2-氯苯丙氨酸的浓度为0.9～1.1mg/ml。所述内标物的加入按照常规方式加入，如0.4ml提取液中加入10ml内标物。优选的，所述研磨采用常规的研磨仪；进一步优选的，研磨过程中，加入的大磁珠直径为5mm，小磁珠直径为3mm。优选的，所述低温静置采用常规的低温方式，如-40℃。优选的，所述高速离心采用本领域常规方式进行，如13000rpm，4℃离心15min。优选的，所述孵育处理具体为：在干燥后的样品中加入甲氧胺盐试剂，孵育25～35min，然后再加入含有0.8～1.2％三甲基氯硅烷的n,o-双(三甲基硅烷基)三氟乙酰胺，在60～80℃的温度下，孵育1～2小时，离心，取上清液。更优选的，所述孵育处理具体为：在干燥后的样品中加入50μl甲氧胺盐试剂(甲氧胺盐酸盐，溶解于吡啶，浓度为20mg/ml)，孵育25～35min，然后再加入含有1.0％三甲基氯硅烷(tmcs)的n,o-双(三甲基硅烷基)三氟乙酰胺(bstfa)，在60～80℃的温度下，孵育1～2小时，离心，取上清液。本发明进一步提出的，步骤2)中，所述gc-ms检测的条件具体为：(1)毛细管柱：agilentdb-5ms(30m×250μm×0.25μm)。(2)进样量0.8～1.2μl。(3)gc-ms检测参数如下：分流模式：splitlessmode隔垫吹扫流速：2.5～3.5ml/min载气：helium柱流速：0.8～1.2ml/min柱箱升温程序：50℃(1min)-310℃(8～12℃/min)，6～10min前进样口温度：270～290℃传输线温度：270～290℃离子源温度：240～260℃电离电压：-65～-75ev质量范围：m/z:50-500扫描速率：18～22spectrapersecond溶剂延迟时间：6～6.5min；优选的，所述gc-ms检测参数如下：分流模式：splitlessmode隔垫吹扫流速：3ml/min载气：helium柱流速：1ml/min柱箱升温程序：50℃(1min)-310℃(10℃/min),8min前进样口温度：280℃传输线温度：280℃离子源温度：250℃电离电压：-70ev质量范围：m/z:50-500扫描速率：20spectrapersecond溶剂延迟时间：6.27min更优选的，质谱检测过程中还包括校验，所述校验具体为：将所述待测样品的混合物作为qc样，先将所述qc样进样一次，然后每分析3～6个待测样品后，再将所述qc样进样一次。通过保证qc样总离子流图的重复性确保实验过程中气相色谱质谱联用仪有较好的稳定性。本发明进一步提出的，步骤3)中，采用xcms处理软件对代谢物原始谱图信息进行提取。所述校验具体为：采用lc-ms数据处理软件xcms对步骤2)所获得的代谢物原始谱图信息，进行保留时间矫正、峰识别、峰提取、峰积分、峰对齐处理。本发明进一步提出的，步骤5)中，卡值标准为t检验(student’st-test)的p值(p-value)小于0.05，并且opls-da模型第一主成分的变量投影重要度(variableimportanceintheprojection,vip)大于1的化合物为转基因玉米和非转基因玉米的差异代谢物。符合卡值标准的化合物即为转双斑萤叶甲mhsnf7基因玉米(tj806)与tj806非转基因玉米的差异代谢物。本发明提供一种优选方案，所述方法包括如下步骤：1)取5～6叶期的玉米嫩叶，采用液氮，研磨至碎末状，获得碎末样品；将所述碎末样品加入提取液中，并加入内标物；在40～50hz的条件下，研磨3～5min；然后在冰水浴中超声4～6min；重复研磨、冰水浴的步骤3次，再低温静置20～40min，高速离心后，弃掉上清液体，保留下层沉淀，将沉淀物干燥处理，在样品中加入甲氧胺盐试剂，孵育25～35min，然后再加入含有0.8～1.2％三甲基氯硅烷的n,o-双(三甲基硅烷基)三氟乙酰胺，在60～80℃的温度下，孵育1～2小时，离心，取上清液，获得含有内标物的待测样品；其中，所述提取液为体积比为3.5～4.5:1的甲醇和水；所述内标物为质量比为1:1.5～2.5核糖醇和l-2-氯苯丙氨酸；2)将所述待测样品进行gc-ms检测，获得所述转基因玉米和非转基因玉米的代谢物原始谱图信息；其中，所述gc-ms检测的条件具体为：(1)毛细管柱：agilentdb-5ms；(2)进样量0.8～1.2μl；(3)gc-ms检测参数如下：分流模式：splitlessmode隔垫吹扫流速：3ml/min载气：helium柱流速：1ml/min柱箱升温程序：50℃(1min)-310℃(10℃/min),8min前进样口温度：280℃传输线温度：280℃离子源温度：250℃电离电压：-70ev质量范围：m/z:50-500扫描速率：20spectrapersecond溶剂延迟时间：6.27min优选的，质谱检测过程中还包括校验，所述校验具体为：将所述待测样品的混合物作为qc样，先将所述qc样进样一次，然后每分析3～6个待测样品后，再将所述qc样进样一次。3)采用lc-ms数据处理软件xcms对步骤2)所获得的代谢物原始谱图信息进行保留时间矫正、峰识别、峰提取、峰积分、峰对齐处理；然后再采用正交偏最小二乘法-判别分析的统计方法进行分析，获得代谢指纹数据；4)根据步骤3)获得代理指纹数据进行建模，获得opls-da模型；5)根据所述opls-da模型，卡值标准为t检验的p值小于0.05，并且opls-da模型第一主成分的变量投影重要度大于1的化合物为转基因玉米和非转基因玉米的差异代谢物。本发明的又一目的在于，提供一种基于上述方法的在转基因玉米安全评价中的应用。本发明建立的一种基于gc-ms获取转基因和非转基因玉米代谢差异的方法。样本前处理过程简单、快速，引入的误差小。与传统方法相比，该分析平台更适合快速色谱分析，可以快速可靠地识别、定量和确认更多化合物；更重要的是，其检测分辨率高，保证谱峰扫描点的完整性，不遗漏任何细节，实现高效可靠的代谢物信息识别，得到的结果更真实可靠。同时，样本预处理过程中加入内标，校正预处理过程和仪器运行过程中响应微小漂移带来的误差。样本序列分析过程中插入由所有样本混合的质量控制(qc)样本，通过对qc样本总离子流图的监控，可以检测分析序列中仪器响应有无漂移，确保了实验数据的可靠性。本发明基于gc-ms筛选出多个对分类有重要贡献的化合物，可以为转基因植物的非预期效应和安全性评价提供依据。附图说明图1为转双斑萤叶甲mhsnf7基因玉米与tj806非转基因玉米叶片样本gc-ms代谢轮廓对比图；图2为opls-da模型得分散点图。具体实施方式以下实施例用于说明本发明，但不用来限制本发明的范围。实施例1本实施例提供一种基于gc-ms获取转基因和非转基因玉米代谢差异的方法，所述方法包括如下步骤：一)待测样品的准备：1)将长至5-6叶期的转双斑萤叶甲mhsnf7基因玉米(tj806)及tj806非转基因玉米嫩叶片收集，并迅速冷冻于液氮中，储存于-80℃冰箱中备用；2)将从低温冰箱中取出的叶片快速地置于液氮预冷的研钵中；倒入少许液氮，用预冷的捣杵进行研磨，直至叶片呈碎末状态；3)称量50±1mg的碎末样品，置于2ml离心管中，将装有样品的离心管置于冰上；在样品离心管中加入一大一小两粒研磨磁珠。大磁珠直径为5mm，小磁珠直径为3mm；4)每个样品管中加入4.5ml提取液(甲醇：水的体积例＝4:1)；5)每个样品管中加入10μl内标混合液(0.5mg/ml核糖醇和1mg/mll-2-氯苯丙氨酸)；6)将样品管置于研磨仪上进行研磨(45hz，4min)；然后在冰水浴中超声5min；重复研磨、冰水浴的步骤3次，再将样品管置于-40℃冰箱中静置半小时；静置结束后，以13000rpm的转速，4℃的温度下，离心15min，弃掉上清液体，保留下层沉淀；7)42℃进行干燥，在干燥后的样品中加入80μl甲氧胺盐试剂(20mg/ml)，80℃孵育30min；8)加入含有1％三甲基氯硅烷(tmcs)的n,o-双(三甲基硅烷基)三氟乙酰胺(bstfa)100μl，70℃孵育1.5小时；离心，取上清液，获得含有内标物的待测样品；二)将所述待测样品进行gc-ms检测，获得所述转基因玉米和非转基因玉米的代谢物原始谱图信息；其中，所述gc-ms检测的条件具体为：(1)毛细管柱：agilentdb-5ms(30m×250μm×0.25μm)；(2)进样量1μl；(3)gc-ms检测参数如下：表1.gc-ms检测参数项目参数分流模式splitlessmode隔垫吹扫流速3ml/min载气helium柱流速1ml/min柱箱升温程序50℃(1min)-310℃(10℃/min),8min前进样口温度280℃传输线温度280℃离子源温度250℃电离电压-70ev质量范围m/z:50-500扫描速率20spectrapersecond溶剂延迟时间3ml/min质谱检测过程中还包括校验，所述校验具体为：将所述待测样品的混合物作为qc样，先将所述qc样进样一次，然后每分析3～6个待测样品后，再将所述qc样进样一次，通过保证qc样总离子流图的重复性确保实验过程中气相色谱质谱联用仪有较好的稳定性；3)采用lc-ms数据处理软件xcms对步骤2)所获得的代谢物原始谱图信息进行保留时间矫正、峰识别、峰提取、峰积分、峰对齐处理；然后再采用正交偏最小二乘法-判别分析的统计方法进行分析，获得代谢指纹数据；4)根据步骤3)获得代理指纹数据进行建模，获得opls-da模型；5)根据所述opls-da模型，卡值标准为t检验的p值小于0.05，并且opls-da模型第一主成分的变量投影重要度大于1的化合物为转基因玉米和非转基因玉米的差异代谢物。符合卡值标准的化合物即为转双斑萤叶甲mhsnf7基因玉米(tj806)与tj806非转基因玉米的差异代谢物；获得7个转双斑萤叶甲mhsnf7基因玉米(tj806)与tj806非转基因玉米的差异代谢物，见表2。表2差异代谢物筛选表idpeaksimilarityrtvipp_value131乙醇胺94410.21512.10.0057480亚麻酸93920.85882.70.0007251α-酮戊二酸93214.25551.70.0430571纤维二糖190124.60122.00.0150478亚油酸88220.7862.30.004695磷酸甲酯8588.999012.20.0016536尿苷284822.72291.30.0260表2为转双斑萤叶甲mhsnf7基因玉米(tj806)与tj806非转基因玉米的差异代谢物。id为物质在本次定性匹配中得到的唯一编号；peak为物质在数据库中的名称；similarity为该物质与质谱检测峰的匹配度打分值；rt为保留时间；vip为来自opls-da模型的vip值；p_value为来自t-test的p值。图1为转双斑萤叶甲mhsnf7基因玉米(tj806)与tj806非转基因玉米叶片样本gc-ms代谢轮廓对比。其中，a为tj806非转基因玉米样本的总离子流图(tic图，上方图)；b为转双斑萤叶甲mhsnf7基因玉米(tj806)样本的总离子流图(tic图，下方图)。两者的差异在于，在同一时间段出现的峰不同(箭头所示)图2为opls-da模型得分散点图。横坐标t[1]p表示第一主成分的预测主成分得分，纵坐标t[1]o表示正交主成分得分，t表示样本在主成分投影值；■为非转基因样本，●为转基因样本，每个点表示一个样品。从opls-da得分图的结果可以看出，两组样本区分非常显著，样本全部处于95％置信区间内。虽然，上文中已经用一般性说明、具体实施方式及试验，对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。当前第1页12