【万能代码分享】一键搞定DEseq2差异分析、热图、火山图、GO分析、KEGG分析、GSEA分析等！

DESeq2是一个比较常用的转录组分析R包，包的使用非常简单，与之前的limma包不一样，DESeq2需要的数据是Row counts矩阵，这点非常重要。所以不管你自己的测序数据，还是GEO下载的数据，你需要明确是不是Row counts矩阵。比如我从GEO下载GSE169758的测序数据（我这里只是举例子），如下操作，打开生物医学之家网站（swyxzj.com）进入GEO官网：

检索GSE号，并且明确知道是测序数据，不是芯片数据：

点击下载raw count：

下载、解压，整理成下面（前面三个对照组，后面3个实验组）：

下面我们可以跑代码了。

二、运行代码：

将raw count矩阵和代码放在同一个文件夹下（RNAseq_Deseq2）：

打开代码，一共是7个步骤：

### 一、DEseq2差异分析+标准化PCA图### 二、ENSG转Symbol----公众号：科研部，参考链接：https://mp.weixin.qq.com/s/l3jgf_GhglZ5pwUNhc1x1g### 三、输出差异基因### 四、作图，差异表达图，热图，火山图### 五、差异基因集的GO分析### 六、差异基因集的KEGG分析### 七、 所有基因集的GSEA富集分析

第一步、差异分析+PCA降维：

如上代码，是3个对照组和3个实验组，所以你如果你的数据是7个，你就修改为7就行。之后，你会生成一个PCA图（感觉分组分的不好，这个数据有点问题）：

同时，生成DEseq2差异分析结果文件：

那我们根据PCA图可以看到有个别样本不行，错综复杂。意思就是这个数据不太好。那我们再换一个我自己的数据重新差异分析看看，数据如下：

读入数据，DEseq2差异分析（唯一要改的就是7个对照组，几个实验组）：

PCA看看分组如何：

分的很开，结果可靠。差异分析结果文件如下（但是是ENSG的，我们后续做GO、kegg、gsea都是需要基因symbol，所以下一步，我们要转换）：

二、ENSG转Symbol :

ENSG转symbol，请仔细看这篇推文（点击查看）。我们测序之后，比对到基因组文件的时候是用的哪个版本的GTF文件，你就需要下载该版本的GTF文件，用ENSGmap进去，得到基因Symbol，那我这里是自己测序的数据，用的是gencode.v39.文件，所以就下载1_gencode.v39.annotation.gtf.gz（怎么下载，请点击查看），放在这个文件夹下：

运行代码（你需要自己下载gtf文件，修改8：21的21，这个代码意思就是说取根据8-21列取每行的中位数，后面好比对到gtf文件时，同名的symbol处理，21的意思是7+n，n是样本数，我是14个样本，所以7+14=21）：

就得到了symbol的DEseq2分析结果文件（有symbol了）：

三、输出差异基因：

我这里取的阈值是logFC为1，矫正后的p值小于0.05，也就是变化倍数大于2，并且有统计学意义的基因作为差异基因，用于后续的分析，打开如下：

我们仔细看看后会发现，我们第二列里面有些也是ENSG号。

我们回到GTF文件看看：

发现有些ENSG基因，它的名字就是ENSG号，数量不多，不影响的。

四、作图，差异表达图，热图，火山图：

首先就是做热图，我们不可能将所有的差异基因都展示出来，所以这里筛选出padj < 0.05和abs(log2FoldChange) >2的基因拿出来做热图展示，就是说明有部分基因变化非常明显：

接着我们画火山图，说明很多基因确实发生了变化（蓝色和红色的基因有很多）：

第五步、差异基因集的GO富集分析：

我们要先加载差异基因文件，然后将其转化为entrezID：

我们发现TEC这个基因有多个entrezID：

这样我们后面用entrezID做GO分析的时候就会报错：错误: near "7006": syntax error。如下：

这个时候，我们就可以随机选取该基因的任何一个entrezID都行。

所以，最终的差异基因GO分析代码是：

这里是将P小于0.05的结果保留下来（万能代码，这一步，你什么都不需要改，后面的KEGG分析也是一样，不需要改）：

第六步、差异基因集的KEGG富集分析：

先输入差异基因，随后kegg富集分析，并保存P小于0.05的结果：

再选择P小于0.05的前10作图，并将KEGGmap到的通路下载下来：

产生结果文件：

第七步、所有基因集的GSEA富集分析：

是用所有基因来做GSEA富集分析，不是用差异基因哦，并且按照logFC排序进行GSEA富集分析，运行到下面代码：

P小于0.05的结果将会被保存起来：

然后把自己想要的通路（hsa号）画出来即可：

每一张都保存起来就行了，就等于一篇论文（你的数据有多新，你就能发多大的文章）：

最终的结果文件：

全套代码一共是300多行，自己多练练就会了，不会写不要紧，看得懂就行了，每一句代码都有注释。代码和演示数据我都上传了：

如下免费获取：

“Raw count的Deseq2包差异分析+热图+火山图+GO分析+KEGG分析+GSEA分析全套万能代码（附代码交流群）”

如下免费下载：

🔽①长按下方二维码关注🔽

②对话框输入关键词：全套代码

更多免费资源：

三、科研绘图类：

21 更多资源在更新中.....

8 PS 2021 (Mac)安装包

四、生信和写作类：

辛苦整理，全文无任何广告！

