奈良女子大集中讲义-bic.kyoto.ppt

奈良女子大集中講義バイオインフォマティクス (7)進化系統樹 阿久津　達也 京都大学　化学研究所 バイオインフォマティクスセンター 講義予定 9月5日 分子生物学概観 分子生物学データベース 配列アラインメント 実習１（データベース検索と配列アラインメント） 9月6日 モチーフ発見 隠れマルコフモデル カーネル法 進化系統樹推定 9月7日 タンパク質立体構造予測 相互作用推定 スケールフリーネットワーク 実習２（構造予測） 進化系統樹推定 系統樹の個数 距離行列法 UPGMA法 近隣結合法 最節約法 進化の確率モデル 最尤法 進化系統樹 種間（もしくは遺伝子間）の進化の関係を表す木 以前は形態的特徴をもとに構成 現在は配列情報をもとに構成 （イメージ図：　実際の進化とは直接、関係無い） 有根系統樹と無根系統樹 節点（頂点）：　木の葉、もしくは、枝分かれ部分 枝（辺）：　種の親子関係を表す 枝長：　進化に要した時間に対応 葉：　現在の種に相当する節点 根：　進化の起点にあたる節点 有根系統樹（根つき木）：　根のある系統樹（下図の左） 無根系統樹（根なし木）：　根のない系統樹（下図の右） 系統樹の同型性 同型な系統樹：本質的に同じつながり方（グラフとして同型）をした系統樹 (a)と(b)：同型 　　(a)と(c)：非同型 (d)と(e)：同型 　　(d)と(f)：非同型 系統樹の個数 葉の個数と節点の個数の関係（葉の個数 ＝ n） 有根系統樹：　2n-1 個の節点、　2n-2 本の枝 無根系統樹：　2n-2 個の節点、　2n-3 本の枝 無根系統樹と有根系統樹の関係 値は無根系統樹の任意の枝に設定可能 　　?　有根系統樹の個数＝(2n-3)×無根系統樹の個数 枝が 2n-3 本の無根系統樹に1本の枝を追加 　　?　2n-3 種類の異なる無根系統樹 #葉 #枝 （無根系統樹） #無根系統樹 #有値系統樹 3 3 1 3 4 5 3 3×5 5 7 3×5 3×5×7 6 9 3×5×7 3×5×7×9 … … … … n 2n-3 (2n-5)!! (2n-3)!! 系統樹の個数（例） 値は無根系統樹の任意の枝に設定可能 　　　?　有根系統樹の個数＝(2n-3)×無根系統樹の個数 枝が 2n-3 本の無根系統樹に1本の枝を追加 　　　?　2n-3 種類の異なる無根系統樹 距離行列法 最初に、種（葉）の間の距離をすべて計算 その距離をもとに系統樹を推定 クラスタリング手法の一種 UPGMA法 距離最短のクラスタをみつけては結合 近隣結合法 近隣となる葉をみつけては結合 アルゴリズム 各配列のみからなるクラスタを作成 クラスタが２個になるまで以下を繰り返す クラスタ間距離が最小のクラスタどうしを併合 新しくできたクラスタと他のクラスタ間の距離を計算 クラスタ間距離 距離の性質 UPGMA法(Unweighted pair group method using arithmetic averages) UPGMA法（例1） 分子時計 分子時計＝進化速度一定性の仮定 枝長＝分子時計により刻まれた時間 分子時計仮説が成立　?　 　「任意の葉までの枝長の和が等しい」が任意節点について成立 　?UPGMA法は系統樹を正しく再構成 　　　　　以下の例で、(a) は仮説を満たさないが、(b)は満たす 近隣結合法（１） UPGMA法では以下の例(a)で、2と3が最初に選ばれたため、正しい系統樹を再構成できなかった 最初に、3と4、もしくは、1と2が選ばれれば良い つまり、兄弟関係（近隣関係）にある葉が選ばれれば良い 近隣結合法では（距離が加法性を満たすという仮定のもとで） 常に近隣関係にある葉を選ぶ 近隣結合法（2） 近隣結合法 距離が加法性を満たせば系統樹を正しく再構成 ただし、計算されるのは無根系統樹 加法性：　節点間の距離　＝　節点間を結ぶパスの枝長の和 アウトグループ 近隣結合法：　無根系統樹しか計算されない 無根系統樹における根の決定法 最も長い葉の間の中点を根とする 対象とする種と離れていることが明らかな種（アウトグループ）を追加して系統樹を作成し、アウトグループに接続する接点を根とする 最節約法（１） 最小置換で各配列を生成する系統樹を計算 ２種類の探索が必要 系統樹の形状（トポロジー)　　?　難しい 内部節点への配列の割り当て　?　動的計画法で計算可能 最節約法（２） Sk(a) : 節点 k に塩基 a を割り当てた時のコスト S(a,b) : a を b に置換するためのコスト 動的計画法により葉から根という順で Sk(a) を計算 配列割り当ては各位置ごとに独立に計算可能 k　が内部節点の場合（i,j