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

登录

该问题已被锁定！

2: 关注

1743: 浏览

转录组定量结果为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 人关注

相关问题

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

stringtie和gffcompare处理转录本的问题: 699 浏览 2 关注 1 回答 0 评论

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

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

转录本坐标转换成基因组坐标: 1098 浏览 4 关注 2 回答 0 评论

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

有关蛋白表达定量的数据库数据打包下载(蛋白组学): 633 浏览 2 关注 1 回答 0 评论

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

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

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

推荐内容

问对突变组和对照组的某基因表达量做T检验: 686 浏览 2 关注 1 回答 0 评论

问请教一个问题，一个转录因子TF只会调控位于同一条链的基因，还是有可能调控反义链的基因？: 894 浏览 2 关注 2 回答 0 评论

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

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

问 log2(fold_change)可以用来做热图吗？遇到想要表达的目的基因log2(fold_change)值为0时，怎么办？: 700 浏览 2 关注 1 回答 0 评论

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

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

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

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

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

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024