差分

このページの2つのバージョン間の差分を表示します。

--- 2022-メタゲノム・edna [2022/11/09 14:36] – [FastQCでシーケンスファイルを開く] suikou
+++ 2022-メタゲノム・edna [2022/11/10 04:34] (現在) – suikou
@@ 行 15: / 行 15: @@
 大雑把なデータ解析の流れとしては、 まずシーケンスデータの品質チェックをFastQCというツールを用いて行います。それから、コンピュータにインストールしたBLASTを使って、すべてのシーケンスデータをデータベースと照合させ、各リードがどの種と相同性があるかを調べます。それからMEGANを使ってBLASTのヒットを集計し、各サンプルに含まれていたバクテリアの種類を網羅的に解析します。
-環境DNAのサンプルは、MEGANの結果で検出された種と、実際に池にいそうな魚なのかの考察や、使用したプライマーの種類(12S全長 or 12S MiFish)の違いが結果に与える影響などを考えてみてください。メタゲノムのほうでは、食品のデータはシャッフルしてお渡ししますので、含まれるバクテリアの種類から、食品サンプルが上記4つの食品のどれであるかを推定します。また、池の水と水槽の水で微生物叢にそれぞれ特徴があるかを考察してください。
 ## サンプルの説明
@@ 行 125: / 行 123: @@
  https://ftp.ncbi.nlm.nih.gov/blast/db/16S_ribosomal_RNA.tar.gz (バクテリアの16S rRNAが約2万種登録されている。)
+```
+#ダウンロードしたファイルを解凍する。
+tar vxf 16S_ribosomal_RNA.tar.gz
+```
 - MitoFish ・・・魚類のミトコンドリアのデータベース。約2万種が登録されている。魚類ミトコンドリア上の16S rRNA、12S rRNAを探す場合に使用する。
@@ 行 219: / 行 222: @@
 ```input_file.fastq```, ```output_file.fasta```は適当に変更すること。
 ### コマンド入力時の便利キー
 |カーソルの上・下キー|前に入力したコマンドを呼び出す。|
-|Ctrl＋Shift＋V|貼り付け|
+|(Win) Ctrl＋Shift＋V or 右クリック, (Mac) Cmd＋V|貼り付け|
 |tabキー|ファイル名・コマンドの自動補完|
 |Ctrl+C|コマンド強制終了(blastを実行中に止めたい場合など)|
@@ 行 285: / 行 287: @@
 {{:pasted:20201026-182323.png}}
 ## D. MEGANによるBLAST結果集計
@@ 行 343: / 行 344: @@
 {{:pasted:20201026-182809.png}}
-．Exportするデータを「taxonPath_to_count」に変更します。(ほかのデータ形式でも勿論可)
+．Exportするデータを「`taxonPath_to_count`」に変更します。(ほかのデータ形式でも勿論可)
 {{:pasted:20201026-182821.png}}
@@ 行 350: / 行 351: @@
 {{:pasted:20201026-182832.png}}
 ## E. データの転送、Excelでの解析
-リモートデスクトップ先で「ファイル管理」を開いて、保存したファイルをコピーしてから、手元のWindowsでファイルエクスプローラーを開いて貼り付けるとファイルを簡単に転送できる。もしくは、メールやGoogle Driveなどで転送するなどしてもよい。
 ExportしたファイルをExcelで開くには、Excelを起動しておき、ExportしたファイルをExcel上にドラッグアンドドロップすれば良いです。
@@ 行 364: / 行 362: @@
 そのほか、「データ」→「フィルター」を使ってみたり、グラフを描いてみたりするのが通常の解析の流れになるかと思います。
 ## F. 明日の内容
@@ 行 390: / 行 387: @@
 検出された魚は三四郎池に棲息していそうな魚かどうか。
-去年の三四郎池のデータとも比較してみる。http://suikou.fs.a.u-tokyo.ac.jp/yosh_data/2021jissyu/2021nanopore.zip (group2-12S-Sanshiroike1.fq) 今年はブラックバスが密かに話題になっているらしいが、去年と比べてブラックバスが増えていそうか。
+去年の三四郎池のデータとも比較してみる。http://suikou.fs.a.u-tokyo.ac.jp/yosh_data/2021jissyu/2021nanopore.zip (group2-12S-Sanshiroike1.fq) 今年は三四郎池のブラックバスが密かに話題になっているらしいが、去年と比べてブラックバスが増えていそうか。
 手法で詳しく説明して欲しい箇所：「電気泳動、DNA精製」（使用したキットはFastGene™ Gel/PCRExtractionキットです。）
@@ 行 398: / 行 395: @@
 今回発酵食品で検出されるバクテリアはほぼ1種類だと思うので、精度の悪いナノポアのリードの精度を向上させる方法を実践してみてください。具体的にはGeneiousでマルチプルアライメントを作成して、コンセンサス配列を作ることで、NCBIのBlastで一致率99%程度のヒットが得られるようになることを確認し、ナノポアのリードはどのような間違いが多いのか考察してみてください。
-加工食品から検出された魚は妥当でしょうか。No Hitのリードを抜き出して、NCBIのBLASTにかけてみると何がヒットしますか？
+加工食品から検出された魚は妥当でしょうか。No Hitのリードを抜き出して、NCBIのBLASTにかけてみると何がヒットしますか？（リードを抜き出す例： `./seqkit grep -rp "1e2f600e-220b-4144-93b0-75ccd2c537de:35:1568:-1:bac-27F-BC05:bac-1492R" ./2022nanopore/y2022-group1-bacteria-16S-Sanshiro.fq`）
 手法で詳しく説明して欲しい箇所：「PCR」（使用したDNAポリメラーゼはrepliQa HiFi ToughMixです。AmpliTaq GoldやEx Taqといった他の酵素と比較して、どういった特徴があるでしょうか。）