转录组定量结果为0 - 做专业的生信问答网站

登录

该问题已被锁定！

2: 关注

3304: 浏览

转录组定量结果为0

转录组

我在服务器上进行转录组分析时（拟南芥），使用hisat2做的比对，用featurecounts进行定量，但是发现有一个基因count和TPM值都为0，但是公司分析结果中这个基因有表达，我使用IGV查看了比对的bam文件，发现是有reads mapping，请问这个原因是？谢谢！

好问题 0 评论收藏举报

查看全部 3 个回答

chatGPT机器人机器人用户来自于: IANA
2023-06-02 18:25

可能的原因有以下几点： 1. 序列覆盖不足：虽然该基因的reads可以成功比对到基因组上，但其覆盖度可能太低，无法达到featurecounts统计时的默认阈值，导致未被计数。 2. 基因注释不完整：featurecounts的定量依赖于基因组注释文件，如果该基因没有被完整地注释到基因组上，可能就无法被featurecounts所识别。 3. 转录本可变剪接：如果该基因存在可变剪接情况，featurecounts默认会将所有转录本的reads数目加起来作为该基因的表达量，而公司的分析结果可能是基于某个特定的转录本或剪接形式计算的表达量。 4. 差异表达分析方法不同：公司分析结果中该基因的表达量可能是基于其他的差异表达分析方法计算得出的，与featurecounts所使用的方法有所不同。针对以上可能原因，可以尝试增加转录组测序深度、更新或补充基因组注释文件、使用更加细致的转录本注释信息、或者尝试其他的差异表达分析方法等来解决该问题。

赞同 0 0评论

关于作者

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

0: 回答

0: 文章

4: 问题

问题动态

发布时间: 2023-06-02 18:22

更新时间: 2023-06-06 12:01

关注人数: 2 人关注

相关问题

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

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

转录组数据样本聚类结果不理想: 1667 浏览 2 关注 1 回答 0 评论

请问无参转录组的GO分析如何进行: 2704 浏览 7 关注 6 回答 0 评论

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

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

测了实验组与对照组的mRNA与miRNA的转录组，miRNA-mRNA互作分析网络怎么处理呢？: 1771 浏览 2 关注 2 回答 0 评论

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

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

绵羊转录组测序后差异表达基因太多了: 2284 浏览 2 关注 2 回答 0 评论

推荐内容

问我测了一批寄生蜂的转录组样品，FPKM值感觉好奇怪，合理吗？: 2124 浏览 2 关注 3 回答 0 评论

问转录组测序问题: 1974 浏览 2 关注 1 回答 0 评论

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

问菌侵染植物的转录组分析中，若菌的一个基因在当前时间没有引起植物基因的变化，过了三个小时后才引起植物基因的变化，如何分析这种现象或如何将这一类基因关联起来，有没有类似的文献推荐？: 1368 浏览 2 关注 1 回答 0 评论

问 HISAT2建立索引: 2252 浏览 1 关注 2 回答 0 评论

问如何根据转录组数据得到新转录本？如何验证一个基因的多个转录本？: 2274 浏览 2 关注 2 回答 0 评论

问转录组和代谢组联合分析的方法有那些？: 1940 浏览 3 关注 2 回答 0 评论

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

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

问聚类数据标准化: 1920 浏览 2 关注 2 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025