ゲノム解析で活躍するコンピューター

Size: px

Start display at page:

Download "ゲノム解析で活躍するコンピューター"

みさきさくいし
7 years ago
Views:

1 ゲノム解析で活躍するコンピューターゲノム情報のためのデータベースと遺伝子機能の解析京都大学化学研究所バイオインフォマティクスセンター五斗進科学カフェ京都第 62 回定例会京都大学 2010/4/10 1

2 本日お話する内容ゲノムとゲノムプロジェクトについてゲノムデータとデータベースについてゲノムデータを使った解析について遺伝子の機能を調べる 2

3 ゲノムとはゲノム (Genome) Gene( 遺伝子 )+ -ome( 全体 ): 遺伝子の総体 Gene( 遺伝子 )+ Chromosome( 染色体 )

4 ゲノムとは個体細胞 DNA の 2 重らせん構造染色体 atggcgacccgcagccctggcgtcgtgattagtgatgatgaaccaggttatgaccttgat ttattttgcatacctaatcattatgctgaggatttggaaagggtgtttattcctcatgga ctaattatggacaggactgaacgtcttgctcgagatgtgatgaaggagatgggaggccat cacattgtagccctctgtgt... ゲノムの全塩基配列

5 DNA とタンパク質 ( セントラルドグマ ) 遺伝子エクソンイントロンゲノム DNA 転写 ACGT: 塩基配列 mrna スプライシング ACGU: 塩基配列成熟 mrna 翻訳タンパク質折り畳み ACDEF : アミノ酸配列機能

6 ゲノム研究の歴史 1900 年前後 1950 年代 1960 年代遺伝法則の発見 ( メンデル ) ゲノム概念の提唱 ( ヴィンクラー ) DNA2 重らせん構造の発見 ( ワトソンクリック ) 遺伝暗号の解読 ( ニーレンバーグら ) 1970 年代 1980 年代 1990 年代 DNA 塩基配列決定法の開発 ΦX174 ファージゲノム配列決定 ( サンガーら ) PCR 法の開発 ( マリス ) ヒトゲノム計画の提案 ( ダルベッコ ) ヒトゲノム計画がスタート

7 ゲノム研究の歴史 1995 年 1996 年 1997 年 1998 年 1999 年 2000 年 2001 年ヘモフィルス菌ゲノムの解読 ( 約 200 万塩基 2000 遺伝子 ) ( 独立生活する生物のゲノム ) 出芽酵母ゲノムの解読 ( 約 1200 万塩基 6000 遺伝子 ) ( 真核生物のゲノム ) 枯草菌ゲノムの解読 ( 約 400 万塩基 4000 遺伝子 ) ( 日本を中心としたグループによる解読 ) 線虫ゲノムの解読 ( 約 9700 万塩基遺伝子 ) ( 多細胞生物のゲノム ) ヒト 22 番染色体ゲノムの解読ショウジョウバエシロイヌナズナゲノムの解読ヒトゲノムの概要配列発表 ( 約 30 億塩基 )

8 ポストゲノムプロジェクトトランスクリプトーム (Transcript + -ome) 転写産物 (RNAのこと) の総体細胞内で実際に mrna として発現している遺伝子プロテオーム (Protein + -ome) タンパク質の総体細胞内で実際にタンパク質として働いている遺伝子メタボローム (Metabolite + -ome) 代謝産物 ( アミノ酸グルコースなど ) の総体代謝系で合成されている化合物 8

9 本日お話する内容ゲノムとゲノムプロジェクトについてゲノムデータとデータベースについてゲノムデータを使った解析について遺伝子の機能を調べる 9

