短読シーケンシング技術を用いた尿細菌の完全ゲノム生成のためのハイブリッド デノボ ゲノムアセンブリ

細菌は、機関審査委員会承認研究19MR0011(UTD)およびSTU 032016-006(UTSW)の一環として、同意した女性から採取された尿から培養された。 1. 改良尿培養 注:すべての培養ステップは、無菌条件下で行う必要があります。すべての器具、ソリューション、メディアを殺菌します。70%エタノールで作業領域を清掃し、ブンゼンバーナーを設置し、炎の近くに慎重に作業して汚染の可能性を減らします。あるいは、クラスIIのバイオセーフティキャビネットは、無菌環境を維持するために使用され得る。病原性の可能性がある微生物への暴露を避けるために、適切な個人用保護具(PPE)を着用してください。 メッキグリセロールストック尿とコロニー分離 室温(RT)でグリセロールストック尿を解凍します。一度解凍し、5sの混合物を混合するためにサンプルを渦を出す。無菌マイクロ遠心チューブで、尿の1:3と1:30希釈液を無菌1xリン酸緩衝塩(PBS)で100 μLの最終体積に調製します。注:グリセロールストック尿は、500μLの希釈されていない尿と500μLの500μLの無菌グリセロールを凍結し、-80°Cで保存することによって調製されます。 37°Cで寒天プレートを15分間使用してください。一般的な尿中細菌属に適した培地の種類や培養条件については 、図1 をご覧ください。めっきする前にピペット加工して希釈した尿をよく混ぜ、希望の寒天板に希釈した尿のプレート100μLをプレートし、滅菌ガラスビーズを使用してサンプルを広げます。1xPBS希釈液のプレート100μLを、成長制御なしとして別のプレート上に取り付ける。注意:一般的な泌尿器病原種(例えば、 大腸菌、クレブシエラ属、エンテロコッカス・フェカリス など)を培養しようとする場合、泌尿器系細菌種の同定が容易となるため、染色体寒天(材料表)を使用することが推奨される(図1)。コリスチンナリジキシック酸(CNA)またはMRS寒天は、非選択的寒天で食潔種を上回る可能性があるグラム陰性泌尿器原体を含むことが知られている尿から潔癖性グラム陽性種(例えば、 ラクトバチルスspp.)を単離するのに有用である。 尿病原菌の場合は24時間、食潔菌では3~5日間、35°Cで所望の大気条件でプレートを反転させる(図1)。 インキュベート期間の後、インキュベーターからプレートを取り出します。各プレートから、独特の色、形態、または血分解パターンを示すコロニーを選びます。 対応する寒天に滅菌ループを使用して細菌コロニーを再ストリークし、望ましい雰囲気の中で2〜5日間反転したプレートをインキュベートし、十分に分離されたコロニーを得る。注:BAPを一次培養に利用する場合、染色原性寒天にコロニーをパッチさせることは、サンプル中の細菌集団の不均一性に関する有用な情報を提供するかもしれない。 液体ブロスおよびグリセロールストッキング細菌分離株の培養 親コロニーの形態に一致する単離されたコロニーが得られたら、単一のコロニーを選び、滅菌接種ループを使用して3mLの液体ブロスに接種する。一般的な尿中微生物叢属の成長をサポートすることができるブロスについては 図1 を参照してください。寒天プレートをパラフィルムで密封し、4°Cで2~4日間保管します。培養物が目に見えて濁るまで、所望の大気条件で液体培養物を1〜5日間インキュベートする。 成長が観察された後、培養物をボルテックスし、2 mLのグリセロール中の500μLの無菌500μLに一晩培養物1mLを加える。シールし、反転によって穏やかに混ぜます。各コロニーに2つのグリセロールストックを用意し(1つはバックアップとして機能します)、-80°Cで保管してください。 2. 16S rRNA遺伝子サンガーシーケンシングによる細菌種の同定 注:微生物の同一性は、代わりに、飛行質量分析(MALDI-TOF)20のマトリックス支援レーザー脱離イオン化時間を使用して確認することができます。 コロニー-ポリメラーゼ連鎖反応(PCR) PCRチューブに25 μLのPCR反応を調製するには、12.5 μLの2x Taqポリメラーゼマスターミックス、0.5 μLの10 μM 8Fプライマー、0.5 μLの10 μM 1492Rプライマー(材料表)、ヌクレアーゼフリー水21.5μL を加えます。注: 複数のサンプルに PCR を実行する場合は、Taq ポリメラーゼ ミックス、プライマー、および無菌ヌクレアーゼフリー水の反応マスター ミックスを作成します。次いで、各PCRチューブに25μLをアリコートした。 コロニーPCRを実行するには、滅菌爪楊枝またはピペットチップを使用して再ストリークからよく隔離されたコロニーをスワイプします。ステップ2.1.1で調製したPCR反応ミックス中のコロニーを再懸濁させる。やさしく混ぜます。2000 x gでクイック スピンでチューブの底にある液体を収集します。注:サンプルに気泡がないことを確認してください。PCR反応ミックスを単独で含有する非鋳型制御(NTC)サンプルを含む。 サンプルチューブをサーモサイクラーに入れ、次のプログラムを実行します: 95 °C 3分間;40サイクル:30sの95°C、30sの51°C、1分30sの72°C;72 °C 10分間10 °Cで保持します。 ゲル抽出と種同定 PCR の実行が完了したら、0.5x トリス-ホレート EDTA (TBE) バッファーで調製した 1% アガロースゲルで PCR 産物を確認します。ゲルを鋳造する前に、臭化エチジウム(EtBr)を加える。次に、少なくとも20 μLのサンプル量を保持するウェルに対して、くしを使用してゲルをキャストします。注意:EtBrは発がん性があると疑われるインターカリング剤です。取り扱い時には必ず手袋とPPEを着用し、機関のガイドラインに従ってEtBrを含む材料を廃棄してください。 ゲルをセットしたら、0.5x TBEバッファーを充填した電気泳動タンクにゲルを入れ、くしを取り除きます。1 kb のラダーを最初のウェルに、10~20 μL の PCR 反応を後続のウェルにロードします。解決されるまで 100 ~ 140 V で実行します。UV光の下でゲルを視覚化し、NTCウェルに存在しない〜1.5キロバイトで明確に定義されたバンドの存在を確認します。注意:紫外線は皮膚や目に有害であり、ゲルを視覚化する際に適切なガードを使用し、適切なPPEを着用してください。注意:コロニーPCRは、いくつかの細菌のために失敗する可能性があります。単離されたgDNAからのPCRを進め、代替オプション22です。 カミソリを使用して〜1.5キロバンドを物品切りし、ゲル切削をきれいなマイクロ遠心チューブに移します。メーカーの指示に従ってゲル抽出プロトコルを進める (資料一覧)。マイクロボリューム分光光度計で精製されたDNAの濃度を測定します。注: 10 ng/μL >濃度が望ましく、A260/280 の間では 1.7~2.0 が許容されます。 各サンプルに対して2つのサンガーシーケンシング反応を準備し、1つは8Fを使用し、もう1つは1492Rプライマーをヌクレアーゼフリーウォーターで使用し、選択したサンガーシーケンシングサービスのガイドラインに従います。 シーケンシングデータを受信したら、DNA配列をNCBI基本局所アライメント検索ツール(BLAST)ウェブサイト(blast.ncbi.nlm.nih.gov/Blast.cgi blastn)にアップロードし、ヌクレオチドBLAST(blastn)を選択し、rRNA/ITSデータベース16SリボソームRNA配列(細菌および古細菌)を選択し、メガブラストプログラムを実行します。分離は、データベースからの参照に対する最高品質のヒットによって識別される可能性があります。注:一部の細菌種は、その16S rRNA配列において高い同一性を示し、この方法だけでは区別できない場合があります。スペシレーションは、同じ属のメンバーを自信を持って区別するためにDNA相同性と生化学的分析を必要とします23. 3. ゲノムDNA(gDNA)の抽出 注:このセクションでは、多様な細菌種からの高品質のゲノムDNAの高収率抽出のために 、材料表 に記載されているgDNA抽出キットで提供される試薬とスピンカラムを利用します。以下に示す推奨される変更と手順を示します。 メーカーの指示に従ってキット試薬を準備します。 適切な滅菌ブロス(図1)で、十分な増殖が観察されるまで、十分な増殖が観察されるまで、図1に示された温度と大気圧で細菌を注入することによって、適切な滅菌ブロス( 図1) に3〜10mL培養を調製する。 インキュベーション後、分光光度計24を用いて培養物の光学濃度を600nm(OD600)で測定する。 一晩培養物を1:10比で希釈して定量化するためのサンプルを準備します。同様に測定のための無菌培養培地のブランクを含める。サンプルの読み取り値から空白の読み取り値を減算し、希釈係数 10 を掛けることで光学密度を計算します。 OD600測定と、種に対するCFU/mL比を事前に確立したOD600を用いて、2 x 109細胞を得るために必要な培養量を算出する。 ペレットに5000 x g で5分間必要な培養量を遠心分離する。上清を吸引し、200 μLの冷たいTEバッファー(手順の開始時に氷上でのプレチル)でペレットを再懸濁します。 5000 x gで 2 分間サンプルを遠心分離します。上清を取り出し、180 μLの酵素リシスバッファー(ELB)でペレットを再懸濁し、20 μLの予熱したRNase A(10 mg/mL)を加えます。グラム陽性菌の効率的なリシスのために、18 μLの変異不反応(25 kU/mL)を加えます。ボルテックスウェル、次いで、2時間回転器上の37°Cでサンプルをインキュベートする。注: グラム陽性菌とグラム陰性菌の両方に対して、メーカーのプロトコルに記載されている ELB を使用することをお勧めします。 製造元の指示に従って進みます。注:必要に応じて、さらに1〜2回溶出手順を繰り返して、追加のgDNA収量を得る。 セクション4で指示された抽出gDNAの品質を評価し、1週間以内に使用する場合は4°CでgDNAを保存します。あるいは、gDNAを-20°Cに保ち、長期保存を行います。 4. 抽出したgDNAの品質の評価 ゲル電気泳動による品質を評価するために、サブセクション2.2に記載されているように1%のアガロースゲルを調製する。きれいなチューブでサンプルを準備する:抽出されたgDNAの1-2 μLとパラフィルム上の2x負荷染料の3 μLを混ぜます。ゲルを一度ロードして実行し、UV ライトの下で視覚化します。注:gDNA抽出の成功は、ゲルの上部にある離散バンドと最小限のスミアによって明らかになります(図2A)。塗り傷はせん断を示す。gDNAバンドが明らかでない場合やスメリングが実質的である場合は、gDNA抽出を繰り返します。RNase AおよびプロテアーゼKのインキュベーション時間を短縮することを検討してください。1.5~3 kb前後の2つのバンドが観察された場合、RNA汚染が示唆される(図2B)。新鮮なRNase Aを準備し、抽出を繰り返します。 マイクロボリューム分光光度計で品質を評価するには、マイクロボリューム分光光度計でgDNA濃度と吸光度比A260/280を測定します。濃度は、1.7~2.0の間で>50 ng/μLおよびA260/280の間で許容されます。注:低gDNA収量は、低入力、高入力、ヌクレアーゼの汚染、不十分なリシスが原因である可能性があります。上記範囲の吸光度比は、RNA汚染を示す。gDNAの品質が悪い場合は抽出を繰り返します。 フルオロメーターによる品質を評価するには、高感度アッセイキットと蛍光計計器を使用してgDNA濃度を定量化するメーカーの指示に従ってください(材料表)。濃度>50 ng/μLが望ましい。 5. ペアエンド次世代短読シーケンシングとライブラリの準備 メモ:短読シーケンスは、異なる読み取り長さと向きで様々な楽器に対して実行することができます。150 bp (300 サイクル) ペアエンドシーケンシングは、細菌の WGS に推奨されます。ライブラリの準備とシーケンシングの両方をコア施設または商業研究所に委託することができます。 メーカーの指示に従ってシーケンスライブラリを準備します (資料一覧)。メーカーの推奨最終ローディングライブラリ濃度に従ってください。ただし、NextSeq インストゥルメントでの読み取り生成を最適化するために、プールされたライブラリを 1.8 pM でロードすることをお勧めします。 オプションですが、バイオアナライザー(材料表)を使用して、プールされたライブラリのフラグメント分布を評価し、フラグメントサイズが平均600 bpであることを確認します。 6. ナノポア MinION シーケンシングライブラリの準備 製造元のプロトコルに従ってシーケンス ライブラリを準備します (資料一覧)。2つのバーコード拡張キットを使用すると、1つのフローセルで最大24個のサンプルを多重化できます。24個のサンプルを多重化する場合は、2つの部分、12個のサンプルでライブラリの準備を行うことをお勧めします。24個のサンプルは、以下に記載されているようにプールすることができる。注:サンプルは、ネイティブバーコードライゲーションを終了すると一晩で4°Cで保存することができます – これは、必要に応じて、プロトコルの停止点を提供します。ライブラリー調製プロトコルのネイティブバーコードライゲーションセクションの終わりには、可能な最大DNA質量(ng)まで各サンプルの正モル量をプールすることが推奨されます。 これを行うには、メーカーの指示に従ってフルオロメーター(材料表)を使用して、バーコードライゲーションに続くすべてのサンプルを定量化します。dsDNA濃度が最も低いサンプルの量を推定し、このサンプルで見つかったdsDNAの合計を計算します。この数を使用して、一緒にプールされる他のすべてのサンプルの等量を決定します。注意:等量計計算はプールされたdsDNAの量を最大にし、大量のプール(>65 μL)を生成するため、プールを集中させるためにクリーンアップが必要です。 dsDNAプールのクリーンアップと濃度 2.5倍の常磁性ビーズ(材料表)をDNAプールに加え、チューブを軽くフリックして内容物を混ぜます。チューブをRTで5分間ローテーターに入れ、サンプルを2000 x g でスピンダウンし、ペレットを磁石に置きます。 250 μL を加えて、作りたての 70% エタノール (ヌクレアーゼを含まない水) を加え、ペレットを乱さないよう注意してください。エタノールを吸引し、エタノールを一度洗浄します。 2回目の吸引の後、サンプルを2000 x g で回転させ、磁石の上に戻します。任意の残留エタノールをピペットオフし、サンプルが約30 sのために乾燥することを可能にする。 マグネットからチューブを取り出し、60~70μLのヌクレアーゼフリー水でペレットを再懸濁します。RTで2分間インキュベートします。ペレットは、エルテが透明になるまで磁石上のサンプルをペレットにし、その後、エルテを除去し、きれいな1.5 mLマイクロ遠心チューブに移します。 フルオロメーターを使用して濃縮プールを定量化し、次にアダプタライゲーションステップに進むアリコートを準備します:65 μL最終体積でサンプルの700 ngを調製します。プールの残りの部分を 4 °C に保持し、最初の実行が終了した後に 2 回目の実行を完了します。 メーカーの指示に従ってアダプターの結紮を進め、フローセルにサンプルをロードします。シーケンス実行を開始します。注:サンプルローディング前のフローセルプライミングポートからの吸気空気および約200 μLの貯蔵バッファー。これは、フローセルのプライミングとサンプルの読み込みが成功するために重要です。フローセルのプライミングポートを介してソリューションを描画および堆積する際に、p1000ピペットとヒントを使用してください。 製造元の指示に従ってライブラリを並び順に指定します。 シーケンス用のオペレーティング・ソフトウェアを開き、 開始をクリックします。実験の名前を入力すると、推奨される命名法には実行日とユーザー名が含まれます。[ キットの選択に進む] をクリックし、使用するライブラリの準備キットとバーコード拡張パックを選択し、[ 実行オプションへ進む] をクリックします。 2 回目の実行に十分なライブラリを準備する場合は、実行の長さを 48 時間に調整します (それ以外の場合は、デフォルトの 72 時間のままにします)。[ ベースコールに進む] をクリックします。 ベースコールオプションの「Config: 高速ベース呼び出し」をチェックし、出力 FASTQ ファイルがバーコードシーケンスからトリミングされ、バーコードに基づいて別々のディレクトリにデマルチプレクスされるように、バーコードが有効に設定されていることを確認します。 [出力に進む ] をクリックします。 出力シーケンスデータを保存する場所を選択します。FASTQ 出力を保存するだけで、FAST5 出力も保存する場合は、約 30 ~ 50 Gb のデータ>予想します。フィルタリングオプションの Qscore: 7 |のチェックを外しますReadlength: セクション 7.2 で説明されているフィルタリングを続行する場合は、フィルタなし、それ以外の場合はオンのままにして 、Readlength を 200 に調整します。 [ セットアップの実行に進む] をクリックし、すべての設定を確認します。設定が正しい場合は、[ スタート]をクリックし、[ 戻る ]をクリックして必要な調整を行います。 必要に応じて、フローセルは、製造業者の指示に従って洗浄され、残りのプールで再ロードされてもよい。最初の実行が完了し、フローセルが洗浄されたら、残りのプールに対して6.2の手順を繰り返します。注:2回目の実行を設定する場合は、48時間を超えるランで以前使用したフローセルに関するメーカーの推奨に従って、バイアス電圧を-250 mVに調整してください。 7. 読み取りの評価と準備 注 : 推奨ディレクトリ構造を 図 4に示します。以下の計算手順に進む前に、 デスクトップ、Long_Reads、Short_Reads、Trimmed_Readsで見つかったディレクトリを作成します。 短い読み取り (図 3)注: 短い読み取りは FASTQ 形式で生成されます。ファイルには、FASTQ ごとに最大 4000 読み取りが含まれています。これらは、多くの場合、圧縮 (.gzアーカイブ) され、複数のファイルに編成されます。プラットフォームによっては、バーコードは通常トリミングされます。zip形式のファイルを受け入れるプログラムもあれば、インポートする前に抽出が必要なものもあります。読み取りは、ゲノムアセンブリ中のデータ精度を確保するために、品質管理(QC)の手順に合格する必要があります。CLCゲノミクスワークベンチが利用できない場合、トリムモマティックなどのQCショート読み取りとトリムに代替プログラムを使用できます。25 またはトリミングのためのトリムガロア(https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/)と読み取り品質を評価するためのFastQC(http://www.bioinformatics.babraham.ac.uk/projects/fastqc/)読み取りの数を平均読み取り長さで乗算し、ゲノムサイズで割ることによって推定される平均短い読み取りカバレッジは、>100xであることが推奨されます。ゲノミクスワークベンチソフトウェア(材料表)を開き、すべてのペアエンドの短読み取りFASTQファイルをインポートします。ペアファイルは自動的に生成されます。 CLC_Dataの下に新しいフォルダを作成するには、上部のツールバーで [新規作成 ]をクリックし、[ フォルダ..]を選択します。をクリックしてファイルを保存します。必要に応じてフォルダに名前を付けるため、推奨される規則ではサンプル ID を使用します。 上部のツールバーで、[ インポート ] ボタンをクリックし、[ イルミナ] を選択します。サンプルに対応する短読み取りファイルに移動し、選択します。ペアリングされた読み取りオプションが選択されていることを確認し、[ 失敗した読み取り] オプションをオフにします。[ 次へ] をクリックし、[ 保存] を選択して、[ 次へ ] をもう一度クリックします。インポートしたファイルを前の手順で作成した新しいフォルダに保存し、[ 完了] をクリックします。 分離のためのすべてのペアファイルのシーケンスリストを作成します。これにより、分析を簡略化するために、読み取りデータが 1 つのファイルに連結されます。 上部のツールバーで、[ 新規 ] ボタンをクリックし、[ シーケンス リスト] を選択します。左側のディレクトリリストで、連結するファイルを選択し、矢印を使用して右側の選択したファイルリストに移動します。[ 次へ] をクリックし、[ 保存] を選択して、[ 次へ ] をもう一度クリックします。シーケンスリストを保存することを選択し、[ 完了]をクリックします。 シーケンスリストが生成されたら、すぐにサンプル ID で名前を変更します。 シーケンス リストでシーケンス読 み込み 用 QC ツールを実行する: この手順では、読み取りが短い NGS によって生成される読み取りの全体的な品質パラメータを評価します。 ツールボックスメニュー(左下ウィンドウ)で 「シーケンス読み込み QC」ツールを検索します。ツールをダブルクリックし、解析するシーケンス リストを選択して [ 次へ] をクリックします。 すべての出力オプションがオンになっていることを確認し、[結果処理] の下の[保存]を選択します。[次へ] をクリックし、出力ファイルを保存するように指定し、[完了] をクリックします。 シーケンスリストで トリム読み込み ツールを実行する:トリムは品質、長さ、およびあいまいさに基づいて行われます。このプロセスでは、シーケンスで使用されるバーコードが、この手順の前にトリミングされていることを前提としています。 ツールボックス (左下) で [読み取り] ツールを検索します。 [Trim Reads]をダブルクリックし、分析するシーケンスリストを選択して[ 次へ]をクリックします。 品質トリミング: 品質スコアの制限を 0.01 に設定し、あいまいなヌクレオチドを 2 のままにします。[ 次へ] をクリックします。注: パラメータはユーザーの裁量で調整できます。推奨される設定です。 [ 自動リードスルー アダプター トリミング ] をオフにします (CLC にインポートする前に、アダプタが読み取りからトリミングされている場合にのみ行います)。[ 次へ ] をクリックし、[ 読み取り時間を次の長さにする]をオンにして、デフォルトの 15 を使用します。 [ 次へ] をクリックし、[ レポートの作成] をオンにして 、[ 保存] をクリックします。[ 次へ ]をクリックし、出力ファイルを保存する場所を指定します。[ 完了 ] をクリックします。 トリムシーケンスリストをエクスポート: 後続のハイブリッドアセンブリと解析はCLCの外部で完了し、トリムされたショート読み取りファイルをエクスポートする必要があります。 左上のディレクトリナビゲーションから、ステップ7.1.4で生成されたトリミングされたファイルを選択し、上部のツールバーで [エクスポート ]をクリックします。エクスポート ファイルの種類として Fastq を選択し、[ 次へ] をクリックします。 [ペア化シーケンスリストを 2 つのファイルにエクスポート]をオンにします。次に、[ 次へ ] をクリックして、ファイルをエクスポートするTrimmed_Readsディレクトリを選択します。[ 完了 ] をクリックします。トリミングされた短読み取りファイルが、拡張子 .fastqの 2 つのファイル (R1 と R2) として正常にエクスポートされたことを確認します。注: トリムされたシーケンス・リストは、通常 CLC によって R1 および R2 として指定される 2 つのファイルにエクスポートする必要があります。ダウンストリームのハイブリッド アセンブリでは、短読データ入力をそのように設定する必要があるため、このことは重要です。 エクスポートしたファイルの名前を変更し、ファイル名にスペースや特殊文字を使用しないでください。わかりやすくするために、推奨される形式はtrimmed_short_fileです。R1.fastq. 長い (最小) 読み取り (図 3)注: ハイブリッド アセンブリの Long (MinION) シーケンス読み取りを準備するための次のパイプラインは、コマンド ラインで実行される NanoFilt および Nanostat プログラム26 を利用します。先に進む前にツールをインストールし、UNIX の基本を理解してこれらのコマンドを実行してください。デフォルトの端末と Bash シェルが推奨されます。一般的なターミナル コマンドと使用方法のレッスン ガイドは、ソフトウェア大工27にあります。以下の手順では、生成されたファイルがバーコード命名法(NB01、NB02など)で命名され、Long_Readsディレクトリに保存されることを前提としています。または、シーケンス実行を設定するときに MinKNOW を使用して読み取りフィルターを実行できます。平均ロング読み取りカバレッジは、>100xであることが推奨されます。推奨される平均読み取り長さは >2000 bp です。したがって、必要な長い読み取りの数は、短い読み取りの数よりも少ないです。 Long_Reads ディレクトリ内の実行で使用される各バーコード (barcode01、barcode02 など) に対して新しいディレクトリを作成します (図 4)。各バーコードに対応するすべての .fastq ファイルを適切なフォルダにコピーします。各実行から各バーコードのすべての .fastq ファイルを結合します。 ターミナルを開き、cdコマンドを使用してLong_Readsディレクトリ内のバーコードディレクトリに移動します: cdデスクトップ/Long_Reads/バーコード01 次のコマンドを実行して、バーコードごとにすべての.fast >qファイルを 1 つの.fastqファイルに連結します。 注: このコマンドは、各 FASTQ ファイルからの読み取りをすべて NB01.fastq という名前の 1 つの大きな 1 つの FASTQ に結合します。 次のコマンドを実行して、サンプルの読み取り品質を評価するために ナノスタット を使用します: 出力をテキストまたは Word ファイルにコピーして、後で参照するために結果を記録します。 MinION 読み取りをフィルタリングするには 、MinION 読み取り < Q < 7、長さ 200 <コマンドを実行して、コマンドを実行して読み取 り |> 読み取り > NB01 _trimmed.fastq.gz ステップ 7.2.6 で生成されたトリミングされたファイルに対して、次のコマンドを実行してNano _trimmed Stat を実行します.gz。 出力をテキストまたは Word ファイルにコピーして結果を記録し、手順 7.2.4 の結果と比較して、フィルタが正常に実行されたことを確認します (表 1)。 シーケンス実行で使用される各バーコードについて、手順 7.2.2 ~ 7.2.8 を繰り返します。注: 手順 7.2.6 で生成された NB01_trimmed.fastq.gz ファイルは、ハイブリッド アセンブリに使用されます。 8. ハイブリッドゲノムアセンブリの生成 注: 次のアセンブリ パイプラインは、セクション 7.1 と 7.2 で準備された短い読み取りと長い読み取りを組み合わせるために Unicycler19、28、29、30を利用します ( 図 3 ) 。ユニサイクラーとその依存関係をインストールし、以下のコマンドを実行します。ステップ 7.1.5 でエクスポートされた短読ファイルは、trimmed_short_file名前が付けられたと見なされます。R1.fastqとtrimmed_short_file。簡単にするために R2.fastq. 短読み取りファイルと長時間読み取りファイルを、Trimmed_Readsという名前の単一のディレクトリに整理します。ディレクトリには、次の情報が含まれている必要があります。 トリムされた長い読み取り用の .fastq.gz ファイル (ステップ 7.2.6 で生成されます)。 2 つの .fastq ファイル (R1 と R2) をトリミングされた短い読み取り (ステップ 7.1.5 で生成されます)。 ターミナルで cd コマンドを使用して読み取りファイルを格納するディレクトリTrimmed_Readsに移動します: cd デスクトップ/Trimmed_Reads 正しいディレクトリに入ったら、2 つの短い読み取りファイルを zip して 、.fastq.gz 形式で次のコマンドを実行して 、gzip trimmed_short_fileを実行します。R1.ファストク R1 と R2 の両方について、ステップ 8.2 を繰り返します。すべての読み取りファイルが .fastq.gz 形式になっていることを確認し、すべてのファイルが同じ分離ファイルと一致することを確認します。 次のコマンドを実行して、ユニサイクラーを使用してハイブリッド アセンブリを開始します。一輪車 -1 trimmed_short_file。R1.fastq.gz -2 trimmed_short_file。R2.fastq.gz -l NB01 _trimmed.fastq.gz -o unicycler_output_directory注: -o は、ユニサイクラー出力が保存されるディレクトリを指定します。事前にディレクトリを生成しないでください。実行時間は、使用するコンピュータの計算能力と、ゲノムサイズと読み取り数によって異なります。これは4時間から1日または2日の間どこでもかかるかもしれません。このプロトコルは、250 Gb RAM、インテル Xeon (R) CPU 2.5 GHz 12 実用的コアと 48 仮想コアを搭載した CentOS Linux 7 マシンで実行されました。また、16 Gb RAM と 2.6 GHz 6 コア プロセッサを搭載したパーソナル コンピュータは、より長い処理時間でこれらのアセンブリを計算できます。 実行が完了したら、一輪.logファイルを確認してエラーがないことを確認します – 生成されたコンティグの数、サイズ、ステータス(完全、不完全)を記録します。 不完全なコンティグが識別された場合 (ユニサイクラーログで不完全と表示される)、ステップ 8.4: –mode 太字のコマンドに次のフラグを追加して、Unicycler を太字モードで再実行します。注: 太字モードでは、アセンブリ中に長い読み取りブリッジに許容される品質しきい値が下がります。これにより、完全なアセンブリが生成される可能性がありますが、アセンブリの品質が低下する可能性があります。必要な場合にのみ、および後で PCR で確認されるコンティグ結合の予備の証拠として、太字モードを使用することをお勧めします。 9. 組み立て品質の評価 メモ:次のプロトコルは、Bandage31と QUAST32を使用する前にセットアップする必要がある 2 つのプログラムを使用します(図 2および図4)。包帯は一度ダウンロードしたインストールを必要とせず、QUASTは基本的なコマンドラインの使用法に精通している必要があります。また、ベンチマークユニバーサルシングルコピーオルソログ(BUSCO)33を使用してゲノムの完全性を評価することをお勧めします。 包帯:ファイルをクリックします。次に、[グラフをロード]を選択し、ステップ8.4でユニサイクラーによって生成されたunicycler_output_directoryに保存されたassembly.gfaファイルを選択します。ロードしたら、左側のツールバーの[グラフの描画]ボタンをクリックし、コンティグ(ノードと呼ばれる)がどのように接続され、アセンブリが完成しているかを評価するために編成されている様子を確認します(図 5)。注: 完全なアセンブリは、両端にリンクされた単一の循環コンティグで表されます(図5A、B)。不完全なアセンブリには、複数のコンティグがリンクされているか、線形である (図 5C)。小さい線形共流は、線形染色体外要素を示す可能性があるため、不完全でない場合があります。また、深さと呼ばれるカバレッジは包帯で注目され、11cyclerで1xに正規化された染色体に対するコンティグの相対的な豊富さを表します。 クアスト ターミナル内で、cdコマンドを使用してユニサイクラー出力を保存するフォルダに移動 unicycler_output_directory Trimmed_Readsします。注: アセンブリが配置されているパスにスペースを入れることができない、つまり、ユニサイクラー出力に通じるディレクトリは、名前にスペースを含めることはできません。または、簡単にアクセスするために、assembly.fasta ファイルをデスクトップにコピーします。 次のコマンドを実行して QUAST を実行quast_output_directoryします。 出力ディレクトリquast_output_directoryで、QUAST によって生成されたレポートを確認します。 10. ゲノムアノテーション 注: 以下のアノテーション パイプラインでは、使用前にインストールする必要があるコマンド ライン ツール Prokka34を使用しています。あるいは、自動化されたGUI K-Base(材料表)を介してProkkaを使用するか、WebサーバーRAST35を介してゲノムにアナプターを付ける。NCBIにゲノムを堆積させる場合、原核生物ゲノムアノテーションパイプライン(PGAP)36を使用して自動的に注釈が付けられます。 ターミナル内で、cdコマンドを使用してユニサイクラー出力を保存するフォルダに移動します(ステップ 9.2.1 を参照)。次に、prokka –プレフィックスsample_ID –outdir prokka_output_directory assembly.fasta コマンドを実行してProkkaを実行します。注: –prefix は指定されたsample_IDに基づいてすべての出力ファイルに名前を付けます。Prokka の出力ディレクトリは事前に作成しないでください。 注釈を確認するには、.tsvテーブルを開くか、生成された.gffファイルをシーケンス分析ソフトウェアにアップロードして、注釈を視覚化して分析します (図 6)。 特定のタイプの注釈は、関心のある遺伝的要因に応じて生成される可能性があります。予備分析37、38 、39、40、41の予備分析のためのゲノム疫学 (www.genomicepidemiology.org/) Web サーバーのユーザー フレンドリー ツールから始めることをお勧めします。CRISPR-casシステムおよびプロファージの検出のための追加ツールが利用可能です (図 3)42,43. 11. データの民主化に関する推奨慣行 可能であれば、生の読み取りデータと組み立てられたゲノムを、NCBIシーケンス読み取りアーカイブ(SRA)やGenbankなどのパブリックリポジトリに保存します。NCBI堆積プロセス中に、ゲノムはPGAPパイプラインを介して自動的にアセプトされます。