今天给大家介绍一个工具FunRich,FunRich (Functional Enrichment analysis tool) 是一个主要用于基因和蛋白质的功能富集和相互作用网络分析的独立的软件工具,我们前面介绍了TCGA数据库的各种数据下载与整理,获得的表达矩阵可以绘制热图,可以进差异分析,差异分析获得的差异表达基因,可以用于后续的KEGG ,GO等分析,我们这里就介绍FunRich工具,可能大家之前都用DAVID这个数据库,不妨试试FunRich这个工具,功能很强大。
往期相关文章:
TCGA数据库:RNA-Seq数据的下载与处理
TCGA数据库:SNP数据的下载整理及其可视化
TCGA数据库:miRNA数据下载与整理
TCGA数据库:ATAC-Seq数据的下载整理及其可视化
TCGA数据库:GDCRNATools包下载数据、处理数据以及差异分析
一文就会TCGA数据库基因表达差异分析
一文解决TCGA数据库临床数据下载与整理
基因表达谱热图绘制
GraphPad Prism绘制生存曲线和R语言批量绘制生存曲线
一个R脚本解决某类功能基因(比如m6A甲基化)临床预后模型分析流程
肿瘤免疫细胞浸润与临床相关性分析
一.软件下载
网址:http://funrich.org/download
![ed4130fed5c78b48ed37aa8b1e98c0f8.png](https://i-blog.csdnimg.cn/blog_migrate/ef11b78b847356215605d43990a1a911.jpeg)
该软件下载后,无需安装,解压即可使用,是不是很方便。
![cf167fc3915b5389cc7df05e299d16dd.png](https://i-blog.csdnimg.cn/blog_migrate/cc9091ea27a523ce883197fc8e017533.jpeg)
二.软件首页
![4a70e5bbc75194d85e33fa5f41ddac7c.png](https://i-blog.csdnimg.cn/blog_migrate/5cfa6d71a095f3f0e29783b8af622961.jpeg)
我们后面一一介绍。
三.导入数据
点击Add dataset ,粘贴我们的基因。我们也可以导入数据集(基因集),可以是我们差异分析获得的差异表达基因【参考文章:TCGA数据库:GDCRNATools包下载数据、处理数据以及差异分析,一文就会TCGA数据库基因表达差异分析】。这里我们以limma包和edgeR包差异的分析结果为例。
![55af3a86d86cf109094fa036d3b4704c.png](https://i-blog.csdnimg.cn/blog_migrate/39fbe19bca435ad4117c2a530d13c51a.jpeg)
点击Apply,我们可以看到,有些基因在数据中没有。
![cabbd02528b9847a55596c426b7f2630.png](https://i-blog.csdnimg.cn/blog_migrate/34dca301a25935228aa4572ea8ffaec7.jpeg)
点击OK后,如下图,我导入了890个基因,数据集里面的基因名和数据库自带的背景数据集基因名可能会不符,只有761个。
![9db7b0e7638251731f05932b81ee805c.png](https://i-blog.csdnimg.cn/blog_migrate/43095f678ba677d3df712141dbcd045f.jpeg)
点击Manage会显示符合的基因列表,可以复制下来。通过Excel中的VLOOK-UP函数与自己原基因列表比较找到不符合的基因,通过其他数据库,找到另外的基因名,再重新导入,也可以忽略。这里就忽略了。我们还可以点击Remove可以删除掉这个数据集,当然,这里就不移除了。
同样的方式,我们导入limma分析的差异基因,这样就有2个数据集啦。
![2632f128341e16bb1580dcd87a8a4fa8.png](https://i-blog.csdnimg.cn/blog_migrate/650448777f0e77036394ab637e8189a2.jpeg)
四.分析
1.韦恩图绘制在Venn Diagram栏中,点击Venn,选择要分析的基因列表,我们选择limma和edgeR数据集。