校准表达量测序差异表达序列文库转录
标题: miRNA差异表达分析时reads数目的校准(TPM)问题
摘要: [miRNA差异表达分析时reads数目的校准(TPM)问题] 在做miRNA的差异表达分析时,很多文章里都要使用TPM(Transcripts per million,公式为:单一miRNA reads数×10^6 总reads数)来校准miRNA的表达量,这样做是为了去除pcr扩增文库时的偏倚吗?然后,如果使用TPM进行校准,那也会产生另外的差异啊我觉得, 关键词:[校准 表达量 测序 差异表达 序列 文库 转录]……
关键词:校准表达量测序差异表达序列文库转录
在做miRNA的差异表达分析时,很多文章里都要使用TPM(Transcripts per million,公式为:单一miRNA reads数×10^6/总reads数)来校准miRNA的表达量,这样做是为了去除pcr扩增文库时的偏倚吗?
然后,如果使用TPM进行校准,那也会产生另外的差异啊我觉得,举个极端的例子:样品A中实际总共表达了10个reads,其中miR1表达了2个;样品B总共表达了100个reads,其中miR1表达了20个,——实际上在两个样品中miR1是显著差异表达的,但经TPM校准后,两个样品中其表达量就都成了20 000,就一样了么,那么校准岂不是产生了原来没有的错误??
我觉得miRNA是短序列,不存在转录组中序列片段长短的问题,那么直接使用reads的count作为其相应miRNA的表达量岂不是更好?
欢迎讨论!
回复
使用TPM(Transcripts per million,公式为:单一miRNA reads数×10^6/总reads数)来校准miRNA的表达量
这其实是个normalization的过程,因为假如你有多个样品(》=2),测序得到的total mapping reads的数目是不一样的,大家必须normalize到同一个水平才方便比较。样品A中实际总共表达了10个reads,其中miR1表达了2个;样品B总共表达了100个reads,其中miR1表达了20个
你说的实际应该是指细胞内真正表达的数量吧,但是实际上通过测序你是拿不到真正意义上的绝对表达量,只是一个相对的表达量。我们在进行2个或者多个样品genes或者miRNA表达差异分析比较的时候,是默认总体RNA表达的一样多的。简单来说就是你例子中sample A和B的RNA总量是一致的,不过A测序测到了10个reads,B测序测的比较多,拿到了100个reads,如果A也达到B的测序深度,其实miR1也就是20个reads,表达量是一致的。
不知道你明白否
回复
哦,明白了,“拿去上机测序的样品总量是一样的”,所以后边才把测得的总量校准到一样!
谢谢你的耐心讲解,终于解了心里纠结很久的困惑!:hand:
相关热词:校准表达量测序差异表达序列文库转录
相关文章
miRNA差异表达分析时reads数目的校准(TPM)问题
1.pH测量注意事项
2.电子天平的使用规程
3.抗氧化的DPPH吸收值出现负的是怎么一回事
4.关于SEM
5.紫外分光光度仪--光度测定
6.发酵罐处理问题
7.有人用流式细胞仪对细菌进行绝对计数的么??
8.求助关于磷脂脂肪酸C19内标和混标的问题
9.恒温培养箱的问题
10.台式实验室pH 型号HANNA-pH211 怎么校准?
11.植物叶片电导率仪使用
12.qPCR出了问题,做实验实在让人不淡定啊
13.关于nanodrop使用 测定相当不准啊
14.MATLAB软件英文词汇
15.QRT-PCR算法
16.用酶标仪测蛋白浓度
17.TBE缓冲液配置问题
18.重组蛋白胞内表达,各位罐上OD最高能到多少?
19.毕赤酵母表达培养基
20.毕赤酵母表达问题
21.iptg诱导大肠表达,iptg加入的时间是什么时候?
22.毕赤酵母摇瓶成絮状
23.做荧光定量PCR的细菌怎么保存啊,要保存时间久一些,又不会影响实验结果,谢谢~
24.关于大肠杆菌(Escherichia coli)表达目标蛋白问题
25.拮抗基因的表达量与拮抗物质的产量
上一篇:高通量测序 序列分析
下一篇:电击转化的问题
校准表达量测序差异表达序列文库转录
标题: miRNA差异表达分析时reads数目的校准(TPM)问题
摘要: [miRNA差异表达分析时reads数目的校准(TPM)问题] 在做miRNA的差异表达分析时,很多文章里都要使用TPM(Transcripts per million,公式为:单一miRNA reads数×10^6 总reads数)来校准miRNA的表达量,这样做是为了去除pcr扩增文库时的偏倚吗?然后,如果使用TPM进行校准,那也会产生另外的差异啊我觉得, 关键词:[校准 表达量 测序 差异表达 序列 文库 转录]……
关键词:校准表达量测序差异表达序列文库转录
在做miRNA的差异表达分析时,很多文章里都要使用TPM(Transcripts per million,公式为:单一miRNA reads数×10^6/总reads数)来校准miRNA的表达量,这样做是为了去除pcr扩增文库时的偏倚吗?
然后,如果使用TPM进行校准,那也会产生另外的差异啊我觉得,举个极端的例子:样品A中实际总共表达了10个reads,其中miR1表达了2个;样品B总共表达了100个reads,其中miR1表达了20个,——实际上在两个样品中miR1是显著差异表达的,但经TPM校准后,两个样品中其表达量就都成了20 000,就一样了么,那么校准岂不是产生了原来没有的错误??
我觉得miRNA是短序列,不存在转录组中序列片段长短的问题,那么直接使用reads的count作为其相应miRNA的表达量岂不是更好?
欢迎讨论!
回复
使用TPM(Transcripts per million,公式为:单一miRNA reads数×10^6/总reads数)来校准miRNA的表达量
这其实是个normalization的过程,因为假如你有多个样品(》=2),测序得到的total mapping reads的数目是不一样的,大家必须normalize到同一个水平才方便比较。样品A中实际总共表达了10个reads,其中miR1表达了2个;样品B总共表达了100个reads,其中miR1表达了20个
你说的实际应该是指细胞内真正表达的数量吧,但是实际上通过测序你是拿不到真正意义上的绝对表达量,只是一个相对的表达量。我们在进行2个或者多个样品genes或者miRNA表达差异分析比较的时候,是默认总体RNA表达的一样多的。简单来说就是你例子中sample A和B的RNA总量是一致的,不过A测序测到了10个reads,B测序测的比较多,拿到了100个reads,如果A也达到B的测序深度,其实miR1也就是20个reads,表达量是一致的。
不知道你明白否
回复
哦,明白了,“拿去上机测序的样品总量是一样的”,所以后边才把测得的总量校准到一样!
谢谢你的耐心讲解,终于解了心里纠结很久的困惑!:hand:
相关热词:校准表达量测序差异表达序列文库转录
相关文章
miRNA差异表达分析时reads数目的校准(TPM)问题
1.pH测量注意事项
2.电子天平的使用规程
3.抗氧化的DPPH吸收值出现负的是怎么一回事
4.关于SEM
5.紫外分光光度仪--光度测定
6.发酵罐处理问题
7.有人用流式细胞仪对细菌进行绝对计数的么??
8.求助关于磷脂脂肪酸C19内标和混标的问题
9.恒温培养箱的问题
10.台式实验室pH 型号HANNA-pH211 怎么校准?
11.植物叶片电导率仪使用
12.qPCR出了问题,做实验实在让人不淡定啊
13.关于nanodrop使用 测定相当不准啊
14.MATLAB软件英文词汇
15.QRT-PCR算法
16.用酶标仪测蛋白浓度
17.TBE缓冲液配置问题
18.重组蛋白胞内表达,各位罐上OD最高能到多少?
19.毕赤酵母表达培养基
20.毕赤酵母表达问题
21.iptg诱导大肠表达,iptg加入的时间是什么时候?
22.毕赤酵母摇瓶成絮状
23.做荧光定量PCR的细菌怎么保存啊,要保存时间久一些,又不会影响实验结果,谢谢~
24.关于大肠杆菌(Escherichia coli)表达目标蛋白问题
25.拮抗基因的表达量与拮抗物质的产量
上一篇:高通量测序 序列分析
下一篇:电击转化的问题