KEGG富集分析弦图绘制

KEGG富集分析弦图绘制

在阅读文献过程中经常可以看到一些比较复杂的弦图,比如说《A simplified herbal formula for the treatment of heart failure: Efficacy, bioactive ingredients, and mechanisms》这篇论文中的弦图(下图所示),首先该图总体想要传递的信息为不同基因所对应的KEGG通路,其次细分下来还包含了不同样本的基因表达量水平,高表达为红色,低表达为蓝色,白色为正常表达量。

attachments-2022-06-9kspXjEH62a96b29b65a3.png

这样的一个复杂的弦图,我们可以通过在线工具Charticulator来实现。Charticulator 是微软研究院2018年推出的一款可视化工具,可以通过交互创作的方式制作图表类型模板。(链接:https://charticulator.com/app/index.html)

Charticulator的用户界面如下图所示,可以看到基本上跟Photoshop等软件类似,包括数据面板(左侧)、Glyph画布(中间上部)、图层面板(中间中部)、可视化属性面板(中间下部)和Chart画布(右侧大面积区域)。

attachments-2022-06-awqsMrxT62a96b53c8014.png

然后有一系列密密麻麻的按钮排列在上侧,从左到右分别是Mark(各种标记)、DDG、nestchart(为了做多层嵌套图表而设计新按钮,在论文中并未涉及,是后期加上的实验性功能)、Guide(辅助线)、PlotSegment(图表分块)、Scaffold(布局脚手架)。与数据绑定的交互一致,这些按钮都遵循拖拽的方式使用。

以上面的弦图为例,简单介绍Charticulator的用法。

一、整理数据   

1.连线信息(基因与信号通路的对应信息),如下图所示

attachments-2022-06-1nJQu1r762a96c4891c52.png

2.节点信息(统计每一通路对应基因数目,以及每一基因对应的通路数目,基因表达水平分类,节点类型),如下图所示

attachments-2022-06-JkUjFUxF62a96c526aad1.png

注:id为基因名或通路名称,total为通路对应基因数目或基因对应的通路数目,C、H7、H14基因在不同样本中的表达量水平(1为中表达,0为低表达,2为高表达,3为通路),type为节点类型(gene或pathway)

二、导入数据

attachments-2022-06-Q3GHROVI62a96c62788ab.png

三、绘制弦图,在绘制弦图前首先捋清楚弦图的各个元素,并排列各个元素的优先顺序,类似PS中的图层。

1.连线links

attachments-2022-06-V6UaLoAX62a96c77dbeb7.pngattachments-2022-06-RH93chwO62a96c8abd2bc.pngattachments-2022-06-wp6xrjqB62a96c9d5069f.png

注:Attribute区Type可调整连线样式

attachments-2022-06-ZjXqN96Y62a96cb1e93eb.png

最终连线样式:

attachments-2022-06-FzFfX9YJ62a96cc2c50f9.png

2.添加基因矩形

由于还要突出不同样本的基因表达量水平,因此优先绘制基因矩形,基因矩形的叠加顺序从下至上依次为外圈、中圈、内圈

(1).外圈

attachments-2022-06-GeB1kELt62a96cd95d9dc.pngattachments-2022-06-9HdD4YPN62a96ced3ddcc.pngattachments-2022-06-fwci2dr462a96d01f2fd3.pngattachments-2022-06-MOyDIIvg62a96d0c4086e.pngattachments-2022-06-1a2rbZyF62a96d19502f8.pngattachments-2022-06-rdmzJ9LK62a96d2504aa8.png

(2).重复1,2步骤,添加中圈矩形

attachments-2022-06-WII4Mw3a62a96d396d1f3.png

重复4,5步骤来填充中圈颜色,可添加矩形边框

(3).内圈同上

效果如下

attachments-2022-06-6UIHb7lO62a96d483322d.png

3.添加通路矩形

attachments-2022-06-dBfCRzKP62a96d53aec01.pngattachments-2022-06-TTKAZaqk62a96d6230611.pngattachments-2022-06-EjXsG3nh62a96d6d6b0a8.png

4.调整links颜色与通路颜色一致

调整后效果如下

attachments-2022-06-NCdgvWGR62a96d7d87395.png

5.添加基因与通路名称

attachments-2022-06-sGe58kzd62a96d89a3af5.pngattachments-2022-06-vvLgSvF262a96d957daf8.pngattachments-2022-06-NDS2W0x762a96d9d83502.pngattachments-2022-06-NxZQvv9462a96dafcb2c6.png

最后调整文本大小,使其相适应

效果如下

attachments-2022-06-vnSosWGX62a96dbfb8e2b.png

注:标题node可通过layers区进行修改删除

attachments-2022-06-ssRi286v62a96dca60035.png

四、结果图导出,可保存为不同格式。

attachments-2022-06-f28z0cBU62a96dd589ab1.png

svg格式可通过AI进行细节修改,比如为连线添加边框线。

最终结果图如下

attachments-2022-06-FTegPCnl62a96de479ed0.png

**************************************************************************

除弦图以外,该工具还可绘制多种图形,以下为一些样式图:

attachments-2022-06-wcRTpr5G62a96e02e7b59.png

延伸阅读

  1. GEO数据库挖掘—WGCNA鉴定骨肉瘤转移相关基因
  2. GEO、TCGA多数据库联合挖掘胰腺导管腺癌预后关键基因
  3. 文献精读-GEO数据挖掘生物信息文章(宫颈癌)
  4. GEO数据挖掘直肠
  5. GEO数据挖掘案例解读-植物篇(拟南芥重金属)
  6. GEO数据如何挖掘?案例解析!
  7. 胃癌免疫侵润预后Signature数据挖掘-糖酵解
  8. IF:8.8|TCGA+ceRNA+肝癌文献解析
  9. IF=8.78|非肿瘤数据挖掘思路(GWAS Catalog)

attachments-2022-06-zi7wOCrr62a96e190a1fa.png

  • 发表于 2022-06-15 13:29
  • 阅读 ( 2810 )
  • 分类:软件工具

你可能感兴趣的文章

相关问题

0 条评论

请先 登录 后评论
你想桃子吗
你想桃子吗

21 篇文章

作家榜 »

  1. omicsgene 702 文章
  2. 安生水 350 文章
  3. Daitoue 167 文章
  4. 生物女学霸 120 文章
  5. xun 82 文章
  6. 红橙子 78 文章
  7. rzx 76 文章
  8. CORNERSTONE 72 文章