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

登录

该问题已被锁定！

2: 关注

1908: 浏览

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

相关问题

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

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

KEGG结果为空？: 1094 浏览 2 关注 3 回答 0 评论

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

转录组组内样本差异大: 524 浏览 2 关注 1 回答 0 评论

在riboseq分析中，如何对没有起始密码子和终止密码子的转录本进行三碱基准确性分析？: 788 浏览 2 关注 1 回答 0 评论

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

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

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

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

推荐内容

问链特异性文库（mRNA/lncRNA/circRNA）如何将RNA类型分开？: 811 浏览 2 关注 1 回答 0 评论

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

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

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

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

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

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

问 WGCNA: 1272 浏览 2 关注 2 回答 0 评论

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

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

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024