病原性細菌の検出には、細菌病原体を迅速に同定する必要がある。このために、通常、病原体は単離され、PCRや全ゲノム配列が行われる。分子タイピングの主な目標の1つは、病原体をクローン群に分類することである。なぜなら、同じ種の系統は宿主に対して大きく異なる効果を持つからである。よく知られている例は、大腸菌O157:H7(Tu、He&Zhou、2014)および大腸菌EC958(Petty et al、2014)のようないくつかの株である。分類するために、MLST(Maiden、2006)またはクローン特異的なマーカーが使用されている(Inouye et al, 2014)。KvarQ(Steiner et al, 2014)、Mykrobe(Bradleyら, 2015)およびSRST2(Inouyeら, 2014)などのWGSデータから直接に関連する突然変異および対立遺伝子を検出することができるいくつかのアプローチが開発されているが、大部分の場合、十分なカバレッジおよび特化したアレルデータベースが必要である(例えばMykrobeはMycobacterium tuberculosisおよびStaphylococcus aureusの同定にのみ使用できる)。そのようなプログラムの使用は、株の同定プロセスを複雑にする。
クローン特異的マーカーを探す代わりに、完全なバクテリアゲノムを参照配列として使用することができる。 Kraken(Wood&Salzberg、2014)またはCLARK(Ounit et al、2015)のようなk-mer長の配列の検出に基づくバクテリア同定プログラムは、RefSeqバクテリアゲノムデータベース全体を使用し、 さらに、それぞれのリードを個別に分類するため、低カバレッジWGSサンプルも処理できる。 Sigma(Ahn、Chai&Pan、2015)のようなアライメントベースのツールと比較して、k-merベースのプログラムは、特に実行時間を考慮すると優れていることが示されている(Lindgreen、Adair&Gardner、2016; Peabody et al、2015 )。 KrakenはNCBIのtaxonomy treeを使用して、リードを別々に識別し、ツリー上の各taxonへのヒットをカウントし、最も多くのヒット数を持つブランチを見つける(kraken紹介)。
StrainSeekerは、異なる分類学レベルで共有されるk-merの数に基づいて、バクテリアの分離株をWGSデータから直接クローンまたはクレードに分類する。NCBI分類法のような既存のtaxonomyシステムに結びついていないリファレンスバクテリアゲノムと系統レベルとの間の系統学的関係を近似するためのガイドツリーを使用するので、大腸菌およびShigella sp. 間のような論争を避けるのに役立つ(ガイドツリーは、ユーザーが提供する必要がある)。
マニュアル
http://bioinfo.ut.ee/strainseeker/index.php?r=site/page&view=manual
webツール
インストール
cent OSに導入した。
依存
- PERL
- R
Builer
- GenomeTester4(GlistMaker、GlistCompare、GlistQuery)
Seeker
- GenomeTester4(GlistMaker、GlistCompare、GlistQuery、GDistribution)
- R scripts for statistical tests
GenomeTester4のインストールは以前紹介しています(リンク)。
本体& helper scirptの ダウンロード
http://bioinfo.ut.ee/strainseeker/index.php?r=site/page&view=downloadable
本体のbuilder.pl、seeker.plと解凍したhelper scirptを同じディレクトリに置いておく。
> perl builder.pl -h
$ perl builder.pl -h
Usage: builder.pl -n <NWK FILE> -d <DIR OF FASTA FILES> -o <USER DEFINED DB NAME> [OPTIONAL PARAMETERS]
Options:
-h, --help - Print this help
-v, --version - Print version of the program
-n, --newick - Guide tree in newick format (same names as fasta files without suffix .fna)
-d, --dir - Directory of fasta files (.fna)
-o, --output - User defined database name
-b, --blacklist - .list file of k-mers unwanted in database (human, plasmids etc)
-w, --word - K-mer length used in database building and later searching (default 32)
-m, --min - Minimal amout of k-mers in node to be considered as subroot (default 250)
-g, --greater - Maximum times child could have more k-mers than parent (default 250)
-t, --threads - Number of cores used
-max - Maximum number of k-mers in one list (default 100000)
> perl seeker.pl -h
$ perl seeker.pl -h
Usage: seeker.pl -d <DB DIR NAME> -i <SAMPLE.fastq> [OPTIONAL PARAMETERS]
Options:
-h, --help - Print this help
-v, --version - Print version of the program
-i, none - Input file (can be multiple, each with own flag)
-o, --output - Output file name (default StrainSeeker_output)
-d, --dir - Path to database directory
-verbose - Print out more of the working process
ラン
1、データベースの構築(300GBくらい空き容量が必要)
perl builder.pl -n refseq_guide_tree.nwk -d strain_fasta_directory -w 32 -b ss_blacklist_w32.list -o my_database
- -n is the guide tree in Newick format, describing the relationships between given strains.
- -d is a directory containing all the .fna files for strains used in the Newick file.
- -b is the path to blacklist (must have the same k-mer length as parameter -w).
- -w is the k-mer length.
- -o user-defined database name.
ここでは既存のデータベースを使う。ダウンロードリンク(リンク)から32-merのDatabase - pre-built from 4,324 NCBI RefSeq strains (w32) か16-merのDatabase - pre-built from 4,324 NCBI RefSeq strains (w16)をダウンロードする。
シーケンスデータは、オーサーらが評価ペーパー(リンク)にしたがって準備したfastqを使う。このリンクからダウンロードできる。
検索する。
perl seeker.pl -i sample.fastq -d ss_db_w32 -o sample_result.txt -verbose
出力の確認(一部)。
> cat sample_result.txt
$ cat sample_result.txt
Sample:sample_result.txt
22.66893% KNOWN Streptococcus_mitis_B6
17.53777% KNOWN Streptococcus_pseudopneumoniae_IS7493
17.02426% KNOWN Streptococcus_pneumoniae_Hungary19A-6
10.73490% RELATED Lactobacillus_helveticus_strain_KLDS18701,Lactobacillus_kefiranofaciens_ZW3,Lactobacillus_helveticus_strain_CAUH18,Lactobacillus_helveticus_H10,Lactobacillus_helveticus_DPC_4571,Lactobacillus_helveticus_CNRZ32,Lactobacillus_helveticus_R0052,Lactobacillus_helveticus_strain_MB2-1,Lactobacillus_helveticus_H9
500MBまでの限定だが、web版もある。
web版(4300のbacteriaデータベース)
http://bioinfo.ut.ee/strainseeker/index.php?r=site/webtool
テスト結果
http://bioinfo.ut.ee/strainseeker/index.php?r=site/results&id=demoresults
http://bioinfo.ut.ee/strainseeker/demo/demoresults/setA1_bacarc_5feb_11min.txt
上のデータを解析すると以下のような結果が得られた。
引用
StrainSeeker: fast identification of bacterial strains from raw sequencing reads using user-provided guide trees.
Roosaare M1, Vaher M1, Kaplinski L1, Möls M1,2, Andreson R1, Lepamets M1, Kõressaar T1, Naaber P3,4, Kõljalg S4,5, Remm M1.
PeerJ. 2017 May 18;5:e3353.