绵羊转录组测序后差异表达基因太多了 - 做专业的生信问答网站

登录

该问题已被锁定！

2: 关注

1915: 浏览

绵羊转录组测序后差异表达基因太多了

转录组差异基因表达分析

一、参考基因组：https://hgdownload.soe.ucsc.edu/goldenPath/oviAri4/bigZips/oviAri4.fa.gz
二、GTF：https://hgdownload.soe.ucsc.edu/goldenPath/oviAri4/bigZips/genes/oviAri4.ncbiRefSeq.gtf.gz
三、hisat2 构建INDEX流程：
1.hisat2_extract_exons.py oviAri4.ncbiRefSeq.gtf > oviAri4.ncbiRefSeq.gtf.exon 2.hisat2_extract_splice_sites.py oviAri4.ncbiRefSeq.gtf > oviAri4.ncbiRefSeq.ss &
3.hisat2-build -p 6 --exon oviAri4.ncbiRefSeq.gtf.exon#提取外显子# --ss oviAri4.ncbiRefSeq.ss#提取剪切位点# oviAri4.fa#参考基因组# oviAri4.fa_ex_ss.gtf > hisat2_build.log 2>&1 &

四、DESeq2得到差异表达基因后，log2foldchange的绝对值> 1 同时padj值小于0.05的筛选条件，共获得了5552个差异表达基因。

想请问一下这么多的差异表达基因可以改筛选条件么？

好问题 0 评论收藏举报

查看全部 2 个回答

孟浩巍超级管理员用户来自于: 北京市
2023-06-14 22:06

你可以画一个MA plot看看，然后再画一个火山图看看这些差异基因的分布。

另外就是，可以增加表达量的筛选，就是ctrl和treatment样品中至少有1个样品是FPKM>1的。或者其他筛选表达量的条件。

赞同 1 0评论

关于作者

: yzx2023 注册会员
这家伙很懒，还没有设置简介

0: 回答

0: 文章

4: 问题

问题动态

发布时间: 2023-06-12 21:56

更新时间: 2023-06-14 22:06

关注人数: 2 人关注

相关问题

二代转录组去除批次效应: 2109 浏览 2 关注 3 回答 0 评论

单细胞转录因子: 1191 浏览 2 关注 1 回答 0 评论

无对照的转录组数据如何寻找差异表达基因: 1593 浏览 2 关注 2 回答 0 评论

如何在差异表达前消除批次效应: 1221 浏览 2 关注 1 回答 0 评论

请问hemmer是否可以使用转录组搜索HMM库？: 1177 浏览 2 关注 1 回答 0 评论

非已知转录因子预测: 1010 浏览 2 关注 1 回答 0 评论

原核生物的转录组测序为什么比真核的贵？: 1687 浏览 2 关注 2 回答 0 评论

如何计算富集分析里的差异倍数: 1234 浏览 2 关注 1 回答 0 评论

转录组: 1636 浏览 2 关注 1 回答 0 评论

麻烦问下，转录组差异分析，我需要筛掉比对率低的bam文件或者外显子比对率的样本吗: 1495 浏览 2 关注 1 回答 0 评论

推荐内容

问用log2(fold_change)数据做热图，遇到想要表达的目的基因log2(fold_change)值为inf或-inf时，怎么办？: 1322 浏览 2 关注 1 回答 0 评论

问请问hisat2+stringtie+deseq2的标准流程: 1907 浏览 2 关注 2 回答 0 评论

问 Bulk-RNAseq多组数据差异表达基因的筛选: 1619 浏览 2 关注 2 回答 0 评论

问进行转录组数据分析时，进行cuffdiff后的输出文件gene_exp.diff中，一个基因出现了两个不同的表达量数据，应该如何处理？: 1015 浏览 2 关注 1 回答 0 评论

问 ciriRNA表达量如何计算？（更）: 1081 浏览 2 关注 1 回答 0 评论

问基因表达: 1156 浏览 2 关注 1 回答 0 评论

问如何利用利用TPM或者FPKM完成DESeq2完成的工作?: 1809 浏览 2 关注 2 回答 0 评论

问普通转录组的跨物种分析: 1602 浏览 1 关注 2 回答 0 评论

问 samtools 有错误: 1126 浏览 2 关注 1 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025