RNA-seq解析（遺伝子発現解析）
1.サービス概要
| 価格 | 40,000円/サンプル ※サンプルが原核生物やSaccharomyces cerevisiaeの場合、rRNA除去が必須です（15,000円/サンプル） ※24サンプル以上で数量割引があります。 |
|---|---|
| データ量 | 2,000万リード/2Gb（サンプルがバクテリア等の場合は、1,000万リード/1Gb） |
| 納期 | 25営業日 ※データ解析不要（生データ納品）の場合は20営業日 |
| 作業内容 | ・送付サンプル（トータルRNA）の品質確認（注1） ・オリゴdTでmRNAを精製、ランダムプライマーでcDNA合成 ・シーケンスライブラリー作製 ・シーケンシング解析（DNBSEQ-G400 100bpペアエンド解析） ・データ解析（注2） hisat2ソフトウェアを用いて参照配列にマッピングいたします。その後、featureCountを用いてカウントを行った後、RPKM/TPM正規化を行います。 ※cufflinksを用いた解析も対応可能です。 |
| 納品物 | 以下のデータをDVDやUSB等の記録媒体に保存して納品します。 ・報告書 ・シーケンス生データ（fastq形式） ・マッピングとカウント結果（エクセルファイル） ・参照配列（fasta形式）と遺伝子アノテーション（gff形式） ・DDBJのデータベース登録に必要なデータ ※データ量が多い場合はHDD納品となり、HDD代として別途10,000円（税別）がかかります。お手持ちのHDDを送付いただき、そちらに保存して納品することも可能です。 ＜HDD納品となる目安＞生データ納品の場合は200Gb、データ解析ありの場合は100Gb |
| ご提供サンプル | トータルRNA 1 μg以上（＞50ng/ul、＞20ul） ※滅菌水に溶解した状態で冷凍便を使用して送付ください。 ※1 μg以上のトータルRNAがご用意できない場合は、ご相談ください。 |
| ご提供情報 | 1.サンプル情報提供用紙→サンプル送付時に同封ください 2.オーダーシート→dna@gikenbio.comへ送信ください 3. データ解析依頼書→ 〃 レファレンス配列(fasta形式)/遺伝子アノテーション情報（gff形式, gtf形式）はCD-RやUSBに入れて送付頂くか、ダウンロードできるURLをご連絡下さい。 De novo transcriptome assemblyの解析をご希望のお客様は、お送り頂かなくても結構です。 |
（注1）rRNAの分解が進んでおり、シーケンシング解析を行ってもデータが得られないと判断した場合は、中止させていただきます（この場合のご請求は行いません）。
（注2）De novo transcriptome assembly 解析の場合、解析方法が通常の解析とは異なります。
2.オプション
| 項目 | 備考 | 価格（税別） | 追加 納期 |
|
|---|---|---|---|---|
| RNA抽出 | QIAGEN社のRNeasy等を使用して、トータルRNAを抽出します。 ＜ご提供サンプル量の目安＞ 細菌：細胞数 1×10^9程度、湿重量 20 mg程度 酵母：細胞数 5×10^7程度、湿重量 25 mg程度 植物（葉）：湿重量 50 mg程度、直径1.5cm 動物細胞：細胞数 1×10^7程度、湿重量 20 mg程度 動物組織：湿重量 20 mg程度 |
20,000円/サンプル | 1週間 | |
| rRNA除去 | siTOOLs Biotech社のriboPOOLを使用して、rRNAを除去します。 | 15,000円/サンプル | なし | |
| デ ータ解析 |
発現比較解析 | <参照ゲノム配列あり> iDEGES正規化（Sun et al., BMC Bioinformatics, 2013）後、edgeR/DESeqを用いて発現変動遺伝子を検出致します。 |
1パターン50,000円 ※追加1パターンごと＋10,000円 |
なし |
| <参照ゲノム配列なし> De novo transcriptome assemblyで得られたcontigに対しRSEMを用いてマッピングを行います。その後、edgeR/DESeqを用いて発現変動遺伝子を検出致します。 |
||||
| Volcano プロット | 発現比較解析結果の発現比とp値で作図します。 | 1パターン10,000円 ※追加1パターンごと＋2,000円 |
なし | |
| 主成分分析（PCA) | 遺伝子発現プロファイルの類似性を座標分布で示します。 | 10,000円 | なし | |
| クラスタリング解析 | 遺伝子発現プロファイルの類似性を階層クラスタリングで示します。 | 10,000円 | なし | |
| 相関図作成 | 遺伝子発現プロファイルの類似性を相関図で示します。 | 10,000円 | なし | |
| ヒートマップ解析 | 遺伝子発現プロファイルを色で表現した可視化グラフを提供します。 | 20,000円 | なし | |
| パスウェイ解析 | KEGG（Kyoto Encyclopedia of Genes and Genomes）を利用したパスウェイ解析です。発現変動解析の結果から同定された遺伝子が、どのような代謝系に関わっているのかを可視化させます。 | 1パターン70,000円 ※追加1パターンごと＋10,000円 |
なし | |
| Gene Ontology (GO) 解析 | 全遺伝子のGO termと発現比較解析の結果から同定された発現変動遺伝子のGO termについて、出現頻度の比較とFisher検定による有意差検定を行います。発現変動遺伝子の特徴的なGO termを推定することができます。 | 1パターン70,000円 ※追加1パターンごと＋10,000円 |
なし | |
| →パスウェイ解析＆GO解析 セット割引 -40,000円 | ||||
| De novo transcriptome assembly | 参照配列がない非モデル生物の遺伝子カタログ作成を、Trinityを用いて行います。 | 100,000円 | なし | |
3.RNA-seq解析（発現解析）の具体例
標準解析
遺伝子発現行列データをエクセルデータで納品します。>>詳しいデータをご覧いただけます
サンプル毎にどの遺伝子に何リードマッピングされたかを示しています。また、RPKM/TPM正規化後の数値も示しています。
発現変動解析（オプション）
1.遺伝子発現比較解析の結果をエクセルデータで納品します。>>詳しいデータをご覧いただけます
===データの見方===
-sample_name: 正規化前のリードカウント
– a.value: log2(平均発現量)
– m.value: log2(発現比)
– p.value: p値
– q.value: q値 (false discovery rate(FDR))
– rank: q値が低い順にランク付けされています。
– estmatedDEG: 「1」は発現変動遺伝子(q<0.05)、「0」は発現変動遺伝子ではないこと(q>0.05)を示します。
2.MAプロットの図をpng形式で納品します。
横軸がlog2(平均発現量(A))、縦軸がlog2(発現比(M))を示しています。
Volcano プロット
Volcano プロットをpng形式で納品します。
Volcano プロット は、MAプロットの親戚の図です。volcano plotは横軸がlog2(発現比)、縦軸がp値を示しています。
主成分分析（オプション）
1.相関行列をエクセルデータで納品します。>>詳しいデータをご覧いただけます
2.主成分分析の結果をpng形式で納品します。
クラスタリング解析（オプション）
クラスタリングの結果をpng形式で納品します。
ヒートマップ解析（オプション）
ヒートマップの図をpng形式で納品します。
パスウェイ解析（オプション）
パスウェイ解析の結果をテキストデータで納品します。エクセルで開くことができます。>>詳しいデータをご覧いただけます
Gene Ontology（GO）解析
Gene Ontology（GO）解析の結果をエクセルデータで納品します。
4.よくある質問
Q1.推奨のRNA抽出キットは？
特にございません。弊社ではQIAGEN社のRNeasyシリーズを使用することが多いです。
Q2.送付の際、バッファーは何でもよい？
滅菌水のほか、Tris や TE 、抽出キットの溶出バッファーでも構いません。但し、塩濃度が高いバッファーの場合、QCの電気泳動がうまくいかない可能性がありますので、避けてください。
Q3.RIN値の目安は？
バイオアナライザーのRIN値で7以上が目安です。送付前にお客様側でも品質確認される際は、参考にしてください。
Q4.発現比較解析などのパターン数って？
例えば、A群・B群・C群があり、A群 VS B群・A群 VS C群の発現比較解析をする場合は2パターン、総当たりで行う場合は3パターンとなります。3群間比較（A群 VS B群 VS C群）も可能ですが、発現比較の発現平均値（A値）と発現比（M値）が算出されないため、MAプロット図の納品はなく、エクセルデータのみの納品となります。
Q5.DNBSEQ-G400の特徴は？
なんといってもコストパフォーマンスがよいことです。低価格で高品質のショートリードが取得できます。
Q6.DNBSEQ-G400の出力データはどんな形式？
出力はfastqファイルで、Q値は最新のイルミナ社製シーケンサーと同様にsanger形式です。
Q7.DNBSEQ-G400の生データのクオリティは？
DNBSEQ-G400(MGI社）とHiSeq（イルミナ社）を用いてシーケンシング解析を行い、取得したリードのQ20およびQ30を比較したところ、2機種のそれらは同等でした。
5.ちらしPDF
6.お問合せ・サンプル送付先
お見積り・お問合せ
お見積りフォーム、お問合せフォーム、メール（dna@gikenbio.com）、または お電話（042-780-8333）でお気軽にご相談ください。
サンプル送付先
移転のため、送付先が変更になりました。ご注意ください。
株式会社生物技研
〒252-0154 神奈川県相模原市緑区長竹657 / TEL 042-780-8333
※平日着でお送りください。発送は元払いでお願いしています。
＜添付書類＞
・サンプル情報提供用紙→サンプル送付時に同封ください。
・オーダーシート→dna@gikenbio.comへ送信ください
・データ解析依頼書→ 〃