10 塩基配列データベース 1970 年代の配列決定技術サンガー法マキサムギルバート法自動 DNAシーケンサーの開発大量の塩基配列の産出 PIRタイプから著者によるサブミットへ三極体制によるデータの収集日本 :DDBJ (DNA DataBank of 遺伝研米国 National Center for Biotechnology Information, National Institute of Health 欧州 European Molecular Biology Laboratory ゲノムプロジェクト

11 GenBank の例 LOCUS X bp DNA linear BCT 12-SEP-1993 DEFINITION E.coli triose phosphate isomerase gene (TPI) (EC ). ACCESSION X00617 VERSION X GI:43111 KEYWORDS glycolysis gluconeogenesis; isomerase. SOURCE Escherichia coli ORGANISM Escherichia coli Bacteria; Proteobacteria; Gammaproteobacteria; Enterobacteriales; Enterobacteriaceae; Escherichia. REFERENCE 1 (bases 1 to 1338) AUTHORS Pichersky,E., Gottlieb,L.D. and Hess,J.F. TITLE Nucleotide sequence of the triose phosphate isomerase gene of Escherichia coli JOURNAL Mol. Gen. Genet. 195 (1-2), (1984) PUBMED COMMENT Data kindly reviewed (30-MAY-1985) by L.D. Gottlieb.

12 FEATURES Location/Qualifiers source /organism="escherichia coli" /mol_type="genomic DNA" /db_xref="taxon:562" CDS /note="unnamed protein product; isomerase" /codon_start=1 /transl_table=11 /protein_id="caa " /db_xref="gi:43112" /db_xref="goa:p0a858" /db_xref="pdb:1tmh" /db_xref="pdb:1tre" /db_xref="uniprotkb/swiss-prot:p0a858" /translation="mrhplvmgnwklngsrhmvhelvsnlrkelagvagcavaiappe MYIDMAKREAEGSHIMLGAQNVNLNLSGAFTGETSAAMLKDIGAQYIIIGHSERRTYH KESDELIAKKFAVLKEQGLTPVLCIGETEAENEAGKTEEVCARQIDAVLKTQGAAAFE GAVIAYEPVWAIGTGKSATPAQAQAVHKFIRDHIAKVDANIAEQVIIQYGGSVNASNA AELFAQPDIDGALVGGASLKADAFAVIVKAAEAAKQA" ORIGIN 1 ctgcaggacg cctactaagg cggcggggaa aaacaaacgt tattacaccg agacagaagg 61 tgcactgcgt tatgttgtcg cggacaacgg cgaaaagggg ctgaccttcg ctgttgaacc 121 aattaagttg gcgctatctg antctcatac tgtttcacag acctgctgcc ctgcggcggc 181 caatcttcct ttattcgctt ataagcgtgg agaattaaaa tgcgacatcc tttagtgatg 241 ggtaactgga aactgaacgg cagccgccac atggttcacg agctggtttc taacctgcgt 301 aaagagctgg caggtgttgc tggctgtgcg gttgcaatcg caccaccgga aatgtatatc ( 中略 ) 1021 ttactttcct taactcttcg ccttaacgca aaatctcaca ctgatgatcc tgaatttcct 1081 cggctgaagc acggttaagc gtcagtagat ttcgttgtgt cgccagcaat acaaatgagt 1141 tatcactctg ccgtaccatc gccagcccgt agcgtcccat atgttcccgc gcctcaggta 1201 cttcttctgc cagcatcata aatgggctgc gttgtaccag ttcgctttcc gttacccgac 1261 gcgcaggtat tcatgcccgc gcaaaccacc tggcagtggc aaccagcggc tgctgatgtt 1321 cgccagattg ttatcgag //

13 塩基配列データベースのサイズエントリーの数塩基の数塩基配列の増加は約 5 年で 10 倍のペース

14 配列データベースの使われ方似たものを探す ( ホモロジー検索 ) ある ( 機能未知の ) 遺伝子の塩基配列を実験で決定その配列をもとにしてデータベース中に似た配列が登録されていないかを探す前提 : 同じ機能を持つ遺伝子は似た配列を持つ a 種分岐 a1 a1 種 1 種 2

15 配列データベースの使われ方大腸菌のある酵素と同じ機能を持つコレラ菌の酵素をアミノ酸配列で比較した結果大腸菌コレラ菌

16 ゲノムプロジェクトが出すデータ個体 atggcgacccgcagccctggcgtcgtgattagtgatgatgaaccaggttatgaccttgat ttattttgcatacctaatcattatgctgaggatttggaaagggtgtttattcctcatgga ctaattatggacaggactgaacgtcttgctcgagatgtgatgaaggagatgggaggccat cacattgtagccctctgtgt... 細胞ゲノムの全塩基配列染色体全遺伝子のカタログ情報個々の遺伝子の機能についての情報ホモロジー検索だけで機能が分かる遺伝子は半分 2/3 程度

17 パスウェイの情報

18 パスウェイの情報

19 パスウェイの情報酵素 / 遺伝子と化合物のネットワークヒトの解糖系 ( 体内に取り込んだ糖を分解して再利用する経路 ) 緑 : ヒトに対応する遺伝子がある酵素白 : ヒトにない酵素

20 本日お話する内容ゲノムとゲノムプロジェクトについてゲノムデータとデータベースについてゲノムデータを使った解析について遺伝子の機能を調べる 20

21 遺伝子の機能予測とゲノムの機能予測 (A) 遺伝子の機能予測 (B) ゲノムの機能予測遺伝子の機能生物の機能立体構造 DB 発現 DB モチーフ DB 相同性検索立体構造予測相互作用 DB パスウェイ解析配列構造解析配列 DB パスウェイ DB 遺伝子配列情報ゲノム情報 ( 遺伝子の集合 ) 21

22 ゲノムの機能予測をした後は機能予測の抜けを探すパスウェイ中で途切れているところ機能未知遺伝子との対応? 様々な情報を比較種に特徴的な機能は何か? 機能未知遺伝子の機能予測パスウェイ ( ネットワーク ) のトポロジーを解析 22

23 機能予測の抜けの例リジン? グルタリル CoA クエン酸回路ゲノム情報から再構築された緑膿菌のリジン分解系 23

24 機能予測の抜けの例リジングルタリル CoA 未知の酵素遺伝子 (missing enzyme) クエン酸回路生化学的な知識による緑膿菌のリジン分解系 24

25 データ統合による知識抽出マイクロアレイ遺伝子発現類似度行列 ( カーネル ) ネットワーク推定酵母 2 ハイブリッド細胞内局在情報系統プロファイル仮定 : 似たパターンを持つ遺伝子同士は機能的に関係している可能性が高い 25

26 カーネル N 個のタンパク質 x 1, x2,..., x N を仮定カーネル K( x, x ) はタンパク質 xとx の類似度 ( 数学的には特徴ベクトルの内積 ) タンパク質セットの類似度行列 K = K( x, x ) ( i, j 1,2,..., N) ij i j = 26

27 ..., ) (, ) ( 2 1 = Φ = Φ x x ) ( ) ( ), ( = = Φ Φ = x x x x K 遺伝子の塩基構成比 x N x x,..., 1, 2 カーネルの例 27

28 カーネルの例 = = : ), ( ), ( ), ( ), ( x x x x x x x x K K K K K カーネル行列一種の類似度行列とみなせる一種の遺伝子類似性ネットワーク = 28

29 多様なデータと表現方法データマイクロアレイ遺伝子発現酵母 2 ハイブリッド ( タンパク質間相互作用 ) 細胞内局在系統プロファイル表現方法数値ベクトルグラフビットベクトルビットベクトル 29

30 多様なデータとデータ型数値ベクトル間の類似度を求める関数 K( x i, x j ) = exp( x i x j 2 2 / σ ) グラフ上の頂点間の類似度を求める関数 K = exp( L) ここで L はグラフのラプラシアン (Kondor, 2002) 30

31 31 多様なデータとデータ型各データにおける遺伝子間の距離をカーネル行列として定義 Kg: ゲノム上での遺伝子間の距離 Ke: 発現パターンの類似度 Kp: 系統プロファイルの類似度カーネルの和を取る K = Kg + Ke + Kp 統合されたカーネル K を用いて遺伝子間の関係を変換教師付き学習

32 教師なしの場合 32

33 機能予測の抜けの例リジングルタリル CoA 未知の酵素遺伝子 (missing enzyme) クエン酸回路生化学的な知識による緑膿菌のリジン分解系 33

34 教師付き学習発現データの類似度行列タンパク質ネットワーク Unknown pathway 34

35 教師付き学習発現データの類似度行列タンパク質ネットワークトレーニングトレーニング Unknown pathway 35

36 教師付き学習発現データの類似度行列タンパク質ネットワークトレーニングテストトレーニングテスト Unknown pathway 36

37 教師付き学習元の空間 x 1 x 3 x 2 : トレーニングセット 37

38 38 教師付き学習元の空間 x 1 x 2 x 3 : 教師なしで直接予測した結果 : トレーニングセット

39 39 教師付き学習元の空間 x 1 x 3 x 2 : 真のネットワーク : トレーニングセット

40 40 教師付き学習ステップ 1: 相互作用するタンパク質ペアが近くにあるような特徴空間に射影元の空間特徴空間 x 1 f f ( x 1 ) x 2 x 3 : 真のネットワーク : トレーニングセット f ( x 3 ) f ( x 2 )

41 41 教師付き学習ステップ 1: 相互作用するタンパク質ペアが近くにあるような特徴空間に射影元の空間特徴空間 x 1 f f ( x 1 ) x 2 x 3 : 真のネットワーク : トレーニングセット f ( x 3 ) f ( x 2 )

42 教師付き学習元の空間特徴空間 x 1 f f ( x 1 ) x 2 x 3 : 真のネットワーク : トレーニングセット f ( x 3 ) f ( x 2 ) : テストセット 42

43 教師付き学習ステップ 2: テストセットに関与するタンパク質間相互作用を予測元の空間特徴空間 x 1 f f ( x 1 ) x 2 x 3 : 真のネットワーク : トレーニングセット f ( x 3 ) f ( x 2 ) : テストセット 43

44 教師付き学習ステップ 2: テストセットに関与するタンパク質間相互作用を予測元の空間特徴空間 x 1 f f ( x 1 ) x 2 x 3 : 真のネットワーク : トレーニングセット f ( x 3 ) f ( x 2 ) : テストセット 44

45 45 アルゴリズム K K 1 2 : 発現データの類似度行列 : ネットワークの類似度行列 ( α 1,α ) 2 = arg max α 1 T K 1 K 2 α 2 ( 1+ λ 1 α T 1 K 12 α) 1/ 2 ( 1 1+ λ 1 α T 2 K 22 α ) 1/ 2 2 f n ( x) = α1 1( x j, x) j= 1 j K データの類似度行列が入力であることが特長

46 性能評価直接予測教師付き学習 46

47 バクテリアの代謝系遺伝子の予測 Missing enzymes 緑膿菌のリジン分解系 47

48 バクテリアの代謝系遺伝子の予測 Lysine degradation of Pseudomonas aeruginosa Glutaryl-CoA dehydrogenase + FAD FADH CO 2 Glutaryl-CoA Crotonoyl-CoA GcdH MVGKASFNWIDPLLLDQQLTEEERMVRDSAYQFAQDKLAPRVLEAFRHEQTDPAIFREMG EVGLLGATIPEQYGGSGLNYVCYGLIAREVERIDSGYRSMMSVQSSLVMVPINEFGTEAQ KQKYLPKLASGEWIGCFGLTEPNHGSDPGSMITRARKVDGGYRLTGSKMWITNSPIADVF VVWAKDDAGDIRGFVLEKGWQGLSAPAIHGKVGLRASITGEIVMDNVFVPEENIFPDVRG LKGPFTCLNSARYGISWGALGAAEACWHTARQYTLDRQQFGRPLAANQLIQKKLADMQTE ITLALQGCLRLGRMKDEGTAAVEITSIMKRNSCGKALDIARMARDMLGGNGISDEFGVAR HLVNLEVVNTYEGTHDVHALILGRAQTGIQAFY Citrate cycle 48

49 バクテリアの代謝系遺伝子の予測 Lysine degradation of Pseudomonas aeruginosa Lysine monooxygenase + O 2 + CO 2 + H 2 O L-lysine 5-aminopentanamide Gene: Unknown Citrate cycle 49

50 バクテリアの代謝系遺伝子の予測 Lysine degradation Pseudomonas aeruginosa????? No sequence data for any other species Gene: Unknown Citrate cycle 50

51 51 予測方法 Known???? 1.Predicting network of 4,225 protein genes in Pseudomonas, and extract lysine degradation pathway 2.Candidate genes are extracted from the genes located between known genes? TCA cycle Known

52 予測に使える情報機能的に関連のあるタンパク質の遺伝子はゲノム上で近い位置にある傾向 (Bork, P. et al., 1998) 機能的に関連のあるタンパク質は同じような進化パターンを持つ傾向 (Pazos, F., 2001; Pellegrini, M. et al, 1999) 52

53 53 系統プロファイル Pellegrini et al. Proc. Natl. Acad. Sci. USA, 96:4285 (1999) オーソログ遺伝子のパターンを分類 E.coli S.cerevisiae B.subtilis H.influenzae 同じパターンを持つ遺伝子は進化的機能的に関連がある

54 54 カーネル ( 類似度の表現 ) ゲノム上での位置 K gen ( x, x ) = exp( d / h) ここで d: 遺伝子 xとx 間の塩基数系統プロファイル K phy ( x, x ) = x x ここで x: 系統プロファイル統合 K int = K gen + K phy

55 55 予測結果 PA0262,260 (hypothetical) PA0265 (dehydrogenase) PA0266 (amino-transferase)

56 56 まとめゲノムをはじめとする多種多様なデータから生命システムを明らかにするにはデータを検索解析技術が重要ハイスループットデータをデータベース化するだけでなく既知の情報もデータベース化して新しいデータと組み合わせて解析できるようにすることが重要

57 京都大学で構築中のデータベース様々な種類のデータを生命現象の総体として立体的に再構築ツールの提供 EGassembler KAAS GENIES KegArray ネットワークの知識 KEGG PATHWAY 高次機能機能の階層分類相互参照用データ KEGG BRITE ツールの提供 e-zyme pathcomp KegArray ゲノムの知識 KEGG GENES/EGENES 化合物の知識 KEGG LIGAND 研究者の知識をゲノムレベルのデータと結びつける KEGG: Kyoto Encyclopedia of Gene and Genomes

58 ご静聴ありがとうございました

KEGG.ppt

KEGG.ppt 1 2 3 4 KEGG: Kyoto Encyclopedia of Genes and Genomes http://www.genome.jp/kegg/kegg2.html http://www.genome.jp/kegg/kegg_ja.html 5 KEGG PATHWAY 生体内(外)の分子間ネットワーク図代謝系 12カテゴリ中間代謝二次代謝薬の代謝全体像制御系 20カテゴリ