miRNA差异表达分析时reads数目的校准(TPM)问题

校准表达量测序差异表达序列文库转录

标题: 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.拮抗基因的表达量与拮抗物质的产量

上一篇:高通量测序 序列分析

下一篇:电击转化的问题


相关内容

  • 全基因组重测序数据分析
  • 全基因组重测序数据分析 1. 简介(Introduction) 通过高通量测序识别发现de novo的somatic 和germ line 突变,结构变异-SNV ,包括重排突变(deletioin, duplication 以及copy number variation)以及SNP 的座位:针对重 ...

  • 肿瘤侵袭转移相关microRNAs的研究进展
  • ·144·中国肺癌杂志2010年2月第13卷第2期Chi n J Lung Can cer, Februar y 2010, Vol.13, No.2 ·综 述· 肿瘤侵袭转移相关microRNAs 的研究进展 王海兵 综述 陈晓峰 王和勇 张慧君 审校 [中图分类号] R734.2 DOI: 10 ...

  • minRNA检测方法
  • microRNA (miRNA )是一种机体内源性表达的单链小分子RNA ,位于基因组非编码区,本身不具有开放阅读框(open reading frame,ORF ),具有高度保守性,时序性和组织特异性.miRNA 广泛存在于各种真核细胞中,不编码任何蛋白质,长度仅为20-24nt .成熟的miRN ...

  • 转录组测序问题集锦
  • 转录组测序问题集锦 转录组是某个物种或者特定细胞类型产生的所有转录本的集合,转录组测序(RNA-seq) 是最近发展起来的利用深度测序技术进行转录分析的方法,可以对全转录组进行系统的研究. Roche GS FLX Titanium .Illumina Solexa GA IIx和AB SOLID ...

  • 全转录组关联分析和前沿热点分析
  • 全转录组是指细胞所能转录出来的所有类型RNA的总和,包括mRNA和非编码RNA(non-coding RNA,ncRNA).目前针对非编码RNA的研究主要集中在small RNA,lncRNA和circRNA.对特定细胞的特定状态下全转录组的研究,可以从RNA层面系统地揭示生物学背后的转录调控规律. ...

  • miR―126调控SKOV3卵巢癌细胞周期的实验研究
  • [摘要] 目的 探讨miR-126对卵巢癌细胞株的细胞周期和形态的影响,评价miR-126靶向调控卵巢癌增殖侵袭的作用. 方法 体外培养人卵巢癌细胞株SKOV3,通过慢病毒包装miR-126并转染卵巢癌SKOV3细胞系,将其分为三组:转染miR-126表达组(SKOV3/LV3-has-miR-12 ...

  • 一种基于大规模并行全基因组测序的单核苷酸多态性检测方法
  • 一种基于大规模并行全基因组测序的单核苷酸多态性 检测方法 技术领域 本发明涉及一种基于新一代大规模并行全基因组测序的单核苷酸多态性的检测方法,属基因工程技术领域. 背景技术 遗传突变影响着生物体表型的变异,罹患疾病的风险,对药物和环境刺激的反应.全基因组连锁分析和定点克隆技术对受单基因影响的孟德尔遗 ...

  • 关于子宫颈癌医学专业的论文摘要样例
  • 关于子宫颈癌医学专业的论文摘要样例 子宫颈癌是全世界妇女发病率和死亡率居第二位的恶性肿瘤[1],且发病率在发展中国家高于发达国家.依照2012年国际癌症研究机构统计全世界的妇女估计有528,000个新发子宫颈癌病例,其中266,000人死于该病,其中约85 %的新发病例和87 %的死亡人数发生在发展 ...

  • miRNA-156基因家族进化及功能分析
  • 本 科 毕 业 论 文 miRNA-156基因家族进化及功能分析 2015年5月20日 独 创 声 明 本人郑重声明:所呈交的毕业论文(设计) ,是本人在指导老师的指导下,独立进行研究工作所取得的成果,成果不存在知识产权争议.尽我所知,除文中已经注明引用的内容外,本论文(设计)不含任何其他个人或集体 ...