生命情報解析 第 4 回 シグナル配列の統計解析 (3)

1. 生命情報解析 第4回シグナル配列の統計解析(3) 慶應義塾大学先端生命科学研究所

2. 確率分布と有意性 (1) P = 0.00243815649926 棒グラフの右側部分の面積の合計が確率、すなわち有意性を表す

3. 確率分布と有意性 (2) 確率 • 確率分布をはっきりさせる • 検定対象の値から右側の面積を求める • “こんなにも大きな値”が出る確率が求まる • 有意性の指標として使う 有意性 確率変数が取る実数 検定対象の値

4. Z Scoreの特徴 • 平均が0、分散(データの散らばり)が1になる • 元の分布が正規分布なら、そのZ Scoreは標準正規分布となる • どんな正規分布でも、Z Scoreに直せば同じ土俵(?)で確率計算ができる • Z Scoreが1.96を超える確率は0.025 1.96 0

5. Z Scoreの計算 サイコロを100回振って、90回”6”の目が出るときのZ Scoreは

6. 塩基の方も… • ここでは簡単のため、１塩基の偏りだけを考える • ゲノム全体の塩基組成を考えて、塩基iが対象となる場所において観測される確率はpiとする • 今、N本の配列のうち、 Ni個について、対象となる位置に塩基iが観測された • この条件では通常、Niは正規分布に従う

7. 頻出塩基の統計的有意性 Z Score = N :解析する配列数 Ni :観測された塩基iの数 Bi:ゲノム中における塩基iの割合 Z Scoreは標準正規分布に従う Z Score > 1.96なら、P < 0.05

8. 複数の塩基の有意性を同時に検定するには? • ゲノム全体の塩基組成をBa=0.3, Bc=0.3, Bg=0.3, Bt=0.1として、与えられた位置における塩基iの個数NiがNa = 40, Nc =40, Ng = 10, Nt = 10のとき、偏りは有意か？ • Z Scoreを４つも計算すると… • ４つも値が出て、取り扱いが煩雑になる • 偶然に高い値を示すものが出やすくなる • Χ2値を使う

9. Χ2値 • n個の互いに独立なZ Score：Z1, Z2, Z3, …, Znがあるとき、 Χn2値 = Z12+Z22+Z32+…+Zn2 • Χn2値は自由度nのΧn2分布に従う

10. Χ2分布

11. Χ2分布に従う値を求める • しかし、Na, Nc, Ng, Ntは互いに独立ではないため、これらのZ Scoreを足しても自由度４のΧ2分布には従わない。そこで… • 塩基がそれぞれBiの頻度で出現するとき、上記Χ2値は自由度３のΧ2分布に従う • Χ2値>12.84ならP < 0.005 但しNは解析対象の配列数

12. Χ2計算の例 • ゲノム全体の塩基組成をBa=0.3, Bc=0.3, Bg=0.3, Bt=0.1として、与えられた位置における塩基iの個数NiがNa = 40, Nc =40, Ng = 10, Nt = 10のとき、偏りは有意か？ • この例では偏りが有意とは言えない。

13. 大腸菌開始コドン周辺の塩基のχ2乗値

14. ２つの数式の関係(1) • 4塩基を2種類に分類して考える • プリン(A,G)、ピリミジン(C,T) • ゲノム中のプリン、ピリミジンの頻度をそれぞれBpur, Bpyrとする • 対象となる位置で観測されたプリン、ピリミジンの頻度をそれぞれPpur, Ppyrとする。但し、Ppur + Ppyr = 1

15. ２つの数式の関係(2) • ２種類の塩基の数をもとに計算したΧ2値は自由度1のΧ2分布に従う • 4種類の塩基の数をもとに計算したΧ2値は自由度3のΧ2分布に従う • 自由度は自由に動ける変数の数を意味する

16. 演習問題 • Ba=0.3,Bc=0.2,Bg=0.2,Bt=0.3として、与えられた位置における塩基iの個数Niが • Na = 50, Nc = 30, Ng = 10, Nt = 10 • Na = 500, Nc = 300, Ng = 100, Nt = 100 • のときの増加情報量、χ2値を求めよ。 log23 ≒1.585、log25 ≒2.322

