差分
このページの2つのバージョン間の差分を表示します。
両方とも前のリビジョン 前のリビジョン 次のリビジョン | 前のリビジョン | ||
2022-メタゲノム・edna [2022/11/09 14:27] – [B. クオリティチェック] suikou | 2022-メタゲノム・edna [2022/11/10 04:34] (現在) – suikou | ||
---|---|---|---|
行 15: | 行 15: | ||
大雑把なデータ解析の流れとしては、 まずシーケンスデータの品質チェックをFastQCというツールを用いて行います。それから、コンピュータにインストールしたBLASTを使って、すべてのシーケンスデータをデータベースと照合させ、各リードがどの種と相同性があるかを調べます。それからMEGANを使ってBLASTのヒットを集計し、各サンプルに含まれていたバクテリアの種類を網羅的に解析します。 | 大雑把なデータ解析の流れとしては、 まずシーケンスデータの品質チェックをFastQCというツールを用いて行います。それから、コンピュータにインストールしたBLASTを使って、すべてのシーケンスデータをデータベースと照合させ、各リードがどの種と相同性があるかを調べます。それからMEGANを使ってBLASTのヒットを集計し、各サンプルに含まれていたバクテリアの種類を網羅的に解析します。 | ||
- | |||
- | 環境DNAのサンプルは、MEGANの結果で検出された種と、実際に池にいそうな魚なのかの考察や、使用したプライマーの種類(12S全長 or 12S MiFish)の違いが結果に与える影響などを考えてみてください。メタゲノムのほうでは、食品のデータはシャッフルしてお渡ししますので、含まれるバクテリアの種類から、食品サンプルが上記4つの食品のどれであるかを推定します。また、池の水と水槽の水で微生物叢にそれぞれ特徴があるかを考察してください。 | ||
## サンプルの説明 | ## サンプルの説明 | ||
行 50: | 行 48: | ||
[[http:// | [[http:// | ||
+ | |||
+ | クオリティスコアQは、エラーの生じる確率 perror から下記のように計算されます。 (出典: https:// | ||
+ | |||
+ | {{: | ||
+ | |||
+ | - クオリティスコアが 10 ならば、シーケンシングエラーが生じる確率は 10.0% であるから、読み取られた塩基の信頼度は 90.0% 。 | ||
+ | - クオリティスコアが 20 ならば、シーケンシングエラーが生じる確率は 1.0% であるから、読み取られた塩基の信頼度は 99.0% 。 | ||
+ | - クオリティスコアが 30 ならば、シーケンシングエラーが生じる確率は 0.1% であるから、読み取られた塩基の信頼度は 99.9% 。 | ||
行 109: | 行 115: | ||
| | ||
- | |||
### 2. データベースのダウンロード | ### 2. データベースのダウンロード | ||
行 118: | 行 123: | ||
| | ||
+ | |||
+ | ``` | ||
+ | # | ||
+ | tar vxf 16S_ribosomal_RNA.tar.gz | ||
+ | ``` | ||
- MitoFish ・・・魚類のミトコンドリアのデータベース。約2万種が登録されている。魚類ミトコンドリア上の16S rRNA、12S rRNAを探す場合に使用する。 | - MitoFish ・・・魚類のミトコンドリアのデータベース。約2万種が登録されている。魚類ミトコンドリア上の16S rRNA、12S rRNAを探す場合に使用する。 | ||
行 145: | 行 155: | ||
## B. クオリティチェック | ## B. クオリティチェック | ||
- | FastQCは主にIllumina用のクオリティチェックツールなので、Nanoporeのデータに対しては適切な評価ができておらず、評価値の〇×は気にしなくてよいです。下記はシーケンスデータのクオリティスコアに関する、平均値等の情報。 | + | FastQCは主にIllumina用のクオリティチェックツールなので、Nanoporeのデータに対しては適切な評価ができておらず、評価値の〇×は気にしなくてよいです。 |
- | Nanoporeのデータだとクオリティスコア10強(精度90%強)となるはずである。ただし、ナノポアのクオリティスコアは実際のスコアとは乖離があるようです。(ナノポアでQ20と出ていても、実際はQ12.5程度。https:// | ||
- | |||
- | リード長の分布が想定される長さになっているかなども確認すること。 (16S全長は1.5 kbp程度、12S全長は1 kbp程度、12S MiFishは250 bp程度) | ||
### FASTQCを実行するには… | ### FASTQCを実行するには… | ||
行 178: | 行 185: | ||
./fastqc | ./fastqc | ||
``` | ``` | ||
- | |||
### FastQCでシーケンスファイルを開く | ### FastQCでシーケンスファイルを開く | ||
{{: | {{: | ||
- | Ctrlキーを押しながらクリックすれば、複数のファイルを一度に選択できるので、必要なファイルを選択して開く。 | + | Ctrlキーを押しながらクリックすれば、複数のファイルを一度に選択できるので、必要なファイルを選択して開く。(ただし一度に複数選択するとメモリーエラーで落ちたりするので、その場合は一つずつ開くこと。) |
{{: | {{: | ||
行 191: | 行 197: | ||
{{: | {{: | ||
+ | Nanoporeのデータだとクオリティスコア10強(精度90%強)となるはずですが、リード全体のクオリティは比較的正しく計算できるみたいだけど、塩基ごとのクオリティはあまり正確ではなく、実際の塩基精度とは乖離があるようです。(ナノポアでQ20と出ていても、実際はQ12.5程度。https:// | ||
+ | |||
+ | リード長の分布が想定される長さになっているかなども確認すること。 (バクテリア16S全長は1.5 kbp程度、魚類ミトコンドリア16Sは600 bp程度、魚類ミトコンドリア12S MiFishは200 bp強) | ||
## C. BLAST | ## C. BLAST | ||
行 213: | 行 222: | ||
```input_file.fastq```, | ```input_file.fastq```, | ||
- | |||
### コマンド入力時の便利キー | ### コマンド入力時の便利キー | ||
|カーソルの上・下キー|前に入力したコマンドを呼び出す。| | |カーソルの上・下キー|前に入力したコマンドを呼び出す。| | ||
- | |Ctrl+Shift+V|貼り付け| | + | |(Win) Ctrl+Shift+V or 右クリック, |
|tabキー|ファイル名・コマンドの自動補完| | |tabキー|ファイル名・コマンドの自動補完| | ||
|Ctrl+C|コマンド強制終了(blastを実行中に止めたい場合など)| | |Ctrl+C|コマンド強制終了(blastを実行中に止めたい場合など)| | ||
行 279: | 行 287: | ||
{{: | {{: | ||
- | |||
## D. MEGANによるBLAST結果集計 | ## D. MEGANによるBLAST結果集計 | ||
行 337: | 行 344: | ||
{{: | {{: | ||
- | 13.Exportするデータを「taxonPath_to_count」に変更します。(ほかのデータ形式でも勿論可) | + | 13.Exportするデータを「`taxonPath_to_count`」に変更します。(ほかのデータ形式でも勿論可) |
{{: | {{: | ||
行 344: | 行 351: | ||
{{: | {{: | ||
- | |||
## E. データの転送、Excelでの解析 | ## E. データの転送、Excelでの解析 | ||
- | |||
- | リモートデスクトップ先で「ファイル管理」を開いて、保存したファイルをコピーしてから、手元のWindowsでファイルエクスプローラーを開いて貼り付けるとファイルを簡単に転送できる。もしくは、メールやGoogle Driveなどで転送するなどしてもよい。 | ||
ExportしたファイルをExcelで開くには、Excelを起動しておき、ExportしたファイルをExcel上にドラッグアンドドロップすれば良いです。 | ExportしたファイルをExcelで開くには、Excelを起動しておき、ExportしたファイルをExcel上にドラッグアンドドロップすれば良いです。 | ||
行 358: | 行 362: | ||
そのほか、「データ」→「フィルター」を使ってみたり、グラフを描いてみたりするのが通常の解析の流れになるかと思います。 | そのほか、「データ」→「フィルター」を使ってみたり、グラフを描いてみたりするのが通常の解析の流れになるかと思います。 | ||
- | |||
## F. 明日の内容 | ## F. 明日の内容 | ||
行 384: | 行 387: | ||
検出された魚は三四郎池に棲息していそうな魚かどうか。 | 検出された魚は三四郎池に棲息していそうな魚かどうか。 | ||
- | 去年の三四郎池のデータとも比較してみる。http:// | + | 去年の三四郎池のデータとも比較してみる。http:// |
手法で詳しく説明して欲しい箇所:「電気泳動、DNA精製」(使用したキットはFastGene™ Gel/ | 手法で詳しく説明して欲しい箇所:「電気泳動、DNA精製」(使用したキットはFastGene™ Gel/ | ||
行 392: | 行 395: | ||
今回発酵食品で検出されるバクテリアはほぼ1種類だと思うので、精度の悪いナノポアのリードの精度を向上させる方法を実践してみてください。具体的にはGeneiousでマルチプルアライメントを作成して、コンセンサス配列を作ることで、NCBIのBlastで一致率99%程度のヒットが得られるようになることを確認し、ナノポアのリードはどのような間違いが多いのか考察してみてください。 | 今回発酵食品で検出されるバクテリアはほぼ1種類だと思うので、精度の悪いナノポアのリードの精度を向上させる方法を実践してみてください。具体的にはGeneiousでマルチプルアライメントを作成して、コンセンサス配列を作ることで、NCBIのBlastで一致率99%程度のヒットが得られるようになることを確認し、ナノポアのリードはどのような間違いが多いのか考察してみてください。 | ||
- | 加工食品から検出された魚は妥当でしょうか。No Hitのリードを抜き出して、NCBIのBLASTにかけてみると何がヒットしますか? | + | 加工食品から検出された魚は妥当でしょうか。No Hitのリードを抜き出して、NCBIのBLASTにかけてみると何がヒットしますか?(リードを抜き出す例: `./seqkit grep -rp " |
手法で詳しく説明して欲しい箇所:「PCR」(使用したDNAポリメラーゼはrepliQa HiFi ToughMixです。AmpliTaq GoldやEx Taqといった他の酵素と比較して、どういった特徴があるでしょうか。) | 手法で詳しく説明して欲しい箇所:「PCR」(使用したDNAポリメラーゼはrepliQa HiFi ToughMixです。AmpliTaq GoldやEx Taqといった他の酵素と比較して、どういった特徴があるでしょうか。) |