绘制包含Fold Change的差异代谢产物HMDB(KEGG)分类图

本文介绍如何使用GraphPad Prism和Adobe Illustrator绘制差异代谢物的HMDB分类图,包括数据准备、初步绘图及精加工步骤。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

前言

前两天画KEGG的分类图师兄路过看到了让帮也画一个,一问是两个分组,于是建议参考这篇 Multi-omics reveals that the rumen microbiome and its metabolome together with the host metabolome contribute to individualized dairy cow performance 画差异代谢物的HMDB图就行,还能标出fold change值,简洁明了

数据准备工作

  • 找到包含差异代谢产物信息的公司结果对应的如Con vs Exp.pos.annotation.csv文件,用 Excel 打开
image-20220112195302841
  • 利用数据透视表,把差异代谢物根据HMDB的SuperClass进行分类,并标出log2FC值,整理成如下(为了之后出图美观,这一步可对每一组的产物根据log2FC值进行排序,选择up值在前还是down值在前;或者把文章中要讨论的关注的点放在前面等)
image-20220112195434325

图片的初步绘制

graphpad prism绘制柱状图

  • 打开 graphpad prism ,选择默认的 Column 模板

    image-20220112200056925

  • 录入数据(可将上面数据透视表的数据在excel中用转置粘贴功能变成横向的,方便graphpad录入)

image-20220112200337510
  • 根据数据集创建Graphs,如下图选择横向柱状图,plot选择均值mean就行,因为每个差异代谢物只对应一个log2FC值
image-20220112200557497
  • 生成后如下图
image-20220112200926731

图片的粗加工

  • 这时graphpad直接生成的图y轴纵坐标连差异代谢物的名字条目都没显示全,我录入了26种差异代谢物,这里只显示了十来个,缺失是因为限制了显示高度

    • 双击y轴差异代谢物的名字,在弹出的对话框中,红框处调整合适的高度大小,绿框是给每个代谢物间加间隔虚线(美观)
    image-20220112201342175
    • 调整之后
    image-20220112201742765
  • 针对x轴和y轴,以及bar的粗细颜色等细节,可以直接导进AI里调整,也可以在graphpad里修,这里举例去掉y轴刻度线:双击y轴,在弹出的对话框选择如下

image-20220112202040654 image-20220112202156874

Adobe Illustrator精加工

  • 将graphpad中的图片导出为pdf格式,打开 Adobe Illustrator,从右上角文件-打开处导入刚刚存的pdf图(pdf格式方便修改,这样导入才能编辑,直接拖入无法编辑)
  • 双击选中y轴和x轴0刻度的分界线按键盘Del键删除,鼠标左键配合shift键多选bar/name调整颜色
  • 选择直线段工具,根据前面数据透视表画竖线,并选择文字工具编辑对应HMDB分类
  • 调整间距,对齐等细节,最终图如下
image-20220112202855478
fold change差异表达基因分析中常用的指标之一,用于衡量基因在不同条件下的表达水平差异fold change表示基因在两个条件之间的表达水平的相对倍数差异。它通过计算两个条件下的平均表达水平之比得出。 例如,如果一个基因在条件A中的表达水平为10,而在条件B中的表达水平为20,则该基因的fold change为20/10=2。这表示在条件B下,该基因的表达水平是在条件A下的2倍。如果fold change大于1,则表示基因在条件B中的表达水平高于条件A;如果fold change小于1,则表示基因在条件B中的表达水平低于条件A。 在差异表达基因分析中,fold change通常与p-value一起使用,以确定哪些基因在不同条件下的表达差异是否显著。p-value用于衡量差异的显著性,而fold change则用于衡量差异的大小。 总结起来,fold change是用于衡量基因在不同条件下表达水平差异的指标,它表示基因在两个条件之间的表达水平的相对倍数差异。它与p-value一起被用来确定差异是否显著。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *3* [如何计算基因表达中的Fold Change值?](https://blog.csdn.net/u013313168/article/details/99309076)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] - *2* [vennt:差异基因表达的动态维恩](https://download.csdn.net/download/weixin_42109178/20258904)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值