7
4

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?

More than 5 years have passed since last update.

Biopyton からClustalW2を使って系統樹を作成する

Last updated at Posted at 2016-11-14

Biopyton からClustalW2を使って系統樹を作成する方法について自分用にメモ.
といっても内容はほとんどBiopython Tutorial and Cookbookに書かれていることを日本語に訳しただけ.

ClustalW2 をインストール

まず,ClustalW2 をダウンロードする.Macの場合には.dmgをマウントして得られるbinファイルを/bin下に置く.

スクリーンショット 2016-11-13 23.59.33.png

系統樹を作成したい株を用意する

次にClustalW2を使用する株のデータを用意する.今回はメタリジウム(Metarhizium)の系統樹をRibosome biogenesis protein YTM1をもとに作成する.ファイルはUniProtからダウンロードした.
使用した株は次の通り.

  • Metarhizium acridum (strain CQMa 102)
  • Metarhizium anisopliae BRIP 53293
  • Metarhizium guizhouense ARSEF 977
  • Metarhizium album ARSEF 1941
  • Metarhizium robertsii
  • Metarhizium rileyi RCEF 4871

Basketに追加した後,FASTA形式でダウンロードする.今回はuniprot-yourlist.fastaという名前で保存した.

BiopythonでClustalW2を実行する

用意した株のデータに対してBiopythonからClustalW2をかける.

from Bio.Align.Applications import ClustalwCommandline

clustalw_cline = ClustalwCommandline("clustalw2", infile="uniprot-yourlist.fasta")
stdout, stderr = clustalw_cline()

するとuniprot-yourlist.alnとuniprot-yourlist.dndの2つのファイルが生成されている.そこでBiopythonのPhyloモジュールを使ってdndファイルの方を読み込み,系統樹を描く.

from Bio import Phylo

tree = Phylo.read("uniprot-yourlist.dnd", "newick")
Phylo.draw(tree)

figure_1.png

draw関数の代わりにdraw_ascii関数を使うと系統樹をアスキーアートで出力してくれる.

  _ tr|E9E7T1|E9E7T1_METAQ
 |
 | , tr|A0A0D9P3B0|A0A0D9P3B0_METAN
 |,|
_||| tr|A0A0A1USL4|A0A0A1USL4_9HYPO
 ||
 || tr|A0A0B4H3C6|A0A0B4H3C6_9HYPO
 |
 |      ______________________________________ tr|A0A0B2X7N3|A0A0B2X7N3_9HYPO
 |_____|
       |______________ tr|A0A167BRY5|A0A167BRY5_9HYPO

Exited with code=0 in 1.1

参考文献

Peter J. A. Cock, Tiago Antao, Jeffrey T. Chang, Brad A. Chapman, Cymon J. Cox, Andrew Dalke, Iddo Friedberg, Thomas Hamelryck, Frank Kauff, Bartek Wilczynski, Michiel J. L. de Hoon: “Biopython: freely available Python tools for computational molecular biology and bioinformatics”. Bioinformatics 25 (11), 1422–1423 (2009). doi:10.1093/bioinformatics/btp163,

Eric Talevich, Brandon M. Invergo, Peter J.A. Cock, Brad A. Chapman: “Bio.Phylo: A unified toolkit for processing, analyzing and visualizing phylogenetic trees in Biopython”. BMC Bioinformatics 13: 209 (2012). doi:10.1186/1471-2105-13-209

7
4
0

Register as a new user and use Qiita more conveniently

  1. You get articles that match your needs
  2. You can efficiently read back useful information
  3. You can use dark theme
What you can do with signing up
7
4

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?