i*****i
发帖数: 154 | 1 马上要开始做了，还有一些疑问，比较初级，望大侠赐教。
我问了Genomic Core，他们说HiSeq2000 的output是180M reads，
我觉得他们说的是单向，pair End 应该是360M， 是不是这样？
另外有些文献讨论多少reads是合适的，目前比较一致的结论好像是30M reads对于人类
基因组就差不太多了。我想知道这里的30M是指单向的，还是pair end的？
对于pair end的数据不了解，希望大家帮助。 |
j*p
发帖数: 411 | 2 1. 180M估计是HiSeq2000 pair-end。pair-end 会产生两个文件，每一个文件大概有
100M reads，每一对read都互相对应。有时候能产生更多的数据，我反正是没见过有
380M这么多的。
2. 如果只想知道普通gene表达，30M reads应该是足够多了。那些做HiSeq的多数是想
要知道例如LncRNA，eRNA之类低表达的基因（需要ChIPseq），或者是想看看有没有不
同alternative splicing， circular RNA等。一般讲数据多，read长，分析后得出的
结果越靠谱。
3. 不同的library construction也会有影响，总的来说，polyA selection会比Ribo 0
，数据看上去要平滑，缺点就是不能够detect没有polyA的gene（废话）。
4. 具体问题具体对待吧。 |
i*****i
发帖数: 154 | 3 谢谢JCP的帮忙。
因为我们的RNA质量不是很高，担心ployA筛选丢掉的东西比较多，目前倾向于用ribo-
zero。 |
s******y
发帖数: 137 | 4 我们最近的hiseq一般都稳定在200M pairs 左右，就是有400M的single end reads.
30M 一般指的是paired end的。不过30M低表达的基因就很难了。
个人认为biological replicates比reads数目更重要，我情愿要30M有3个重复也不要
100M没有重复的结果来做数据分析。
RNA质量你的RIN score是多少，很多地方小于7或8就不收了。å |
