代谢组学技术服务-非靶向代谢组学
样本名称 | 最低量 | 一般量 |
动物和临床各类组织样本 | 5-10mg | 50mg |
血样(血清、血浆和全血) | 10μl | 50μl |
尿样 | 10μl | 50μl |
粪便和肠道内容物 | 10mg | 25mg |
体液样本(脑脊液和唾液等) | 5-10μl | 25μl |
植物组织样本(根、茎、叶、花和果实等) | 10mg | 50mg |
细胞和微生物菌体 | 1*10^5cells | 1*10^6cells |
培养基和发酵液 | 10μl | 50μl |
1.数据检查
对所有样本的总离子流色谱图(TIC)色谱图进行可视化检查,如图1:
2.数据预处理:
讲仪器检测得到的原始数据转化为通用格式,然后对质谱数据进行峰匹配,峰对齐和保留时间校正,得到去卷积的质谱数据。然后,对来自于同一个物质的各个峰数据进行归属分析。
3.多维统计分析:
对数据进行后处理,将处理后的数据导入到SIMCA软件进行多维统计分析,在软件中首先进行Pareto格式化(par scaling)和平均中心化(mean-centering)处理,然后再进行PCA、PLS-DA和OPLS-DA等多维统计分析。
● 3.1整体PCA分析和组件PLS-DA
采用SIMCA软件对整体样本进行PAC分析(主成分分析),用于解释和分析各组样本之间的代谢差异。如图2:整体PCA分析
对样本进行PLS-DA分析(偏最小二乘法判别分析),用于分析组间的差异代谢物信息的显著性,并进行Permutation test (置换检验),结合模式的Q2和R2Y,验证模型的可靠性。如图3:组间样本PLS-DA分析和Permutation test
● 3.2 组件OPLS-DA分析和volcano plot分析
为消除无关噪音信息及准确获得两组样本间的显著性差异代谢物信息,我们采用组间OPLS-DA(正交偏最小二乘法判别分析)进行分析,进而获得OPLS-DA模型和相关值(VIP值),作为下一步进行差异物定性筛选的参考。并可对OPLS-DA分析的结果进行Loading Plot 分析和S-plot 分析,另外使用R语言平台进行Volcano Plot (火山图)分析。如图4:组件样本OPLS-DA分析和volcano plot分析、S-plot分析和Volcano Plot分析
4.单维统计分析
对数据进行单维统计分析,dui对数据分别jin'xing进行Shapiro Wilk's test、Welch's Test和Wilcoxon Mann-Whitney test(U test)。呈正态分布的变量则采用Wilk's test结果,而呈非正态分布的变量则采用Wilcoxon Mann-Whitney test结果,综合得到各变量在各对比组之间的显著性分析结果(p-value)
5.差异代谢物定性
多维统计分析(VIP>1)结合单维统计(p-value<0.05)寻找差异性表达代谢物,
采用软件人工结合的方式,使用RT和特征M/Z(GC-MS)或精确分子量和二级质谱(LC-MS)与数据库进行逐一对比分析,差异性代谢物的定性方法为:搜索自建的标准物质数据库、Fiehn GC/MS Metabolomics RTL Library 、Golm Metabolome Database、Metlin、HMDB、KEGG、Lipid和NIST等商业数据库。差异物列表示例如表1:组间样本的差异性代谢物示例
Metabolites | p-value | VIP | FC(A/B) | HMDB | KEGG | Pathway(KEGG) |
pyruvic acid | 2.16E-03 | 1.71 | -0.51 | HMDB00243 | C00022 | Glycolysis/Gluconegenesis;Citrate cycle(TCA cycle);Pentose phosphate pathway |
glucose | 4.11E-02 | 1.82 | 2.65 | HMDB00122 | C00031 | Glycolysis/Gluconegenesis;Pentose phosphate pathway;Galactose metabolism |
gluconic acid | 1.52E-02 | 1.63 | 0.68 | HMDB00625 | C00257 | Pentose phosphate pathway |
mannitol | 4.11E-02 | 1.45 | 1.47 | HMDB00765 | C00392 | Fructose and mannose metabolsm |
dulcitol | 2.16E-03 | 1.74 | 3.72 | HMDB00107 | C01697 | Galactose metabolism |
galactonic acid | 2.16E-03 | 1.96 | 2.24 | HMDB00565 | C00880 | Galactose metabolism |
ethanolamine | 2.16E-03 | 2.10 | 1.00 | HMDB00149 | C00189 | Glycerophospholopid metabolism |
6.相关性分析
●6.1 Pearson Correlation 分析
为了表征各差异性代谢物之间的(浓度)相关性,我们会对这些物质的定量信息经行Pearson Correlation分析。如图5:差异性代谢物的相关性矩阵图
●6.2热图分析:
为了表示差异物之间的聚类关系,我们会对这些物质的定量信息进行热图分析,如图6:差异性代谢物的热图
7.代谢通路分析
●7.1代谢通路进行归类分析
我们采用KEGG数据库对每个差异代谢物所属的代谢通路进行归类分析,如图7:差异性代谢物所属KEGG代谢通路实例
●7.2metametabonanalyst pathway分析
我们用软件metametabonanalyst对差异性代谢物进行pathway analysis,metabolome view如图8所示:
Pathway views如表2所示,total表示该途径所含的代谢物总数,hits表示该途径含有差异性代谢物数,-log(p)表示图8纵坐标值,impact表示图8横坐标值。参数如表2组间样本的pathway view所示:
Pathyway | Total | Expected | Hits | Raw p | -log(p) | Holm adjust | FDR | Impact |
Pantothenate and CoA biosynthesis | 16 | 0.39402 | 2 | 0.05672 | 2.8695 | 1 | 1 | 0 |
Zeatin biosynthesis | 16 | 0.39402 | 2 | 0.05672 | 2.8695 | 1 | 1 | 0 |
Butanoate metabolism | 20 | 0.49252 | 2 | 0.08457 | 2.4702 | 1 | 1 | 0 |
Alanine,aspartate and glutamate metabolism | 21 | 0.51715 | 2 | 0.09208 | 2.3851 | 1 | 1 | 0 |
C5-Branched dibasic acid metabolism | 4 | 0.098505 | 1 | 0.09504 | 2.3534 | 1 | 1 | 0 |
Citrate cycle(TCA cycle) | 20 | 0.49252 | 1 | 0.39526 | 0.9282 | 1 | 1 | 0 |
Pyruvate metabolism | 20 | 0.49252 | 1 | 0.39526 | 0.9282 | 1 | 1 | 0.1148 |
Purine metabolism | 55 | 1.3544 | 2 | 0.39526 | 0.9244 | 1 | 1 | 0.065 |
单位名称: |
详细地址:
上海市嘉定区宝安公路4997号国际医疗产业园1号楼
|
qq:
|
官网: |
联系电话: |
Email: |