刘小泽写于19.1.25
最近要开始一批无参转录组学习,我喜欢在学东西时先学一些知识储备,现在不太喜欢拿来代码直接去跑,可能对背后的原理性东西想要做更多了解了。比如这个结果是怎么得到的?怎么看懂结果?我怎么才可以根据自己的需要获得部分结果。这些都需要很长的时间去积累,而流程,只要调试好就可以去跑,但即使不报错的结果也未必是对的,这一点需要注意
先上文章
综述类:
- A survey of best practices for RNA-seq data analysis
https://genomebiology.biomedcentral.com/articles/10.1186/s13059-016-0881-8 - De novo transcript sequence reconstruction from RNA-Seq: reference generation and analysis with Trinity
https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3875132/ - RNAseq data analysis from data carpentry
- Great introduction of RNA-seq from sample preparation to data analysis
RNA-seqlopedia - RNA-seq differential expression & pathway analysis with Sailfish, DESeq2, GAGE, and Pathview
细节类:
关于样本数量:How many biological replicates are needed in an RNA-seq
experiment and which differential expression tool should
you use?
https://www.ncbi.nlm.nih.gov/pubmed/27022035后续标准化、定量、差异分析:
Genome-wide repressive capacity of promoter DNA methylation is revealed through epigenomic manipulation不比对直接定量:sleuth 、Kallisto 、Salmon: the sucessor of Salfish ~ Accurate, Versatile and Ultrafast Quantification from RNA-seq Data using Lightweight-Alignment.
转录组分析前期
这部分内容是测序公司经常说的内容,主要考虑以下三点
生物学重复的数量:至少三个以上的生物学重复
-
文库类型:
- mRNA片段如何获取?
第一种:真核生物成熟的mRNA一般带有polyA尾巴,一般方式就是对具有PolyA尾巴的片段进行捕获。这种方法对mRNA的完整性要求比较高,如果样本发生了降解,那么就会丢失一定的转录组信息
第二种:去除Total RNA中占比很高(超过90%)的rRNA,剩下的RNA中就有mRNA(约1-2%)。样本降解的影响比第一种要小,但需要更高的测序数据量因此成本会高【同时由于原核的mRNA没有polyA,所以只能用去除rRNA的方式】 - 是否需要链特异性文库?
因为现在虽然是叫RNA-seq,但其实还是需要RNA反转录成cDNA然后扩增进行测序,最后是测得还是ATCG。
第一种:利用普通文库。这种方法会同时测到模板链以及反向互补链的信息,不能判断原来的mRNA是什么方向(正链还是负链),既然确定链的来源有偏差,那么也对转录本的定量产生了影响
第二种:利用链特异性文库。这种方法可以将建库过程产生的反向互补序列的文库直接去除掉,可以确定mRNA的来源链,转录本定量的准确性会更高
- mRNA片段如何获取?
-
测序深度:
主要考虑目标转录组的复杂度。低丰度的转录本的定量需要更高的测序深度,但是过高的测序深度也会带来转录本噪声。一般的测序结果中利用饱和曲线可以较好的评估测序深度是否合适
先研究研究流程、文章,最主要用的还是基于比对进行定量的trinity以及不基于比对进行定量的sleuth、kallisto,下一步开始介绍Trinity的使用,这可是一个非常大的套件