(57)【要約】【課題】ヘテロ型の「ハプロタイプ」の組に由来する「ディプロタイプ」相互の区別を可能とする、あるいは、区別不能な「ディプロタイプ」を最小限にするための、変異セットの選別法の提供。【解決手段】識別すべき「ハプロタイプ」の種類ＮHに対して、変異セットの多型部位数Ｎを（ＮH−１）≦２N-1を満たす範囲に設定し、「ハプロタイプ」相互の識別が可能な「変異セット」候補を複数選別する第１の工程を具えている。さらに、第１の工程で選別される、複数の「変異セット」候補について、「ディプロタイプ」の表現型をも識別可能か否かの検証を行う第２の工程を設けている。この第１の工程、第２の工程の二つの関門を通過したものを、「最小数」の多型部位で構成される「変異セット」として、選定する方法としている。【選択図】なし【特許請求の範囲】
【請求項１】
ハプロタイプの組み合わせであるディプロタイプを識別する変異セットの選別法において、
ハプロタイプ相互の識別が可能な最小の変異セットを複数選別する第一の工程と、
前記第一の工程で選別された複数の変異セット候補から、ディプロタイプの表現型をも識別可能な最小の変異セットを選別する第二の工程とを有する
ことを特徴とする変異セットの選別法。
【請求項２】
Ａ（アデニン）、Ｇ（グアニン）、Ｔ（チミン）、Ｃ（シトシン）の塩基配列で構成される複数の個体のＤＮＡ配列についてアライメントし、
アライメントした結果、１％以上の頻度で見られる変異を多型と定義し、
定義された多型部位を複数箇所、個々の個体を識別する識別子として、該識別子を用いて、前記変異セットが定義されている
ことを特徴とする請求項１に記載の変異セットの選別法。
【請求項３】
前記変異セットは、
個々の個体が有するハプロタイプを識別する識別子として、
定義された多型部位の複数箇所について、野生型アレル、変異型アレルを識別可能な文字表記を行って、ハプロタイプ配列の形態で定義されている
ことを特徴とする請求項２に記載の変異セットの選別法。
【請求項４】
前記変異セットは、
多型と定義される変異として、
一塩基多型、インサーション、デリーション、リピート配列を含むセットである
ことを特徴とする請求項２または３に記載の変異セットの選別法。
【請求項５】
前記ハプロタイプ相互の識別が可能な最小の変異セットを複数選別する第一の工程では、
ハプロタイプ配列の形態で定義されている個体数分の変異セットを、並列することで、各個体のハプロタイプ配列を行に、定義された多型部位の複数箇所の各変異を列として、行列を構築した上で、
特定の個体のハプロタイプ配列を識別するために、ハプロタイプの情報量として必要十分な最小数の変異セットを複数抽出し、
該ハプロタイプの情報量として必要十分な最小数の変異セット複数を、ディプロタイプの表現型をも識別可能な最小の変異セットの候補として、選別する
ことを特徴とする請求項３に記載の変異セット選別方法。
【請求項６】
前記ハプロタイプ相互の識別が可能な最小の変異セットを複数選別する第一の工程において、
特定の個体のハプロタイプ配列を識別するために、ハプロタイプの情報量として必要十分な最小数の変異セット複数の抽出は、下記のステップ（ａ）〜（ｉ）を含む操作によって行われる
ことを特徴とする請求項５に記載の変異セット選別方法。
（ａ）複数の変異セットの候補を多型情報データベースから抽出し、変異部位だけで構成されたハプロタイプ配列を取得するステップ；
（ｂ）識別したい複数の個体の配列行に対して、アライメントされた変異セットの列について、列方向に野生型か変異型かを表現するパターンが同一、もしくは反転しているセットが他にある場合には、その列のセットを除外するステップ；
（ｃ）識別したい複数の個体の配列行に対して、一つ参照行を決め、各変異セットの列方向で参照行と異なる要素がある場合に、消去し参照行と識別できることを表示するステップ；
（ｄ）順に全ての個体の配列行に対して、参照行として、ステップ（ｃ）を繰り返し、全ての配列行の組合せについて比較し、消去演算するステップ；
（ｅ）全ての個体の配列行について、行方向に消去された要素数をカウントし、最も少ない配列行の消去された列を変異セットの第１の候補とするステップ；
（ｆ）第１の候補の変異列を固定し、他の異なる列と組合せ、全ての組合せで行方向にスクリーニングして、消去した要素を少なくとも１つ含む場合には、その行を消去し他の行と識別できることを表示し、第１の候補が複数ある場合には、その各々に対して処理を実行するステップ；
（ｇ）全ての個体の配列行について、行方向に消去された要素数をカウントし、最も少ない配列行の消去された列セットを、変異セットの第２の候補とし、第２の候補は、第１の候補を継承するステップ；
（ｈ）第２の候補の変異列セットを固定し、他の異なる列と組合せ、ステップ（ｆ）からステップ（ｇ）の操作を繰り返すステップ；
（ｉ）候補の変異列セットにおいて、列方向の要素が全て消去されると、一個体の特定配列を識別するための必要十分な最小数の変異セットとして決定するステップ。
【請求項７】
前記第一の工程で選別された複数の変異セット候補から、ディプロタイプの表現型をも識別可能な最小の変異セットを選別する第二の工程では、
特定の個体のハプロタイプ配列を識別するために、ハプロタイプの情報量として必要十分な最小数の変異セット複数として抽出された、複数の変異セット候補から、個人のディプロタイプが重複しないように、最小数の最適な変異セットを抽出する
ことを特徴とする請求項６に記載の変異セット選別方法。
【請求項８】
前記第一の工程で選別された複数の変異セット候補から、ディプロタイプの表現型をも識別可能な最小の変異セットを選別する第二の工程においては、
ディプロタイプの表現型をも識別可能な最小の変異セットの選別は、
前記の最小数変異セットの候補に対して、全て個体の配列行の組合せについてディプロタイプを演算する下記の（ｊ）〜（ｌ）のステップを含む操作によって行われる
ことを特徴とする請求項７記載の変異セット選別方法。
（ｊ）野生型アレル、変異型アレルを識別し、同一列の要素を第１の演算方法で処理することでディプロタイプの要素として、野生型ホモ、変異型ホモ、ヘテロを出力するステップ。
（ｋ）前記のディプロタイプの組合せ全てに対して、同一変異列の要素を第２の演算方法で処理することで同一のディプロタイプが複数あるか否かを検証するステップ。
（ｌ）・同一ディプロタイプが複数ある場合には、当該の変異セットを候補から除外し、別箇の変異セットでステップ（ｊ）〜ステップ（ｌ）を繰り返す；
・全てのディプロタイプが単独で存在する場合には、その変異セットを最適化された候補とする；
あるいは
・全ての最小数変異セットについて候補がない場合には、変異セットの列数を増やして、前記工程１のステップ（ａ）からステップ（ｉ）、ならびにステップ（ｊ）からステップ（ｌ）の処理を繰り返す；
前記の三種の操作のいずれかを行うステップ。
【請求項９】
前記ステップ（ｊ）の第１の演算方法の処理において、
前記の変異セットの候補に対して、野生型アレル、変異型アレルを識別可能な文字表記したハプロタイプ配列について、各配列の異なる全ての組合せの間で和演算処理をし、その演算結果により、野生型ホモ、変異型ホモ、ヘテロを識別し、それをディプロタイプ配列とする
ことを特徴とする請求項８に記載の変異セット選別方法。
【請求項１０】
前記ステップ（ｋ）の第２の演算方法の処理において、
前記のディプロタイプ配列について異なる全ての組合せの間で減演算処理をし、その演算結果がすべて０になる組合せについて、異なるハプロタイプ配列のセットで同一のディプロタイプ配列を表現する組み合わせが存在すると判断する
ことを特徴とする請求項９に記載の変異セット選別方法。
【請求項１１】
前記ステップ（ｋ）の第２の演算方法の処理において、
前記ディプロタイプ配列について異なる全ての組合せに間で比較演算処理をし、その演算結果が全て同値である組合せについて、異なるハプロタイプ配列のセットで同一のディプロタイプ配列を表現する組み合わせが存在すると判断する
ことを特徴とする請求項９に記載の変異セット選別方法。
【請求項１２】
異なる全てのハプロタイプの組合せが異なるディプロタイプ型を示すように選択された変異セット、もしくは、異なるハプロタイプの組合せによって同じディプロタイプ型を示す組数が最小数になるように選択された変異セットを選別する
ことを特徴とする請求項８〜１１のいずれか一項に記載の変異セット選別方法。
【請求項１３】
ハプロタイプ相互の識別、ならびに、ディプロタイプの表現型をも識別可能な最小の変異セットに対して、該最小の変異セットを構成する、多型部位の複数箇所について、それぞれの塩基の種別を特定するために利用されるプローブ・セットを設計する方法であって、
該プローブ・セットの設計工程は、下記のステップＡとステップＢを含む
ことを特徴とするプローブ・セットの設計方法。
ステップＡ：
請求項１〜１１のいずれか一項に記載の変異セットの選別方法を適用して、最小の変異セットを選定する。
ステップＢ：
選定された最小の変異セットを構成する各多型部位について、該多型部位をその塩基配列中に含ように、プローブの塩基配列を選択する。
【請求項１４】
ハプロタイプ相互の識別、ならびに、ディプロタイプの表現型をも識別可能な最小の変異セットに対して、該最小の変異セットを構成する、多型部位の複数箇所について、それぞれの塩基の種別を特定するために利用されるプライマー・セットを設計する方法であって、
該プライマー・セットの設計工程は、下記のステップＡとステップＢを含む
ことを特徴とするプローブ・セットの設計方法。
ステップＡ：
請求項１〜１１のいずれか一項に記載の変異セットの選別方法を適用して、最小の変異セットを選定する。
ステップＢ：
選定された最小の変異セットを構成する各多型部位について、該多型部位の塩基の種類を特定可能なように、プライマーの塩基配列を選択する。
【請求項１５】
請求項１〜１１のいずれか一項に記載される変異セットの選別方法を適用して、最小の変異セットを選定するための、一連の数値演算処理、ならびに、該数値演算処理結果を用いる論理演算処理をコンピュータシステム上で実行するためのプログラムを、該コンピュータシステムにおいて実行可能なプログラム言語により表記したソフトウエアを記録媒体上に電子的に収納してなる
ことを特徴とする、ソフトウエア収納記録媒体。
【請求項１６】
ハプロタイプの組み合わせであるディプロタイプを識別する変異セットを選別する作業を実施するための、コンピュータシステムであって、
該システムは、
ハプロタイプ相互の識別が可能な最小の変異セットを複数選別する、第一の演算処理を実施する第一の工程と、
前記第一の工程で選別された複数の変異セット候補から、ディプロタイプの表現型をも識別可能な最小の変異セットを選別する、第二の演算処理を実施する第二の工程とを具えている
ことを特徴とする変異セットの選別システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の配列情報の中から特定配列を特徴づけている要素を抽出する情報処理方法に関する。特に、遺伝子多型情報の組み合わせ（以降、ハプロタイプと呼ぶ）を構成する変異セットを最適化する情報処理技術に関する。
【背景技術】
【０００２】
分類学上は同一種に分類される生物群においても、個体間の比較を行うと、その遺伝子に起因する多様性を示す。具体的には、特定の遺伝子の塩基配列を個体間で比較すると、塩基配列の相違が離散的に存在する。特に、遺伝子の塩基配列中の特定部位の塩基が、個体群全体において、ある確率分布で相違している場合がある。この一塩基の変異は、「一塩基多型（ＳＮＰ）」と称されている。一方、ゲノムＤＮＡ上には、多数の遺伝子が存在しているが、「一塩基多型」を示す遺伝子が複数個存在している。さらに、一つの遺伝子上に、「一塩基多型（ＳＮＰ）」を示す部位が、複数箇所存在することも少なくない。
【０００３】
遺伝子の塩基配列中に「一塩基多型（ＳＮＰ）」が存在する場合、該遺伝子の発現により産生される産物、例えば、タンパク質、ペプチド性生理物質、ＲＮＡに差違が生じることもある。あるいは、産生される産物自体には、差違は存在しないが、その産生量に差違が生じる場合もある。その結果として、個体を示す体質に相違が生じる場合も知られている。すなわち、個体の示す体質に影響する遺伝的要因として、遺伝子の多型（ＳＮＰｓ）が着目されている。個体の遺伝子多型を知ることで、例えば、薬剤の代謝に関連する情報から直接疾患要因を特定する情報まで、種々の情報が得られる。その目的で、個体から採取されたゲノムＤＮＡ試料を利用して、その個体が有する遺伝子の多型を検出し、特定するための、遺伝子検査薬、キット、装置が開発されている。
【０００４】
一方、既に判明している遺伝子の多型（ＳＮＰｓ）の総数は、膨大な数に達しつつある。また、薬剤代謝、疾患罹患の危険性など、個体の体質的要因の決定は、代謝酵素タンパク質群や疾患関連タンパク質群など、複数種のタンパク質が関与していることも少なくない。個々の遺伝子、例えば、代謝酵素遺伝子や、疾患関連遺伝子においても、それぞれの遺伝子中には、多数の多型が存在して、この「複数個の多型の組み合わせ」に依って、その遺伝子産物に顕著な差違が生じることが少なくない。すなわち、特定の遺伝子から発現されるタンパク質のアミノ酸配列を対比すると、複数種の変異体が存在しており、これらの変異体をコードする遺伝子は、特定の「複数個の多型の組み合わせ」を有することが多い。この「複数個の多型の組み合わせ」を構成する個々の多型は、「連鎖」しており、一つの「ハプロタイプ」を構成している。「ハプロタイプ」を構成する複数個の多型に関しては、その「連鎖マップ」に従って、系統的に多型を検出し、特定することが重要となっている。
【０００５】
従来、遺伝子の多型情報を抽出する方法としては、先ず、多数の個体の遺伝子の塩基配列をアライメントし、変異が見出される部位（多型部位）を特定する。次いで、該多型部位に変異を有する遺伝子に関連した表現系と、検出された多型との相関を検証した上で、該遺伝子のアレル（対立遺伝子）群を特定する。特定されたアレル（対立遺伝子）の種類が増えてくると、該遺伝子中に存在することが確認された多型部位の総数も増えてくる。しかしながら、各アレル（対立遺伝子）で見出される、複数個の特定の多型は、全部もしくは一部「連鎖」しており、個々のアレルを判別する際に、全ての多型部位を調査する必要がない場合も少なくない。
【０００６】
この点に着目し、対象の遺伝子において、既に確認され、また、その塩基の出現頻度も確定している、多数の多型部位の中で、特定の「ハプロタイプ」の決定に必須な「最小数」の多型部位を含む変異セットを選択する手法が開発されている。例えば、特許文献１に開示される「配列を同定するための変異要素のセット抽出方法」を、該遺伝子の塩基配列に応用すると、特定の「ハプロタイプ」の決定に必須な「最小数」の多型部位のセット（変異セット）が取得される。
【特許文献１】特開２００５−１９０４２７号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
特許文献１に開示される手法により取得される、「最小数」の多型部位のセット（変異セット）は、評価対象の核酸分子が、該遺伝子の既知の「ハプロタイプ」のいずれに属するかを決定する目的としては、十分なものである。一方、個体のゲノムが、２個一組の相同染色体ｎ組からなる２倍体型の構成である場合、一組の相同染色体上の、相同の遺伝子座に、一対のアレル（対立遺伝子）が存在している。すなわち、対象の遺伝子Ｃには、複数種（ＮH種）のアレル（対立遺伝子）｛Ｃｉ：ｉ＝１〜ＮH｝が存在しており、各個体は、ホモ型の（Ｃi、Ｃi）、あるいは、ヘテロ型の（Ｃi、Ｃj）の組み合わせを有することになる。ホモ型の（Ｃi、Ｃi）の場合には、前記「最小数」の多型部位のセット（変異セット）を用いて、その「ハプロタイプ」を特定することが可能である。一方、ヘテロ型の（Ｃi、Ｃj）の場合には、その組み合わせによっては、「最小数」の多型部位のセット（変異セット）により判定される多型の種類を比較する限り、同一の「ディプロタイプ」を表現することがある。従って、そのような場合には、前記「最小数」の多型部位のセット（変異セット）を用いて、検出される多型情報のみからは、該個体の「ハプロタイプ」、すなわち、（Ｃi、Ｃj）の組み合わせを確定できない。その際には、対象とする個体の属する家系、集団（人種、民族）における各アレルの出現頻度（確率分布）を考慮して、該「ディプロタイプ」を表現する複数種の「ハプロタイプ」の組み合わせ（Ｃi、Ｃj）間で、その確率を統計的に推定する。例えば、各個体の有する「ハプロタイプ」の組み合わせ（Ｃi、Ｃj）は、両親由来の２種の配偶子に由来することを考慮して、両家系における各アレルの出現頻度（確率分布）に基づき、各（Ｃi、Ｃj）の組み合わせが出現する確率を統計的に推定する。その推定結果において、最も高い確率を示す（Ｃi、Ｃj）の組み合わせを、対象の個体は有すると推断していた。
【０００８】
一方、対象とする遺伝子において、確認されている多型部位の総数を、ＮSNPとし、その多型部位（ｃk：ｋ＝１〜ＮSNP）における、塩基の出現パターン（ｃ1，…，ｃk，…，ｃNSNP）を考える。各多型部位（ｃk：ｋ＝１〜ＮSNP）に出現する塩基の種類が２である場合、可能な出現パターン（ｃ1，…，ｃk，…，ｃNSNP）の総数は、２の（ＮSNP）乗となる。しかしながら、実際に、アレル（対立遺伝子）として、確認されている、出現パターン（ｃ1，…，ｃk，…，ｃNSNP）の種類は、前記の最大数（２NSNP）よりも、格段に少ない数（ＮH）である。この点を利用して、「最小数」の多型部位のセット（変異セット）、すなわち、ＮSNPより遥かに小さな「最小数」：Ｎmin個の多型部位のみで、ＮH種類のアレル（対立遺伝子）｛Ｃi：ｉ＝１〜ＮH｝の区別を行っている。すなわち、ＮH＜２Nminとなっている。
【０００９】
さらに、ホモ型の（Ｃi、Ｃi）、あるいは、ヘテロ型の（Ｃi、Ｃj）の組み合わせの総数は、｛ＮH＋1/2ＮH×（ＮH−１）｝＝1/2ＮH×（ＮH＋１）である。従って、少なくとも、｛1/2ＮH×（ＮH＋１）｝＜1/2（２Nmin）×（２Nmin＋１）となっている。1/2（２Nmin）×（２Nmin＋１）＜２NSNPの条件が満たされる場合、変異セットを構成する、多型部位の数：Ｎsetを増し、その多型部位：ｃkの組み合わせを適正に選別すると、ヘテロ型（Ｃi、Ｃj）を区別することが可能となる場合もある。
【００１０】
しかしながら、ヘテロ型（Ｃi、Ｃj）の「ハプロタイプ」の組に由来する「ディプロタイプ」相互の区別を可能とする、あるいは、区別不能な「ディプロタイプ」を最小限にするための、変異セットの選別法は、未だ具現化されていなかった。
【００１１】
本発明は、前記課題を解決するものである。本発明の目的は、ヘテロ型（Ｃi、Ｃj）の「ハプロタイプ」の組に由来する「ディプロタイプ」相互の区別を可能とする、あるいは、区別不能な「ディプロタイプ」を最小限にするための、変異セットの選別法を提供することにある。
【課題を解決するための手段】
【００１２】
本発明では、各「ハプロタイプ」が、複数箇所の多型部位において、観測される塩基配列が、野生型アレル（Ｗ）か、変異型アレル（Ｍ）かの差違によって、相互に識別できる点に着目する。例えば、各「ハプロタイプ」：Ｃi：ｉ＝１〜ＮHは、複数箇所の多型部位：ｃk：ｋ＝１〜ＮSNPにおける塩基配列の区分：ｃk（ｉ）：ｋ＝１〜ＮSNPを示す、行ベクトル（…，ｃk（ｉ），…）を利用して、表記できる。その際、各「ハプロタイプ」：Ｃiを示す行ベクトル（…，ｃk（ｉ），…）は、該多型部位：ｃk：ｋ＝１〜ＮSNPの総数、「ＮSNP次元のベクトル空間」上で、互いに相違する行ベクトル群を構成している。一方、二つの「ハプロタイプ」の組み合わせ（Ｃi，Ｃj）に起因する「ディプロタイプ」の表現型も、「ＮSNP次元のベクトル空間」上で、行ベクトル（…，｛ｃk（ｉ）＋ｃk（ｊ）｝，…）と表記できる。
【００１３】
この表記法を採用する際、各「ハプロタイプ」、ならびに「ハプロタイプ」の組（Ｃi，Ｃj）に起因する「ディプロタイプ」の表現型の識別も可能な、「最小数」の多型部位で構成される「変異セット」は、下記の手順で選別できることを見出した。
【００１４】
すなわち、本発明にかかる変異セットの選別法では、
ハプロタイプの組み合わせであるディプロタイプを識別する変異セットの選別法において、下記の二つの工程を有することを特徴としている。
【００１５】
・ハプロタイプ相互の識別が可能な最小の変異セットを複数選別する第一の工程。
【００１６】
・前記第一の工程で選別された複数の変異セット候補から、ディプロタイプの表現型をも識別可能な最小の変異セットを選別する第二の工程。
【発明の効果】
【００１７】
本発明にかかる方法により選別される「最小数」の多型部位で構成される「変異セット」は、複数の「ハプロタイプ」相互の識別に加えて、相同染色体上に一組の「ハプロタイプ」を有する個体が示す「ディプロタイプ」相互の識別も可能である。すなわち、個体が示す「ディプロタイプ」相互の識別により、該「ディプロタイプ」の表現を与える、一組の「ハプロタイプ」の組み合わせを特定することが可能である。有性生殖による繁殖を行う生物では、各個体の有する遺伝子型は、二つの配偶子の由来する相同染色体上に存在するアレル（対立遺伝子）の組み合わせに因って決定される。該相同の遺伝子座に存在するアレル（対立遺伝子）は、父系と母系の二つの系統に由来している。この二つのアレル（対立遺伝子）に起因する発現がなされると、例えば、ｍＲＮＡ中には、両者の塩基配列を有する２種のｍＲＮＡ分子が共存している「ディプロタイプ」の表現型を示す。該遺伝子中に存在する多型部位の塩基配列を判定すると、「ディプロタイプ」の表現型においては、一組の「ハプロタイプ」が示す多型部位の情報が総合された状態で観測される。該「最小数」の多型部位で構成される「変異セット」は、一組の「ハプロタイプ」の組み合わせの相違に起因する、「ディプロタイプ」の表現型の相違を識別することで、その個体における「ハプロタイプ」の組み合わせを特定することを可能にする。加えて、一組の「ハプロタイプ」の組み合わせの相違に起因する、「ディプロタイプ」の表現型の相違を識別する上で、測定すべき多型部位の個数を「最小数」とする利点を有している。この「最小数」の多型部位で構成される「変異セット」を採用すると、各「ハプロタイプ」が示す多型部位を検出用のプローブ・セットや、プライマー・セットの設計も、より効率的に行うことが可能になる。
【００１８】
本発明にかかる方法では、「ディプロタイプ」の表現型の識別を可能とする、測定すべき多型部位の個数の下限を、識別すべき「ハプロタイプ」の種類数より決定する。そして、この下限条件を満足する、複数の「ハプロタイプ」相互の識別可能な、「最小数」の多型部位で構成される「変異セット」を複数の候補として選別し、その候補の中から、「ディプロタイプ」の表現型の相違を識別することが可能なものを選定する。従って、複数の「ハプロタイプ」相互の識別は勿論、一組の「ハプロタイプ」の組み合わせの相違に起因する、「ディプロタイプ」の表現型の相違をも識別可能な「最小数」の多型部位で構成される「変異セット」を効率的に選別できる。
【発明を実施するための最良の形態】
【００１９】
本発明にかかる、「ハプロタイプ」のみでなく、「ディプロタイプ」の識別も可能である、「最小数」の多型部位で構成される「変異セット」を選別する方法を以下に詳しく説明する。加えて、該選別方法を適用して、選別される「最小数」の多型部位で構成される「変異セット」に対して、その多型部位の塩基を判定するためのプローブ・セットの設計方法、あるいは、塩基配列解析用プライマー・セットの設計方法を併せて説明する。
【００２０】
また、該「最小数」の多型部位で構成される「変異セット」を選別する方法に基づく、「変異セット」選別用ソフトウエア、該ソフトウエアを実装し、「変異セット」選別の操作を行うシステムに関しても、詳しい説明を行う。同じく、前記プローブ・セットの設計方法に基づく、プローブ・セット設計用ソフトウエア、該ソフトウエアを実装し、プローブ・セット設計の操作を行うシステムに関しても、詳しい説明を行う。前記塩基配列解析用プライマー・セットの設計方法に基づく、プライマー・セット設計用ソフトウエア、該ソフトウエアを実装し、プライマー・セット設計の操作を行うシステムに関しても、詳しい説明を行う。
【００２１】
まず、本発明の変異セットの選別法では、
ハプロタイプの組み合わせであるディプロタイプを識別する変異セットの選別法において、下記の二つの工程を有することを特徴としている。
【００２２】
・ハプロタイプ相互の識別が可能な最小の変異セットを複数選別する第一の工程。
【００２３】
・前記第一の工程で選別された複数の変異セット候補から、ディプロタイプの表現型をも識別可能な最小の変異セットを選別する第二の工程。
【００２４】
その際、「変異セット」として、下記のように定義される変異セットを用いることが好ましい。まず、Ａ（アデニン）、Ｇ（グアニン）、Ｔ（チミン）、Ｃ（シトシン）の塩基配列で構成される複数の個体のＤＮＡ配列についてアライメントし、アライメントした結果、１％以上の頻度で見られる変異を多型と定義する。次に、定義された多型部位を複数箇所、個々の個体を識別する識別子として、選択し、この識別子を用いて、変異セットを定義する。
【００２５】
特に、前記変異セットは、個々の個体が有するハプロタイプを識別する識別子として、定義された多型部位の複数箇所について、野生型アレル、変異型アレルを識別可能な文字表記を行って、ハプロタイプ配列の形態で定義することが好ましい。なお、前記変異セットは、多型と定義される変異として、一塩基多型、インサーション、デリーション、リピート配列を含むセットとすることができる。
【００２６】
一方、前記ハプロタイプ相互の識別が可能な最小の変異セットを複数選別する第一の工程では、下記の操作を行う形態とすることが好ましい。
【００２７】
・ハプロタイプ配列の形態で定義されている個体数分の変異セットを、並列することで、各個体のハプロタイプ配列を行に、定義された多型部位の複数箇所の各変異を列として、行列を構築する。
【００２８】
・その上で、特定の個体のハプロタイプ配列を識別するために、ハプロタイプの情報量として必要十分な最小数の変異セットを複数抽出する。
【００２９】
・該ハプロタイプの情報量として必要十分な最小数の変異セット複数を、ディプロタイプの表現型をも識別可能な最小の変異セットの候補として、選別する。
【００３０】
例えば、前記る第一の工程において、特定の個体のハプロタイプ配列を識別するために、ハプロタイプの情報量として必要十分な最小数の変異セット複数の抽出は、下記のステップ（ａ）〜（ｉ）を含む操作によって行われることが好ましい。
【００３１】
ステップ（ａ）：複数の変異セットの候補を多型情報データベースから抽出し、変異部位だけで構成されたハプロタイプ配列を取得する。
【００３２】
ステップ（ｂ）：識別したい複数の個体の配列行に対して、アライメントされた変異セットの列について、列方向に野生型か変異型かを表現するパターンが同一、もしくは反転しているセットが他にある場合には、その列のセットを除外する。
【００３３】
ステップ（ｃ）：識別したい複数の個体の配列行に対して、一つ参照行を決め、各変異セットの列方向で参照行と異なる要素がある場合に、消去し参照行と識別できることを表示する。
【００３４】
ステップ（ｄ）：順に全ての個体の配列行に対して、参照行として、ステップ（ｃ）を繰り返し、全ての配列行の組合せについて比較し、消去演算する。
【００３５】
ステップ（ｅ）：全ての個体の配列行について、行方向に消去された要素数をカウントし、最も少ない配列行の消去された列を変異セットの第１の候補とする。
【００３６】
ステップ（ｆ）：第１の候補の変異列を固定し、他の異なる列と組合せ、全ての組合せで行方向にスクリーニングして、消去した要素を少なくとも１つ含む場合には、その行を消去し他の行と識別できることを表示する。逆に、第１の候補が複数ある場合には、その各々に対して処理を実行する。
【００３７】
ステップ（ｇ）：全ての個体の配列行について、行方向に消去された要素数をカウントし、最も少ない配列行の消去された列セットを、変異セットの第２の候補とする。第２の候補は、第１の候補を継承する。
【００３８】
ステップ（ｈ）：第２の候補の変異列セットを固定し、他の異なる列と組合せ、ステップ（ｆ）からステップ（ｇ）の操作を繰り返す。
【００３９】
ステップ（ｉ）：候補の変異列セットにおいて、列方向の要素が全て消去されると、一個体の特定配列を識別するための必要十分な最小数の変異セットとして決定する。
【００４０】
また、前記第一の工程で選別された複数の変異セット候補から、ディプロタイプの表現型をも識別可能な最小の変異セットを選別する第二の工程では、下記の操作を行う形態とすることが好ましい。すなわち、特定の個体のハプロタイプ配列を識別するために、ハプロタイプの情報量として必要十分な最小数の変異セット複数として抽出された、複数の変異セット候補から、個人のディプロタイプが重複しないように、最小数の最適な変異セットを抽出する。
【００４１】
例えば、前記第二の工程においては、ディプロタイプの表現型をも識別可能な最小の変異セットの選別は、以下の手順で行うことが好ましい。すなわち、前記の最小数変異セットの候補に対して、全て個体の配列行の組合せについてディプロタイプを演算する下記の（ｊ）〜（ｌ）のステップを含む操作によって行うことが好ましい。
【００４２】
ステップ（ｊ）：野生型アレル、変異型アレルを識別し、同一列の要素を第１の演算方法で処理することでディプロタイプの要素として、野生型ホモ、変異型ホモ、ヘテロを出力する。
【００４３】
ステップ（ｋ）：前記のディプロタイプの組合せ全てに対して、同一変異列の要素を第２の演算方法で処理することで同一のディプロタイプが複数あるか否かを検証する。
【００４４】
ステップ（ｌ）：下記の三種の操作のいずれかを行う。
【００４５】
・同一ディプロタイプが複数ある場合には、当該の変異セットを候補から除外し、別箇の変異セットでステップ（ｊ）〜ステップ（ｌ）を繰り返す。
【００４６】
・全てのディプロタイプが単独で存在する場合には、その変異セットを最適化された候補とする。
【００４７】
あるいは
・全ての最小数変異セットについて候補がない場合には、変異セットの列数を増やして、前記工程１のステップ（ａ）からステップ（ｉ）、ならびにステップ（ｊ）からステップ（ｌ）の処理を繰り返す。
【００４８】
例えば、前記ステップ（ｊ）の第１の演算方法の処理において、
前記の変異セットの候補に対して、野生型アレル、変異型アレルを識別可能な文字表記したハプロタイプ配列について、各配列の異なる全ての組合せの間で和演算処理を行う。その演算結果により、野生型ホモ、変異型ホモ、ヘテロを識別し、それをディプロタイプ配列とすることが好ましい。
【００４９】
また、前記ステップ（ｋ）の第２の演算方法の処理において、
前記のディプロタイプ配列について異なる全ての組合せの間で減演算処理を行う。その演算結果がすべて０になる組合せについて、異なるハプロタイプ配列のセットで同一のディプロタイプ配列を表現する組み合わせが存在すると判断することができる。
【００５０】
あるいは、前記ステップ（ｋ）の第２の演算方法の処理において、
前記ディプロタイプ配列について異なる全ての組合せに間で比較演算処理を行う。その演算結果が全て同値である組合せについて、異なるハプロタイプ配列のセットで同一のディプロタイプ配列を表現する組み合わせが存在すると判断することもできる。
【００５１】
最終的に、異なる全てのハプロタイプの組合せが異なるディプロタイプ型を示すように選択された変異セットを選別することが望ましい。もしくは、異なるハプロタイプの組合せによって同じディプロタイプ型を示す組数が最小数になるように選択された変異セットを選別することもできる。
【００５２】
また、本発明では、ハプロタイプ相互の識別、ならびに、ディプロタイプの表現型をも識別可能な最小の変異セットに対して、その変異検出用のプローブ・セットを設計する方法を提供する。すなわち、該最小の変異セットを構成する、多型部位の複数箇所について、それぞれの塩基の種別を特定するために利用されるプローブ・セットを設計する方法を提供する。本発明にかかるプローブ・セットの設計方法における、該プローブ・セットの設計工程は、下記のステップＡとステップＢを含む。
【００５３】
ステップＡ：
上で述べた本発明の変異セットの選別方法を適用して、最小の変異セットを選定する。
【００５４】
ステップＢ：
選定された最小の変異セットを構成する各多型部位について、該多型部位をその塩基配列中に含ように、プローブの塩基配列を選択する。
【００５５】
さらに、本発明では、ハプロタイプ相互の識別、ならびに、ディプロタイプの表現型をも識別可能な最小の変異セットに対して、その変異検出用のプライマー・セットを設計する方法を提供する。すなわち、該最小の変異セットを構成する、多型部位の複数箇所について、それぞれの塩基の種別を特定するために利用されるプライマー・セットを設計する方法を提供する。本発明にかかるプライマー・セットの設計方法における、該プライマー・セットの設計工程は、下記のステップＡとステップＢを含む。
【００５６】
ステップＡ：
上で述べた本発明の変異セットの選別方法を適用して、最小の変異セットを選定する。
【００５７】
ステップＢ：
選定された最小の変異セットを構成する各多型部位について、該多型部位の塩基の種類を特定可能なように、プライマーの塩基配列を選択する。
【００５８】
なお、本発明にかかる変異セットの選別方法は、一連の数値演算処理、ならびに、該数値演算処理結果を用いる論理演算処理を行うことで、最小の変異セットを特定する形態として実施できる。これらの演算処理を含む操作は、コンピュータシステム上で実行するため、専用のソフトウエアとして、記録媒体上に記録する形態にすることが望ましい。
【００５９】
従って、本発明は、前記の「専用のソフトウエア」を収納する、ソフトウエア収納記録媒体の発明をも提供する。
【００６０】
すなわち、本発明のソフトウエア収納記録媒体は、下記のような構成上の特徴を有する。
【００６１】
・上述する本発明の変異セットの選別方法を適用して、最小の変異セットを選定するためのものである。
【００６２】
・その方法における、一連の数値演算処理、ならびに、該数値演算処理結果を用いる論理演算処理をコンピュータシステム上で実行するためのプログラムを、主要な要素としている。
【００６３】
・そのプログラムを、該コンピュータシステムにおいて実行可能なプログラム言語により表記したソフトウエアを記録媒体上に電子的に収納してなるものである。
【００６４】
また、本発明の変異セットの選別方法は、上述するように、「専用のソフトウエア」を用いて、コンピュータシステム上において、その設計操作を実施することが一般的である。従って、本発明は、本発明の変異セットの選別方法を適用して、最小の変異セットを選定する操作を実施する用途に適するコンピュータシステムの発明をも提供する。
【００６５】
すなわち、本発明にかかる変異セットの選別システムは、下記の構成上の特徴を有する。
【００６６】
・ハプロタイプの組み合わせであるディプロタイプを識別する変異セットを選別する作業を実施するための、コンピュータシステムである。
【００６７】
・該システムは、下記の二つの工程を具えている。
【００６８】
ハプロタイプ相互の識別が可能な最小の変異セットを複数選別する、第一の演算処理を実施する第一の工程。
【００６９】
前記第一の工程で選別された複数の変異セット候補から、ディプロタイプの表現型をも識別可能な最小の変異セットを選別する、第二の演算処理を実施する第二の工程。
【００７０】
なお、本発明にかかる変異セットの選別システムにおいても、上記の本発明の変異セットの選別方法における好適な構成を選択すると、好ましい形態となる。
【００７１】
例えば、本発明の一態様では、各「ハプロタイプ」が、複数箇所の多型部位において、観測される塩基配列が、野生型アレル（Ｗ）か、変異型アレル（Ｍ）かの差違によって、相互に識別できる点に着目する。すなわち、各「ハプロタイプ」：Ｃi：ｉ＝１〜ＮHは、複数箇所の多型部位：ｃk：ｋ＝１〜ＮSNPにおける塩基配列の区分：ｃk（ｉ）：ｋ＝１〜ＮSNPを示す、行ベクトル（…，ｃk（ｉ），…）を利用して、表記できる。その際、各「ハプロタイプ」：Ｃiを示す行ベクトル（…，ｃk（ｉ），…）は、該多型部位：ｃk：ｋ＝１〜ＮSNPの総数、「ＮSNP次元のベクトル空間」上で、互いに相違する行ベクトル群を構成している。一方、二つの「ハプロタイプ」の組み合わせ（Ｃi，Ｃj）に起因する「ディプロタイプ」の表現型も、「ＮSNP次元のベクトル空間」上で、行ベクトル（…，｛ｃk（ｉ）＋ｃk（ｊ）｝，…）と表記できる。
【００７２】
この表記法を採用する際、各「ハプロタイプ」、ならびに「ハプロタイプ」の組（Ｃi，Ｃj）に起因する「ディプロタイプ」の表現型の識別も可能な、「最小数」の多型部位で構成される「変異セット」は、下記の手順で選別できることを見出した。
【００７３】
まず、「最小数」：Ｎの多型部位で構成される「変異セット」により、各「ハプロタイプ」：Ｃi、ならびに「ハプロタイプ」の組（Ｃi，Ｃj）に起因する「ディプロタイプ」の表現型を行ベクトルで表記することを考える。このＮ個の要素からなる行ベクトルは、基となる「ＮSNP次元のベクトル空間」上で定義される行ベクトル（…，ｃk（ｉ），…）を、「Ｎ次元のベクトル空間」に投影したものに相当する。その「Ｎ次元のベクトル空間」に投影した際、各「ハプロタイプ」：Ｃiを示す行ベクトルが、同一にならないことが、「ハプロタイプ」相互の識別が可能な「変異セット」の必要条件（ａ）となる。同様に、「Ｎ次元のベクトル空間」に投影した際、「ディプロタイプ」の表現型を表記する行ベクトルが、同一にならないことが、「ディプロタイプ」相互の識別が可能な「変異セット」の必要条件（ｂ）となる。
【００７４】
必要条件（ａ）は、「Ｎ次元のベクトル空間」において、「ハプロタイプ」の行ベクトルに可能な選択肢の総数：２N通りが、識別すべき「ハプロタイプ」の総数ＮHより多い：ＮH≦２Nと表現できる。一方、必要条件（ｂ）は、少なくとも、「Ｎ次元のベクトル空間」において、「ハプロタイプ」の行ベクトル中に、互いに逆ベクトルとなる組が、二つ以上存在していない条件とすることもできる。すなわち、「Ｎ次元のベクトル空間」において、「ハプロタイプ」の行ベクトルに可能な選択肢は、１／２・（２N）の互いに逆ベクトルとなる「ベクトル対」で構成されている。従って、必要条件（ｂ）は、「ハプロタイプ」：Ｃiを示すＮ次元のベクトル中に、互いに逆ベクトルとなる組が、二つ以上存在していない条件：（ＮH−１）≦１／２・２Nと表現できる。この必要条件（ｂ）：（ＮH−１）≦２N-1は、勿論、必要条件（ａ）：ＮH≦２Nをも満たしている。
【００７５】
従って、必要条件（ｂ）：（ＮH−１）≦２N-1を満たす、最小のＮを先ず決定する。次いで、基となる「ＮSNP次元のベクトル空間」上で定義される行ベクトル（…，ｃk（ｉ），…，ｃNSNP（ｉ））を、「Ｎ次元のベクトル空間」に投影した行ベクトル（…，ｃk'（ｉ），…，ｃN（ｉ））を考慮する。この各「ハプロタイプ」：Ｃiを示す行ベクトル（…，ｃk'（ｉ），…，ｃN（ｉ））が、互いに、相違したものとなるようなＮの多型部位で構成される「変異セット」（ｃk'）を複数選別する。選別される、「最小数」：Ｎの多型部位で構成される複数の「変異セット」を候補として、「ディプロタイプ」の表現型を表記する、Ｎ次元の行ベクトル（…，｛ｃk'（ｉ）＋ｃk'（ｊ）｝，…）を考慮する。その際、このＮ次元の行ベクトル（…，｛ｃk'（ｉ）＋ｃk'（ｊ）｝，…）相互が、同一にならないことを基準に、「最小数」：Ｎの多型部位で構成される「変異セット」を選定する。
【００７６】
本発明の一態様の「最小数」の多型部位で構成される「変異セット」の選別方法は、上記の必要条件（ｂ）：（ＮH−１）≦２N-1を満たす範囲で、「ハプロタイプ」相互の識別が可能な「変異セット」候補を複数選別する第１の工程を具えている。さらに、第１の工程で選別される、複数の「変異セット」候補について、「ディプロタイプ」の表現型をも識別可能か否かの検証を行う第２の工程を設けている。この第１の工程、第２の工程の二つの関門を通過したものを、「最小数」の多型部位で構成される「変異セット」として、選定する方法としている。
【００７７】
また、本発明の一態様にかかる、「最小数」の多型部位で構成される「変異セット」を選別する方法は、具体的には、その構成、ならびに、その好ましい形態は、下記のように記述することができる。
【００７８】
この態様においては、本発明にかかる「最小数」の多型部位で構成される「変異セット」を選別する方法は、下記の構成を有することが好ましい。
【００７９】
この態様の「変異セット」の選別方法は、
個体の有する相同染色体上、相同遺伝子座を占めるアレル（対立遺伝子）に関して、該個体が属する集団中において、複数種のアレル（対立遺伝子）の存在が確認されている際に適用する。
【００８０】
その複数種のアレル（対立遺伝子）各々に対してＡ（アデニン）、Ｇ（グアニン）、Ｔ（チミン）、Ｃ（シトシン）の塩基で構成される塩基配列に関して、相互の塩基配列をアライメントする。該塩基配列のアライメントによって対比した際、見出される複数の変異箇所中、前記集団中における出現頻度の確率が１％以上である複数の箇所を選択し、該遺伝子の多型部位と定義する。
【００８１】
選択された該遺伝子の多型部位の総数をＮsnp0と表記する。
【００８２】
前記複数種のアレル（対立遺伝子）相互の区別を行う指標として、該Ｎsnp0個の多型部位における塩基の種別を用い、各多型部位の塩基として、該集団中における出現頻度の最大の塩基を「野生型アレル（Ｗ）」、それ以外の塩基を「変異型アレル（Ｍ）」と定義する。
【００８３】
前記Ｎsnp0個の多型部位における塩基の種別の２値化処理を行って、前記複数種のアレル（対立遺伝子）を識別する識別子として、次の行ベクトルを考える。すなわち、前記Ｎsnp0個の多型部位：ｃk（ｋ＝１〜Ｎsnp0）における、２値化処理された塩基の種別で構成される、Ｎsnp0個の要素：ｃk（ｉ）（ｋ＝１〜Ｎsnp0）からなる、Ｎsnp0次元の行ベクトルを定義する。
【００８４】
前記複数種のアレル（対立遺伝子）を、該Ｎsnp0次元の行ベクトルに基づき分類し、該Ｎsnp0次元の行ベクトルが同一となるアレル（対立遺伝子）の群を複数種定義し、この定義された複数種のアレル（対立遺伝子）の群を該遺伝子の「ハプロタイプ」と称する。
【００８５】
定義される、該遺伝子の複数種の「ハプロタイプ」の総数をＮHと表記し、該ＮH種の「ハプロタイプ」をＣi：i＝１〜ＮHと表記する。
【００８６】
各個体中に存在する相同遺伝子座を占める、二つのアレル（対立遺伝子）の組み合わせに起因して、その「ハプロタイプ」の組み合わせ（Ｃi，Ｃj）に由来する、前記Ｎsnp0個の多型部位における、塩基の種別；ＷとＭの組み合わせの表現型を考える。そして、前記Ｎsnp0個の多型部位における、塩基の種別；ＷとＭの組み合わせの表現型を、各個体の「ディプロタイプ」：（Ｃi，Ｃj）と定義する。
【００８７】
本発明では、各個体の有する、二つのアレル（対立遺伝子）の組み合わせを分類する上で利用可能な、各「ハプロタイプ」相互の識別、同時に、その組み合わせである「ディプロタイプ」：（Ｃi，Ｃj）の相違の識別を行うことを目標とする。更に、各「ハプロタイプ」相互の識別と、「ディプロタイプ」：（Ｃi，Ｃj）の相違の識別が可能な、前記Ｎsnp0個の多型部位：ｃk（ｋ＝１〜Ｎsnp0）中から選択される、「最小数」の多型部位で構成される「変異セット」を選定する方法が、本発明で提供される。
【００８８】
その際、該選定方法は、下記の手順を有する。
【００８９】
前記「最小数」の多型部位で構成される「変異セット」を構成する多型部位の個数をＮと表記する際、
条件：（ＮH−１）≦２N-1を満たすＮの最小値：Ｎminを算定する。
【００９０】
その上で、前記Ｎsnp0個の多型部位：ｃk（ｋ＝１〜Ｎsnp0）中から、該Ｎmin個以上の多型部位を選択して、構成される「変異セット」の集合の中から、「最小数」の多型部位で構成される「変異セット」を選定する。
【００９１】
該「最小数」の多型部位で構成される「変異セット」を選定する工程は、下記のステップ（ａ）〜（ｉ）の操作を含んでなる第１の工程と、ステップ（ｊ）〜（ｌ）の操作を含んでなる第２の工程を含む。
【００９２】
（第１の工程）：各「ハプロタイプ」相互の識別を行うことが可能な、Ｎmin個以上の多型部位で構成される「変異セット」を、「最小数」の多型部位で構成される「変異セット」の候補として抽出する、下記のステップ（ａ）〜（ｉ）の操作を含んでなる工程；
ステップ（ａ）：
前記集団中において、その存在が確認されている複数種のアレル（対立遺伝子）の塩基配列の情報を、該遺伝子の多型情報を収納するデータベースから抽出する。抽出された塩基配列の情報に基づき、前記Ｎsnp0個の多型部位：ｃk（ｋ＝１〜Ｎsnp0）における、２値化処理された塩基の種別で構成される、Ｎsnp0個の要素：ｃk（ｉ）（ｋ＝１〜Ｎsnp0）からなる、Ｎsnp0次元の行ベクトルを作成する。次いで、作成されるＮsnp0次元の行ベクトルの集合中から、互いに、相違するＮsnp0次元の行ベクトルを示す、該ＮH種の「ハプロタイプ」：Ｃi：i＝１〜ＮHを表記する、Ｎsnp0次元の行ベクトルＮH種を選別する。
【００９３】
ステップ（ｂ）：
前記Ｎsnp0次元の行ベクトルＮH種を用いて、ＮH行×Ｎsnp0列の行列を形成し、その行列の各列より、Ｎsnp0個のＮH次元の列ベクトルを作成する。各列ベクトル間の対比を行い、相互に同一の列ベクトルの群、ならびに、相互に逆ベクトルの関係にある列ベクトルの群を選択する。選択される列ベクトルの群中の、一つの列ベクトルを選択し、該群に含まれる他の列ベクトルを除去し、前記ＮH行×Ｎsnp0列の行列から、各列ベクトルが相互に独立している列ベクトルにより構成される部分行列を作成する。該部分行列の列の総数をＮsnpと表記し、該ＮH行×Ｎsnp列の行列における、各列を与えている多型部位：ｃk（ｋ＝１〜Ｎsnp）を、該Ｎmin個以上の多型部位を選択して、構成される「変異セット」の可能な要素として選択する。
【００９４】
ステップ（ｃ）：
ＮH行×Ｎsnp列の行列において、各「ハプロタイプ」：Ｃiを表記する行ベクトル（ｃk（ｉ））を、参照行として選択する。他の「ハプロタイプ」：Ｃjを表記する行ベクトル（ｃk（ｊ））に対して、前記参照行の行ベクトル（ｃk（ｉ））と、各要素ｃk（ｊ）とｃk（ｉ）の比較を行う。比較する要素ｃk（ｊ）とｃk（ｉ）が、異なる場合、両者間の区別が可能として、ｃk（ｊ）を「−」の符号に置換する。この操作を、前記参照行の「ハプロタイプ」：Ｃiを表記する行ベクトル（ｃk（ｉ））以外の、（ＮH−１）個の行ベクトルに施す。この対比操作の結果を、ＮH行×Ｎsnp列の行列として、表記し、「ハプロタイプ」：Ｃiを表記する行ベクトル（ｃk（ｉ））を参照行とする、第１の対比結果の行列とする。
【００９５】
ステップ（ｄ）：
引き続き、各「ハプロタイプ」：Ｃiを表記する行ベクトル（ｃk（ｉ））を参照行として、同じ対比操作を行う。最終的に、各「ハプロタイプ」：Ｃiを表記する行ベクトル（ｃk（ｉ））を参照行とする、第１の対比結果の行列を、計ＮH個作成する。
【００９６】
ステップ（ｅ）：
作成される計ＮH個の第１の対比結果の行列を参照し、前記の対比操作が施された後の行、ＮH×（ＮH−１）のうち、「−」の符号へ置き換えが、最小数である行を選択する。この選択された行において、「−」の符号へ置き換えがなされている要素（列）を特定する。この特定される要素（列）を与えている多型部位を特定する。この特定される多型部位を、「最小数」の多型部位で構成される「変異セット」を構成する多型部位の「第１の候補」：ｃk1とする。
【００９７】
ステップ（ｆ）：
次いで、「第１の候補」：ｃk1と、他の多型部位：ｃk：ｃk≠ｃk1との組み合わせ：（ｃk1，ｃk）として、計（Ｎsnp−１）組を考える。
【００９８】
ＮH行×Ｎsnp列の行列において、各「ハプロタイプ」：Ｃiを表記する行ベクトルを、参照行として選択する。他の「ハプロタイプ」：Ｃjを表記する行ベクトルに対して、前記参照行の行ベクトルと、各要素対（ｃk1（ｉ），ｃk（ｉ））と（ｃk1（ｊ），ｃk（ｊ））の比較を行う。比較する要素対（ｃk1（ｉ），ｃk（ｉ））と（ｃk1（ｊ），ｃk（ｊ））が、異なる場合、両者間の区別が可能として、（ｃk1（ｊ），ｃk（ｊ））を「−」の符号に置換する。この操作を、前記参照行の「ハプロタイプ」：Ｃiを表記する行ベクトル以外の、（ＮH−１）個の行ベクトルに施す。この対比操作の結果を、列要素として、計（Ｎsnp−１）個の（ｃk1，ｃk）を有する、ＮH行×（Ｎsnp−１）列の行列として、表記する。このＮH行×（Ｎsnp−１）列の行列を、「ハプロタイプ」：Ｃiを表記する行ベクトルを参照行とする、第２の対比結果の行列とする。
【００９９】
引き続き、各「ハプロタイプ」：Ｃiを表記する行ベクトルを参照行として、同じ対比操作を行う。最終的に、各「ハプロタイプ」：Ｃiを表記する行ベクトルを参照行とする、第２の対比結果の行列を、計ＮH個作成する。
【０１００】
ステップ（ｇ）：
作成される計ＮH個の第２の対比結果の行列を参照し、前記の対比操作が施された後の行ＮH×（ＮH−１）のうち、「−」の符号へ置き換えが、最小数である行を選択する。この選択された行において、「−」の符号へ置き換えがなされている要素対（ｃk1，ｃk）を特定する。この特定される要素対（ｃk1，ｃk）を与えている、他の多型部位：ｃkを特定する。この特定される、他の多型部位：ｃkを、「最小数」の多型部位で構成される「変異セット」を構成する多型部位の「第２の候補」：ｃk2とする。
【０１０１】
ステップ（ｈ）：
次いで、「第１の候補」：ｃk1、「第２の候補」：ｃk2と他の多型部位：ｃkとの組み合わせ：（ｃk1，ｃk2，ｃk）として、計（Ｎsnp−２）組を考える。
【０１０２】
前記ステップ（ｆ）の対比操作に準じて、対比操作を進め、その結果を、列要素として、計（Ｎsnp−２）個の（ｃk1，ｃk2，ｃk）を有する、ＮH行×（Ｎsnp−２）列の行列として、表記する。このＮH行×（Ｎsnp−２）列の行列を、「ハプロタイプ」：Ｃiを表記する行ベクトルを参照行とする、第３の対比結果の行列とする。
【０１０３】
引き続き、各「ハプロタイプ」：Ｃiを表記する行ベクトルを参照行として、同じ対比操作を行う。最終的に、各「ハプロタイプ」：Ｃiを表記する行ベクトルを参照行とする、第３の対比結果の行列を、計ＮH個作成する。
【０１０４】
作成される計ＮH個の第３の対比結果の行列を参照し、前記の対比操作が施された後の行、ＮH×（ＮH−１）のうち、「−」の符号へ置き換えが、最小数である行を選択する。この選択された行において、「−」の符号へ置き換えがなされている要素組（ｃk1，ｃk2，ｃk）を特定する。この特定される要素組（ｃk1，ｃk2，ｃk）を与えている、他の多型部位：ｃkを特定する。この特定される、他の多型部位：ｃkを、「最小数」の多型部位で構成される「変異セット」を構成する多型部位の「第３の候補」：ｃk3とする。
【０１０５】
なお、前記「第２の候補」：ｃk2が複数選別されている場合、各（ｃk1，ｃk2）について、同様の対比操作を行う。
【０１０６】
ステップ（ｉ）：
選別される多型部位の総数が、上記のＮmin個に達するまで、前記ステップ（ｈ）に準じて、多型部位の数を増した組を考え、対比操作を進める。
【０１０７】
選別される多型部位の総数が、上記のＮmin個に達すると、その対比操作の結果を示す対比結果の行列中には、列要素（ｃk1，…，ｃkNmin）は、全て「−」に書き換えられるものが存在する。この（ｃk1，…，ｃkNmin）の組は、各「ハプロタイプ」相互の識別を行うことが可能な、Ｎmin個の多型部位で構成される「変異セット」として、選別される。このステップ（ｉ）において、選別されたＮmin個の多型部位で構成される複数個の「変異セット」は、下記の第２の工程における、「ディプロタイプ」：（Ｃi，Ｃj）の相違の識別を行うことが可能な「変異セット」の候補とされる。
【０１０８】
（第２の工程）：第１の工程で選別される、Ｎmin個の多型部位で構成される「変異セット」の候補から、「ディプロタイプ」：（Ｃi，Ｃj）の相違の識別も可能な「変異セット」を選択する、下記のステップ（ｊ）〜（ｌ）の操作を含んでなる工程；
ステップ（ｊ）：
Ｎmin個の多型部位で構成される「変異セット」の候補：（ｃk1，…，ｃkNmin）の一つについて、各「ハプロタイプ」Ｃiを表記する、Ｎmin次元の行ベクトルを定義する。各「ハプロタイプ」の組み合わせ（Ｃi，Ｃj）に由来する、「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトルを作成する。このＮmin次元の行ベクトルは、「ハプロタイプ」の組み合わせ（Ｃi，Ｃj）における、各要素の組み合わせに基づき、「野生型ホモ」、「変異型ホモ」、「ヘテロ」の３種の指標のいずれかにより、表記する。
【０１０９】
「野生型ホモ」は、その多型部位において、二つの「ハプロタイプ」（Ｃi，Ｃj）が、ともに、野生型アレル（Ｗ）である組み合わせ（Ｗ＋Ｗ）を意味する。
【０１１０】
「変異型ホモ」は、その多型部位において、二つの「ハプロタイプ」（Ｃi，Ｃj）が、ともに、変異型アレル（Ｍ）である組み合わせ（Ｍ＋Ｍ）を意味する。
【０１１１】
「ヘテロ」は、その多型部位において、二つの「ハプロタイプ」（Ｃi，Ｃj）が、一方は野生型アレル（Ｗ）、他方は、変異型アレル（Ｍ）である組み合わせ（Ｗ＋Ｍ）を意味する。
【０１１２】
「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトルは、各「ハプロタイプ」Ｃiを表記する、Ｎmin次元の行ベクトルを用いて、第１の演算処理を施すことで、その要素を算定する。
【０１１３】
ステップ（ｋ）：
各「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトル相互を比較し、同一のＮmin次元の行ベクトルで表記される、複数の「ディプロタイプ」：（Ｃi，Ｃj）が存在するか、否かを判定する。
【０１１４】
この各「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトル相互の比較は、二つのＮmin次元の行ベクトルに対して、第２の演算処理を施すことで、同一か否かを判定する。
【０１１５】
ステップ（ｌ）：
同一のＮmin次元の行ベクトルで表記される、複数の「ディプロタイプ」：（Ｃi，Ｃj）が存在する場合、その「変異セット」：（ｃk1，…，ｃkNmin）は、「ディプロタイプ」：（Ｃi，Ｃj）の識別も可能な「変異セット」ではないと判定する。
【０１１６】
同一のＮmin次元の行ベクトルで表記される、複数の「ディプロタイプ」：（Ｃi，Ｃj）が存在しない場合、その「変異セット」：（ｃk1，…，ｃkNmin）は、「ディプロタイプ」：（Ｃi，Ｃj）の識別も可能な「変異セット」であると判定する。
【０１１７】
他のＮmin個の多型部位で構成される「変異セット」の候補：（ｃk1，…，ｃkNmin）に関して、前記ステップ（ｉ）〜（ｌ）を繰り返す。
【０１１８】
最終的に、ステップ（ｌ）において、「ディプロタイプ」：（Ｃi，Ｃj）の識別も可能な「変異セット」であると判定された「変異セット」：（ｃk1，…，ｃkNmin）を、「最小数」の多型部位で構成される「変異セット」と選定する。
【０１１９】
その際、好ましくは、第１の演算手段は、
各「ハプロタイプ」Ｃiを表記する、Ｎmin次元の行ベクトルに基づき、「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトルを、下記のように定義する。
【０１２０】
すなわち、「ハプロタイプ」Ｃiを表記する、Ｎmin次元の行ベクトルと、「ハプロタイプ」Ｃjを表記する、Ｎmin次元の行ベクトルとの「和」によって、これを定義しＮmin次元の行ベクトルとする。
【０１２１】
また、第２の演算手段は、下記の手段で構成することができる。
【０１２２】
比較される二つの「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトル間の「差ベクトル」を求める。
【０１２３】
該「差ベクトル」が、零ベクトルである場合、同一のＮmin次元の行ベクトルであると判定する。
【０１２４】
該「差ベクトル」が、零ベクトルでない場合、同一のＮmin次元の行ベクトルではないと判定する。
【０１２５】
あるいは、第２の演算手段は、下記の手段で構成することもできる。
【０１２６】
比較される二つの「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトル間で、各要素を相互に比較する。
【０１２７】
全ての要素が同一である場合、同一のＮmin次元の行ベクトルであると判定する。
【０１２８】
同一でない要素が存在する場合、同一のＮmin次元の行ベクトルではないと判定する。
【０１２９】
上記態様では、本発明の方法においては、
前記遺伝子の多型部位における変異とは、一塩基多型、インサーション、デリーション、リピート配列のいずれかに起因する変異であることができる。
【０１３０】
また、本発明にかかる、プローブ・セットの設計方法の一態様では、下記の構成を有することが好ましい。
【０１３１】
この態様のプローブ・セットの設計方法は、
個体の有する相同染色体上、相同遺伝子座を占めるアレル（対立遺伝子）に関して、該個体が属する集団中において、複数種のアレル（対立遺伝子）の存在が確認されている際に適用する。
【０１３２】
その各々複数種のアレル（対立遺伝子）のＡ（アデニン）、Ｇ（グアニン）、Ｔ（チミン）、Ｃ（シトシン）の塩基で構成される塩基配列に関して、相互の塩基配列をアライメントする。この塩基配列のアライメントによって対比した際、見出される複数の変異箇所中、前記集団中における出現頻度の確率が、１％以上である複数の箇所を選択し、該遺伝子の多型部位と定義する。
【０１３３】
選択された該遺伝子の多型部位の総数をＮsnp0と表記する。
【０１３４】
前記複数種のアレル（対立遺伝子）相互の区別を行う指標として、該Ｎsnp0個の多型部位における塩基の種別を用い、各多型部位の塩基として、該集団中における出現頻度の最大の塩基を「野生型アレル（Ｗ）」、それ以外の塩基を「変異型アレル（Ｍ）」と定義する。
【０１３５】
前記Ｎsnp0個の多型部位における塩基の種別の２値化処理を行って、前記複数種のアレル（対立遺伝子）を識別する識別子として、次の行ベクトルを考える。これは、、前記Ｎsnp0個の多型部位：ｃk（ｋ＝１〜Ｎsnp0）における、２値化処理された塩基の種別で構成される、Ｎsnp0個の要素：ｃk（ｉ）（ｋ＝１〜Ｎsnp0）からなる、Ｎsnp0次元の行ベクトルと定義する。
【０１３６】
前記複数種のアレル（対立遺伝子）を、該Ｎsnp0次元の行ベクトルに基づき分類し、該Ｎsnp0次元の行ベクトルが同一となるアレル（対立遺伝子）の群を複数種定義し、この定義された複数種のアレル（対立遺伝子）の群を該遺伝子の「ハプロタイプ」と称する。
【０１３７】
定義される、該遺伝子の複数種の「ハプロタイプ」の総数をＮHと表記し、該ＮH種の「ハプロタイプ」をＣi：i＝１〜ＮHと表記する。
【０１３８】
各個体中に存在する相同遺伝子座を占める、二つのアレル（対立遺伝子）の組み合わせに起因して、その「ハプロタイプ」の組み合わせ（Ｃi，Ｃj）に由来する、前記Ｎsnp0個の多型部位における、塩基の種別；ＷとＭの組み合わせの表現型を考える。そして、前記Ｎsnp0個の多型部位における、塩基の種別；ＷとＭの組み合わせの表現型を、各個体の「ディプロタイプ」：（Ｃi，Ｃj）と定義する。
【０１３９】
本発明は、各個体の有する、二つのアレル（対立遺伝子）の組み合わせを分類する上で利用可能な、各「ハプロタイプ」相互の識別、同時に、その組み合わせである「ディプロタイプ」：（Ｃi，Ｃj）の相違の識別を行うことを目標とする。更に、「ハプロタイプ」相互の識別と、「ディプロタイプ」：（Ｃi，Ｃj）の相違の識別が可能な、前記Ｎsnp0個の多型部位：ｃk（ｋ＝１〜Ｎsnp0）中から選択される、「最小数」の多型部位で構成される「変異セット」を選定する。
【０１４０】
その上で、該「最小数」の多型部位で構成される「変異セット」における、各多型部位における塩基の種別を特定するために利用されるプローブ・セットを設計する方法である。
【０１４１】
本発明のプローブ・セットの設計方法において、該プローブ・セットの設計工程は、下記のステップＡとステップＢを含む。
【０１４２】
ステップＡ：
上述する本発明の変異セットの選別方法を適用して、「最小数」の多型部位で構成される「変異セット」を選定する。
【０１４３】
ステップＢ：
選定された「最小数」の多型部位で構成される「変異セット」の各多型部位について、該多型部位をその塩基配列中に含むように、プローブの塩基配列を選択する。
【０１４４】
また、本発明の一態様にかかる、プライマー・セットの設計方法は、下記の構成を有することが好ましい。
【０１４５】
この態様のプライマー・セットの設計方法は、
個体の有する相同染色体上、相同遺伝子座を占めるアレル（対立遺伝子）に関して、該個体が属する集団中において、複数種のアレル（対立遺伝子）の存在が確認されている際に適用する。
【０１４６】
その各々複数種のアレル（対立遺伝子）のＡ（アデニン）、Ｇ（グアニン）、Ｔ（チミン）、Ｃ（シトシン）の塩基で構成される塩基配列に関して、相互の塩基配列をアライメントする。この塩基配列のアライメントによって対比した際、見出される複数の変異箇所中、前記集団中における出現頻度の確率が、１％以上である複数の箇所を選択し、該遺伝子の多型部位と定義する。
【０１４７】
選択された該遺伝子の多型部位の総数をＮsnp0と表記する。
【０１４８】
前記複数種のアレル（対立遺伝子）相互の区別を行う指標として、該Ｎsnp0個の多型部位における塩基の種別を用い、各多型部位の塩基として、該集団中における出現頻度の最大の塩基を「野生型アレル（Ｗ）」、それ以外の塩基を「変異型アレル（Ｍ）」と定義する。
【０１４９】
前記Ｎsnp0個の多型部位における塩基の種別の２値化処理を行って、前記複数種のアレル（対立遺伝子）を識別する識別子として、次の行ベクトルを考える。すなわち、前記Ｎsnp0個の多型部位：ｃk（ｋ＝１〜Ｎsnp0）における、２値化処理された塩基の種別で構成される、Ｎsnp0個の要素：ｃk（ｉ）（ｋ＝１〜Ｎsnp0）からなる、Ｎsnp0次元の行ベクトルを定義する。
【０１５０】
前記複数種のアレル（対立遺伝子）を、該Ｎsnp0次元の行ベクトルに基づき分類し、該Ｎsnp0次元の行ベクトルが同一となる複数種のアレル（対立遺伝子）の群を定義し、この定義された複数種のアレル（対立遺伝子）の群を該遺伝子の「ハプロタイプ」と称する。
【０１５１】
定義される、該遺伝子の複数種の「ハプロタイプ」の総数をＮHと表記し、該ＮH種の「ハプロタイプ」をＣi：i＝１〜ＮHと表記する。
【０１５２】
各個体中に存在する相同遺伝子座を占める、二つのアレル（対立遺伝子）の組み合わせに起因して、その「ハプロタイプ」の組み合わせ（Ｃi，Ｃj）に由来する、前記Ｎsnp0個の多型部位における、塩基の種別；ＷとＭの組み合わせの表現型を考える。そして、前記Ｎsnp0個の多型部位における、塩基の種別；ＷとＭの組み合わせの表現型を、各個体の「ディプロタイプ」：（Ｃi，Ｃj）と定義する。
【０１５３】
各個体の有する、二つのアレル（対立遺伝子）の組み合わせを分類する上で利用可能な、各「ハプロタイプ」相互の識別、同時に、その組み合わせである「ディプロタイプ」：（Ｃi，Ｃj）の相違の識別を行うことを目標とする。「ハプロタイプ」相互の識別と、「ディプロタイプ」：（Ｃi，Ｃj）の相違の識別が可能な、前記Ｎsnp0個の多型部位：ｃk（ｋ＝１〜Ｎsnp0）中から選択される、「最小数」の多型部位で構成される「変異セット」を選定する。
【０１５４】
その上で、該「最小数」の多型部位で構成される「変異セット」における、各多型部位における塩基の種別を特定するために利用されるプライマー・セットを設計する方法である。
【０１５５】
この態様のプライマー・セットの設計方法において、該プライマー・セットの設計工程は、下記のステップＡとステップＢを含む。
【０１５６】
ステップＡ：
上述する本発明の一態様の変異セットの選別方法を適用して、「最小数」の多型部位で構成される「変異セット」を選定する。
【０１５７】
ステップＢ：
選定された「最小数」の多型部位で構成される「変異セット」の各多型部位について、該多型部位の塩基配列の解析用プライマーに適する、プライマーの塩基配列を選択する。
【０１５８】
なお、本発明の一態様の変異セットの選別方法は、一連の数値演算処理、ならびに、該数値演算処理結果を用いる論理演算処理を行うことで、「最小数」の多型部位で構成される「変異セット」を特定する形態として実施できる。これらの演算処理を含む操作は、コンピュータシステム上で実行するため、専用のソフトウエアとして、記録媒体上に記録する形態にすることが望ましい。
【０１５９】
従って、この態様は、前記の「専用のソフトウエア」を収納する、ソフトウエア収納記録媒体の発明をも提供する。
【０１６０】
すなわち、この態様のソフトウエア収納記録媒体は、下記のような構成上の特徴を有する。
【０１６１】
・上述する本発明の一態様の変異セットの選別方法を適用して、「最小数」の多型部位で構成される「変異セット」を選定するためのものである。
【０１６２】
・その方法における、一連の数値演算処理、ならびに、該数値演算処理結果を用いる論理演算処理をコンピュータシステム上で実行するためのプログラムを、主要な要素としている。
【０１６３】
・そのプログラムを、該コンピュータシステムにおいて実行可能なプログラム言語により表記したソフトウエアを記録媒体上に電子的に収納してなるものである。
【０１６４】
また、本発明の一態様の変異セットの選別方法は、上述するように、「専用のソフトウエア」を用いて、コンピュータシステム上において、その設計操作を実施することが一般的である。従って、本発明は、本発明の変異セットの選別方法を適用して、「最小数」の多型部位で構成される「変異セット」を選定する操作を実施する用途に適するコンピュータシステムの発明をも提供する。
【０１６５】
すなわち、本発明の一態様にかかる「最小数」の多型部位で構成される「変異セット」の選定を実行するためのコンピュータシステムは、下記の構成上の特徴を有する。
【０１６６】
この態様の「変異セット」選定用コンピュータシステムは、
個体の有する相同染色体上、相同遺伝子座を占めるアレル（対立遺伝子）に関して、該個体が属する集団中において、複数種のアレル（対立遺伝子）の存在が確認されている際に利用される。
【０１６７】
その各々複数種のアレル（対立遺伝子）のＡ（アデニン）、Ｇ（グアニン）、Ｔ（チミン）、Ｃ（シトシン）の塩基で構成される塩基配列に関して、相互の塩基配列をアライメントする。該塩基配列のアライメントによって対比した際、見出される複数の変異箇所中、前記集団中における出現頻度の確率が、１％以上である複数の箇所を選択し、該遺伝子の多型部位と定義する。
【０１６８】
選択された該遺伝子の多型部位の総数をＮsnp0と表記する。
【０１６９】
前記複数種のアレル（対立遺伝子）相互の区別を行う指標として、該Ｎsnp0個の多型部位における塩基の種別を用い、各多型部位の塩基として、該集団中における出現頻度の最大の塩基を「野生型アレル（Ｗ）」、それ以外の塩基を「変異型アレル（Ｍ）」と定義する。
【０１７０】
前記Ｎsnp0個の多型部位における塩基の種別の２値化処理を行って、前記複数種のアレル（対立遺伝子）を識別する識別子として、次の行ベクトルを考える。すなわち、前記Ｎsnp0個の多型部位：ｃk（ｋ＝１〜Ｎsnp0）における、２値化処理された塩基の種別で構成される、Ｎsnp0個の要素：ｃk（ｉ）（ｋ＝１〜Ｎsnp0）からなる、Ｎsnp0次元の行ベクトルを定義する。
【０１７１】
前記複数種のアレル（対立遺伝子）を、該Ｎsnp0次元の行ベクトルに基づき分類し、該Ｎsnp0次元の行ベクトルが同一となる複数種のアレル（対立遺伝子）の群を定義し、この定義された複数種のアレル（対立遺伝子）の群を該遺伝子の「ハプロタイプ」と称する。
【０１７２】
定義される、該遺伝子の複数種の「ハプロタイプ」の総数をＮHと表記し、該ＮH種の「ハプロタイプ」をＣi：i＝１〜ＮHと表記する。
【０１７３】
各個体中に存在する相同遺伝子座を占める、二つのアレル（対立遺伝子）の組み合わせに起因して、その「ハプロタイプ」の組み合わせ（Ｃi，Ｃj）に由来する、前記Ｎsnp0個の多型部位における、塩基の種別；ＷとＭの組み合わせの表現型を考える。そして、前記Ｎsnp0個の多型部位における、塩基の種別；ＷとＭの組み合わせの表現型を、各個体の「ディプロタイプ」：（Ｃi，Ｃj）と定義する。
【０１７４】
本発明は、各個体の有する、二つのアレル（対立遺伝子）の組み合わせを分類する上で利用可能な、各「ハプロタイプ」相互の識別、同時に、その組み合わせである「ディプロタイプ」：（Ｃi，Ｃj）の相違の識別を行うことを目標とし、各「ハプロタイプ」相互の識別と、「ディプロタイプ」：（Ｃi，Ｃj）の相違の識別が可能な、前記Ｎsnp0個の多型部位：ｃk（ｋ＝１〜Ｎsnp0）中から選択される、「最小数」の多型部位で構成される「変異セット」を選定する。この「最小数」の多型部位で構成される「変異セット」を選定する作業を実施するための、コンピュータシステムである。
【０１７５】
その際、該システムは、
まず、前記「最小数」の多型部位で構成される「変異セット」を構成する多型部位の個数をＮと表記する際、条件：（ＮH−１）≦２N-1を満たすＮの最小値：Ｎminを算定する演算処理を実施する。
【０１７６】
次いで、前記Ｎsnp0個の多型部位：ｃk（ｋ＝１〜Ｎsnp0）中から、該Ｎmin個以上の多型部位を選択して、構成される「変異セット」の集合の中から、「最小数」の多型部位で構成される「変異セット」を選定する演算処理を実施するシステムである。
【０１７７】
該「最小数」の多型部位で構成される「変異セット」を選定する演算処理工程は、下記のステップ（ａ）〜（ｉ）の演算操作を含んでなる第１の工程と、ステップ（ｊ）〜（ｌ）の演算操作を含んでなる第２の工程を含む。
【０１７８】
（第１の工程）：各「ハプロタイプ」相互の識別を行うことが可能な、Ｎmin個以上の多型部位で構成される「変異セット」を、「最小数」の多型部位で構成される「変異セット」の候補として抽出する、下記のステップ（ａ）〜（ｉ）の操作を含んでなる工程；
ステップ（ａ）：
前記集団中において、その存在が確認されている複数種のアレル（対立遺伝子）の塩基配列の情報を、該遺伝子の多型情報を収納するデータベースから抽出する。抽出された塩基配列の情報に基づき、前記Ｎsnp0個の多型部位：ｃk（ｋ＝１〜Ｎsnp0）における、２値化処理された塩基の種別で構成される、Ｎsnp0個の要素：ｃk（ｉ）（ｋ＝１〜Ｎsnp0）からなる、Ｎsnp0次元の行ベクトルを作成する。次いで、作成されるＮsnp0次元の行ベクトルの集合中から、互いに、相違するＮsnp0次元の行ベクトルを示す、該ＮH種の「ハプロタイプ」：Ｃi：i＝１〜ＮHを表記する、Ｎsnp0次元の行ベクトルＮH種を選別する。
【０１７９】
ステップ（ｂ）：
前記Ｎsnp0次元の行ベクトルＮH種を用いて、ＮH行×Ｎsnp0列の行列を形成し。その行列の各列より、Ｎsnp0個のＮH次元の列ベクトルを作成する。各列ベクトル間の対比を行い、相互に同一の列ベクトルの群、ならびに、相互に逆ベクトルの関係にある列ベクトルの群を選択する。選択される列ベクトルの群中の中から、一つの列ベクトルを選択し、該群に含まれる他の列ベクトルを除去し、前記ＮH行×Ｎsnp0列の行列から、各列ベクトルが相互に独立している列ベクトルにより構成されている部分行列を作成する。該部分行列の列の総数をＮsnpと表記し、該ＮH行×Ｎsnp列の行列における、各列を与えている多型部位：ｃk（ｋ＝１〜Ｎsnp）を、該Ｎmin個以上の多型部位を選択して、構成される「変異セット」の可能な要素として選択する。
【０１８０】
ステップ（ｃ）：
ＮH行×Ｎsnp列の行列において、各「ハプロタイプ」：Ｃiを表記する行ベクトル（ｃk（ｉ））を、参照行として選択する。他の「ハプロタイプ」：Ｃjを表記する行ベクトル（ｃk（ｊ））に対して、前記参照行の行ベクトル（ｃk（ｉ））と、各要素ｃk（ｊ）とｃk（ｉ）の比較を行う。比較する要素ｃk（ｊ）とｃk（ｉ）が、異なる場合、両者間の区別が可能として、ｃk（ｊ）を「−」の符号に置換する。この操作を、前記参照行の「ハプロタイプ」：Ｃiを表記する行ベクトル（ｃk（ｉ））以外の、（ＮH−１）個の行ベクトルに施す。この対比操作の結果を、ＮH行×Ｎsnp列の行列として、表記し、「ハプロタイプ」：Ｃiを表記する行ベクトル（ｃk（ｉ））を参照行とする、第１の対比結果の行列とする。
【０１８１】
ステップ（ｄ）：
引き続き、各「ハプロタイプ」：Ｃiを表記する行ベクトル（ｃk（ｉ））を参照行として、同じ対比操作を行う。最終的に、各「ハプロタイプ」：Ｃiを表記する行ベクトル（ｃk（ｉ））を参照行とする、第１の対比結果の行列を、計ＮH個作成する。
【０１８２】
ステップ（ｅ）：
作成される計ＮH個の第１の対比結果の行列を参照し、前記の対比操作が施された後の行、ＮH×（ＮH−１）のうち、「−」の符号へ置き換えが、最小数である行を選択する。この選択された行において、「−」の符号へ置き換えがなされている要素（列）を特定する。この特定される要素（列）を与えている多型部位を特定する。この特定される多型部位を、「最小数」の多型部位で構成される「変異セット」を構成する多型部位の「第１の候補」：ｃk1とする。
【０１８３】
ステップ（ｆ）：
次いで、「第１の候補」：ｃk1と、他の多型部位：ｃk：ｃk≠ｃk1との組み合わせ：（ｃk1，ｃk）として、計（Ｎsnp−１）組を考える。
【０１８４】
ＮH行×Ｎsnp列の行列において、各「ハプロタイプ」：Ｃiを表記する行ベクトルを、参照行として選択する。他の「ハプロタイプ」：Ｃjを表記する行ベクトルに対して、前記参照行の行ベクトルと、各要素対（ｃk1（ｉ），ｃk（ｉ））と（ｃk1（ｊ），ｃk（ｊ））の比較を行う。比較する要素対（ｃk1（ｉ），ｃk（ｉ））と（ｃk1（ｊ），ｃk（ｊ））が、異なる場合、両者間の区別が可能として、（ｃk1（ｊ），ｃk（ｊ））を「−」の符号に置換する。この操作を、前記参照行の「ハプロタイプ」：Ｃiを表記する行ベクトル以外の、（ＮH−１）個の行ベクトルに施す。この対比操作の結果を、列要素として、計（Ｎsnp−１）個の（ｃk1，ｃk）を有する、ＮH行×（Ｎsnp−１）列の行列として、表記する。このＮH行×（Ｎsnp−１）列の行列を、「ハプロタイプ」：Ｃiを表記する行ベクトルを参照行とする、第２の対比結果の行列とする。
【０１８５】
引き続き、各「ハプロタイプ」：Ｃiを表記する行ベクトルを参照行として、同じ対比操作を行う。最終的に、各「ハプロタイプ」：Ｃiを表記する行ベクトルを参照行とする、第２の対比結果の行列を、計ＮH個作成する。
【０１８６】
ステップ（ｇ）：
作成される計ＮH個の第２の対比結果の行列を参照し、前記の対比操作が施された後の行ＮH×（ＮH−１）のうち、「−」の符号へ置き換えが、最小数である行を選択する。この選択された行において、「−」の符号へ置き換えがなされている要素対（ｃk1，ｃk）を特定する。この特定される要素対（ｃk1，ｃk）を与えている、他の多型部位：ｃkを特定する。この特定される、他の多型部位：ｃkを、「最小数」の多型部位で構成される「変異セット」を構成する多型部位の「第２の候補」：ｃk2とする。
【０１８７】
ステップ（ｈ）：
次いで、「第１の候補」：ｃk1、「第２の候補」：ｃk2と他の多型部位：ｃkとの組み合わせ：（ｃk1，ｃk2，ｃk）として、計（Ｎsnp−２）組を考える。
【０１８８】
前記ステップ（ｆ）の対比操作に準じて、対比操作を進め、その結果を、列要素として、計（Ｎsnp−２）個の（ｃk1，ｃk2，ｃk）を有する、ＮH行×（Ｎsnp−２）列の行列として、表記する。このＮH行×（Ｎsnp−２）列の行列を、「ハプロタイプ」：Ｃiを表記する行ベクトルを参照行とする、第３の対比結果の行列とする。
【０１８９】
引き続き、各「ハプロタイプ」：Ｃiを表記する行ベクトルを参照行として、同じ対比操作を行う。最終的に、各「ハプロタイプ」：Ｃiを表記する行ベクトルを参照行とする、第３の対比結果の行列を、計ＮH個作成する。
【０１９０】
作成される計ＮH個の第３の対比結果の行列を参照し、前記の対比操作が施された後の行、ＮH×（ＮH−１）のうち、「−」の符号へ置き換えが、最小数である行を選択する。この選択された行において、「−」の符号へ置き換えがなされている要素組（ｃk1，ｃk2，ｃk）を特定する。この特定される要素組（ｃk1，ｃk2，ｃk）を与えている、他の多型部位：ｃkを特定する。この特定される、他の多型部位：ｃkを、「最小数」の多型部位で構成される「変異セット」を構成する多型部位の「第３の候補」：ｃk3とする。
【０１９１】
なお、前記「第２の候補」：ｃk2が複数選別されている場合、各（ｃk1，ｃk2）について、同様の対比操作を行う。
【０１９２】
ステップ（ｉ）：
選別される多型部位の総数が、上記のＮmin個に達するまで、前記ステップ（ｈ）に準じて、多型部位の数を増した組を考え、対比操作を進める。
【０１９３】
選別される多型部位の総数が、上記のＮmin個に達すると、その対比操作の結果を示す対比結果の行列中には、列要素（ｃk1，…，ｃkNmin）は、全て「−」に書き換えられるものが存在する。この（ｃk1，…，ｃkNmin）の組は、各「ハプロタイプ」相互の識別を行うことが可能な、Ｎmin個の多型部位で構成される「変異セット」として、選別される。このステップ（ｉ）において、選別されたＮmin個の多型部位で構成される複数個の「変異セット」は、下記の第２の工程における、「ディプロタイプ」：（Ｃi，Ｃj）の相違の識別を行うことが可能な「変異セット」の候補とされる。
【０１９４】
（第２の工程）：第１の工程で選別される、Ｎmin個の多型部位で構成される「変異セット」の候補から、「ディプロタイプ」：（Ｃi，Ｃj）の相違の識別も可能な「変異セット」を選択する、下記のステップ（ｊ）〜（ｌ）の操作を含んでなる工程；
ステップ（ｊ）：
Ｎmin個の多型部位で構成される「変異セット」の候補：（ｃk1，…，ｃkNmin）の一つについて、各「ハプロタイプ」Ｃiを表記する、Ｎmin次元の行ベクトルを定義する。各「ハプロタイプ」の組み合わせ（Ｃi，Ｃj）に由来する、「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトルを作成する。このＮmin次元の行ベクトルは、「ハプロタイプ」の組み合わせ（Ｃi，Ｃj）における、各要素の組み合わせに基づき、「野生型ホモ」、「変異型ホモ」、「ヘテロ」の３種の指標のいずれかにより、表記する。
【０１９５】
「野生型ホモ」は、その多型部位において、二つの「ハプロタイプ」（Ｃi，Ｃj）が、ともに、野生型アレル（Ｗ）である組み合わせ（Ｗ＋Ｗ）を意味する。
【０１９６】
「変異型ホモ」は、その多型部位において、二つの「ハプロタイプ」（Ｃi，Ｃj）が、ともに、変異型アレル（Ｍ）である組み合わせ（Ｍ＋Ｍ）を意味する。
【０１９７】
「ヘテロ」は、その多型部位において、二つの「ハプロタイプ」（Ｃi，Ｃj）が、一方は野生型アレル（Ｗ）、他方は、変異型アレル（Ｍ）である組み合わせ（Ｗ＋Ｍ）を意味する。
【０１９８】
「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトルは、各「ハプロタイプ」Ｃiを表記する、Ｎmin次元の行ベクトルを用いて、第１の演算処理を施すことで、その要素を算定する。
【０１９９】
ステップ（ｋ）：
各「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトル相互を比較し、同一のＮmin次元の行ベクトルで表記される、複数の「ディプロタイプ」：（Ｃi，Ｃj）が存在するか、否かを判定する。
【０２００】
この各「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトル相互の比較は、二つのＮmin次元の行ベクトルに対して、第２の演算処理を施すことで、同一か否かを判定する。
【０２０１】
ステップ（ｌ）：
同一のＮmin次元の行ベクトルで表記される、複数の「ディプロタイプ」：（Ｃi，Ｃj）が存在する場合、その「変異セット」：（ｃk1，…，ｃkNmin）は、「ディプロタイプ」：（Ｃi，Ｃj）の識別も可能な「変異セット」ではないと判定する。
【０２０２】
同一のＮmin次元の行ベクトルで表記される、複数の「ディプロタイプ」：（Ｃi，Ｃj）が存在しない場合、その「変異セット」：（ｃk1，…，ｃkNmin）は、「ディプロタイプ」：（Ｃi，Ｃj）の識別も可能な「変異セット」であると判定する。
【０２０３】
他のＮmin個の多型部位で構成される「変異セット」の候補：（ｃk1，…，ｃkNmin）に関して、前記ステップ（ｉ）〜（ｌ）を繰り返す。
【０２０４】
最終的に、ステップ（ｌ）において、「ディプロタイプ」：（Ｃi，Ｃj）の識別も可能な「変異セット」であると判定された「変異セット」：（ｃk1，…，ｃkNmin）を、「最小数」の多型部位で構成される「変異セット」と選定する。
【０２０５】
その際、好ましくは、第１の演算手段は、下記の操作で構成する。
【０２０６】
各「ハプロタイプ」Ｃiを表記する、Ｎmin次元の行ベクトルに基づき、「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトルを、下記のように定義する。
【０２０７】
すなわち、「ハプロタイプ」Ｃiを表記する、Ｎmin次元の行ベクトルと、「ハプロタイプ」Ｃjを表記する、Ｎmin次元の行ベクトルとの「和」によって、定義されるＮmin次元の行ベクトルとする。
【０２０８】
また、第２の演算手段は、下記の操作で構成することができる。
【０２０９】
比較される二つの「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトル間の「差ベクトル」を求める。
【０２１０】
該「差ベクトル」が、零ベクトルである場合、同一のＮmin次元の行ベクトルであると判定する。
【０２１１】
該「差ベクトル」が、零ベクトルでない場合、同一のＮmin次元の行ベクトルではないと判定する。
【０２１２】
あるいは、第２の演算手段は、下記の操作で構成することもできる。
【０２１３】
比較される二つの「ディプロタイプ」：（Ｃi，Ｃj）を表記する、Ｎmin次元の行ベクトル間で、各要素を相互に比較する。
【０２１４】
全ての要素が同一である場合、同一のＮmin次元の行ベクトルであると判定する。
【０２１５】
同一でない要素が存在する場合、同一のＮmin次元の行ベクトルではないと判定する
本発明のシステムにおいては、
前記遺伝子の多型部位における変異は、一塩基多型、インサーション、デリーション、リピート配列のいずれかに起因する変異であることができる。
【０２１６】
一つの態様においては、本発明にかかるコンピュータシステムは、そのハードウエアの構成は、複数の機器例えば、ホストコンピュータ、インタフェース機器、ネットワーク機器のような統合システムであってもよい。あるいは、前記のシステムの演算機能、入出力機能などを具える、一体型のコンピュータであってもよい。
【０２１７】
一方、本発明にかかるソフトウエア収納記録媒体では、記憶媒体として、例えば、フロッピーディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、ＤＶＤ、磁気テープ、不揮発性のメモリカード、ＲＯＭなどが利用できる。
【０２１８】
なお、実際のコンピュータシステムにおける、一連の処理・操作の実行は、前記コンピュータシステムの記録媒体中に収納されるソフトウエアの読み込みによって行われる。すなわち、ソフトウエアの実装とは、実際に、ソフトウエアの読み込みがなされた状態、あるいは、ソフトウエアの読み込みが可能な状態に、ソフトウエア収納記録媒体をシステムに装着された状態を意味する。ソフトウエアの主要な要素のプログラムは、一連の処理ステップの進行を規定するものであり、その際、利用される各種の演算処理自体は、システム側に設けられていてもよい。勿論、該プログラム中のサブ・セットとして、一体化されていてもよく、プログラム自体に、演算処理を内在する形態であってもよい。
【０２１９】
なお、プログラムにより実行される演算処理は、下記の実施形態に記載する各ステップの作業内容を、利用するシステムで実行可能な演算処理に適する形態に定式化したものとする。
【０２２０】
以下に、具体例に基づき、本発明をより詳しく説明する。なお、下記の具体例は、本発明の最良の実施形態の一例であるが、本発明は、かかる態様に限定されるものではない。
【０２２１】
〔第１の実施形態〕
図１〜図８を参照して、本発明にかかる第１の実施形態を、以下に説明する。
【０２２２】
８つの多型部位中から、Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６（Ｃi：ｉ＝１〜６）を識別でき、同時に、ホモ型（Ｃi、Ｃi）、ヘテロ型（Ｃi、Ｃj）の組み合わせの「ディプロタイプ」相互の区別も可能な「最小数の変異セット」を抽出する例を説明する。
【０２２３】
（ステップａ−ｂ）
この実施形態では、各「ハプロタイプ」：Ｃiは、８つの多型部位における、変異の種別：ｃk（ｉ）（ｋ＝１〜８）を用いて、Ｃi＝（ｃ1（ｉ），…，ｃk（ｉ），…，ｃ8（ｉ））と表記する。また、全体の集団における出現頻度確率を評価した場合、各多型部位ｃk（ｋ＝１〜８）に関して、最も出現頻度の高い「塩基種類」を、野生型アレル（Ｗ）の塩基と定義する。それ以外の塩基を、変異型アレル（Ｍ）と定義する。「ハプロタイプ」：Ｃiの多型部位ｃk（ｋ＝１〜８）が、野生型アレル（Ｗ）の場合、ｃk（ｉ）＝−１と、変異型アレル（Ｍ）の場合、ｃk（ｉ）＝＋１と表記する。
【０２２４】
Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６（Ｃi：ｉ＝１〜６）に関して、確認されている多型部位の総数はＮSNP、各Ｃiは、（ｃ1（ｉ），…，ｃk（ｉ），…，ｃNSNP（ｉ））となっている。その際、各多型部位ｃk（ｋ＝１〜ＮSNP）における列（ｃk（１），…，ｃk（ｉ），…，ｃk（６））を比較し、同じ列ベクトルとなるものがある場合、その複数の多型部位ｃkの内一つのみを選択する。すなわち、確認されている多型部位において、同じ列ベクトルを示す複数の多型部位に関しては、それぞれ、その一つを予め選別している。
【０２２５】
加えて、各多型部位ｃk（ｋ＝１〜ＮSNP）における列（ｃk（１），…，ｃk（ｉ），…，ｃk（６））を比較し、相互に、要素が反転した列ベクトルとなるものがある場合も、その複数の多型部位ｃkの内一つのみを選択している。この二つ列は、互いに相違はしているが、各要素を−１⇔＋１と変換すると、同じものとなるため、互いに、「独立」するものとは見做せない。従って、その一方のみを、予め選別している。
【０２２６】
以上の操作によって、「等価」と見做せる列ベクトルのうち、一つを残すことによって、検討の対象となる多型部位ｃkを限定している。ここでは、選別後、ｋ＝１〜８に限定されている。
【０２２７】
図１に、Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６を、それぞれＣi＝（ｃ1（ｉ），…，ｃk（ｉ），…，ｃ8（ｉ））と表記し、列（Ｃ1，…，Ｃ6）を、６行×８列の行例〔ｃk（ｉ）〕として、示している。その際、各列（ｃk（１），…，ｃk（ｉ），…，ｃk（６））は、互いに相違しており、また、各行（ｃ1（ｉ），…，ｃk（ｉ），…，ｃ8（ｉ））も、互いに相違している。
【０２２８】
（ステップｃ−ｉ）
Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６を示す行ベクトルＣi＝（ｃ1（ｉ），…，ｃk（ｉ），…，ｃ8（ｉ））：ｉ＝１〜６は、８次元空間において、互いに区別可能なベクトルと見做せる。
【０２２９】
８つの多型部位ｃk（ｋ＝１〜８）のうち、Ｎ個の多型部位ｃkを選択して、Ｎ次元空間を想定する。Ｎ個の多型部位ｃkの選択は、8ＣN通りあり、それぞれの組み合わせを、Ｎｊ：ｊ＝１〜8ＣNと表記する。
【０２３０】
組み合わせ：Ｎｊにおける、Ｎ次元空間上において、Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６を表記する行ベクトルＣiNjを考える。この行ベクトルＣiNjは、前記組み合わせ：Ｎｊ：（ｃkNj：k＝１〜Ｎ）に従って、Ｎ次元の行ベクトル（ｃ1Nj（ｉ），…，ｃkNj（ｉ），…，ｃNNj（ｉ））と表記される。一方、Ｎ次元の行ベクトルにおいて、各要素を、＋１、−１のいずれかとすると、異なる行ベクトルは、２N通りしかない。
【０２３１】
この２N通りの行ベクトル中に、Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６に対応して、互いに異なる６つの行ベクトルが含まれるためには、少なくとも、６＜２Nの条件を満足する必要がある。換言すると、この事例では、Ｎは、３以上である必要がある。
【０２３２】
なお、２N通りの行ベクトルは、その要素は、＋１、−１のいずれかであり、全ての要素について、−１⇔＋１と変換すると、同じものとなる、互いに逆ベクトルの関係にある、ベクトル組、２N-1組で構成されている。
【０２３３】
一方、ヘテロ型（Ｃi、Ｃj）組み合わせの「ディプロタイプ」を、行ベクトル（｛ｃ1Nj（ｉ）＋ｃ1Nj（ｊ）｝，…，｛ｃkNj（ｉ）＋ｃkNj（ｊ）｝，…，｛ｃNNj（ｉ）＋ｃNNj（ｊ）｝）と表記する。その際、ヘテロ型（Ｃi、Ｃj）組み合わせが、互いに逆ベクトルの関係にある、ベクトル組であると、｛ｃkNj（ｉ）＋ｃkNj（ｊ）｝＝０となる。従って、少なくとも、この互いに逆ベクトルの関係にある、ベクトル組に相当するヘテロ型（Ｃi、Ｃj）の組み合わせの「ディプロタイプ」は、区別不能となる。
【０２３４】
換言するならば、ヘテロ型（Ｃi、Ｃj）の組み合わせの「ディプロタイプ」の全てを区別可能とする必要条件は、下記のようになる。各「ハプロタイプ」：Ｃiが、互いに逆ベクトルの関係にある、ベクトル組、２N-1組から、独立に選択されたものとなっていることである。この事例では、互いに逆ベクトルの関係にある、ベクトル組、２N-1組から、６つの組を選択できることがその条件となる。従って、少なくとも、６＜２N-1の条件を満足する必要がある。換言すると、この事例では、Ｎは、４以上である必要がある。
【０２３５】
実際に、６つの「ハプロタイプ」：Ｃi：ｉ＝１〜６について、「最小数」の多型部位を含む変異セット、すなわち、３個の多型部位からなるセット（ｃkNj：ｉ＝１〜３）を考える。その場合、互いに逆ベクトルの関係にある、ベクトル組は、２2組＝４組しかない。従って、少なくとも、４つの「ハプロタイプ」：Ｃiは、互いに逆ベクトルの関係にある、ベクトル組となっている。すなわち、ヘテロ型（Ｃi、Ｃj）組み合わせの「ディプロタイプ」のうち、二つは、（｛ｃ1Nj（ｉ）＋ｃ1Nj（ｊ）｝，｛ｃ2Nj（ｉ）＋ｃ2Nj（ｊ）｝，｛ｃ3Nj（ｉ）＋ｃ3Nj（ｊ）｝）＝（０，０，０）となっている。
【０２３６】
仮に、４つの「ハプロタイプ」：Ｃi：ｉ＝１〜３，６について、「最小数」の多型部位を含む変異セット、すなわち、３個の多型部位からなるセット（ｃkNj：ｉ＝１〜３，６）を考える。その場合、互いに逆ベクトルの関係にある、ベクトル組は、２2組＝４組であるから、４つの「ハプロタイプ」：Ｃiは、互いに逆ベクトルの関係にないように選択することが可能となる。一般に、ＮH種の「ハプロタイプ」：Ｃiを判別でき、同時に、｛1/2ＮH×（ＮH−１）｝種のヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の区別も可能な「最小数」の多型部位を含む変異セットの数：Ｎは、ＮH≦２N-1の条件を満たす必要がある。単に、ＮH種の「ハプロタイプ」：Ｃiを判別できる、「最小数」の多型部位を含む変異セットの数：Ｎminは、場合によっては、ＮH≦２Nminの条件を満たすのみで十分である。
【０２３７】
以下に、単に６種の「ハプロタイプ」：Ｃiを判別できる、「最小数」の多型部位を含む変異セット、すなわち、多型部位の数：Ｎmin＝３の変異セットを選別する工程を先ず説明する。さらに、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の区別も可能な「最小数」の多型部位を含む変異セットの候補となる、多型部位の数：Ｎ＝４の変異セットを選別する工程を説明する。
【０２３８】
６種の「ハプロタイプ」：Ｃi：ｉ＝１〜６を識別する上では、８個の多型部位ｃk（ｋ＝１〜８）全てを判定要素とするのは、確かに冗長である。そのため、有効な「最小数」：Ｎminの多型部位を抽出し、最小数の変異セットの組を選別する。ここでは、選別方法として、特開２００５−１９０４２７号公報に開示される「配列を同定するための変異要素のセット抽出方法」を応用する。
【０２３９】
（ステップｃ，ｄ）
６行×８列の行例〔ｃk（ｉ）〕に示す、Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６のうち、「参照行」を一つだけ固定して、列方向に「参照行」の要素と異なる場合には識別可の記号として−に置換する。「参照行」と同値である場合には、そのままの値を保存する。図２は、「参照行」として、Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６を、それぞれ選び、前記の比較演算処理を行った結果を、６通りの「６行×８列の行例」として、図示したものである。各「６行×８列の行例」中、黒枠で囲まれた行が「参照行」である。なお、この比較演算処理を施す６行×８列の行例〔ｃk（ｉ）〕を、最上部に示す。この６行×８列の行例〔ｃk（ｉ）〕は、図１と同じ行列である。
【０２４０】
（ステップｅ）
比較演算処理を行った結果のうち、Ｈａｐｌｏｔｙｐｅ３を「参照行」とした結果に着目する。Ｈａｐｌｏｔｙｐｅ５の行に対する比較演算処理結果では、多型部位：ｃ5（多型５）のみが、「識別可（−）」となっており、それ以外の７つの要素：ｃk（５）：ｋ＝１〜４、６〜８は、「消去」されていない。すなわち、少なくとも、Ｈａｐｌｏｔｙｐｅ３とＨａｐｌｏｔｙｐｅ５との間の識別を行う上では、多型部位：ｃ5（多型５）の判定が不可欠であることが判る。勿論、Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６の識別を行う際には、Ｈａｐｌｏｔｙｐｅ３とＨａｐｌｏｔｙｐｅ５との間の識別は必須である。換言するならば、Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６の識別を行う上でも、多型部位：ｃ5（多型５）の判定が不可欠であることになる。従って、「最小数」の多型部位を含む変異セットにおける、多型部位の「第１の候補」として、多型部位：ｃ5（多型５）が選択される。
【０２４１】
（ステップｆ）
次に、多型部位：ｃ5（多型５）と、他の一つの多型部位を組み合わせ、（ｃ5，ｃk）：ｋ＝１〜４、６〜８を考慮する。６行×７列の行例〔（ｃ5（ｉ），ｃk（ｉ））〕を考え、同様に、「参照行」を一つだけ固定して、列方向に「参照行」の要素と異なる場合には識別可の記号として−に置換する。「参照行」と同値である場合には、そのままの値を保存する。すなわち、（ｃ5（ｉ），ｃk（ｉ））相互を比較して、少なくとも、一方の値が相違していれば、「識別可」とする。（ｃ5（ｉ），ｃk（ｉ））が同じ値の場合、ｃk（ｉ）の値を表記する。図３に、前記（ｃ5（ｉ），ｃk（ｉ））相互の比較演算処理結果を、５通りの「６行×７列の行例」として、図示する。列の表示：１５、２５、３５、４５、５６、５７、５８は、おのおの、多型１と５、多型２と５、多型３と５、多型４と５、多型５と６、多型５と７、多型５と８の組み合わせを示す。
【０２４２】
（ステップｇ）
比較演算処理を行った結果のうち、Ｈａｐｌｏｔｙｐｅ３を「参照行」とした結果に着目する。Ｈａｐｌｏｔｙｐｅ４の行に対する比較演算処理結果では、多型４と５、多型５と７、多型５と８の組み合わせは、「識別可（−）」となっている。それ以外の４つの組み合わせは、「消去」されていない。
【０２４３】
すなわち、Ｈａｐｌｏｔｙｐｅ３とＨａｐｌｏｔｙｐｅ４との間の識別を行う上では、多型４と５（ｃ4，ｃ5）、多型５と７（ｃ5，ｃ7）、多型５と８（ｃ5，ｃ8）の組み合わせ判定のいずれかが不可欠であることが判る。勿論、Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６の識別を行う際には、Ｈａｐｌｏｔｙｐｅ３とＨａｐｌｏｔｙｐｅ４との間の識別は必須である。換言するならば、Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６の識別を行う上でも、多型４と５（ｃ4，ｃ5）、多型５と７（ｃ5，ｃ7）、多型５と８（ｃ5，ｃ8）の組み合わせ判定のいずれかの判定が不可欠であることになる。従って、「最小数」の多型部位を含む変異セットにおける、多型部位の「第２の候補」として、多型４と５（ｃ4，ｃ5）、多型５と７（ｃ5，ｃ7）、多型５と８（ｃ5，ｃ8）の組み合わせが選択される。
【０２４４】
（ステップｈ）
次に、多型４と５（ｃ4，ｃ5）と、他の一つの多型部位を組み合わせ、（ｃ4，ｃ5，ｃk）：ｋ＝１〜３、６〜８を考慮する。６行×６列の行例〔（ｃ4（ｉ），ｃ5（ｉ），ｃk（ｉ））〕を考え、同様に、「参照行」を一つだけ固定して、列方向に「参照行」の要素と異なる場合には識別可の記号として−に置換する。「参照行」と同値である場合には、そのままの値を保存する。すなわち、（ｃ4（ｉ），ｃ5（ｉ），ｃk（ｉ））相互を比較して、少なくとも、一方の値が相違していれば、「識別可」とする。（ｃ4（ｉ），ｃ5（ｉ），ｃk（ｉ））が同じ値の場合、ｃk（ｉ）の値を表記する。このとき、多型３、４と５（ｃ3，ｃ4，ｃ5）、多型６、４と５（ｃ4，ｃ5，ｃ6）の組み合わせでは、Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６の識別が可能であることが判る。
【０２４５】
同様に、多型５と８（ｃ5，ｃ8）と、他の一つの多型部位を組み合わせ、（ｃ5，ｃ8，ｃk）：ｋ＝１〜４、６、７を考慮する。６行×６列の行例〔（ｃ5（ｉ），ｃ8（ｉ），ｃk（ｉ））〕を考え、同様に、「参照行」を一つだけ固定して、列方向に「参照行」の要素と異なる場合には識別可の記号として−に置換する。「参照行」と同値である場合には、そのままの値を保存する。すなわち、（ｃ5（ｉ），ｃ8（ｉ），ｃk（ｉ））相互を比較して、少なくとも、一方の値が相違していれば、「識別可」とする。（ｃ5（ｉ），ｃ8（ｉ），ｃk（ｉ））が同じ値の場合、ｃk（ｉ）の値を表記する。このとき、多型３、５と８（ｃ3，ｃ5，ｃ8）、多型６、５と８（ｃ5，ｃ6，ｃ8）の組み合わせでは、Ｈａｐｌｏｔｙｐｅ１〜Ｈａｐｌｏｔｙｐｅ６の識別が可能であることが判る。
【０２４６】
また、多型５と７（ｃ5，ｃ7）と、他の一つの多型部位を組み合わせ、（ｃ5，ｃ7，ｃk）：ｋ＝１〜４、６、８を考慮する。６行×６列の行例〔（ｃ5（ｉ），ｃ7（ｉ），ｃk（ｉ））〕を考え、同様に、「参照行」を一つだけ固定して、列方向に「参照行」の要素と異なる場合には識別可の記号として−に置換する。「参照行」と同値である場合には、そのままの値を保存する。すなわち、（ｃ5（ｉ），ｃ7（ｉ），ｃk（ｉ））相互を比較して、少なくとも、一方の値が相違していれば、「識別可」とする。（ｃ5（ｉ），ｃ7（ｉ），ｃk（ｉ））が同じ値の場合、ｃk（ｉ）の値を表記する。いずれの組み合わせでも、まだ、「ハプロタイプ」相互の識別できていない組み合わせが、少なくとも一つ残っている。
【０２４７】
以上の比較演算処理を行った結果を、図４に併せて示す。
【０２４８】
（ステップｉ）
以上のステップｃ〜ｈの結果、全ての「ハプロタイプ」：Ｃi：ｉ＝１〜６を識別することが可能な、「最小数」Ｎmin＝３の多型部位の組み合わせで構成される「変異セット」が選別されている。選別された多型部位の組み合わせは、（ｃ3，ｃ4，ｃ5）、（ｃ4，ｃ5，ｃ6）、（ｃ3，ｃ5，ｃ8）、（ｃ5，ｃ6，ｃ8）の４種のセットである。
【０２４９】
上記のヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の区別も可能な「最小数」の多型部位を含む変異セットが選択されるための必要条件；ＮH≦２N-1は、ＮH＝６では、Ｎ＝４以上である。従って、この「最小数」Ｎmin＝３の多型部位の組み合わせで構成される「変異セット」は、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の区別も可能な「最小数」の多型部位を含む変異セットでは無いことを、実際に検証する。
【０２５０】
（ステップｊ〜ｌ）
（ｃ3，ｃ4，ｃ5）、（ｃ4，ｃ5，ｃ6）、（ｃ3，ｃ5，ｃ8）、（ｃ5，ｃ6，ｃ8）の４種のセットについて、「ハプロタイプ」：Ｃi：ｉ＝１〜６を表記する行ベクトルを、図５に示す。これらの行ベクトルを、相互に比較すると、互いに、逆ベクトルとなっているベクトル組が、２組存在している。
【０２５１】
ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表記する行ベクトルとして、（ｃ3，ｃ4，ｃ5）のセットでは、（｛ｃ3（ｉ）＋ｃ3（ｊ）｝，｛ｃ4（ｉ）＋ｃ4（ｊ）｝，｛ｃ5（ｉ）＋ｃ5（ｊ）｝）を定義する。（ｃ4，ｃ5，ｃ6）のセットでは、（｛ｃ4（ｉ）＋ｃ4（ｊ）｝，｛ｃ5（ｉ）＋ｃ5（ｊ）｝，｛ｃ6（ｉ）＋ｃ6（ｊ）｝）を定義する。（ｃ3，ｃ5，ｃ8）のセットでは、（｛ｃ3（ｉ）＋ｃ3（ｊ）｝，｛ｃ5（ｉ）＋ｃ5（ｊ）｝，｛ｃ8（ｉ）＋ｃ8（ｊ）｝）を定義する。（ｃ5，ｃ6，ｃ8）のセットでは、（｛ｃ5（ｉ）＋ｃ5（ｊ）｝，｛ｃ6（ｉ）＋ｃ6（ｊ）｝，｛ｃ8（ｉ）＋ｃ8（ｊ）｝）を定義する。
【０２５２】
図５には、前記の定義に従って、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表記する行ベクトルを併せて示す。ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」に関して、その「ハプロタイプ」の組み合わせを、例えば、「Ｈａｐｌｏｔｙｐｅ１とＨａｐｌｏｔｙｐｅ２のディプロタイプ」は、「Ｈａｐｌｏｔｙｐｅ１＋２」のように表記している。「＋２」は変異型アレルホモ、「−２」は野生型アレルホモ、「０」はヘテロを意味する。
【０２５３】
図５に示すヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表記する行ベクトルを参照すると、いずれの「変異セット」においても、（０，０，０）を示す「ディプロタイプ」が二つ存在している。少なくとも、この二つの「ディプロタイプ」間の区別が不能であることが判る。
【０２５４】
その他にも、区別不能な複数の「ディプロタイプ」が存在することを確認する。その目的で、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を示す行ベクトル中に、「０」を二つ以上含むものを抽出する。抽出した結果を、図６に示す。この「０」を二つ以上含むヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の相互比較を行う。すなわち、二つのヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表記する行ベクトル間の「差ベクトル」を求める。すなわち、各要素間の差を算出し、行ベクトルの形式で表記する。二つのヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表記する行ベクトルが同じ場合、求められる「差ベクトル」は、（０，０，０）となる。図６の下段には、この二つのヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表記する行ベクトル間の「差ベクトル」を示す。
【０２５５】
（ｃ3，ｃ4，ｃ5）、（ｃ4，ｃ5，ｃ6）、（ｃ3，ｃ5，ｃ8）、（ｃ5，ｃ6，ｃ8）の４種のセットの何れにおいても、（０，０，０）となる「差ベクトル」が三つ存在している。
【０２５６】
変異セット（ｃ3，ｃ4，ｃ5）においては、下記のヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の間で、相互の区別が不能となっている。
【０２５７】
・「Ｈａｐｌｏｔｙｐｅ１とＨａｐｌｏｔｙｐｅ３のディプロタイプ」と「Ｈａｐｌｏｔｙｐｅ５とＨａｐｌｏｔｙｐｅ６のディプロタイプ」；（０，−２，０）
・「Ｈａｐｌｏｔｙｐｅ１とＨａｐｌｏｔｙｐｅ４のディプロタイプ」と「Ｈａｐｌｏｔｙｐｅ２とＨａｐｌｏｔｙｐｅ５のディプロタイプ」；（０，０，０）
・「Ｈａｐｌｏｔｙｐｅ２とＨａｐｌｏｔｙｐｅ３のディプロタイプ」と「Ｈａｐｌｏｔｙｐｅ４とＨａｐｌｏｔｙｐｅ６のディプロタイプ」；（０，０，−２）
変異セット（ｃ4，ｃ5，ｃ6）においては、下記のヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の間で、相互の区別が不能となっている。
【０２５８】
・「Ｈａｐｌｏｔｙｐｅ１とＨａｐｌｏｔｙｐｅ３のディプロタイプ」と「Ｈａｐｌｏｔｙｐｅ５とＨａｐｌｏｔｙｐｅ６のディプロタイプ」；（−２，０，０）
・「Ｈａｐｌｏｔｙｐｅ１とＨａｐｌｏｔｙｐｅ４のディプロタイプ」と「Ｈａｐｌｏｔｙｐｅ２とＨａｐｌｏｔｙｐｅ５のディプロタイプ」；（０，０，０）
・「Ｈａｐｌｏｔｙｐｅ２とＨａｐｌｏｔｙｐｅ３のディプロタイプ」と「Ｈａｐｌｏｔｙｐｅ４とＨａｐｌｏｔｙｐｅ６のディプロタイプ」；（０，−２，０）
変異セット（ｃ3，ｃ5，ｃ8）においては、下記のヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の間で、相互の区別が不能となっている。
【０２５９】
・「Ｈａｐｌｏｔｙｐｅ１とＨａｐｌｏｔｙｐｅ３のディプロタイプ」と「Ｈａｐｌｏｔｙｐｅ２とＨａｐｌｏｔｙｐｅ５のディプロタイプ」；（０，０，−２）
・「Ｈａｐｌｏｔｙｐｅ１とＨａｐｌｏｔｙｐｅ４のディプロタイプ」と「Ｈａｐｌｏｔｙｐｅ５とＨａｐｌｏｔｙｐｅ６のディプロタイプ」；（０，０，０）
・「Ｈａｐｌｏｔｙｐｅ２とＨａｐｌｏｔｙｐｅ４のディプロタイプ」と「Ｈａｐｌｏｔｙｐｅ３とＨａｐｌｏｔｙｐｅ６のディプロタイプ」；（０，−２，０）
変異セット（ｃ5，ｃ6，ｃ8）においては、下記のヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の間で、相互の区別が不能となっている。
【０２６０】
・「Ｈａｐｌｏｔｙｐｅ１とＨａｐｌｏｔｙｐｅ３のディプロタイプ」と「Ｈａｐｌｏｔｙｐｅ２とＨａｐｌｏｔｙｐｅ５のディプロタイプ」；（０，０，−２）
・「Ｈａｐｌｏｔｙｐｅ１とＨａｐｌｏｔｙｐｅ４のディプロタイプ」と「Ｈａｐｌｏｔｙｐｅ５とＨａｐｌｏｔｙｐｅ６のディプロタイプ」；（０，０，０）
・「Ｈａｐｌｏｔｙｐｅ２とＨａｐｌｏｔｙｐｅ４のディプロタイプ」と「Ｈａｐｌｏｔｙｐｅ３とＨａｐｌｏｔｙｐｅ６のディプロタイプ」；（−２，０，０）
以上のように、「ハプロタイプ」の識別が可能な「最小数」の多型部位を含む変異セットに関して、上記の必要条件；ＮH≦２N-1を満たさない場合、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」間に区別不能な組み合わせが残ることが検証される。
【０２６１】
従って、必要条件；ＮH≦２N-1；ＮH＝６では、Ｎ＝４以上を満たすように、４つの多型部位の組み合わせで構成される、「ハプロタイプ」の識別が可能な「変異セット」を新たに見つける。
【０２６２】
その際、（ｃ3，ｃ4，ｃ5）、（ｃ4，ｃ5，ｃ6）、（ｃ3，ｃ5，ｃ8）、（ｃ5，ｃ6，ｃ8）の４種のサブ・セット自体でも、「ハプロタイプ」の識別が可能である。従って、さらに一つの多型部位を付加した、（ｃ3，ｃ4，ｃ5，ｃk）、（ｃ4，ｃ5，ｃ6，ｃk）、（ｃ3，ｃ5，ｃ8，ｃk）、（ｃ5，ｃ6，ｃ8，ｃk）型の「変異セット」も、「ハプロタイプ」の識別が可能な「変異セット」である。
【０２６３】
（ｃ3，ｃ4，ｃ5，ｃk）、（ｃ4，ｃ5，ｃ6，ｃk）、（ｃ3，ｃ5，ｃ8，ｃk）、（ｃ5，ｃ6，ｃ8，ｃk）型の「変異セット」として、下記の１６種類が存在する。
【０２６４】
・（ｃ3，ｃ4，ｃ5，ｃk）型の「変異セット」
（１，３，４，５）、（２，３，４，５）、（３，４，５，６）、（３，４，５，７）、（３，４，５，８）；５種類
・（ｃ4，ｃ5，ｃ6，ｃk）型の「変異セット」
（１，４，５，６）、（２，４，５，６）、（４，５，６，７）、（４，５，６，８）；（５−１）種類
・（ｃ3，ｃ5，ｃ8，ｃk）型の「変異セット」
（１，３，５，８）、（２，３，５，８）、（３，５，６，８）、（３，５，７，８）：（５−１）種類
・（ｃ5，ｃ6，ｃ8，ｃk）型の「変異セット」
（１，５，６，８）、（２，５，６，８）、（５，６，７，８）：（５−２）種類
図４に示す、（ｃ3，ｃ4，ｃ5）、（ｃ4，ｃ5，ｃ6）、（ｃ3，ｃ5，ｃ8）、（ｃ5，ｃ6，ｃ8）以外のサブ・セット；合計｛（３×６）−４｝のサブ・セットについて、一つの多型部位を付加したセットを考える。上記と同様の演算処理を行って、得られる行列型の演算結果に基づき、「ハプロタイプ」の識別が可能な「変異セット」をさらに探索する。「ハプロタイプ」の識別が可能な「変異セット」として、下記の「変異セット」１３種が、さらに選択される。
【０２６５】
・（ｃ1，ｃ4，ｃ5，ｃk）型の「変異セット」
（１，２，４，５）、（１，４，５，８）；２種類
・（ｃ2，ｃ4，ｃ5，ｃk）型の「変異セット」
（２，４，５，７）；１種類
・（ｃ4，ｃ5，ｃ7，ｃk）型の「変異セット」
（４，５，７，８）；１種類
・（ｃ1，ｃ5，ｃ7，ｃk）型の「変異セット」
（１，２，５，７）、（１，３，５，７）、（１，５，６，７）、（１，５，７，８）；４種類
・（ｃ2，ｃ5，ｃ7，ｃk）型の「変異セット」
（２，３，５，７）、（２，５，６，７）、（２，５，７，８）；３種類
・（ｃ3，ｃ5，ｃ7，ｃk）型の「変異セット」
（３，５，７，８）；１種類
・（ｃ1，ｃ5，ｃ8，ｃk）型の「変異セット」
（１，２，５，８）；１種類
従って、１６種類＋１３種類の「変異セット」が、Ｎ＝４の多型部位の組み合わせで構成される、「ハプロタイプ」の識別が可能な「変異セット」として、選別される。
【０２６６】
この２９種の「変異セット」について、各「変異セット」における、６種の「ハプロタイプ」：Ｃi：ｉ＝１〜６を表記する行ベクトルを作成する。
【０２６７】
例えば、変異セット（ｃ1，ｃ4，ｃ5，ｃ8）における、６種の「ハプロタイプ」：Ｃi：ｉ＝１〜６を表記する行ベクトル（ｃ1（ｉ），ｃ4（ｉ），ｃ5（ｉ），ｃ8（ｉ））を、図７に示す。また、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表す行ベクトル（｛ｃ1（ｉ）＋ｃ1（ｊ）｝，｛ｃ4（ｉ）＋ｃ4（ｊ）｝，｛ｃ5（ｉ）＋ｃ5（ｊ）｝，｛ｃ8（ｉ）＋ｃ8（ｊ）｝）を定義する。このヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表す行ベクトルを、図７の下段に示す。次いで、二つのヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表す行ベクトルの「差ベクトル」を算出する。算出される「差ベクトル」を、図８に示す。
【０２６８】
まず、図７に示す、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表す行ベクトル中に、（０，０，０，０）となるものは、一つしか存在していない。また、図８に示される「差ベクトル」中に、（０，０，０，０）となるものは存在していない。すなわち、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表す行ベクトルが、同一となるものは存在していない。換言すると、この変異セット（ｃ1，ｃ4，ｃ5，ｃ8）を用いる際には、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の全てが識別可能であることが判る。
【０２６９】
同様の検証操作を、残る２８種類の「変異セット」に関しても、実施して、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の全てを識別可能な「変異セット」を選択する。以上の検証操作の結果、下記する合計５種の「変異セット」が、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の全てを識別可能である「変異セット」として、選択された。
【０２７０】
「変異セット」（１，４，５，８）、（４，５，７，８）、（１，５，７，８）、（３，４，５，８）、（４，５，６，８）
なお、この５種の「変異セット」は、「ハプロタイプ」の全てを識別可能であるので、ホモ型（Ｃi、Ｃi）「ディプロタイプ」も識別可能である。
【０２７１】
上記の例では、必要条件；（ＮH−１）＜ＮH≦２N-1を満たしている、４つの多型部位の組み合わせで構成される「変異セット」中に、「ハプロタイプ」のみでなく、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の識別も可能である「変異セット」が見出されている。従って、前記の５種の「変異セット」は、「ハプロタイプ」のみでなく、「ディプロタイプ」の識別も可能である、「最小数」の多型部位で構成される「変異セット」となっている。
【０２７２】
上記の例では、「ディプロタイプ」：（Ｃi、Ｃj）を表す行ベクトルを、「ハプロタイプ」ＣiとＣjを表す行ベクトルの「和」：（Ｃi＋Ｃj）を示す行ベクトルとして、定義している。この「和ベクトル」：（Ｃi＋Ｃj）の総数は、「ハプロタイプ」の種類ＮHに対して、ＮH＋1/2ＮH・（ＮH−１）＝1/2ＮH・（ＮH＋１）となっている。一方、Ｎ次元で、各要素が、「＋２」、「−２」、「０」のいずれである行ベクトルの種類の総数は、３Nである。「ディプロタイプ」：（Ｃi、Ｃj）を表す行ベクトルの相違を識別できるためには、1/2ＮH・（ＮH＋１）≦３Nの条件を満たすことが必要である。例えば、ＮH＝６の場合、1/2ＮH・（ＮH＋１）＝２１となり、Ｎ＝３の場合、３N＝２７となり、前記の条件を満たしてはいる。しかし、実際には、Ｎ＝３の場合には、「ディプロタイプ」：（Ｃi、Ｃj）を表す行ベクトルの相違を識別できる「変異セット」は存在しない。その点からも、上記の必要条件；（ＮH−１）＜ＮH≦２N-1を採用することが適切であることが示される。
【０２７３】
上記の実施形態においては、表現される「ディプロタイプ」が、同一か否かを判定する手段として、「ディプロタイプ」を表記する行ベクトル間で、「差ベクトル」を算出し、算出される「差ベクトル」が、零ベクトルか否かを判定基準としている。この二つの行ベクトルが同一か否かを判定する手段としては、他の比較演算処理を採用することもできる。例えば、二つの行ベクトルＡ，Ａ’の長さ（ｎｏｒｍ（Ａ）、ｎｏｒｍ（Ａ’））をそれぞれ算出する。次いで、二つの行ベクトルの内積：Ａ・Ａ’を算出し、ｎｏｒｍ（Ａ）＝ｎｏｒｍ（Ａ’）、Ａ・Ａ’＝ｎｏｒｍ（Ａ）・ｎｏｒｍ（Ａ’）の場合、二つの行ベクトルが同一であると判定する手段を採用することもできる。
【０２７４】
上記の実施形態では、各多型部位における塩基種類（変異）を、「野生型アレル（Ｗ）」と「変異型アレル（Ｍ）」の二つの区分に分類して、「２値化」を行っている。「野生型アレル（Ｗ）」を「−１」、「変異型アレル（Ｍ）」を「＋１」と表記している。また、表現される「ディプロタイプ」においては、「野生型アレル（Ｗ）」のホモ、「変異型アレル（Ｍ）」のホモ、ヘテロの三つの区分となり、それに対応する「三つの値」を利用している。すなわち、「野生型アレル（Ｗ）」のホモ（Ｗ＋Ｗ）を「−２」、「変異型アレル（Ｍ）」のホモ（Ｍ＋Ｍ）を「＋２」、ヘテロ（Ｗ＋Ｍ）を「０」と表記している。
【０２７５】
上記の「ハプロタイプ」を表記する行ベクトル間における要素ｃk（ｉ）相互の比較演算処理が可能である限り、「野生型アレル（Ｗ）」と「変異型アレル（Ｍ）」の「２値化」表記法は、上記の形式以外を採用することができる。その際、「ディプロタイプ」を表記する行ベクトルの要素｛ｃk（ｉ）＋ｃk（ｊ）｝に関して、「三つの値」を与えるように、「和演算」を適宜定義することが可能である。
【０２７６】
例えば、「野生型アレル（Ｗ）」を「Ａ」、「変異型アレル（Ｍ）」を「Ａ*」として、表記することができる。その際、「Ａ」は、絶対値｜Ａ｜＝１の複素数とし、「Ａ*」は、その共役的複素数とする。一方、要素｛ｃk（ｉ）＋ｃk（ｊ）｝を与える、「和演算」を、｛Ａ＋Ａ｝≡Ａ・Ａ，｛Ａ*＋Ａ*｝≡Ａ*・Ａ*，｛Ａ＋Ａ*｝≡Ａ・Ａ*＝１（但し、Ａ・Ａ≠Ａ*・Ａ*）と、「複素数の乗算」の形式で定義することも可能である。
【０２７７】
従来の「ハプロタイプ」の全てを識別可能な「最小数」の多型部位で構成される「変異セット」を選定する方法は、「ハプロタイプ」のみでなく、「ディプロタイプ」の識別も可能であるための必要条件：ＮH≦２N-1を考慮していないものであった。従って、前記の必要条件：ＮH≦２N-1を満足しない「多型部位の個数：Ｎ」、上記の例では、Ｎ＝３の条件で、「ハプロタイプ」の全てを識別可能な「最小数」の多型部位で構成される「変異セット」を選択している。このＮ＝３の条件で選択される「変異セット」は、必要条件：ＮH≦２N-1を満足してないため、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」のうち、幾つかに関しては、その組み合わせ（Ｃi、Ｃj）の識別ができないものであった。
【０２７８】
本発明では、「ハプロタイプ」のみでなく、「ディプロタイプ」の識別も可能である、「最小数」の多型部位で構成される「変異セット」を選択するため、前記の必要条件：ＮH≦２N-1を満足する「変異セット」を対象として、選別操作を行っている。その際、必要条件：ＮH≦２N-1を満足する「変異セット」中から、「ハプロタイプ」の全てを識別可能な「最小数」の多型部位で構成される複数個の「変異セット」を、候補として選択する。次いで、この「変異セット」候補について、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の全てを識別可能か否かの検証操作を施し、実際に、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の全てを識別可能なものを選別している。
【０２７９】
本発明の「最小数」の多型部位で構成される「変異セット」の選別法では、各多型部位における、塩基種類を、「野生型アレル（Ｗ）」と「変異型アレル（Ｍ）」の二つの区分に分類して、「２値化」を行っている。従って、「一塩基多型」のみでなく、インサーション型、デリーション型の変異であっても、「野生型アレル（Ｗ）」と「変異型アレル（Ｍ）」の二つの区分に分類して、「２値化」可能なものも対象とできる。
【０２８０】
一方、「ハプロタイプ」と、「ハプロタイプ」対で表記される「ディプロタイプ」を対象とするため、一般に、一組の相同染色体上の、相同の遺伝子座に、一対のアレル（対立遺伝子）が存在している系を対象とする。すなわち、識別の対象となる「ハプロタイプ」対が、一組の相同染色体上の、相同の遺伝子座に、一対のアレル（対立遺伝子）として存在する限り、個体の生物種、属、株に制限はない。通常、一つの遺伝子内に存在する複数個の多型部位によって決定される「ハプロタイプ」を対象とする。しかしながら、複数の遺伝子からなる「遺伝子群」を対象として、この「遺伝子群」中に散在している複数個の多型部位によって決定される「ハプロタイプ」を対象とすることもできる。
【０２８１】
本発明では、既に「塩基配列」が判明している複数種の「ハプロタイプ」を対象として、その塩基配列中に存在する「多型部位」が、「野生型アレル（Ｗ）」と「変異型アレル（Ｍ）」のいずれかを測定することを前提としている。本発明にかかる方法は、この測定すべき「多型部位」の種類を選択する際に利用される。すなわち、各個体が有している「ハプロタイプ」対に対して、該遺伝子の塩基配列中に存在している複数個の「多型部位」のうち、測定すべき「最小数」の「多型部位」の種類（位置）を選択する際に利用される。各「多型部位」の種類（位置）に応じて、その部位の塩基（変異）が「野生型アレル（Ｗ）」と「変異型アレル（Ｍ）」のいずれかを判定するため、それぞれ、専用のプローブまたは、配列解析用のプライマーを作製する。その際、作製すべきプローブまたは、プライマーの種類の総数を、必要最小限に留める上では、本発明にかかる選別方法を応用して、「最小数」の多型部位で構成される「変異セット」を決定することが有効である。勿論、「最小数」の多型部位で構成される「変異セット」に対応する「変異測定」を行うため、測定に利用する検出デバイスのコスト・パフォーマンスを向上させる目的でも有効である。
【産業上の利用可能性】
【０２８２】
本発明にかかる「塩基配列を識別する変異セット」の選別法は、その塩基配列中に変異を有する遺伝子の変異分析、また、該遺伝子の変異分析を応用した診断の分野において、該遺伝子の変異分析に利用される核酸プローブ、プライマーを設計する際利用できる。
【図面の簡単な説明】
【０２８３】
【図１】本発明の第１の実施形態において、識別の対象とする「ハプロタイプ」６種と、「変異セット」を構成する際、その要素の候補となる、８箇所の多型部位を説明する図である。「ハプロタイプ」６種が示す該８箇所の多型部位における塩基の種別；変異型アレル（Ｍ）、野生型アレル（Ｗ）を、行列形式で表記している。
【図２】本発明の第１の実施形態において、「ハプロタイプ」６種を識別可能な、「最小数」の多型部位で構成される「変異セット」を選定する工程中、該「変異セット」に含まれる多型部位の「第１の候補」（ｃk1）を選択するステップｃ〜ｅを説明する図である。「ハプロタイプ」６種を表記する行ベクトルの一つを「参照行」とし、他の行ベクトルと、要素｛ｃk｝が同じか、異なるかの比較演算を行い、異なると判断される結果を、「−」で示す。
【図３】本発明の第１の実施形態において、「ハプロタイプ」６種を識別可能な、「最小数」の多型部位で構成される「変異セット」を選定する工程中、該「変異セット」に含まれる多型部位の「第２の候補」（ｃk2）を選択するステップｆ〜ｇを説明する図である。「第１の候補」（ｃk1）と、他の多型部位（ｃk）との組み合わせ（ｃk1，ｃk）について、その要素対｛ｃk1，ｃk｝を考慮する。「ハプロタイプ」６種を表記する行ベクトルの一つを「参照行」とし、他の行ベクトルと、要素対｛ｃk1，ｃk｝が同じか、異なるかの比較演算を行い、異なると判断される結果を、「−」で示す。
【図４】本発明の第１の実施形態において、「ハプロタイプ」６種を識別可能な、「最小数」の多型部位で構成される「変異セット」を選定する工程中、該「変異セット」に含まれる多型部位の「第２の候補」（ｃk2）を選択するステップｈを説明する図である。「第１の候補」（ｃk1）、「第２の候補」（ｃk2）と、他の多型部位（ｃk）との組み合わせ（ｃk1，ｃk2，ｃk）について、その要素ベクトル｛ｃk1，ｃk2，ｃk｝を考慮する。「ハプロタイプ」６種を表記する行ベクトルの一つを「参照行」とし、他の行ベクトルと、要素ベクトル｛ｃk1，ｃk2，ｃk｝が同じか、異なるかの比較演算を行い、異なると判断される結果を、「−」で示す。
【図５】本発明の第１の実施形態において、３箇所の多型部位で構成される「変異セット」中から選別された、「ハプロタイプ」６種を識別可能な「変異セット」を示す図である。また、該「変異セット」において、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表記する行ベクトルの要素｛ｃk（ｉ）＋ｃk（ｊ）｝を併せて示す。
【図６】本発明の第１の実施形態において、「ハプロタイプ」６種を識別可能な、３箇所の多型部位で構成される「変異セット」において、表現が等しくなるヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を抽出する（ステップｋ）の操作を示す図である。上段に、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表記する行ベクトル中、「０」と表記される要素を２以上含むものを選別した結果をしめす。下段には、上段に示すヘテロ型（Ｃi、Ｃj）「ディプロタイプ」において、二つのヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表記する行ベクトルの「差ベクトル」の算出結果を示す。「差ベクトル」が、零ベクトルとなり、二つのヘテロ型（Ｃi、Ｃj）「ディプロタイプ」が同一の表現を示す組み合わせを、網掛け表示で示す。
【図７】本発明の第１の実施形態において、４箇所の多型部位で構成される「変異セット」中から選別された、「ハプロタイプ」６種を識別可能な「変異セット」の一例：（ｃ1，ｃ4，ｃ5，ｃ8）を示す図である。また、該「変異セット」において、ヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表記する行ベクトルの要素｛ｃk（ｉ）＋ｃk（ｊ）｝を併せて示す。
【図８】本発明の第１の実施形態において、「ハプロタイプ」６種を識別可能な、４箇所の多型部位で構成される「変異セット」において、表現が等しくなるヘテロ型（Ｃi、Ｃj）「ディプロタイプ」の有無を確認する（ステップｋ）の操作を示す図である。二つのヘテロ型（Ｃi、Ｃj）「ディプロタイプ」を表記する行ベクトルの「差ベクトル」の算出結果を示す。
【図１】
2008165375.tif 000002
【図２】
2008165375.tif 000003
【図３】
2008165375.tif 000004
【図４】
2008165375.tif 000005
【図５】
2008165375.tif 000006
【図６】
2008165375.tif 000007
【図７】
2008165375.tif 000008
【図８】
2008165375.tif 000009