PEDを用いた麹菌の多型解析例

• 図をクリックすると拡大されます。
  
• SNP以外に、挿入・欠失変異も正確に検出できる特徴があります。
• 変異部分を増幅するプライマー対の配列も得られます。

Linux OS (Ubuntu)のインストール

• PEDはUnix系のOSで動作するため、ここでは動作が確認されているUbuntuをインストールします。
  Ubuntu Desktop 日本語 Remixより日本語版のISOイメージのファイルをダウンロードします。
• Ubuntu 20.04 LTSインストールガイド【スクリーンショットつき解説】が詳しいので参考にしてインストールしてください。
  PEDは解析時間をできるだけ短くするため、CPUの限界近くまで負荷がかかりますのでほかのOSとの共存は避けた方がいいです。
  使わなくなった古いコンピュータなどにインストールして専用にするのがお勧めです。
  マルチスレッドで動作しますので、4あるいは8コアのCPUを用いると早く処理できます。

PEDのセットアップと解析例

• 必要なソフトをダウンロードとインストールを行います。

  $ sudo apt update
  $ sudo apt upgrade
  $ sudo apt install curl
  $ sudo apt install git
  $ sudo apt install sra-toolkit
  $ git clone https://github.com/akiomiyao/ped.git

• Ubuntuのバージョンによっては、sra-toolkitのインストールに失敗することがあります。
  また最近のバージョンのfastq-dumpは使用前にvdb-configコマンドで実行環境を設定する必要があります。
  設定方法は、[fastq-dumpの設定]のページに記載しました。
• NCBIのSRAアーカイブからGenome sequencing of industrial Asperigillus oryzae koji strainsのページのRunsを開いてfastq形式のファイル(DRR163504)をダウンロードします。
  

  $ cd ped
  $ git pull
  $ perl download.pl accession=DRR163504

  git pullで最新のスクリプトに更新されます。git pullは毎回行う必要はありませんが、プログラムが更新されている場合はpullしてください。
  pedのディレクトリの中に、DRR163504の名称のディレクトリができて、その中のreadというディレクトリの中にfastqファイルがダウンロードされます。
  
• DDBJのDRAからもダウンロードできます。
  DDBJからのDRR163504のダウンロードページ
  2行目のFASTQをクリックしてbz2形式で圧縮されたファイルをダウンロードします。ダウンロードしたファイルを DRR163504/read のディレクトリに置きます。解凍する必要はありません。
  同様の構造の任意の名前のディレクトリを作れば手持ちのデータの解析ができます。
  
• DRR163504の配列を麹菌の参照配列RIB40にマップして多型を検出します。

  $ perl ped.pl target=DRR163504,ref=RIB40

  麹菌の参照配列は設定済みなので自動的にダウンロードして解析用データを作成して一連の作業が自動的に進みます。
  解析が無事終わるとDRR163504のディレクトリの中に、DRR163504.vcfという結果をまとめたファイルができます。
  各座位に対してプライマーが作成可能な場合は、プライマーペアの配列データも出力されます。
• 注意点 IVGでは長い欠失の情報をうまく表示できない場合があるので、80塩基以上の欠失に関しては、vcfファイルに出力しないように設定しています。
  大きな構造変異を調べたい場合は、targetのディレクトリの中に出力された、svの拡張子のついたファイルを直接ご参照ください。
  

Dockerでの解析

• Dockerを用いれば、必要な時に設定済み実行環境を取り込んで解析できます。
  グーグルクラウド(GCP)やアマゾンウェブサービス(AWS)などのクラウドでの実行も可能です。
• Dockerのセットアップ
  

  $ sudo apt install docker
  $ sudo apt install docker.io

• 解析例
  

  $ sudo docker pull akiomiyao/ped
  
  $ sudo docker run -w /ped -v `pwd`:/work akiomiyao/ped perl download.pl accession=DRR163504,wd=/work
  $ sudo docker run -w /ped -v `pwd`:/work akiomiyao/ped perl check_length.pl target=DRR163504,wd=/work
  $ sudo docker run -w /ped -v `pwd`:/work akiomiyao/ped perl ped.pl target=DRR163504,ref=RIB40,wd=/work

  上記のコマンドをコピーしてターミナルウインドウにペーストして実行してください。
  

IVGブラウザによる結果の確認

• Integrative Genomics Viewer (IVG) はBroad Instituteで作成されたゲノムブラウザです。
  http://software.broadinstitute.org/software/igv/download からダウンロードできます。
  RIB40.fastaをクリックしてダウンロードします。
  IVGを開いてGenomesのプルダウンメニューよりLoad Genome from Fileを選んで、RIB40.fastaをロードします。
  Aspergillus_oryzae.ASM18445v3.47.gff3 をクリックしてダウンロードして、メインの表示ボックスにドラッグします。
  DRR163504.vcf(クリックしてダウンロードできます)
  をメインの表示ボックスにドラッグすれば、Aspergillus oryzae TK-1の多型データが取り込まれます。
  各染色体の番号をクリックすると、染色体ごとの多型情報が表示されます。
  興味のある部分をクリックすると、どんどん拡大されます。
  多型をマウスで指すと詳細情報が表示されます。
  多型の種類や多型頻度、また、プライマー対の配列情報と増幅産物のサイズも表示されます。
  ページのトップの画像が表示例です。
  多型の部分を増幅するプライマー配列も表示されます。至適配列が見つからない部位では表示されません。

PEDに関する資料

• スクリプト: https://github.com/akiomiyao/ped/
• 論文: Polymorphic edge detection (PED): two efficient methods of polymorphism detection from next-generation sequencing data (2019) BMC Bioinformatics 20(1):362

お問い合わせ