快捷搜索:  汽车  科技

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)下面以GO富集弦图为例来解读其所展示的含义:GO plot可以得到包括气泡图、弦图(chord plot)在内的高级GO富集可视化图形,科学家只需几行代码就可以生成有洞察力的图形!柱形图或点图只能突出展示富集率而不能展示蛋白和通路的关系以及指示通路功能强弱的变化。如何对GO分析结果进行全面而详细的可视化仍然是一种挑战!2015 年Wencke Walter等【1】在Bioinformatics发表了一个基于ggplot2的R包——GO plot(http://wencke.github.io/),文章发表以来引用次数高达127次。

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(1)

解螺旋公众号·陪伴你科研的第2152天

Cell阶段总结:合成蛋白水解靶向嵌合体的发展历程

基因本体(Gene ontology,GO)分析在生物信息学领域中广泛使用,包括细胞组分(Cellular component)、分子功能(Molecular function)与生物过程(Biological process)。

文献中通常使用柱形图或点图来展示GO富集分析结果:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(2)

柱形图或点图只能突出展示富集率而不能展示蛋白和通路的关系以及指示通路功能强弱的变化。

如何对GO分析结果进行全面而详细的可视化仍然是一种挑战!

2015 年Wencke Walter等【1】在Bioinformatics发表了一个基于ggplot2的R包——GO plot(http://wencke.github.io/),文章发表以来引用次数高达127次。

GO plot可以得到包括气泡图、弦图(chord plot)在内的高级GO富集可视化图形,科学家只需几行代码就可以生成有洞察力的图形!

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(3)

下面以GO富集弦图为例来解读其所展示的含义:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(4)

GO富集弦图即显示了参与GO term的基因有哪些;

还按照变化倍数对上调/下调蛋白进行了排序;

还能够显示到底是哪些蛋白参与了何种GO term;

并将GO富集分析的term按照强弱进行了排序。

GO富集弦图可以反应全面详细的数据信息,提高文章图片档次!

那么问题来了,不会R语言的小白该如何进行GO分析?如何绘制GO富集弦图呢?

今天半夏就给大家分享无代码GO分析并作弦图。

1

Metascape进行富集分析

Metascape网址:http://metascape.org/gp/

其界面如下:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(5)

Metascape可提交单个或多个基因序列,可提交xlsx、csv与txt格式文件,也可以直接进行粘贴。我们选择直接粘贴基因序列:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(6)

待分析完成后打开Analysis reporte page,此报告页面可在72小时内访问:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(7)

可以得到的结果:

(1)富集条形图(Bar graph of enriched terms)

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(8)

图片可以保存也可以点击PDF图标可保存为pdf文件。

(2)输入基因列表的统计数据

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(9)

(3)Enriched term网络

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(10)

可导出Cytoscape(官网http://www.cytoscape.org/)可以编辑的CYS文件、PDF文件与在线打开,点击Web在线打开网络:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(11)

GO network显示模式有三种:Color by cluster,即根据聚类的簇显示颜色、Color by p value,即颜色由p值决定与Default默认单一颜色:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(12)

Layout的模式有3种:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(13)

鼠标点击节点可以显示节点信息:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(14)

鼠标选中节点可拖动节点至任意位置:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(15)

Enriched term网络右侧工具分别为:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(16)

选中1个或多个节点,点击Show/Hide Table可显示节点信息:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(17)

从表格中可以得到很多信息,例如Fgg、Kng2与Fga参与了GO Biological process中的negative regulation of blood coagulation,其LogP约为-3.97:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(18)

(4)蛋白质相互作用富集分析

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(19)

对于每个给定的基因列表,利用BioGrid6数据库进行了蛋白质-蛋白质相互作用富集分析。

得到的网络包含与列表中至少一个其他成员形成物理相互作用的蛋白质子集。

如果网络包含3到500个蛋白质,则应用MCODE7算法来识别紧密连接的网络。

导出分析结果:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(20)

选择All in One Zip File:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(21)

打开metascape_result Excel文件可得到基因注释与富集情况,可得到哪些基因参与了何种生物过程,还能根据p值对富集生物过程强弱的变化进行排序:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(22)

2

无代码GO富集弦图绘制

工具是微生信网站,步骤:

1. 打开GO富集弦图绘图模块:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(23)

2.按照右侧绘图即数据格式说明准备数据:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(24)

metascape_result文件中找到GO富集注释信息:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(25)

调整数据为右侧要求格式,即第一列是基因,最后一列是logFC,其余列是GO term,1表示该基因在这个GO term中,0表示不在:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(26)

输入数据,设置图片宽、高与文字大小,排序选择logFC:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(27)

点击提交得到炫目多彩的GO富集弦图:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(28)

在生成的图上点击鼠标右键,在弹出的窗口中选择另存为(直接是svg格式文件,文件后缀不用改),若不行请换浏览器,推荐chrome。

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(29)

如何修改图片,以满足杂志需求?

默认Arial字体,可下载Inkscape软件,或者Acrobat IIustrator (AI)、PS软件,打开下载的svg文件,修改文字,布局等后,导出高分辨率出版级矢量图,例如pdf。

3

扩展

Cytoscape进行富集分析

Cytoscape官网:http://www.cytoscape.org/。

下载安装Cytoscape后Apps -> App manager安装BiNGO插件。

Cytoscape -> Apps -> BiNGO打开插件:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(30)

输入数据,设置参数:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(31)

得到分析结果:

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(32)

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(33)

4

其他富集分析

1. David (https://david.ncifcrf.gov/home.jsp)

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(34)

2. Gene Ontology Consortium (http://geneontology.org/)

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(35)

3. Webgenstal (http://www.webgestalt.org/)

go富集分析结果怎么可视化(无代码绘制生信美图之GO富集弦图)(36)

弦图除了可以用来绘制GO富集分析结果还可以用来绘制KEGG pathway结果。

今天给大家分享无代码GO分析并作炫目多彩的弦图就到此为止了,希望对大家有所帮助!

参考文献:

  1. Walter W Sánchez-Cabo F Ricote M. GOplot: an R package for visually combining expression data with functional analysis. Bioinformatics. 2015;31(17):2912–2914. doi:10.1093/bioinformatics/btv300
  2. Zhou Y Zhou B Pache L et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523.

10天领悟3分SCI

点击报名

解螺旋-中国临床医生科研成长平台

猜您喜欢: