3' 末端配列 A seq2 とライブラリの準備

キャプチャと mRNA の 3' 端のシーケンスは、mRNA プロセシングの研究と遺伝子の発現定量をことができます。真核生物 Mrna がポリ a の尾のためビードを固定化したオリゴ チミジン (cDNA 合成の主なことができます oligo(dT)) 分子と合計セル lysates から効率的に精製することができます。ただし、このアプローチには、2 つ欠点があります。最初、トラン スクリプト内部にある A の伸張できる cDNA 合成、スプリアス多サイトの結果の主要なも。第二に、均質多ストレッチは、トラン スクリプト識別のための情報されていない別のシーケンスの特定課題を提起します。ポリ a を逆転写など、これらの制限を回避するために様々 なアプローチが提案されているツインテール RNase H 消化 (3 P seq ¹)、続いて 20 で終わるカスタム配列のプライマーの使用 Ts (2 P seq ²) の選択基準RNase H 消化 (3' 読み取り³)、そしてヘアピン (A seq ⁴) で 3' アダプターを含むオリゴ dT プライマーの使用によって続いて CU₅T₄₅プライマーで 50 以上のヌクレオチドのポリ a の尾を持つ RNA 断片。

最近開発された A seq2 方法⁵は多と特に発生するアダプターの自己 ligation によって生成されるダイマーの割合を最小限に抑えるために同時にシーケンスをバイパスすることを目的とするときのアダプターのモル濃度挿入濃度を上回る。この問題は、両方のアダプターが、seq2、3' アダプターが逆転写後の RNA のフラグメントと Cdna の 5' 末端に 5' アダプターの 5' 端に結紮がようにポリヌクレオチド終了の同じ型に組み合わされてときに排除できます。方法は当社従来 - seq - シーケンスが、5'-に - 3 年より便利な ' により、正確に方向制御 RNA 断片化-、ポリ a サイトの識別の高精度を維持しながら。典型的なサンプルのシーケンス処理された読み取りの約 80% はゲノムに一意にマップし、以上 20,000 多サイト クラスター、注釈付き 3' UTRs と重複の 70% 以上の識別に 。

簡単に言えば、A seq2 プロトコルは mRNA の断片化と逆の補数 3' 5' RNA のフラグメントの端にアダプターの ligation から始まります。ポリ (A)-含む Rna は、リバース 3' 末端、位置 4 dU および 5' 端、ビオチンでアンカー ヌクレオチドを含む 25 のヌクレオチド (nt) 長いランダムプライマー プライマーと転写磁気ストレプトアビジン ビーズに cDNA のバインディングが可能です。、ビオチンを含むプライマーのほとんどはウラシル DNA グリコシラーゼ (UDG) と DNA グリコシラーゼ リアーゼ エンドヌクレアーゼ VIII を含むユーザー酵素ミックスによるデュで胸の谷間で cDNA から削除されます。胸の谷間に残るポリ a テールの位置をマークするこの反応 5' アダプター、および 3 つの Ts 左の結紮のためそのまま終了を残します。受信者の 5' 端に結紮 5' と 3' の両方のアダプターが添付されているので、アダプターの二量体は生成されません。4 ヌクレオチド ランダム-mers 読み取りの初めに導入された最新のシーケンシングの楽器のクラスター分解能を可能し、検出と PCR 増幅アーチファクト除去のため分子意識 (UMI) としても利用できます。UMI のサイズは、他の研究⁶ではさらに増加することができます。プロトコルすべて無作為化の四量体で、5' 端始まるによって PCR 増幅デジタルアーティファクトの補正で 3 診断 Ts を持っている読み取りの 3 ts 処理に続いて始まる mRNA 3' 端に相補的な逆が読み取りを生成します。行政、3' アダプター配列除去を悪用し、逆に補完。A 豊富なサイトは内部でランダムプライマー プライミングから由来している可能性があります読み取りまた計算識別、破棄されます。偽のサイトは一般的に 18 の特徴の一つを欠いているとする必要があります保存された多信号 〜 21 ヌクレオチド上流明らか胸の谷間サイト⁷の。

プロトコル細胞培養と一晩をカウントされません約 8 h の実践的な時間が必要です。解析ソフトウェアで高精度な多サイト id を関連付けられている読みます。ポリ a サイトからクラスターが作成この注釈付きの遺伝子とこれらの原稿 (コントロール siRNA と si HNRNPC 処理細胞の 2 つの生物学的複製) 84% 重複と 3' UTR、75% の重複、86% どちらかとさらに強調表示 4 のサンプルに基づいて、3' UTR またはターミナル エクソン。3' 端レプリケートのサンプルでの発現のピアソンの相関係数は 0.92 と、0.9 以上の値は、通常、メソッドで取得します。したがって、A seq2、非常に再現性のある結果を与える便利な方法です。

1 ですセル成長および mRNA 分離

1. ~ 80% 合流でよくあたり 10 の ⁶ セル × 1 6 ウェル プレートの実験的設計に応じて細胞を成長。
。
2. は、成長培地を削除し、セルで 1 回洗浄リン酸緩衝生理食塩水。直接 mRNA 分離キットから 1 mL の溶解バッファーを追加することによって、プレート上のセルを溶解させます。転送粘性 1 mL ピペット チップを 15 mL プラスチック チューブに溶解。ゴム製へらを使用してプレート表面から細胞材料を完全にデタッチします
。
3. せんライセート含む粘性 DNA ライセート、もはや粘性までいくつか上下にピストンの動き活発で 23 G 注射針に接続されている 1 mL 注射器。注射針をチューブから溶解液を取り出すように底の中心にポイントします
。
4. は、注射器を使用して 1.5 mL チューブにライセート転送します。20,000 × g で 5 分間、残骸を削除する 4 ° C をスピンします。プロトコルを通して DNA 低いバインド 1.5 mL バイアルを使用します
。
5. 遠心分離機の実行中は、換散バッファーの 500 μ L と磁気ラックに再懸濁オリゴ (dT) ₂₅ 磁気ビーズの 300 μ L を洗ってください。ラックに 2-3 回チューブをミックスします。ソリューションが明確だった後、バッファーを削除します。手順 1.4 から明確な上澄みを収集し、ビーズに追加します。再懸濁し、ホイールの回転で 10 分間のチューブを配置
6. は、磁気ラックにチューブを置きます。2 分追加 0.8 mL バッファー A 後 mRNA 分離キットから透明な液体を削除します。2-3 回、ラックに 180 ° 度の管の電源を入れます。この洗濯、手順 A. のバッファーをもう一度
7. 手順 1.6 0.8 mL のバッファー B にビードを 2 回洗浄しなさい
。 ビーズ、バインドされた mRNA を溶出する
8. は、33 μ H ₂ O を追加し、ビードを再停止しなさい。5 分加熱ブロック 75 ° C まで加熱します。チューブ、1 の s および場所をすぐにスピン磁気ラックにそれら。上清を新しいチューブに転送します。さらに使用するまで-80 ° c のサンプルを格納できます
。 アルカリ加水分解、33 バッファー追加 66 μ
9. μ L mRNA (ステップ 1.8) 混合し、加熱ブロック 95 ° C で 5 分間まったく熱します。すぐに氷のチューブをチル
10. 分離 RNA、RNA のクリーンアップ キットレンズ
    。 メモ: ボリュームを確認します。それは、100 μ L をする必要があります。
    1. 追加 350 μ L RLT バッファー キット 250 から μ L エタノール。列と 30 のスピンに負荷室温 (RT) で 8,000 の x g で s。キットから 500 μ L RPE バッファーで洗浄します。500 μ L 80% エタノールで洗浄します。列を乾燥する 20,000 × g で 5 分間スピンします。36 μ L H ₂ O 列を追加、列を破棄し、溶出液を保存する 20,000 x g に 1 分のためにスピンします
    。

2. 5 ' リン酸化と DNase 処理終了

1. 追加 5 μ L ポリヌクレオチド バッファー、5 μ L 10 mM ATP 1 μ L リボヌクレアーゼ阻害剤 1 μ L DNase と 2 μ L ポリヌクレオチド キナーゼサンプルし、, 37 ° C で 30 分必要に応じて 1.1 ボリュームを混合することによりプロトコルを通してマスター反応混合物の準備のため n × (n = サンプル数) の各コンポーネント
。
2. バッファーを変更し、次の手順でポリ a 付加しないようにスピン列に ATP を削除します。
   1. Prespin スピン - 735 x g で 1 分間に列は列を新しい 1.5 mL バイアルに転送し、列にキナーゼ反応を読み込みます。スピン 735 x g に 2 分の列列を破棄収集反応チューブを氷の場所や-80 デパート ° C

3。ブロック 3 ' 三リン酸コルジセピンで終わる

注: ブロック、3 が不可欠です ' 以降ライゲーション反応 3 で自分の concatemerization を避けるために RNA のフラグメントの端 '、(既にブロックされていない終了。3 の付加によって扱われる加水分解後の繰返し) リン酸 ' dATP (三リン酸コルジセピン) チェーン ターミネーター塩基ポリ a ポリメラーゼの助けを借りて。ここでは、0.5 mg/mL の濃度で表現され、 ⁸ で説明されているように精製酵母ポリ a ポリメラーゼ (yPAP) が使用されました。酵母や 大腸菌 PAP 両方 3 を追加するためほぼ同じ活動がある ' dATP をし、商業的に購入することができます (資料の表を参照してください).

1. 追加 13.5 μ L 5 倍濃縮多ポリメラーゼ反応バッファー、10 mM 3 の 2 μ L ' dATP、1 μ L RNase 阻害剤との反応に 1 μ L のポリ a ポリメラーゼ ステップ 2.2.1。ミックスと各反応を 30 分追加 32.5 μ L H ₂ O の 37 ° C で 1 s. 加温のためのスピン。1.10.1 のステップのように RNA を浄化します。溶出が 14 μ L H ₂ o ・ RNA

4。逆に 3 の結紮 ' 5 アダプター ' RNA のフラグメントの末尾

1. が 6 μ L 追加 3 μ l 10 x T4 RNA 結紮バッファー、3 μ l の 10 mM ATP にボリュームを減らすために 10 分間真空濃縮、反応。、15 μ L 0.1 mM 逆の補数 3 1 μ、1 μ L RNase 阻害剤 PEG­ 8000 ' アダプター " revRA3 " (材料表参照) 1 μ L 高濃度 RNA リガーゼ 1、ミックスと。
2. では、1,000 rpm で断続的な混合加熱ミキサー上で 16 時間 24 の ° C で反応を孵化させなさい。各反作用に 70 μ L H ₂ O を加え、混ぜます。1.10.1 のステップのように RNA を浄化します。14 μ L H ₂ o. サンプルを-80 ° C でこの時点で格納できる RNA を溶出します
。

5。逆に転写 (RT)

1. 場所 11 μ。 転送に 200 μ l の PCR 反応に、量を削減する 3 分の真空濃縮で溶出チューブ。1 μ L 0.05 mM RT プライマーを追加 " バイオ デュ dT25 "。PCR サーマルサイクラーの 70 ° C で 5 分間加熱し、5 分の RT のまま
2. は、1 μ L の 10 mM dNTPs、4 μ L 5 x 逆転写酵素バッファー、1 μ L 0.1 M DTT、1 μ L RNase 阻害剤と 1 μ L 逆転写酵素を追加します。ミックスし、PCR サーマルサイクラーの 80 ° c 10 分 55 ° C まで、10 分の反応の熱します。氷の上や長期保管-80 ° C で維持します
。

6。ウラシル DNA グリコシラーゼ酵素ミックス消化

1. ピペット 100 μ L 1.5 mL バイアルにストレプトアビジン ビーズ 800 μ L ビオチン結合バッファーで再懸濁し、磁気ラックに置き。2-3 回チューブを反転します。クリア時のバッファーを削除します。洗浄手順を繰り返します。200 μ L ビオチン結合バッファーのビードを再停止しなさい
。
2. は、ビーズ ソリューションへの逆転写反応を追加し、ホイールの回転に 4 ° C で 20 分を孵化させなさい。2 倍のビオチン結合バッファーでビーズを踏む磁気ラック 6.1 および 2 x 10 のバッファーで洗浄。50 μ L の 10 のバッファーのビードを再懸濁します 2 μ L ウラシル DNA グリコシラーゼ酵素ミックスを追加して、断続的な混合ミキサーで 37 ° C で 1 時間インキュベートします
。
3. 追加 50 μ L H ₂ O、リボヌクレアーゼ H の 11 μ L バッファー、1 μ L RNase H、反応します。磁気ラックに 20 分位チューブの 37 ° C で孵化させなさい、新しいチューブに劈開の cDNA を含有する液体を転送
4. 劈開の cDNA を浄化します。
   1. 胸の谷間反応する PCR 精製キットの PB のバッファーの追加 550 μ L。10 μ L の 3 M 酢酸ナトリウム、pH が 5.2 pH を下げるために追加します。反応をロードします。最小限の溶出スピン列と 1 分 17,000 x g でスピン
   2. 列および 17,000 x g で 1 分間にスピンを追加 750 μ L バッファー PE を流れを破棄します。乾燥する 1 分 17,000 x g で列をスピンします。1.5 mL のバイアルに列を転送、16 μ L H ₂ O および 17,000 x g で 1 分間でスピン 7 μ L のボリュームに集中する 8 分の真空濃縮で反応の場所を追加します
   。

7。5 の結紮 ' 5 アダプター ' cDNA の端

1. 分離 cDNA T4 RNA リガーゼ 1 バッファー、3 μ L 10 mM ATP 15 μ L x 3 μ 10 を追加 PEG­-8000、1 μ L 50 μ M " revDA5 " オリゴ、、1 μ L 高濃度 T4 RNA リガーゼ 1。H ₂ O 各反応に 24 ° C、20 h. 追加 70 μ L インキュベートします。サンプルは-20 ° C でこの時点で保存できます
。

8。PCR 増幅のライブラリ、サイズ選択をパイロット

1. パイロットの反応指数段階の増幅をライブラリに到達するサイクルの PCR の最適な数を決定します。
   1. ピペット 25 μ L の DNA ポリメラーゼのミックス、20 μ L の ligation の反作用、2 μ L H ₂ O、1.5 μ 10 μ M 前方 PCR プライマー (RP1)、1.5 μ 10 μ M 逆 PCR インデックス プライマー 200 μ L の PCR チューブにします
   。
   2. は、次のプログラムで、サーマルサイクラーを実行: 95 ° C、20 の 20 のサイクルが続く 3 分 s 98 ° C、20 s 67 ° C と 30 s 72 ° c. を集める 7 μ 因数 cycler から直接 6、8、10、12、14、16、18 のサイクル後。バッファー (50% のグリセロール、0.05% キシレンシアノール) x 1 μ 10 を追加します。注: 場合はバーコードを結合するときに多重化を使用してサプライヤーの推奨事項に従ってください
   。
   3. 、1:10、緑の蛍光染料の 00 希釈を含む 1 x TBE バッファーに 2% の agarose のゲルの小さなスロット別々 の製品。
      1. ロード因数 2% の agarose のゲル、ゲルのドキュメンテーション システムの 100 ボルト 15 分 PCR の製品の可視化移行のためのゲルを実行します
      。