觉得有用的话，您就点个在看、点赞!

http://mp.weixin.qq.com/s?__biz=MzI4NjYwMDQ5Nw==&mid=2247492350&idx=2&sn=2ffa3e38a946102a9b16742ff40a8cfb

科研大亨

请给我发信息，我24小时在线，领取资源请发有效截图哦！

最新文章

2024年度中国博士后科学基金第76批面上资助拟资助人员名单公示（医药类）

手把手教你WB蛋白定量、拉直，用于SCI投稿！

TCPA数据库蛋白质组学数据挖掘发SCI实操课，免费送！

15套GEO/TCGA数据挖掘发SCI实操课免费送！

药物靶点孟德尔随机化分发SCI实战课免费送

手把手教你孟德尔随机化分析发一篇Nature:

8899元NHANES数据库挖掘发SCI实操课免费送！

资源共享VIP会员群，点击进入！

中国顶尖科学家排名发布：没见钟南山、张文宏，但施一公，管轶、李兰娟等学者上榜

【万能代码分享】一键搞定DEseq2差异分析、热图、火山图、GO分析、KEGG分析、GSEA分析等！

8899元NHANES数据库挖掘临床数据发SCI实操课免费送！

手把手教你用ChatGPT分析大量文献，找科研思路：

手把手教你鼠标点点点复现这篇SCI（零代码）：SND1基因在肿瘤中的意义！换一个基因/疾病，模仿操作可发SCI：

药物靶点孟德尔随机化分发SCI实战课免费送

手把手教你孟德尔随机化分析发一篇Nature:

资源共享VIP会员群，点击进入！

2293名博士博士后，获经费资助！

手把手教你组图投SCI：如何将Western Blot (WB) 条带对齐和水平校正？

GraphPad绘制柱形图，自动配色、自动添加统计学P值：

非肿瘤自噬基因纯生信分析发SCI实战课免费送

手把手教你分析双硫死亡发SCI，不做实验：

Graphpad+SPSS统计作SCI图课程免费送！

单细胞+孟德尔随机化分析实操课免费送

资源共享VIP会员群，点击进入！

“拼妈式”学术不端？浙江大学最高层次奖学金本科生获得者与其教授母亲同发论文遭质疑！校方回应

手把手教你组图投SCI：如何将Western Blot (WB) 条带对齐和水平校正？

GraphPad绘制柱形图，自动配色、自动添加统计学P值：

非肿瘤自噬基因纯生信分析发SCI实战课免费送

手把手教你分析双硫死亡发SCI，不做实验：

Graphpad+SPSS统计作SCI图课程免费送！

单细胞+孟德尔随机化分析实操课免费送

资源共享VIP会员群，点击进入！

重磅！2024年度“全球高被引科学家”名单公布（内附完整名单）：

非肿瘤自噬基因纯生信分析发SCI实战课免费送

手把手教你分析双硫死亡发SCI，不做实验：

Graphpad+SPSS统计作SCI图课程免费送！

单细胞+孟德尔随机化分析实操课免费送

【实操】7688元临床预测模型构建&机器学习实操课免费送！！！

8899元NHANES数据库挖掘发SCI实操课免费送！

资源共享VIP会员群，点击进入！

重磅！2024年度“全球高被引科学家”名单公布（内附完整名单）：

非肿瘤自噬基因纯生信分析发SCI实战课免费送

手把手教你分析双硫死亡发SCI，不做实验：

Graphpad+SPSS统计作SCI图课程免费送！

单细胞+孟德尔随机化分析实操课免费送

【实操】7688元临床预测模型构建&机器学习实操课免费送！！！

8899元NHANES数据库挖掘发SCI实操课免费送！

资源共享VIP会员群，点击进入！

非肿瘤自噬基因纯生信分析发SCI实战课免费送

手把手教你孟德尔随机化分析发一篇Nature:

分类

时事

民生

政务

教育

文化

科技

财富

体娱

健康

情感

旅行

百科

职场

楼市

企业

乐活

学术

汽车

时尚

创业

美食

幽默

美体

文摘

原创标签

时事社会财经军事教育体育科技汽车科学房产搞笑综艺明星音乐动漫游戏时尚健康旅游美食生活摄影宠物职场育儿情感小说曲艺文化历史三农文学娱乐电影视频图片新闻宗教电视剧纪录片广告创意壁纸头像心灵鸡汤星座命理教育培训艺术文化金融财经健康医疗美妆时尚餐饮美食母婴育儿社会新闻工业农业时事政治星座占卜幽默笑话独立短篇连载作品文化历史科技互联网

发布位置

广东北京山东江苏河南浙江山西福建河北上海四川陕西湖南安徽湖北内蒙古江西云南广西甘肃辽宁黑龙江贵州新疆重庆吉林天津海南青海宁夏西藏香港澳门台湾美国加拿大澳大利亚日本新加坡英国西班牙新西兰韩国泰国法国德国意大利缅甸菲律宾马来西亚越南荷兰柬埔寨俄罗斯巴西智利卢森堡芬兰瑞典比利时瑞士土耳其斐济挪威朝鲜尼日利亚阿根廷匈牙利爱尔兰印度老挝葡萄牙乌克兰印度尼西亚哈萨克斯坦塔吉克斯坦希腊南非蒙古奥地利肯尼亚加纳丹麦津巴布韦埃及坦桑尼亚捷克阿联酋安哥拉