ログイン 概要 よくある質問

strand specific RNA-Seqデータ解析におけるtophatの利用について質問させてください。 現在SOLiD™ Whole Transcriptome Analysis Kitで得られたstrand specific RNA-Seqのデータを解析しております。

tophatとcufflinksを用いて発現量について観ようとしているのですが、 tophatのlibrary-typeをtophatのマニュアル通りに --library-type fr-secondstrand と指定すると その後cuffcompareをかけたときのClass Codesが cが3765 eが122 xが10668 とあまりにもxが多い結果となっています。

--library-type fr-firststrand で行ってみると cが9739 eが6870 xが33 となります。

この結果から本当に --library-type fr-secondstrand をつかっていいものかと思っています。

tophatのマニュアルにをよく読むと fr-firststrandはdUTP, NSR, NNSR fr-secondstrand はLigation, Standard SOLiD に適していると書かれています。

私は一般的なRNA ligationを用いたstrand specificなRNA-Seqはfirst strandを dUTPやNNSRはsecond strandをシーケンスすると思っていたのでが違うのでしょうか? それともtophatのマニュアルがまちがえているのでしょうか。

以上長文となってしまいましたが質問させていただければと思います。

質問日 Apr 30 '11 at 02:16

snk_u's gravatar image

snk_u
1112

edited May 18 '11 at 22:29

mn3's gravatar image

mn3 ♦♦
5154922


直接の回答になるわけではないですが、同様の質問がSEQanswersにもありました。 http://seqanswers.com/forums/showthread.php?t=10810 どうもTopHatの挙動が微妙なようですね。

回答日 Apr 30 '11 at 02:44

gaou's gravatar image

gaou ♦♦
22125

edited May 18 '11 at 22:29

mn3's gravatar image

mn3 ♦♦
5154922

夜分遅くに回答ありがとうございます。

やはりTopHatのssRNA-Seqの解析は少し微妙といったところなのでしょうか...

回答日 Apr 30 '11 at 02:53

snk_u's gravatar image

snk_u
1112

あなたの回答
プレビューをトグルする

この質問をフォローする

By Email:

Once you sign in you will be able to subscribe for any updates here

By RSS:

回答

回答とコメント

タグ:

×47
×6
×2

質問日: Apr 30 '11 at 02:16

閲覧数: 27,424 回

最終更新日: May 18 '11 at 22:29

powered by OSQA