最(大)節約法: Maximum Parsimony Method

準備:PHYLIPのダウンロード

予備知識確認テスト

授業を開始する前に、次の問いに答えてもらいます。

Table 1.

ex1.jpg

系統樹の表記法(ニューイックフォーマット, Newic Formst)

下の系統樹は (M(N(A(B C)))) と表記する。

ex3.jpg

最節約法:6月19日授業課題

1.ホームページからデータファイル(&ref(): File not found: "rbcL_plant.fst" at page "授業/H22/系統学特論/05";)をダウンロード。BioEditで開き、NEXUSファイルに変換。
Download(&ref(): File not found: "rbcL_plant.fst" at page "授業/H22/系統学特論/05";) and open it by BioEdit. Save the aligned data as a NEXSUS file.

2.PAUP* DEMOでNEXUSファイルを開き、Heuristic Search法で解析する。得られた最短系統樹をファイルに保存する。
Open the above NEXUS file with PAUP* DEMO, and analyze it by "Heuristic Search" method. Two shortest trees will be obtained. Save them in a tree file.

ファイル名は  学籍番号(半角英数文字のみ).emf
file name  "Your_ID.emf"

梶田までメールの添付書類として送信。件名は「6月19日課題」にすること。
Send the file to Kajita by e-mail.

提出期限:6月28日
Deadline: 28 June.

6月19日授業の課題はここまでにしましたが、できるなら3以下をやって提出してもらっても良いです。
-------------------

3.次に、Bootstrap解析(Heuristic search, 1000回〔どうしてもできなければ100回でよい〕)を行う。得られた結果の表は保存しておく。〔ここまでの操作にはPAUP Blockを用いても良い〕
Do bootstrap analysis for 1,000 replicates. Keep the obtained bootstrap values of each clades for later use.

4.メモ帳等のテキストエディタで2で作った最短系統樹のファイルPedic.treを開き、2つ目の系統樹のクレードに、3の解析結果を見ながらBootstrap値を記入して保存する。

Open the tree file (TEXT file format) saved in the strep 2, and add bootstrap values for each clades.

5.上で保存したファイルをTreeViewで開き、 2つ目の系統樹をPhylogramモードで、Bootstrap値と共に表示させる(show internal edge label)。内群が単系統群になるようにルートをつける。
Open the above file with TreeView, and display the second tree in Phylogram view with boot strap values ("show internal edge label"). Define a root and make the ingroup taxa as an monophyletic group.

6.TreeViewで、File > Save as graphic を選び、enhanced metafileの形式で2つ目の系統樹のみを保存する。このファイルをmメールの添付書類として梶田宛に提出する。
Choose "Save as graphic" in the File menu, and save only the second tree as "enhanced metafile" format. Send the file by e-mail to Kajita as an attachment file.

ファイル名は 学籍番号(半角英数文字のみ).emf
The file name should be "Your Student ID.emf"


授業/H17/系統学特論/PAUPブロック

最節約法の基礎

最節約法を基本から勉強したいという人には、次の参考書(なんと無料!)をお勧めします。1コマの講義で節約法の全貌と、ソフトウェアを使った系統推定を理解するのは、初めての人には大変だと思います。以下の文献よ読んで、復習してください。

The Compleat Cladist (日本語訳が「系統分類学入門」というタイトルで、文一総合出版から出ている)http://taxonomy.zoology.gla.ac.uk/teaching/CompleatCladist.pdf

Basics of Cladistic Analysis (上よりもさらに簡単な内容)http://www.gwu.edu/~clade/faculty/lipscomb/Cladistics.pdf

BioEditでNEXUSファイルを作る

1. 例題ファイルfilePedic_align.fstをBioEditで開く

2. BioEdit: File > Export > Sequence Alignment

スライド0.jpg
スライド1.JPG

PAUP*を使う

1. File > Open : Pedic.nex を Edit モードで開く

スライド2.JPG
スライド3.JPG

NEXUSファイルの中身。[ ] に挟まれた部分はコメント。

スライド4.JPG

〔参考〕FASTA形式からNEXUS形式には簡単に変更できる

2. データをメモリに読み込む File > Execute

スライド5.JPG

3. Branch and Boundで解析:BandB

スライド6.JPG
スライド7.JPG

4. 画面に系統樹を表示: showtree 1

スライド8.JPG

5. ファイルに系統樹を保存: savetrees file=Pedic.tre brlen=yes from=1 to=2

6. TreeView>で系統樹を表示: File > Open

スライド9.JPG

7. PAUP*でBootstrap解析: bootstrap nrep=1000 search=heuristic brlens=yes

スライド10.JPG

8. 系統樹を保存: savetrees file=Pedic_boot.tre savebootp=nodelabels from=1 to=1

スライド11.JPG

9. ディスプレーバッファーを編集: Edit > Edit Display Buffer

スライド12.JPG

10. TreeViewでBootstrap値つき系統樹を開く

スライド13.JPG

11. 最節約系統樹にBootstrap値をのせる:メモ帳で編集

スライド14.JPG

12. TreeViewでBootstrap値つき最節約系統樹を開く

スライド15.JPG

13. PAUP Blockをファイルに埋め込んで、解析を自動化する

スライド17.JPG