Loading [MathJax]/jax/input/TeX/config.js
前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >专栏 >任何新手小白都可以使用这个差异基因分析神器:DESeq2 | 零代码,无需编程

任何新手小白都可以使用这个差异基因分析神器:DESeq2 | 零代码,无需编程

作者头像
简说基因
发布于 2025-03-28 08:33:12
发布于 2025-03-28 08:33:12
1970
举报
文章被收录于专栏:简说基因简说基因

差异基因分析软件,有生物学重复的情况下,DESeq2 大概是首选。

本文旨在带领小伙伴们进行零代码差异基因分析,任何人都可以简单上手。

闲话少叙,我们先从上传数据开始。

上传数据

新建历史(可选)

  • • 新建一个历史记录,并命名为:DESeq2

上传数据

  • • 点击左上角的 “上传”按钮
  • • 弹出上传页面:可以上传本地文件,或者粘贴网址自动下载网络文件

以上传本地文件为例:

  • • 全选所有文件,并点击 “打开”
  • • 点击 “开始”
  • • 开始上传后,右侧历史面板中,黄色显示正在上传,绿色显示上传完成

进行分析

进入网站:usegalaxy.cn,在左侧工具搜索框中输入:deseq2,结果中第一个便是。

名词解释

因子:即实验设计,可以简单地理解为实验分组。因子名称,即是分组名称,比如:药物疗效。

因子水平:可简单理解为处理组、未处理组,或者其他分组情况。

DESeq2 可支持多因子,多水平的实验设计。

应用举例

例 1:处理组 vs 未处理组
参数设置
  • • 因子名称:Treatment
  • • 因子水平 1:treated(要注意,感兴趣的分组要作为第 1 个水平,后续差异结果是水平 1 vs 水平 2,所以这个顺序非常重要)
  • • 因子水平 2:untreated
  • • 计数文件:这里以来自于 featureCounts,Stringtie 等软件的原始 Count 值为例,每一个样本对应一个 Count 文件。在计数文件框中选择每一组对应的 count 文件。
  • • 其他参数保持默认,点击 “运行工具”,等待结果
查看结果
  • • 差异表格

可以看到 log2FC,p-value, p-adj 值了。

  • • PCA

可以看到,treated 和 untreated 在 PC1 上很好地分开了,说明处理与否对结果有影响。

  • • 热图

当然 DESeq2 还可以输出其他结果,可以在参数设置界面选择。

例 2:双端测序 vs 单端测序
参数设置(如图)
查看结果
  • • PCA

可以看到,PE 和 SE 在 PC2 上能很好地区分,说明测序方式对结果有影响。

例 3:处理组 vs 未处理组 + 双端测序 vs 单端测序

同时考虑处理与测序两种情况。参数设置:设置完一个因子后,可点击 “Insert Factor” 增加一个因子。

查看结果
  • • PCA

可以看到,处理方式在 PC1 方向上能很好地区分,测序技术在 PC2 上能很好地区分。说明处理与否对结果的影响最大,测序技术的影响排第二。

好了,DESeq2 进行基因水平的差异分析就介绍完了。有不明白的小伙伴可以到群里,或者知识星球提问。

