(57)【要約】【課題】ハプロタイプブロックおよび遺伝的構造の影響を考慮することによって、高精度な診断支援を行なうシステムを提供する。【解決手段】ハプロタイプブロック推定手段１３によってハプロタイプブロックの位置を推定し、ハプロタイプブロックごとに解析を行なうことによって、個体のハプロタイプパターンを高精度に推定する。遺伝的構造推定手段１５によって個体のハプロタイプパターンによるクラスタリングを行ない、集団をいくつかの亜集団に分割することによって、集団に存在する遺伝的構造の影響を除去する。遺伝的構造情報データベース１６と診療情報データベース１１とを用いて臨床情報と遺伝子情報との関連性を解析することによって、高精度な診断支援知識を提供する。臨床情報と遺伝子情報との関連性の解析によって得られた診断支援知識を基礎として、罹患危険度算出手段１９によって所定の個体が疾患に罹患する危険度を算出する。【選択図】図１
2005276022.tif 000002
【特許請求の範囲】
【請求項１】
複数の個体の診療情報を格納する診療情報データベースと、集団の遺伝子多型に関する情報を格納する遺伝子多型情報データベースと、前記遺伝子多型情報データベースの情報を基礎として前記集団のハプロタイプブロックおよび前記各ハプロタイプブロック内におけるハプロタイプ頻度を推定するハプロタイプブロック推定処理プログラムと、前記集団の推定された前記各ハプロタイプブロック内におけるハプロタイプパターンおよび前記ハプロタイプ頻度を格納するハプロタイプ情報データベースと、前記ハプロタイプ情報データベースの情報を基礎として前記集団に存在する遺伝的構造を推定して、前記集団を複数の亜集団に分割する遺伝的構造推定処理プログラムと、分割された前記亜集団ごとの前記ハプロタイプ情報および前記各個体の前記各亜集団への帰属度情報を格納する遺伝的構造情報データベースと、前記診療情報データベースおよび前記遺伝的構造情報データベースの情報を基礎として被診断者のハプロタイプと形質との関連性を解析する関連性解析処理プログラムと、前記関連性解析処理プログラムによって得られた情報を格納する診断支援知識データベースと、前記診断支援知識データベースの情報を基礎として所定の個体が疾患に罹患する危険度を算出する罹患危険度算出処理プログラムとを有することを特徴とする診断支援システム。
【請求項２】
前記遺伝的構造推定処理プログラムは、前記各ハプロタイプブロック内に存在するハプロタイプ間に定義される距離にもとづくクラスタリングを行なう処理と、前記クラスタリングによって得られた前記亜集団ごとに前記ハプロタイプパターンおよび前記ハプロタイプ頻度を求める処理と、適切な前記亜集団の数を判定する処理と、得られた前記亜集団に対する前記各個体の帰属度を求める処理を行なう請求項１に記載の診断支援システム。
【請求項３】
前記距離が、各ハプロタイプ間における組み換えおよび突然変異の起こりやすさにより定義される請求項２に記載の診断支援システム。
【請求項４】
遺伝子多型に関する情報を格納する遺伝子多型情報データベースの情報を基礎としてハプロタイプブロックおよび各ハプロタイプブロック内におけるハプロタイプ頻度を推定するステップ、推定された前記各ハプロタイプブロック内におけるハプロタイプパターンおよびハプロタイプ頻度をハプロタイプ情報データベースに格納するステップと、前記ハプロタイプ情報データベースの情報を基礎として集団に存在する遺伝的構造を推定して、前記集団を複数の亜集団に分割する遺伝的構造を推定するステップと、分割された前記亜集団ごとの前記ハプロタイプ情報および前記各個体の前記各亜集団への帰属度情報を遺伝的構造情報データベースに格納するステップと、複数の個体の診療情報を格納する診療情報データベースと前記遺伝的構造情報データベースの情報を基礎としてハプロタイプと形質との関連性を解析する関連性解析ステップと、前記関連性解析ステップで得られた情報を診断支援知識データベースに格納するステップと、前記診断支援知識データベースの情報を基礎として所定の個体が疾患に罹患する危険度を算出する罹患危険度算出ステップとよりなることを特徴とする診断支援方法。
【請求項５】
前記遺伝的構造を推定するステップは、前記各ハプロタイプブロック内に存在するハプロタイプ間に定義される距離にもとづくクラスタリングを行なう処理と、前記クラスタリングによって得られた前記亜集団ごとに前記ハプロタイプパターンおよび前記ハプロタイプ頻度を求める処理と、適切な前記亜集団の数を判定する処理と、得られた前記亜集団に対する前記各個体の帰属度を求める処理を行なうことよりなる請求項４に記載の診断支援方法。
【請求項６】
請求項５に記載の診断支援方法において、前記距離が、各ハプロタイプ間における組み換えおよび突然変異の起こりやすさにより定義されることを特徴とする診断支援方法。
【請求項７】
複数の個体の診療情報を格納する診療情報データベースと、遺伝子多型に関する情報を格納する遺伝子多型情報データベースと、前記遺伝子多型情報データベースの情報を基礎としてハプロタイプブロックおよび前記各ハプロタイプブロック内におけるハプロタイプ頻度を推定するハプロタイプブロック推定処理プログラムと、推定された前記各ハプロタイプブロック内におけるハプロタイプパターンおよび前記ハプロタイプ頻度を格納するハプロタイプ情報データベースと、前記ハプロタイプ情報データベースの情報を基礎として集団に存在する遺伝的構造を推定して、前記集団を複数の亜集団に分割する遺伝的構造推定処理プログラムと、分割された前記亜集団ごとの前記ハプロタイプ情報および前記各個体の前記各亜集団への帰属度情報を格納する遺伝的構造情報データベースと、前記診療情報データベースおよび前記遺伝的構造情報データベースの情報を基礎としてハプロタイプと形質との関連性を解析する関連性解析処理プログラムと、前記関連性解析処理プログラムによって得られた情報を格納する診断支援知識データベースと、前記診断支援知識データベースの情報を基礎として所定の個体が疾患に罹患する危険度を算出する罹患危険度算出処理プログラムとを有する診断支援システムに接続して受けることができるサ診断支援ービスであって、診断支援ービスを受けるものは、被診断者の個体から受け取った前記所定の個体遺伝子型データおよび形質データを診断支援システムに送信して、診断支援システムは、前記集団に存在する遺伝的構造に関する情報と、前記所定の個体の前記各亜集団への帰属度と、前記所定の個体が疾患に罹患する危険度とを算出して前記診断支援ービスを受けるものに提供することを特徴とする診断支援サービス。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、臨床情報と遺伝子情報との関連性を解析し、臨床上有用な情報を抽出して提示する診断支援システムおよび診断支援方法に関する。
【背景技術】
【０００２】
ヒトゲノムプロジェクトは配列決定をほぼ終了し、ポストシークエンスの時代に突入した。今後は、蓄積された膨大な遺伝子情報の医学への有効利用が期待されている。遺伝子と疾患との関連性の解明が進むと、個人の遺伝子型を基礎として疾患の発症リスクを予測することが可能になり、個人の遺伝的体質に応じた疾患の予防、早期発見、治療を行なうことが可能になる。これらを実現するためには、臨床情報と遺伝子情報との関連性の解析が必要である。
【０００３】
臨床情報と遺伝子情報との関連性の解析において強力な手法の一つに、遺伝統計学的解析法がある。遺伝統計学的解析法は、個人の遺伝子情報と疾患の有無とをデータとして、統計学を用いて疾患に関連する遺伝子を探索する方法であり、機序が未知である疾患の関連遺伝子も発見できる可能性があるため、次第に重要性を増している。遺伝統計学的解析法は、複数の遺伝子座（染色体上の遺伝子の位置）間の連鎖を利用して特定の形質に関係する遺伝子領域を探索する技術である。形質とは、個体レベルで観察される各種の形態的特徴のことで、疾患罹患の有無、身長、目や髪の色等が形質である。連鎖とは、「２つの異なる形質はそれぞれ分離独立して遺伝する」というメンデルの独立法則の例外である。
【０００４】
ある２つの形質を規定する遺伝子座が染色体上で近い位置に存在しているとき、それらの遺伝子は分離独立せずに、つながったまま親から子へ遺伝する。この状態を、２つの遺伝子座が連鎖しているという。減数分裂の際に、両親から伝わった１対の染色体間に部分的な交換が起こり、子に伝える遺伝子の組み合わせが両親由来のものと異なる場合がある。この現象を組み換えという。
【０００５】
１回の減数分裂においてある２つの遺伝子座間に組み換えが起こる確率を組み換え割合という。２つの遺伝子座間の距離が近いほど、組み換え割合は小さい。すなわち、連鎖する可能性が高い。遺伝統計学的解析法では、組み換えの情報を基礎として染色体上に網羅された遺伝子多型（一塩基多型（ＳＮＰｓ）やマイクロサテライト等）と疾患関連遺伝子との連鎖の有無を検定することによって、疾患関連遺伝子座を絞り込む。
【０００６】
遺伝統計学的解析法は現在までにいくつかの手法が報告されている。単一遺伝子疾患については、大家系のデータを用いたパラメトリック連鎖解析によって、これまでに数多くの原因遺伝子が同定されてきた。今後の疾患原因遺伝子探索の研究においては、複数の遺伝的要因と環境要因とによって発症する多因子疾患（complex disease）の原因遺伝子の探索が主流になるものと考えられる。当初は、多因子疾患の原因遺伝子についても、多数の小家系のデータを用いたノンパラメトリック連鎖解析（罹患同胞対解析）によってその同定が可能と考えられていた。しかし、一般に浸透率（発症する確率）の低い多因子疾患の原因遺伝子を直接同定することは困難な場合が多い。最近では、その検出力の高さと解析のしやすさから、疾患集団と正常集団において着目する遺伝子多型の対立遺伝子（アレル）頻度を比較する相関解析（関連解析ともいう）が注目されている。
【０００７】
従来の相関解析では、真に形質と関係している遺伝子を見落としたり、目的とする形質とまったく関係がない遺伝子を誤って選択したりする可能性が比較的高かった。一般に、前者は偽陰性の問題、後者は偽陽性の問題として取り扱われる問題である。解析結果に偽陰性や偽陽性が生じる理由としては、単一の遺伝子多型もしくは狭い範囲の遺伝子多型によるハプロタイプのみを用いて遺伝子と形質との関係を解析していること、ハプロタイプを用いた解析を行なう際にハプロタイプブロックを考慮していないこと、対象とする集団に存在する多様性（これを遺伝的構造と呼ぶことにする）を考慮していないこと、等が挙げられる。
【０００８】
ハプロタイプとは、連鎖している複数の座位における同じ親由来のアレルの組み合わせのことをいう。染色体上の近い距離に存在する複数の座位におけるアレルは、世代交代における組み換えの影響を受けずに連鎖した状態で次の世代へと伝達される。その結果、何世代にもおよぶ世代交代を経た後も、近い距離に存在する複数の座位にはお互いに相関関係が見られる。この状態を連鎖不平衡と呼ぶ。近年、例えば、非特許文献１（Gabriel SB et al.: The Structure of Haplotype Blocks in the Human Genome, Science, Vol.296, pp.2225‐2229, 2002）等によって、ゲノム上には連鎖不平衡が比較的強い状態で保たれたハプロタイプブロックとよばれる部分と、組み換えが高頻度で起こるために座位間の連鎖不平衡を弱めるホットスポットとよばれる部分とが交互に存在することが報告されている。
【０００９】
この事実は、ハプロタイプブロックの位置を正確に推定することができれば、ハプロタイプブロック内の数個の座位の遺伝子型を測定するだけで正確なハプロタイプパターンを決定することが可能であることを意味している。また、この事実は同時に、ホットスポットを跨ぐような複数の座位を用いて解析を行なった場合には、遺伝学的には意味がない偽陽性の結果が多く出てしまうことを意味している。
【００１０】
一般に、相関解析を行なう際には、注目する形質に応じて対象とする集団を群分けすることが多い。ある集団の中から多数の患者（case）と対照者（control）とをサンプルし、着目するアレルの頻度を患者群と対照者群とで比較し、アレル頻度に有意な差がみられる多型の座位を検出する症例対照研究（case−control study）が最も有名である。症例対照研究においては、患者の集団と対照者の集団とが注目する形質以外は完全にマッチした集団であることが前提となっている。
【００１１】
しかし、この前提は常に成立しているわけではない。対象とする集団に遺伝的構造が存在する場合は特に問題となる。患者群と対照者群とを遺伝的に異なる、全く別の集団からサンプリングしてしまった場合などには、遺伝的構造が解析結果に大きく影響を与える。簡単な例を挙げて集団の遺伝的構造の影響を説明する。例えば、アメリカで鎌形赤血球症の患者群と対照者群を集めようとすると、患者群にはアフリカ由来の人々が多く含まれ、対照者群にはヨーロッパ由来の人々が多く含まれるはずである。遺伝的構造の影響を考えずにこの２つの集団を比較すると、本来アフリカ人とヨーロッパ人とでアレル頻度に差がある多くの座位が、鎌形赤血球症の原因座位として検出されてしまう。このように、集団の遺伝的構造は解析結果に多くの偽陽性を生じさせる。また、集団の遺伝的構造は解析結果に偽陽性だけではなく偽陰性を生じさせることもある。
【００１２】
【非特許文献１】Gabriel SB et al.: The Structure of Haplotype Blocks in the Human Genome, Science, Vol.296, pp.2225‐2229, 2002
【発明の開示】
【発明が解決しようとする課題】
【００１３】
上述したように、相関解析を行なう際に、対象とする集団に存在するハプロタイプブロックの影響や遺伝的構造の影響を考慮しなかった場合、解析の際に多くの偽陰性や偽陽性が生じ、解析結果に多大な影響を与えるという問題があった。そこで、本発明では、ハプロタイプブロックおよび遺伝的構造の影響を考慮することによって、高精度な診断支援を行なうシステムを提供することを目的とする。
【課題を解決するための手段】
【００１４】
本発明の診断支援システムおよび診断支援方法は、ハプロタイプブロック推定手段によって遺伝子多型情報を基礎として組み換えの位置を推定してハプロタイプブロックの位置を推定し、ハプロタイプブロックごとに解析を行なうことによって、個体のハプロタイプパターンを高精度に推定する。推定されたハプロタイプ頻度情報や個体のハプロタイプパターン情報はハプロタイプ情報データベースに格納される。また、遺伝的構造推定手段によって個体のハプロタイプパターンによるクラスタリングを行ない、集団をいくつかの亜集団に分割することによって、集団に存在する遺伝的構造の影響を除去し、臨床情報と遺伝子情報との関連性を高精度に解析することを可能にする。遺伝的構造推定手段によって得られた結果は遺伝的構造情報データベースに格納され、遺伝的構造情報データベースと診療情報データベースとを用いて臨床情報と遺伝子情報との関連性を解析することによって、高精度な診断支援知識の提供が可能となる。臨床情報と遺伝子情報との関連性の解析によって得られた診断支援知識は診断支援知識データベースに格納され、罹患危険度算出手段によって診断支援知識データベースの情報を基礎として所定の個体が疾患に罹患する危険度を算出する。
【発明の効果】
【００１５】
本発明の診断支援システムおよび診断支援方法は、ハプロタイプブロック推定アルゴリズムによって組み換えの位置を推定してハプロタイプブロックの位置を推定し、ハプロタイプブロックごとに解析を行なうことによって、個体のハプロタイプパターンを高精度に推定することを可能にする。また、遺伝的構造推定アルゴリズムによって個体のハプロタイプパターンによるクラスタリングを行ない、集団をいくつかの亜集団に分割することによって、集団に存在する遺伝的構造の影響を除去し、臨床情報と遺伝子情報との関連性を高精度に解析することを可能にする。
【発明を実施するための最良の形態】
【００１６】
図１は、本発明の診断支援システムの構成例を示す図である。本発明の診断支援システム１１１は、いわゆるパソコン等の電子計算機を主体として構成される。システムバス５に処理装置１、メモリー２、入力装置３、表示装置４および外部記憶装置１０が接続される。外部記憶装置１０内には、複数の個体（被診断者）の診療情報を格納する診療情報データベース１１、複数の個体（被診断者）の遺伝子多型に関する情報を格納する遺伝子多型情報データベース１２、該遺伝子多型情報データベース１２の情報を基礎としてハプロタイプブロックの位置を推定し、ハプロタイプブロックごとに集団のハプロタイプ頻度および個体のハプロタイプパターンを推定して得られたハプロタイプブロックごとに集団のハプロタイプ頻度情報と個体のハプロタイプパターンを格納するハプロタイプ情報データベース１４、該ハプロタイプ情報データベース１４の情報を基礎として集団の遺伝的構造を推定し、ハプロタイプブロックごとに個体のハプロタイプパターンによるクラスタリングを行ない、集団をいくつかの亜集団に分割するとともに、各個体の各亜集団への帰属度を推定して得られた、分割された亜集団ごとのハプロタイプ情報および各個体の各亜集団への帰属度情報を格納する遺伝的構造情報データベース１６、前記診療情報データベース１１および遺伝的構造情報データベース１６の情報を基礎として、亜集団のハプロタイプブロックごとに個体のハプロタイプパターンと形質との関連性を解析し、疾患に罹患する危険度を算出する関連性解析によって得られた知識を格納する診断支援知識データベース１８と、前記遺伝子多型情報データベース１２の情報から前記ハプロタイプ情報データベース１４の情報を導出するためのハプロタイプブロック推定処理プログラム１３、前記ハプロタイプ情報データベース１４の情報から前記遺伝的構造情報データベース１６の情報を導出する遺伝的構造推定処理プログラム１５、前記診療情報データベース１１および前記遺伝的構造情報データベース１６の情報から前記診断支援知識データベース１８の情報を導出する関連性解析処理プログラム１７、前記診断支援知識データベース１８の情報を基礎として所定の個体が疾患に罹患する危険度を算出する罹患危険度算出処理プログラム１９が内蔵される。もちろん、これらの他に、電子計算機としての機能を果たすために必要とされるデータベースおよび処理プログラムが備えられる。
【００１７】
ここで、上述のデータベースは集団のデータが扱われるものであり、診断支援知識データベース１８の情報は、その集団に対して有効なものである。また、これらのデータベースの内容は診断を受けた人のデータの蓄積により、より充実したものとなる。
【００１８】
本発明の診断支援システムは、ハプロタイプブロック推定処理プログラム１３によって遺伝子多型情報を基礎として組み換えの位置を推定してハプロタイプブロックの位置を推定し、ハプロタイプブロックごとに解析を行なうことによって、個体のハプロタイプパターンを高精度に推定する。推定されたハプロタイプ頻度情報や個体のハプロタイプパターン情報はハプロタイプ情報データベース１４に格納される。また、遺伝的構造推定手段１５によって個体のハプロタイプパターンによるクラスタリングを行ない、集団をいくつかの亜集団に分割することによって、集団に存在する遺伝的構造の影響を除去し、臨床情報と遺伝子情報との関連性を高精度に解析することを可能にする。遺伝的構造推定処理プログラム１５によって得られた結果は遺伝的構造情報データベース１６に格納され、遺伝的構造情報データベース１６と診療情報データベース１１とを用いて臨床情報と遺伝子情報との関連性を解析することによって、高精度な診断支援知識の提供が可能となる。臨床情報と遺伝子情報との関連性の解析によって得られた診断支援知識は診断支援知識データベース１８に格納され、罹患危険度算出処理プログラム１９によって診断支援知識データベース１８の情報を基礎として所定の個体が疾患に罹患する危険度を算出する。
【００１９】
診療情報データベース１１には、個体の氏名、住所、生年月日、家族構成等の基本データや、個体の既往歴、家族歴、主訴、所見、検査結果、生活習慣、症状経過、治療経過、薬剤の処方に関する情報等の臨床データや、インフォームドコンセントに関するデータ等を格納する。遺伝子多型情報データベース１２には、多型に関する基本情報（位置、測定方法、多型タイプ（ＳＮＰ、ＳＴＲＰ等）、アレル頻度等）や、個体の遺伝子多型測定結果（塩基配列パターン、ホモ、ヘテロ等）や、検査に用いた検体の識別情報や、保存状態等の検体管理データ等を格納する。
【００２０】
次に、ハプロタイプブロック推定処理プログラム１３について説明する。先に述べたように、ハプロタイプブロック内では比較的強い状態で連鎖不平衡が保たれている。また、例えば、先に述べた非特許文献１に示されるように、ハプロタイプブロック内ではハプロタイプの多様性が比較的小さいことも知られている。ハプロタイプブロックの位置を推定するためには、ゲノム上のある領域における連鎖不平衡の強さを定義する必要がある。
【００２１】
一般に、連鎖不平衡の強さは２つの座位間の連鎖不平衡係数Ｄ’を用いて表されることが多い。本発明では、例えば、ある領域における複数の座位の連鎖不平衡係数が次式の条件を満たすような場合に、その領域をハプロタイプブロックとして定義する。
ｍｉｎ（｜Ｄ’｜）＞０．８
推定したハプロタイプブロックについて、各ハプロタイプブロック内における集団のハプロタイプ頻度および個体のハプロタイプパターンを推定する。個体が持つ２つのハプロタイプの組み合わせのことをディプロタイプ形と呼ぶことにする。遺伝子型データから個体のディプロタイプ形を推定するための手法は現在までにいくつか提案されている。代表的なものとしては、例えば、文献：Excoffier L & Slatkin M: Maximum‐likelihood estimation of molecular haplotype frequencies in a diploid population, Mol Biol Evol, Vol. 12, pp. 921-927, 1995に示すようなＥＭアルゴリズムを用いた手法や、文献：Stephens M et al.: A new statistical method for haplotype reconstruction from population data, Am J Hum Genet, Vol. 68, pp. 978‐989, 2001に示すようなＰＨＡＳＥ法等がある。
【００２２】
ＥＭアルゴリズムを用いて母集団のハプロタイプ頻度と個体のディプロタイプ形とを推定する方法について以下で説明する。いま、ｎ個の個体からなるサンプル集団を考える。この集団において、連鎖する複数のマーカー座位におけるハプロタイプを考え、その母集団における頻度をＦ＝（Ｆ１，Ｆ２，・・・，ＦＭ）とする。Ｍは可能なハプロタイプの総数である。例えば、マーカー座位が全てＳＮＰ座位である場合、座位数をＬとするとＭ＝２Ｌである。各個体の、連鎖する複数のマーカー座位における遺伝子型の観察データをＧ＝（Ｇ１，Ｇ２，・・・，Ｇｎ）とする。多くの場合、Ｇｉは不完全データである。したがって、Ｇｉに対応するディプロタイプ形は１つに定まらないことが多い。このような場合は、可能なディプロタイプ形の上の確率分布（これをディプロタイプ分布と呼ぶ）を定義する。個体ｉ（ｉ＝１，２，・・・，ｎ）について、Ｇｉに対応するディプロタイプ形をＤｉｊ（ｊ＝１，２，・・・，ｍｉ）とする。ここで、ｍｉはＧｉに対して可能なディプロタイプの数であり、ｍｉの最大値はＭである。
【００２３】
図２は、母集団のハプロタイプ頻度と個体のディプロタイプ形とを推定するハプロタイプブロック推定処理プログラム１３の例を示す図である。
【００２４】
ステップ２１：まず、可能なＭ個のハプロタイプ（それぞれＨ１，Ｈ２，・・・，ＨＭとする）に対して、ハプロタイプ頻度の初期値Ｆ（０）を与える。ハプロタイプ頻度の合計は１である。
【００２５】
次に、ｔ＝０，１，２，・・・について、以下のステップ２２〜ステップ２５によってＦ（ｔ）からＦ（ｔ＋１）を計算する。
【００２６】
ステップ２２：各ディプロタイプ形Ｄｉｊは２つのハプロタイプＨｌ，Ｈｍによって構成されている。ただし、１≦ｌ≦Ｍ、１≦ｍ≦Ｍである。母集団のハプロタイプ頻度Ｆ（ｔ）が与えられているとき、Ｄｉｊが得られる確率は式（１）の通りである。
【００２７】
【数１】
2005276022.tif 000003
【００２８】
したがって、遺伝子型の観察データＧｉのもとでの、個体ｉのディプロタイプ形がＤｉｊである事後確率Ｐｒ（Ｄｉｊ｜Ｇｉ）は、ベイズの定理より、式（２）となる。
【００２９】
【数２】
2005276022.tif 000004
【００３０】
これを全てのｊ（ｊ＝１，２，・・・，ｍｉ）について計算すれば個体ｉのディプロタイプ分布が定まる。これをサンプル集団における全ての個体について適用する。
【００３１】
ステップ２３：個体のディプロタイプ分布が定まると、サンプル集団における全ての個体のディプロタイプ分布から母集団のハプロタイプ頻度の期待値を計算することができる。母集団のハプロタイプ頻度の期待値は、式（３）となる。
【００３２】
【数３】
2005276022.tif 000005
【００３３】
ここで、ＮＤｊｋｉはディプロタイプ形Ｄｊｋの中に含まれるＨｉの数（すなわち０，１，２のいずれか）である。
【００３４】
ステップ２４：このとき、全体の尤度は、個体ごとの全てのディプロタイプ形の尤度を結合し、さらに全ての個体の尤度を結合することによって、式（４）で表すことができる。
【００３５】
【数４】
2005276022.tif 000006
【００３６】
ステップ２５：Ｆ（ｔ＋１）＝Ｅ［Ｆ（ｔ）］としてＦを更新する。Ｌ（Ｆ）の値が収束したか否か判定する。Ｌ（Ｆ（ｔ＋１））−Ｌ（Ｆ（ｔ））＜βを満足すれば収束としてステップ２６に進み、満足しなければ、ステップ２２に戻り、ステップ２５まで繰り返す。ここでβは閾値である。
【００３７】
ステップ２６：収束した時点でのＥ［Ｆ］＝Ｆ（ＥＭ）を母集団におけるハプロタイプ頻度の最尤推定値とし、このときのＰｒ（Ｄ｜Ｇ）を母集団におけるハプロタイプ頻度の最尤推定値のもとでの個体のディプロタイプ分布とする。
【００３８】
ハプロタイプ情報データベース１４には、前述したように、遺伝子多型情報データベース１２の情報を基礎として、ハプロタイプブロックの位置を推定し、ハプロタイプブロックごとに集団のハプロタイプ頻度および個体のハプロタイプパターンを推定して得られたハプロタイプブロックごとに集団のハプロタイプ頻度情報と個体のハプロタイプパターンを格納するとともに、ハプロタイプブロックの設定に必要な基本情報と、各ハプロタイプブロック内におけるハプロタイプパターンおよびハプロタイプ頻度情報とを格納する。
【００３９】
図３は、ハプロタイプブロックの設定に必要な基本情報の格納データ例を示す図である。例えば、遺伝子ＧＥＮＥ＿１については、ＳＮＰ多型である多型ＰＯＬ＿１および多型ＰＯＬ＿２とＳＴＲＰ多型であるＰＯＬ＿３がテーブルに登録されており、ＰＯＬ＿１，ＰＯＬ＿２およびＰＯＬ＿３がハプロタイプブロックＨＢ＿１を構成していることを示している。図３に示したデータ以外にも、例えば、ハプロタイプブロックの長さ、ハプロタイプブロックを構成する多型の選択基準（アレル頻度やアミノ酸変異の有無等）、連鎖不平衡係数、ハプロタイプブロックを構成する多型が存在する遺伝子の位置等を格納してもよい。
【００４０】
図４は、各ハプロタイプブロック内におけるハプロタイプパターンおよびハプロタイプ頻度情報の格納例を示す図である。例えば、ハプロタイプブロックＨＢ＿１内には、ハプロタイプＨＴ＿１，ハプロタイプＨＴ＿２，ハプロタイプＨＴ＿３およびハプロタイプＨＴ＿４の４つのハプロタイプが存在し、各ハプロタイプの母集団における頻度はそれぞれ０．５０，０．２８，０．１５および０．０７であることを示している。
【００４１】
図５は、個体ごとのハプロタイプパターンの格納例を示す図である。例えば、個体ＰＥＲＳＯＮ＿１はハプロタイプブロックＨＢ＿１についてはハプロタイプＨＴ＿１を２つ有しており（２つのハプロタイプＨＴ＿１から構成されるディプロタイプ形を有しており）、そのディプロタイプ形を有している確率が１．００であることを示している。同様に、個体ＰＥＲＳＯＮ＿１はハプロタイプブロックＨＢ＿２については２つのハプロタイプＨＴ＿５から構成されるディプロタイプ形（確率０．９５）またはハプロタイプＨＴ＿５およびハプロタイプＨＴ＿６から構成されるディプロタイプ形（確率０．０５）を有しており、ハプロタイプブロックＨＢ＿ｍについては２つのハプロタイプＨＴ＿Ｙから構成されるディプロタイプ形（確率１．００）を有していることを示している。
【００４２】
次に、遺伝的構造推定処理プログラム１５について説明する。本発明では、集団の遺伝的構造を推定するために、個体のハプロタイプパターンによるクラスタリングを行ない、集団をいくつかの亜集団に分割する。本発明では、各ハプロタイプ間における突然変異および組み換えの起こりやすさによって決定される距離を新しく定義し、この距離を用いて個体のクラスタリングを行なう。以下では、本発明のクラスタリング手法について述べる。
【００４３】
図６は、あるハプロタイプブロック内においてハプロタイプ１〜ハプロタイプ５に示す５つのハプロタイプが観察された例を説明する図である。各ハプロタイプ間の距離を計算するために、まず図６に示すようなハプロタイプ進化系統樹を作成する。ハプロタイプ進化系統樹を作成する方法は、例えば、文献：McPeek MS & Strahs A:Assessment of linkage disequilibrium by the decay of haplotype sharing、 with application to fine‐scale genetic mapping、Am J Hum Genet、Vol. 65, pp. 858‐875, 1999に示される方法等、現在までにいくつか報告されている。
【００４４】
本発明では、進化系統樹の枝（edge）が１回の突然変異もしくは１回の組み換えによる進化を表すように進化系統樹を作成する。図６のハプロタイプ１からハプロタイプ５への進化のように、１回の突然変異もしくは１回の組み換えだけでは進化を表すことができない場合は、実際には観察されていない補助的なハプロタイプを挿入して進化系統樹を作成する。図６のハプロタイプ６は、この補助的なハプロタイプの例である。
【００４５】
次に、作成した系統樹のそれぞれの枝について、その進化が組み換えによるものなのか突然変異によるものなのかを決定する。例えば、図６において、ハプロタイプ１からハプロタイプ４への進化は組み換えによるものと考えられるが、ハプロタイプ１からハプロタイプ２への進化やハプロタイプ１からハプロタイプ３への進化は、突然変異による進化と組み換えによる進化の両方が考えられる。
【００４６】
あるハプロタイプＨＳが別のハプロタイプＨＴへと進化したときの尤度は、式（５）で表される。
【００４７】
【数５】
2005276022.tif 000007
【００４８】
ここで、ｍｕｔ．は突然変異（mutation）を、ｒｅｃ．は組み換え（recombination）を示す。（式）５は、あるハプロタイプＨＳが別のハプロタイプＨＴへと進化したときの尤度はその進化が突然変異によるものであると仮定した場合の尤度と組み換えによるものであると仮定した場合の尤度との和で表されることを示している。ここで、ある座位ｊにおける突然変異率をγｊ、ハプロタイプにおけるｋ番目のギャップの組み換え割合をθとすると、Ｐｒ（ｍｕｔ．｜ｍｕｔ．  ｏｒ  ｒｅｃ．）＝Ａ／（Ａ＋Ｂ）であり、また、Ｐｒ（ｒｅｃ．｜ｍｕｔ．  ｏｒ  ｒｅｃ．）＝Ｂ／（Ａ＋Ｂ）である。ただし、Ａは式（６）、Ｂは式（７）に示す通りである。
【００４９】
【数６】
2005276022.tif 000008
【００５０】
【数７】
2005276022.tif 000009
【００５１】
図６におけるハプロタイプ１からハプロタイプ４への進化のように、ハプロタイプを構成する多型が２つ以上の座位で異なる場合は、その進化が組み換えによるものであることが明らかであり、Ｐｒ（ＨＴ｜ＨＳ、ｍｕｔ．）＝０である。組み換えによる進化の場合は、例えば、図６のハプロタイプ１からハプロタイプ４への進化の場合、ハプロタイプ１およびハプロタイプ４の右側に共通する部分ハプロタイプＧＣＣＣＴＣＴＡＴ上のどのギャップ（両端を含む）において組み換えが起こっても、見かけ上は同じハプロタイプが形成される。そこで、ＨＳとＨＴとが、ｋ０番目のギャップまでは見かけ上同じアレルから構成されており（これをＩＢＳ（identical by state）という）、それ以降の部分で異なっているとすると、組み換えによる進化の場合の尤度は式（８）と表される。
【００５２】
【数８】
2005276022.tif 000010
【００５３】
いま、ＨＳがＬ個の座位によって構成されていることとし、ＨＳのうちの座位ｍ，ｍ＋１，・・・，ｎの部分で構成される部分ハプロタイプをＨＳ｛ｍ：ｎ｝と表すことにする。ＨＴについても同様に表すことにすると、式（９）となる。
【００５４】
【数９】
2005276022.tif 000011
【００５５】
ここで、ある２つのハプロタイプがＩＢＤ（identical by descent）であるとは、同祖由来のアレルを共有していることを示す。また、ある２つのハプロタイプが見かけ上はＩＢＳであっても実際にはＩＢＤである場合もあるため、これをＩＢＳ＊と表すことにする。
【００５６】
ベイズの定理を適用すると、式（１０）となる。
【００５７】
【数１０】
2005276022.tif 000012
【００５８】
ここで、式（１１）と仮定することができ、
【００５９】
【数１１】
2005276022.tif 000013
【００６０】
式（１２）はＨＴ｛１：ｋ｝の頻度であるから、式（１０）の値は容易に計算することができる。
【００６１】
【数１２】
2005276022.tif 000014
【００６２】
本発明では、式（５）で表される尤度を各ハプロタイプ間の距離として新しく定義し、この距離を用いて個体のクラスタリングを行なうこととした。したがって、ｋ番目のハプロタイプブロックについてＨｋａｋ、Ｈｋｂｋのハプロタイプをもつ個体とＨｋｃｋ、Ｈｋｄｋのハプロタイプをもつ個体との距離ｄｋを、式（１３）のように定義する。
【００６３】
【数１３】
2005276022.tif 000015
【００６４】
ハプロタイプブロックの数をｍとすると、２つの個体間の距離ｄは、全てのハプロタイプブロックにおける距離を結合して、式（１４）となる。
【００６５】
【数１４】
2005276022.tif 000016
【００６６】
次に、個体の帰属度の推定方法、すなわち、遺伝的構造推定処理プログラム１５について説明する。本発明では、各個体が、先に述べたクラスタリング手法によって生成された亜集団のうちのどの亜集団に属するかという情報を、個体の帰属度として定義する。
【００６７】
図７は、個体の帰属度を推定する遺伝的構造推定処理プログラム１５を示す図である。
【００６８】
ステップ７１：図６を参照して説明した方法によって、各ハプロタイプごとにハプロタイプ間の距離を決定する。
【００６９】
ステップ７２：ハプロタイプ間の距離にもとづくクラスタリングを行なう。
【００７０】
ステップ７３：ステップ７２の結果から、ｎ個の個体からなる集団がＮ個の亜集団に分割されたとする。このとき、ある個体ｉがある亜集団ｊに分類されているとすると、個体ｉの亜集団ｊへの帰属度は１００％であり、個体ｉの亜集団ｊ以外の亜集団への帰属度は０％である。ハプロタイプブロックの数をｍとすると、全体の尤度は式（１５）と表すことができる。
【００７１】
【数１５】
2005276022.tif 000017
【００７２】
ここで、Ｐｒ（Ｄ｜Ｇ）は個体の最尤ディプロタイプ分布であり、式（１６）はある亜集団ｊのｋ番目のハプロタイプブロックにおける個体ｉの最尤ディプロタイプ分布を示す。
【００７３】
【数１６】
2005276022.tif 000018
【００７４】
ステップ７４：Ｌ（Ｎ）の値が収束したか否か判定する。Ｌ（Ｎｋ−１）−Ｌ（Ｎｋ）＜βを満足すれば収束としてステップ７５に進み、満足しなければ、ステップ７１に戻り、ステップ７４まで繰り返す。ここでβは閾値である。
また、式（１７）は個体ｉの亜集団ｊへの帰属度である。
【００７５】
【数１７】
2005276022.tif 000019
【００７６】
ステップ７５：式（１５）で表される尤度が最大となるときのＮが亜集団数の最尤推定値である。この最尤推定値をパラメータとして採用する。
【００７７】
ステップ７６：式（１５）で表される尤度を基礎として個体の各亜集団への帰属度を計算する。例えば、Ｎ＿｛ｋ｝個の亜集団があり、次の連結ステップで亜集団Ｎ＿｛ｌ｝と亜集団Ｎ＿｛ｌ＋１｝とが連結されてＮ＿｛ｋ−１｝個の亜集団が形成されるとすると、このステップにおいて尤度に変化がなく、かつこのときに尤度が最大となる場合には、亜集団Ｎ＿｛ｌ｝および亜集団Ｎ＿｛ｌ＋１｝に分類されている全ての個体について、亜集団Ｎ＿｛ｌ｝および亜集団Ｎ＿｛ｌ＋１｝への帰属度をそれぞれ５０％ずつとする。
【００７８】
遺伝的構造情報データベース１６には、先にも述べたように、各亜集団におけるハプロタイプパターンおよびハプロタイプ頻度情報と、個体ごとの各亜集団への帰属度情報とを格納する。
【００７９】
図８は、各亜集団におけるハプロタイプパターンおよびハプロタイプ頻度情報の格納例を示す図である。例えば、亜集団ＳＵＢＰＯＰ＿１および亜集団ＳＵＢＰＯＰ＿２内にハプロタイプブロックＨＢ＿１、ＨＢ＿２がある。ここで、亜集団ＳＵＢＰＯＰ＿１内にはハプロタイプＨＴ＿１，ハプロタイプＨＴ＿２，ハプロタイプＨＴ＿３およびハプロタイプＨＴ＿４の４つのハプロタイプが存在し、亜集団ＳＵＢＰＯＰ＿２内には別のハプロタイプＨＴ＿７，ハプロタイプＨＴ＿８，ハプロタイプＨＴ＿９の３つのハプロタイプが存在することを示している。
【００８０】
一方、図４を参照して分かるように、例えば、ハプロタイプブロックＨＢ＿１内には、ハプロタイプＨＴ＿１，ハプロタイプＨＴ＿２，ハプロタイプＨＴ＿３およびハプロタイプＨＴ＿４の４つのハプロタイプが存在し、各ハプロタイプの母集団における頻度はそれぞれ０．５０，０．２８，０．１５および０．０７であることを示している。また、ハプロタイプブロックＨＢ＿１内に、別のハプロタイプＨＴ＿７，ハプロタイプＨＴ＿８およびハプロタイプＨＴ＿９の３つのハプロタイプが存在し、各ハプロタイプの母集団における頻度はそれぞれ０．３４，０．３３および０．３３であることを示している。
【００８１】
図９は、個体ごとの各亜集団への帰属度情報の格納例を示す図である。例えば、個体ＰＥＲＳＯＮ＿１は亜集団ＳＵＢＰＯＰ＿１への帰属度は１．００（１００％と百分率で表記してもよい）であり、個体ＰＥＲＳＯＮ＿２は亜集団ＳＵＢＰＯＰ＿１への帰属度が０．５０（５０％）、亜集団ＳＵＢＰＯＰ＿３への帰属度が０．５０（５０％）であることを示している。
【００８２】
次に、関連性解析処理プログラム１７によって、診療情報データベース１１および遺伝的構造情報データベース１６の情報を基礎として、各亜集団のハプロタイプブロックごとに個体のハプロタイプパターンと形質との関連性を解析する手順について説明する。関連性解析処理プログラム１７は、特定のハプロタイプを所有する個体の群と所有しない個体の群との間の形質を比較して（例えば、疾患の発症の有無を比較して）両群間のオッズ比等を計算し、特定のハプロタイプを所有する個体の群が特定のハプロタイプを所有しない個体の群と比較してどの程度疾患を発症するリスクが高まるかどうかを推定する。
【００８３】
本発明では、例えば、特定のハプロタイプを所有する個体の群の、特定のハプロタイプを所有しない個体の群に対する疾患発症のオッズ比をハプロタイプ相対リスクとして定義する。多くの場合、特定のハプロタイプの所有の有無、疾患の発症の有無（臨床イベントの有無や薬剤の副作用の有無等でもよい）によって２×２分割表を作成し、この２×２分割表の独立性の検定（χ２乗検定やＦｉｓｈｅｒの直接確率法を用いる）によって特定のハプロタイプの所有の有無が疾患の発症の有無に与える影響を計算する。形質がいくつかのカテゴリに分割できないような場合は、ｔ検定やＷｉｌｃｏｘｏｎテスト等を実施し、特定のハプロタイプを所有する個体の群と所有しない個体の群との形質の差を比較してもよい。
【００８４】
関連性解析処理プログラム１７によって得られた知識は、診断支援知識データベース１８に格納される。
【００８５】
図１０は、診断支援知識データベース１８の記述例を示す図である。各亜集団におけるハプロタイプ相対リスク情報の格納例を示している。ハプロタイプ相対リスクは、疾患の発症の有無や臨床イベントの有無、検査結果の正常・異常、薬剤の副作用の有無等、様々な臨床データについて定義することが可能であり、ここでは、心疾患、糖尿病、疾患Ｘの発症の有無に対する各亜集団ごとのハプロタイプ相対リスク情報の格納例を示している。例えば、ハプロタイプＨＴ＿１は亜集団ＳＵＢＰＯＰ＿１内では心疾患に対する相対リスクが１．５０であり、糖尿病、疾患Ｘに対する相対リスクがそれぞれ１．３５，１．００であることを示している。また、同時に、ハプロタイプＨＴ＿１は亜集団ＳＵＢＰＯＰ＿２内では心疾患に対する相対リスクが２．００であり、糖尿病、疾患Ｘに対する相対リスクがそれぞれ１．８９，１．００と変化することを示している。
【００８６】
罹患危険度算出処理プログラム１９は、遺伝的構造情報データベース１６および診断支援知識データベース１８を参照して所定の個体が疾患に罹患する危険度を算出する。個体ｉがある疾患に罹患する危険度Ｒｉは、ハプロタイプブロックの数をｍ、集団内に存在する亜集団の数をＮ、亜集団ｊのハプロタイプブロックｋにおける個体ｉのハプロタイプ相対リスクをｒｉｊｋとすると、式（１８）と表すことができる。
【００８７】
【数１８】
2005276022.tif 000020
【００８８】
図１１は、外部診療機関１１２から接続パス３１，３２およびインターネット３０を介して本発明の診断支援システム１１１にアクセスし、本発明の診断支援システム１１１を利用した診断支援を受ける場合のシステム例を示す図である。外部診療機関１１２も、いわゆるパソコン等の電子計算機を備え、システムバス５に処理装置１、メモリー２、入力装置３、表示装置４および外部記憶装置１０が接続される。ただし、外部診療機関１１２は、本発明のように、大規模な母集団のデータを扱うものではないので、複数の個体（被診断者）の診療情報を格納する診療情報データベース１１３および複数の個体（被診断者）の遺伝子多型に関する情報を格納する遺伝子多型情報データベース１１４は小規模のもので良い。単に、被診断者の診断に際して、個別に、本発明の診断支援システム１１１を利用した診断支援を受けることだけであれば、診療情報データベース１１３および遺伝子多型情報データベース１１４は無くても良い。尤も、本発明の診断支援システム１１１は、これを利用する外部診療機関１１２が、被診断者のデータを収集してこれを提供してくれることで、データを充実させてシステムをより完全なものにしていくことが望ましい。外部診療機関１１２が本発明の診断支援システム１１１を利用した診断支援を受ける場合は、外部診療機関１１２は診療情報データベース１１３および遺伝子多型情報データベース１１４から個体の遺伝子データおよび形質データを抽出して本発明の診断支援システム１１１へ送付する。外部診療機関１１２が診療情報データベース１１３および遺伝子多型情報データベース１１４を持たないときは、これらの情報を入力装置３から入力して、本発明の診断支援システム１１１に送付するものとすれば良い。本発明の診断支援システム１１１は、これらのデータを基礎に、算出された疾患に対する罹患危険度情報、遺伝的構造情報、個体の各亜集団への帰属度情報等を依頼元の外部診療機関１１２へ提供する。計算機の処理フローは特に説明するまでもない。
【図面の簡単な説明】
【００８９】
【図１】本発明の診断支援システムの構成例を示す図である。
【図２】母集団のハプロタイプ頻度と個体のディプロタイプ形とを推定するハプロタイプブロック推定処理プログラム１３の例を示す図である。
【図３】ハプロタイプブロックの設定に必要な基本情報の格納データ例を示す図である。
【図４】各ハプロタイプブロック内におけるハプロタイプパターンおよびハプロタイプ頻度情報の格納例を示す図である。
【図５】個体ごとのハプロタイプパターンの格納例を示す図である。
【図６】あるハプロタイプブロック内においてハプロタイプ１〜ハプロタイプ５に示す５つのハプロタイプが観察された例を説明する図である。
【図７】個体の帰属度を推定する遺伝的構造推定処理プログラム１５を示す図である。
【図８】各亜集団におけるハプロタイプパターンおよびハプロタイプ頻度情報の格納例を示す図である。
【図９】個体ごとの各亜集団への帰属度情報の格納例を示す図である。
【図１０】診断支援知識データベース１８の記述例を示す図である。
【図１１】外部診療機関１１２から接続パス３１，３２およびインターネット３０を介して本発明の診断支援システム１１１にアクセスし、本発明の診断支援システム１１１を利用した診断支援を受ける場合のシステム例を示す図である。
【符号の説明】
【００９０】
１…処理装置、２…メモリー、３…入力装置、４…表示装置、５…システムバス、１０…外部記憶装置、１１…診療情報データベース、１２…遺伝子多型情報データベース、１３…ハプロタイプブロック推定処理プログラム、１４…ハプロタイプ情報データベース、１５…遺伝的構造推定処理プログラム、１６…遺伝的構造情報データベース、１７…関連性解析処理プログラム、１８…診断支援知識データベース、１９…罹患危険度算出処理プログラム、２１…ハプロタイプ頻度初期値設定ステップ、２２…ディプロタイプ分布計算ステップ、２３…尤度計算ステップ、２４…ハプロタイプ頻度期待値計算ステップ、２５…ハプロタイプ頻度・ディプロタイプ分布最尤推定値採用ステップ、７１…ハプロタイプ間距離決定ステップ、７２…クラスタリング実行ステップ、７３…尤度計算ステップ、７４…パラメータ採用ステップ、７５…帰属度計算ステップ、１１１…診断支援システム、１１２…外部診療機関、１１３…外部診療の診療情報データベース、１１４…外部診療機関の遺伝子多型情報データベース。
【図１】
2005276022.tif 000021
【図２】
2005276022.tif 000022
【図３】
2005276022.tif 000023
【図４】
2005276022.tif 000024
【図５】
2005276022.tif 000025
【図６】
2005276022.tif 000026
【図７】
2005276022.tif 000027
【図８】
2005276022.tif 000028
【図９】
2005276022.tif 000029
【図１０】
2005276022.tif 000030
【図１１】
2005276022.tif 000031