2. パイロットの反応 ( 図 2) のボリューム 2 倍に大規模な PCR の反作用のためのパイロットの反応指数関数的増幅の先頭サイクル数を使用して。
   1. 大規模な PCR の反作用の集中、PCR 精製キットの最初の反応を脱塩し、1 x TBE バッファー内の 2% アガロースゲル上ワイド スロット上の製品を分離します
   。
3. 200-350 nt DNA を含むゲルのスライスをカット製品。最大 30 分間常温カオトロ ピック バッファーのゲルを溶かします。ゲル抽出キットとゲルのスライスから DNA を抽出します。A 豊富な DNA ⁹ のバインディングでバイアスを防ぐために 50 ° C に加熱しないでください
。 シーケンスの
4. 送信します
   。 注: 通常、50 サイクル シングル読み取り (SR50) は十分です (例えば、については https://www.illumina.com/technology/next-generation-sequencing.html を参照).

9。データ処理

注: gitlab リポジトリ (https://git.scicore.unibas.ch/zavolan_public/A-seq2-processing) で使用可能なソフトウェア (fastq 形式) での結果のシーケンス データが処理されます。分析には 4 つの主な手順が含まれています: (1) ダウンロード git リポジトリ、仮想環境のインストール (2) (3) 特定のパラメーターの設定構成ファイルと (4) で起動による分析 ‘ snakemake ’ ¹⁰. ステップ 4 で実行全体の分析 1 つだけコマンドが必要です。分析の詳細なステップバイ ステップの説明は、gitlab リポジトリ内の README ファイルで見つけることができます、簡単な説明は以下。すべての個別の処理ステップが公開されているツールは、いずれかの外部ソースからの実行あるいは社内準備します。計算のパイプラインは、snakemake パッケージの利用できる ¹⁰ anaconda ベース ¹¹ python 3 仮想環境に依存します。それは Unix ライクなオペレーティング システムを持つマシン上で実行し、インストールされている CentOS 6.5 オペレーティング システムと Linux 環境で 40 GB 使用可能な RAM 試験しました。ソフトウェアの依存関係は、仮想環境内で自動的に制御されます。次の公に利用可能なソフトウェア ツールが必要でありそれにより環境と共にインストール: snakemake (v3.9.1) ¹⁰, fastx ツールキット (v0.0.14) ¹² 星 (v2.5.2a) ¹³、cutadapt (v1.12) ¹⁴, samtools (v1.3.1) ^{14 , 15}, bedtools (v2.26.0) ^{16 , 17}.

1. Cdna を読み取りからデータ前処理
   注: シーケンスの深さは実行の間に異なる場合があり、楽器によって 1 つのサンプルからのデータを複数のシーケンスのファイルに分割できます。この場合は、次の手順で使用されている 1 つの入力ファイルに 1 つのサンプルに対応するファイルを連結します。
   1. Fastq から fasta 形式にファイルを変換します
   。
   2. (読み取りの位置 5、6、7 で 3 thymidines) の正しい構造を持つエキスを読み取ります
      。 メモ: 上記実験的プロトコルに従って準備が正しく読み取りが必要構造 (5 から ' 終わり): 4 ヌクレオチド バーコード - 3 thymidines - 逆に成績証明書 3 の補数 ' 終わり
   。
   3. は、シーケンスの説明の行で開始四量体に関する情報を格納します
      。 注: 後における増幅デジタルアーティファクトの補正を容易にするユニークな分子識別子 (UMI) として、テトラマー
   。
   4. 読んでから最初の 7 つのヌクレオチドを削除 ' s 5 ' 終了します
   。
   5. と同じ読み取りの 1 つだけのコピーを保つことによって増幅アイテム挿入のシーケンスと海を修正します
   。
   6. 3 の部分を削除する ' エンド アダプター シーケンスとその逆の補数シーケンスに一致します。のみ、最小長さの読み取りを続行 (既定: 15 nt).
      注: 元の mRNA フラグメントと、3 サイクル シーケンス数の長さに応じて '、読み取りの終了は 3 の部分を含めることができます ' アダプターは、この手順で削除されます
   。
2. 次の条件を満たすすべての読み取りを抽出: 最大 2 不明なヌクレオチド (' N ')、として、最大 80% と読むない A. 最後のヌクレオチドこれらの読み取り、解析で使用される十分な品質であることと見なされます
。
3. 読み取りをスプライシングの読み取りを処理し、BAM 形式で出力ファイルを生成するツールをゲノムにマップします。
   1. スターの場合は、読み取りする必要がありますマップするゲノムのインデックス ファイルを作成します。人間のゲノムのこのステップは 35 GB のメモリ (RAM) を必要があります
   。
   2. は、読み取りをゲノムにマップします
      。 注: (星固有の注意事項) ソフト クリップは 3 のマッピングを強制するために無効になって ' それぞれの終わりを読んでこれはヌクレオチド胸の谷間サイトのすぐ上流
   。
4. BAM をベッド ファイルに変換します。読み取りは、複数の場所にマップする場合これら最低の編集距離だけを維持します
   。 注: 特定の場所でマップを読むのコピー数はスコアとして使用されます。複数の場所にマップされる読み取り、読み取りのマップ先の場所の 1/数と等しい重量の各場所で分数カウントされます
。 可能性が高い順序の誤りによって異なる
5. 崩壊を読み取ります。同じ場所に 2 つの異なる読み取りマップする場合 (マッピングの開始と終了位置が同じ) と同じ海を共有、PCR 重複としてそれらを考慮して 1 つしか保存します
。
6. 推論すべての個々 の mRNA 3 ' 処理サイトを終了します
   。 注: 個々 の読みは 3 で証拠を提供する ' その最後の 4 つのヌクレオチドがエラーなしゲノムにマップされるときに終了します。先の位置、3 ' 読み取りのマップの端、胸の谷間のサイトとして格納されます
。
7. 検出 3 ' 内部プライミングから由来するかもしれないサイトを終了します。10 ゲノムで胸の谷間サイトの下流に nt を満たすとき内部のプライミング加工品としてサイトを定義する次の条件の 1 つ: として六つ以上が含まれています、6 年連続を含むまたは次テトラマーのいずれかから始まります: AAAA, AGAA, 佐賀, AAAG.
8. 3 の個々 のテーブルを生成 ' ベッド形式の処理サイトを終了します
。
9. 識別がポリ a サイト クラスターを個別調整
   。 注: ここで説明する手順に従って前の発行 ⁵ で導入された手順に従います。
   1. 個々 の 3 を集めることによって開始 ' 研究のすべてのサンプルで得られた処理サイトを終了します
   。
   2. 9 の地域で知られている多信号 ⁷ に注釈を付ける各個別 3 周りヌクレオチド ' エンド処理サイト
   。
   3. 識別多サイト各サンプルの背景の上次のように表されます。
      1. は、現在のサンプル内の元の式でサイトを並べ替えます。上から下、ゲノムのあらかじめ定義された距離内にある場合より高いランク付けされたサイトとランクの低いサイトを関連付けるサイトのリストを走査 (既定: 25 nt アップまたはダウン ストリーム) 上位サイトから
         。 注: 上位サイトに関連付けられたすべての低ランキングのサイト定義式はこれらのサイトのすべてを文書化する読み取りの数クラスター
      。
      2. 式でこれらのクラスターを並べ替えるし、最高から最低式に式しきい値 c で注釈付きポリ a を持つクラスターの割合信号定義済みしきい値 (下のドロップを決定するクラスターのリストを走査既定値: 90%).
      3. サイトのカットオフ以下の任意のクラスターからを破棄します
      。
   4. クラスター近接 3 ' のサンプルの間で得られるサイトを終了します
      。 注: 並べ替え 3 ' 終了のサンプルをサポートの数し、正規化の合計によってサイトをまず処理カウントを読む (あたりの読み取り数百万 (RPM)) のサンプルの間で。上から順に、上位サイトへの間隔で定義済みの上限よりも大きくない場合、低ランクのサイトをより高いランク サイトに関連付けるリスト内の移動 (デフォルト: 12 nt)。いずれかの構成する 3 ' エンド サイト注釈付き多信号と重複または多信号の直接内部のプライミングを検出する詳細な検査に対応するクラスターをマークする下流
   。
   5. マージ多サイト クラスター
      。 注: クラスターが推定内部プライミング候補としてマークされると、それが 2 つのクラスターを共有、多信号に応じて下流クラスターにマージまたはクラスターで最もダウン ストリーム サイト最低に位置する多信号がある場合は保持上流の距離 (既定: 15 nt)。場合に密接に間隔をあけられたクラスターを結合する最後に、: 同じ多込むを共有 (i) または (ii) 結果のクラスターのスパンが最大値を超えない (既定: 25 nt).
   6. すべての 3 からカウントを読む正規化総ベッド ファイル形式でクラスターを保存 ' スコアと各クラスター内のサイトを終了します
   。

ポリ (A)-アルカリ加水分解による断片化の培養細胞から分離された RNA を含むと Cdna のランダムプライマー プライマーが付いているトランスクリプションを逆で作られました。結果として得られる cDNA はストレプトアビジン ビーズに固定化した、ウラシル特定切除反応でデュ開裂、5 ' アダプターを結紮、3' 端に裂かれたフラグメントの挿入されました。図 1は、実験のグラフィカルな概要を示しています。

Hela 細胞 HEK293 細胞の 10 の6セルされたプロシージャの終わりに蛋白質のコーディングの遺伝子の大半のためのポリ a サイトを特定するのに足りる。ただし、他の種類の細胞や組織の識別された多サイト数実験で使用されているセルの数と彩度をテストする必要があります増加します。パイロットの PCR の代表の結果のステップし、DNA のフラグメントのシーケンスの前にサンプルの分析は、図 2に示します。

シーケンサーから得られるゲノムにマップする準備ができている品質チェック、アダプターによりトリミングされた読み取りで終わる fastq ファイルから始まって計算の分析の前処理のステップを図 3に示します。図 4は、mRNA 3' 末端処理特定のサンプルで特定されたサイトのカタログに対応するゲノムと末尾に読み取りのマッピングで始まる解析手順を示します。複数のサンプルを分析するときの追加手順が個々 のサンプルで発見されたサイトを処理して、3' 末端に一致し、サンプルの間で彼らの余剰を報告する実行されます。これらの手順は、図 5のとおりです。

したがって、サンプル化されている、利用可能な処理パイプラインを通じて (fastq 形式) のファイルを読み取り結果のシーケンスの分析は簡単です。サンプルについての情報を構成ファイルに追加すると、パイプラインの実行は出力ファイルの 2 つの主なタイプになります: 1) ベッド-ファイルとすべて 3' 終了 (は例えば「個々 のサンプルの識別処理サイトsample1.3pSites.noIP.bed.gz「)、2) 研究のすべてのサンプルの間ですべてのポリ a サイト クラスター (clusters.merged.bed) とベッド ファイル。出力では、すべての読み取り (例: "sample1 それぞれの個々 のサンプルからのゲノム座標も含まれています。STAR_out/Aligned.sortedByCoord.out.bam") その後 IGV16のようなゲノムのブラウザーで表示できます。読み取りのプロファイルの外観検査は一般的にポリ a ゲノムの研究で実施された特定の摂動によって発生する変更サイトで配信の最初の一見を提供します。たとえば、図 6の HNRNPC 蛋白質のノックダウンする特定の遺伝子の応答が表示されます。

これらのゲノム広い分布の概要は (表 1) を設けています。具体的には、"カウント/annotation_overlap"ディレクトリ内の出力ファイルを含む特定の注釈機能と重複サイトの分数 (gtf ファイルから入力として提供; 注釈: 3' UTR、ターミナルのエクソン、エクソン、イントロン、遺伝子間)。最後に、各サンプルは、個別の処理ステップの結果は、("sample1.summary.tsv"など) をも保存されます。これの数が含まれています: ゲノムに一意に対応する各サンプル、5' 末端の予想される構造の読み取り、読み取りにおける PCR 重複を崩壊した後で、読み取り、高品質 9.2 のステップで定義されている条件に従って読む読み取り(シーケンス エラーに起因するそれらを崩壊した後手順 9.5 を参照)、複数のマッピングを読み取り (シーケンス エラーに起因するそれらを崩壊した後手順 9.5 を参照)、生 (クラスター化されていない) 3' 端の各サンプルでは、raw 3' 末端処理サイトのサイトを処理候補内部プライミングすることがなくユニークな 3' 終了処理サイト内部プライミング候補と最終的なポリ a サイトのクラスターのセットなしすべてのサンプルから。

図 1: A seq2 プロトコルの主な手順です。個々 のステップは、図の左側に示されます。挿入の RNA のフラグメントは、逆のトランスクリプション; 後の cDNA の赤く緑の線として描かれています。アダプターは、明るい青またはオレンジに着色されています。この図の拡大版を表示するのにはここをクリックしてください。

図 2: パイロット PCR および最終製品プロファイル。PCR の反作用から (、) の因数は異なるサイクルで収集された、2% アガロースゲル上区切られます。左に数字は、ヌクレオチド DNA の梯子のそれぞれのバンドのサイズを示します。この実験で 12 サイクル (*) は、大規模な PCR の反作用のため選ばれました。サイズの後サンプルの例 (b) 選択は、明らかに平均サイズは約 280 のヌクレオチドのフラグメント サイズ アナライザーで実行します。[FU] 左に数字は相対的な信号強度を示します。この図の拡大版を表示するのにはここをクリックしてください。

図 3: 読み取りをシーケンシングの前処理のアウトライン。シーケンサー機器関連ソフトウェアによって生成される読み取り fastq ファイルは、対応するゲノムにマップする高品質の読み取りを識別するために処理されます。「データ処理」に記載されているプロトコル個々 手順へのリンクが、パイプライン個々 ステップの入力/出力仕様を示します。この図の拡大版を表示するのにはここをクリックしてください。

図 4: シーケンスの概要読んで処理、個別 3' 末端の世代にゲノムにマッピングの手順から処理サイト。図は i へのリンクと、パイプラインの個々 のステップの入力/出力仕様「データ処理」に記載されているプロトコル個々 手順を実行。ユーザーに配信される主な出力ファイルは、太字で示されます。この図の拡大版を表示するのにはここをクリックしてください。

図 5: 共同規制 3' 末端配列サイトのクラスターを生成する撮影されている手順の概要です。「データ処理」に記載されているプロトコル個々 手順へのリンクが、パイプライン個々 ステップの入力/出力仕様を示します。メイン出力ファイルは、太字で示されます。この図の拡大版を表示するのにはここをクリックしてください。

図 6: 3' のプロファイルの例の結果終了 IGV 16ゲノム ブラウザーに表示される NUP214 遺伝子のターミナルのエクソンに沿って処理読み取ります。A seq2 読み取り HEK 293 細胞、コントロール siRNA または HNRNPC siRNA との 2 つのサンプルの調製。解析パイプラインによって注釈が付けられる多サイト記載されている読み取り IGV のゲノムのブラウザーへの入力として使用された BAM 形式で保存されました。MRNA 3' にマップ読み取りピークの 3' 端 Ensembl で注釈が終了します。プロファイルは、HNRNPC ノック ダウン時に長い 3' UTR アイソの使用の増加を示します。この図の拡大版を表示するのにはここをクリックしてください。

表 1: 解析パイプラインの出力例。個々 のステップで得られた読み取りの概要。

著者が明らかに何もありません。