本文参与 腾讯云自媒体同步曝光计划,分享自微信公众号。
原始发表:2025-03-27,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 简说基因 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
暂无评论
推荐阅读
编辑精选文章
换一批
哈佛课程笔记-理解DEseq2差异分析原理
纵轴是基因,横轴是count数,这个图揭示了某个样本中所有基因的count数分布。
生信技能树
2022/03/03
7.6K0
哈佛课程笔记-理解DEseq2差异分析原理
RNA-seq入门实战(四):差异分析前的准备——数据检查
连续两次求贤令:曾经我给你带来了十万用户,但现在祝你倒闭,以及 生信技能树知识整理实习生招募,让我走大运结识了几位优秀小伙伴!大家开始根据我的ngs组学视频进行一系列公共数据集分析实战,其中几个小伙伴让我非常惊喜,不需要怎么沟通和指导,就默默的完成了一个实战!
生信技能树
2022/07/26
2.6K1
RNA-seq入门实战(四):差异分析前的准备——数据检查
RNA-seq入门实战(五):差异分析——DESeq2 edgeR limma的使用与比较
连续两次求贤令:曾经我给你带来了十万用户,但现在祝你倒闭,以及 生信技能树知识整理实习生招募,让我走大运结识了几位优秀小伙伴!大家开始根据我的ngs组学视频进行一系列公共数据集分析实战,其中几个小伙伴让我非常惊喜,不需要怎么沟通和指导,就默默的完成了一个实战!
生信技能树
2022/07/26
18.2K0
RNA-seq入门实战(五):差异分析——DESeq2 edgeR limma的使用与比较
DESeq2差异表达分析
在前文scRNA-seq marker identification(二),我们我们提到了差异分析,下面我们来详细了解下
生信技能树jimmy
2020/12/24
6.1K0
表达差异基因分析
1安装BiocManage,再安装DESeq2包 > # <差异基因分析> > # 1.判断是否有BiocManager包,若不存在则安装 > options(repos=structure(c(CRAN="https://mirrors.tuna.tsinghua.edu.cn/CRAN/"))) #设置清华镜像,加速下载 > > if (!requireNamespace("BiocManager", quietly = TRUE)) + install.packages("BiocManager
爱学习的小明明
2020/09/22
1.4K0
DESeq2差异基因分析和批次效应移除
差异基因鉴定 基因表达标准化 不同样品的测序量会有差异,最简单的标准化方式是计算 counts per million (CPM),即原始reads count除以总reads数乘以1,000,000。 这种计算方式的缺点是容易受到极高表达且在不同样品中存在差异表达的基因的影响;这些基因的打开或关闭会影响到细胞中总的分子数目,可能导致这些基因标准化之后就不存在表达差异了,而原本没有差异的基因标准化之后却有差异了。 RPKM、FPKM和TPM是CPM按照基因或转录本长度归一化后的表达,也会受到这一影响。 ca
生信宝典
2018/06/26
6.8K0
pseudobulks单细胞差异基因分析
之前绘制过FindMarkers/FindAllmarkers差异分析后的单细胞差异基因火山图,除了FindMarkers/FindAllmarkers这种方法以外,pseudobulks是另一种单细胞差异基因分析的方法,这次就来学习和整理一下。
凑齐六个字吧
2024/08/15
4810
pseudobulks单细胞差异基因分析
基因芯片数据分析(六):DESeq2包的基本原理
DESeq2是另外一个分析差异基因的R包,它的功能很多,使用也比较复杂。我们在前面提到过,RPKM,FPKM与TPM是常用的用于均一化不同的样本reads数的方法,不过DESeq2和edgeR并不使用前面的三种方法,因为在对文库进行均一化时,存在两个问题,如下所示:
DoubleHelix
2019/12/13
14K1
你确定你的差异基因找对了吗?
PS:如果你的转录组实验分析报告没有这三张图,就把我们生信技能树的这篇教程甩在他脸上,让他瞧瞧,学习下转录组数据分析。
生信技能树
2019/09/27
1.5K0
你确定你的差异基因找对了吗?
DEApp | RNA-seq差异表达分析工具
在之前介绍[[RNA-seq相关内容介绍]]的视频当中,作者提到了一个用来分析 RNA-seq 差异表达分析的工具。DEApp: https://yanli.shinyapps.io/DEApp/ 。 所以今天就来简单的介绍一下这个差异表达分析的工具。
医学数据库百科
2022/02/08
1.8K0
DEApp | RNA-seq差异表达分析工具
差异分析|DESeq2完成配对样本的差异分析
前段时间拿到一个RNA-seq测序数据(病人的癌和癌旁样本,共5对)及公司做的差异分析结果(1200+差异基因),公司告知用的是配对样本的DESeq分析。
生信补给站
2021/03/03
7.2K6
差异分析|DESeq2完成配对样本的差异分析
(DESeq2) Why are some p values set to NA?
在上一期奇怪的转录组差异表达矩阵之实验分组中,我们谈到DESeq2输出NA的问题,这周我们仍使用上周 GSE126548-分组差异并不大,这个数据集来进行分析
生信菜鸟团
2023/09/09
3.5K0
(DESeq2) Why are some p values set to NA?
跟着存档教程动手学RNAseq分析(五):DESeq2基因水平差异表达分析
跟着存档教程动手学RNAseq分析(四):使用DESeq2进行DE分析的QC方法
王诗翔呀
2022/12/30
2.9K0
跟着存档教程动手学RNAseq分析(五):DESeq2基因水平差异表达分析
对转录组测序的counts矩阵去除批次效应
相信大家听完了我B站的RNA-seq分析流程后,对这个数据的应用方向都不陌生。最近连续收到好几个求助,都是关于转录组测序的counts矩阵去除批次效应,值得写推文解答一下咯!
生信技能树
2020/02/20
6.5K0
RNA-seq 差异分析的细节详解 (5)
本系列[1]将开展全新的转录组分析专栏,主要针对使用DESeq2时可能出现的问题和方法进行展开。
数据科学工厂
2024/12/30
1750
RNA-seq 差异分析的细节详解 (5)
RNA-seq 差异分析的点点滴滴(3)
本系列[1]将开展全新的转录组分析专栏,主要针对使用DESeq2时可能出现的问题和方法进行展开。
数据科学工厂
2024/12/30
1200
RNA-seq 差异分析的点点滴滴(3)
比较不同流程(limma/voom,edgeR,DESeq2 )差异分析的区别
距离第一次听说生信已经十几年了,现在是邋遢大叔重新开始学代码,精力确实已不像从前,各位入坑还是要乘早。后来约莫在5年前,课题组当时有个RNA-Seq数据,lab meeting时听瑞典小哥在汇报DEGs筛选,当时感觉好是神奇。其实陆陆续续也有过学习的念头,但在对自己的各种纵容下,想法又逐渐隐没。直到2月前,机缘巧合参加了生信技能树培训,才进一步强化了自己学习生信技术的信念。
生信技能树
2021/01/05
5.2K0
limma/voom,edgeR,DESeq2分析注意事项,差异分析表达矩阵与分组信息
给粉丝朋友们带来了很多理解上的挑战,所以我们开辟专栏慢慢介绍其中的一些概念性的问题,上一期: 箱线图的生物学含义
生信技能树
2019/08/22
14.2K0
DESeq2差异表达分析(二)
DESeq2工作流程的下一步是QC,它包括样本级和基因级的步骤,对计数数据执行QC检查,以帮助我们确保样本/重复 看起来很好。
生信技能树jimmy
2020/12/24
6.6K0
DESeq2差异分析及VST变换的探索
关于TCGA的差异分析之前介绍过,不过略微有些不够完整,而且主要是演示的TCGAbiolinks这个包,对于DEseq2介绍的不够,所以今天专门说一下使用DEseq2进行差异分析。
医学和生信笔记
2023/02/14
3.6K3
DESeq2差异分析及VST变换的探索
推荐阅读
相关推荐
哈佛课程笔记-理解DEseq2差异分析原理
更多 >
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档