17. 演習問題 解答 • (1), (2)ともに • (1) • (2)

18. 特定のシグナル配列の存在頻度 • 様々な塩基配列の偏りを調べるのではなく、特定のシグナル配列の存在頻度を調べたい (ex. SD配列 “AGG”) • 最も単純なのは、頻度＝あるシグナル配列が観測される配列数÷解析対象の配列数

19. 大腸菌開始コドン周辺の”AGG”の頻度

20. 頻出塩基配列パターンの統計的有意性 Z Score = Nobs:パターン観測数 N:解析する配列数 p:ゲノム中におけるパターンの割合 パターンの出現頻度がpのとき、Z Scoreは標準正規分布に従う Z Score > 1.96なら、P < 0.05

21. 大腸菌開始コドン周辺の”AGG”のZ-Score

22. 翻訳開始シグナル抽出結果 Escherichia coli 16S rRNA 3- terminal: gcggttggatcacctcctta3 Expected SD Sequence: 5taaggaggtgatccaaccgc Pat. Z-Sc. Pos. agga 94.97 -11 ggag 82.94 -10 aagg 58.15 -11 gagg 53.08 -11 gaga 42.23 -9

23. シグナル配列出現の評価 • 塩基の偏り • 偏りの程度 … エントロピー、増加情報量 • 偏りの有意性 … Χ2値 • 配列パターン • 出現の程度 … 頻度 • 出現の有意性 … Z Score

24. 転写 AUG RNA 翻訳 分子レベルの生命現象の根幹～ セントラルドグマ ～ ATG TAA DNA UAA タンパク質 機能

25. RNAレベルで機能する分子 • tRNA • rRNA • Other non-coding RNA • Translational regulation by mRNA

26. tRNA

27. ACC tRNA ACGAGUACA UGCUCAUGUUGG

28. rRNA

29. リボソーム Methionine fMet-tRNAf 16S rRNA AUUCCUCC mRNA AUG AGGAGG 開始コドン Shine-Dalgarno sequence 16S rRNAの3‘末端はShine-Dalgarno配列と 対合する

30. 翻訳での遺伝子の発現制御 Fe Ferritin gene 5’ 3’

31. 二次構造による終止コドンの読み飛ばし 二次構造 mRNA UAA 通常の長さのタンパク質 Steneberg, P. 2001 リードスルーによってできた長いタンパク質

32. Function of readthrough product is stronger hdc 1034 2981 4274 UAA UAA AUG hdc gene is expressed in tracheoles in larvae of D. melanogaster Possibility of Regulation by readthrough? Long product Short product Branching of lumens are inhibited strongly Branching of lumens are inhibited weakly Steneberg and Samakovlis, 2001

33. 転写 ＡUG ＡTG mRNA cDNA 逆転写 cDNA配列を用いた転写産物の収集 遺伝子 ＡＴＧ TAA ＤＮＡ UAA UAA

34. マウスcDNA配列の網羅的収集

35. コード領域を持たないcDNA? ゲノム cDNA Numata et al. 2003

36. さらにコード領域を持たない多数のcDNA?

37. タンパク質をコードしないcDNA配列が多くあるタンパク質をコードしないcDNA配列が多くある cDNA ゲノム ゲノムの62.5%をカバー 多くのRNAは翻訳されなくても機能を持つ？ 非翻訳RNAが多量に存在？

38. 多数の非翻訳RNAの存在が予想されて いるものの、ほとんどは機能未知

39. RNAの二次構造予測 • 一本鎖RNAはDNAに比べ、自由な構造を取ることが可能 • RNAが機能する上で立体構造が重要になってくる • 二次構造は、どの塩基とどの塩基が結合しているかを表す • 一次配列から二次構造を予測しよう！

40. tRNAの二次構造予測の例 GenBank tRNA配列 http://www.genome.ad.jp/dbget-bin/www_bget?gb:ECOCPTGG http://www.genome.ad.jp/dbget-bin/www_bget?gb:ECOPHER Zukerのmfold http://www.bioinfo.rpi.edu/applications/mfold/old/rna/form1.cgi