无需写代码的高颜值富集分析神器

富集分析是了解一个基因集功能倾向性的一个方式,在组学研究领域应用广泛。常见的有基于差异基因的Over-representation分析,也就是常说的GO、KEGG富集分析和Functional class scoring分析,如GSEA。这两种富集分析算法不同,但可以都支持同样的注释集,如GO、KEGG或其他类型的注释。基本原则只有一个:基因集的基因名字与注释集的基因名字能匹配。剩下的就是了解下原理去操作了。

下面三篇文章建议仔细多读几遍,尤其是原理解释部分。对以后的分析事半功倍。做了两年多的培训的体验是:学不会是因为没时间。时间用到了,自然就会了,培训现场的进度就是最好的支持证据。我们的培训希望在3天时间告诉大家需要3个月到半年可以掌握的知识,多读几遍也就都会了 (跑题了)。

今天给大家介绍的是Cytoscape的插件——ClueGO。ClueGO能以网络图的形式来展示GO富集结果。

除了在线工具和R代码,今天带来的是网络绘制工具Cytoscape的一个插件ClueGO,支持200多个物种,支持数据库实时更新,支持多个注释数据集,支持多集合富集结果比较和网络图展示。配合Cytoscape的网络操作,可以做更多探索。下面看看怎么安装和使用吧。

ClueGO安装

官网下载

(http://apps.cytoscape.org/apps/cluego)
40c8fcfa9df78346da641004e4ec48d1.png

Cytoscape的APP Manager中下载

16c8c4e73d262bbe714f1de59f2622ab.png
67aa40ef70fa8165d5b7e1ef113d66f5.png

ClueGO激活

ClueGO第一次使用需要按照提示申请license,申请网址为:http://www.ici.upmc.fr/cluego/cluegoLicense.shtml。点击Request a license key,填写基本信息之后提交审核,审核时间可能需要几天。拿到license之后点击App-ClueGO,在弹出的提示框中输入license即可激活ClueGO。

简单使用ClueGO

1. 上传数据

可以粘贴基因的idsymbol等或从本地上传文件 (富集分析的基因集只需要基因名字),我们以ClueGO的示例文件“CluePediaExampleFile_Th1RelatedGenes.txt”为例,并选择其对应物种和自动识别Gene id类型选项。
330ae1a2d0c722a43c43a9c0ab7e9f93.png

2.选择富集数据集和富集标准

默认是BP(Biological Process),我们增加常用的CC(Cellular Component)和MF(Molecular Function)
07b001f129ca8e02c94d7068ed04c3f8.png
然后勾选以下选项,目的是只显示p值<=0.05的富集结果。具体选择的富集显著性P value,可自行调整。
1e7c8d1da1d4f1e49698ad92dc2541e3.png

3.富集分析和结果展示

点击start开始进行富集分析
b376647348bcdb1e66cdbd1a8b695a34.png
运行结束后,会跳出一个简述富集结果信息的提示框,包括物种基因名字匹配的注释信息重叠度冗余度等信息。
3e00091da8708a4d520c23ca66168341.png
Cytoscape的网络视图中展示的是富集结果的网络图。每个节点代表性富集通路节点的连线表示通路之间的共有的基因数颜色表示该节点的富集情况分类(隶属于哪个功能组,颜色也可以和表格形式的的ClueGO结果对应)。
9647f8e9d26c025ceedd83930c704e05.png

结果导出为矢量图
e4771e968ae4a78217ac4098f221969b.png

结果保存为表格或Excel格式
a1177648ac24c0f83d66516842586be4.png

ClueGO优势

1. 物种丰富,并支持手动下载

本地有人类和小鼠可选,也可根据需要下载,支持近200个物种。


6aed48a23f8bbd59f04745d58f0165ef.png

2. 物种注释手动更新

d9d5c37ea907fc3f6b3d869bf91122e5.png
(前面推荐的GOEAST是每周更新,这个可以认为是实时更新。R中的注释包相比来说更新的慢了些,需要社区发布新版本或自己整理。)

3. 多种网络布局和自由网络设计

多种网络布局类型选择
1aa79006ff1535a3f58af7eba3e7b49c.png
style选项可对网络的颜色,节点,网络等进行修改
51ace41ff00b08d0c2c58aae5c7a86af.png

4. 多个数据集比较

同时导入2组或多组基因列表,可以进行比较分析。
35de09ef10babf91debf40a9b743486b.png
富集结果:
首先跳出了一个简述富集结果信息的提示框,和单组分析类似,包括物种基因名字匹配的注释信息重叠度冗余度等信息,但是额外统计了了两组共有的基因的富集结果。
8dd5f9d780d47b32d02d3961a4227932.png
Cytoscape的网络视图中展示的是富集结果的网络图。

每个节点代表性富集通路节点的连线表示通路之间的共有的基因数颜色表示该节点的富集情况分类,而两组单独的,综合的,共有的富集通路等均可在ClueGO results详细看到。
8425d4a3a4cbcc10bc539fe641fd0b6b.png
27ef4a64750c667da7f4774a40ab3f6d.png

5. 链接外源数据库

看到感兴趣的基因、通路,可以直接选择需要的外源数据库查询获得更多信息。
486050e56d10e1d6846de147f8f39301.png

更多Cytoscape的使用

1. Cytoscape教程1
2. Cytoscape之操作界面介绍
3. 新出炉的Cytoscape视频教程
4. Cytoscape制作带bar图和pie图节点的网络图
5. Cytoscape: MCODE增强包的网络模块化分析
6. 生信宝典之傻瓜式 (五) 文献挖掘查找指定基因调控网络
7. 没钱买KEGG怎么办?REACTOME开源通路更强大

学习津贴

单篇留言点赞数的第一位(点赞数至少为8)可获得我们赠送的在线基础课的9折优惠券

越留言,越幸运。

主编会在每周选择一位最有深度的留言,评论者可获得我们赠送的任意一门在线课程的9折优惠券(偷偷告诉你,这个任意是由你选择哦)。

高颜值免费在线绘图

c9c7e9cfeb7baf262380e0b1311c79e4.png

猜你喜欢

10000+:菌群分析 宝宝与猫狗 梅毒狂想曲 提DNA发NatureCell专刊肠道指挥大脑

系列教程:微生物组入门 Biostar 微生物组  宏基因组

专业技能:学术图表 高分文章 生信宝典 不可或缺的人

一文读懂:宏基因组 寄生虫益处 进化树

必备技能:提问 搜索  Endnote

文献阅读 热心肠 SemanticScholar Geenmedical

扩增子分析:图表解读 分析流程 统计绘图

16S功能预测   PICRUSt  FAPROTAX  Bugbase Tax4Fun

在线工具:16S预测培养基 生信绘图

科研经验:云笔记  云协作 公众号

编程模板: Shell  R Perl

生物科普:  肠道细菌 人体上的生命生命大跃进  细胞暗战 人体奥秘  

写在后面

为鼓励读者交流、快速解决科研困难,我们建立了“宏基因组”专业讨论群,目前己有国内外5000+ 一线科研人员加入。参与讨论,获得专业解答,欢迎分享此文至朋友圈,并扫码加主编好友带你入群,务必备注“姓名-单位-研究方向-职称/年级”。PI请明示身份,另有海内外微生物相关PI群供大佬合作交流。技术问题寻求帮助,首先阅读《如何优雅的提问》

befc13fa90b8941ec1644c86fd5cc5bc.png

学习16S扩增子、宏基因组科研思路和分析实战,关注“宏基因组”

点击阅读原文,跳转最新文章目录阅读

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值