(57)【特許請求の範囲】
【請求項１】
未知の核酸断片とハイブリダイゼーション反応させて遺伝子解析を行うためのプローブに用いられる塩基配列を設計するプローブ設計装置を用いてプローブを設計するプローブ設計方法であって、
生成手段が、ターゲットの塩基配列に基づいて得られる複数の部分塩基配列をノードの経路上に配置したツリーを生成する生成工程と、前記ツリーにおいて、各階層は部分塩基配列における各塩基位置と１対１に対応し、各ノードは部分塩基配列における所定位置の塩基に対応し、
評価手段が、前記ツリーから選択されたノードに関して、該選択されたノードと接続する前記ツリー上の経路に存在する各ノードに配置された部分塩基配列の数に基づいて、該選択されたノードが表す部分塩基配列の各位置における特異性を算出し、算出された複数の特異性とそれらの位置とに基づいてプローブとしての利用可能性を評価する評価工程と、
決定手段が、前記評価工程による評価結果に基づいてプローブとして使用すべき部分塩基配列を決定する決定工程とを備えることを特徴とするプローブ設計方法。
【請求項２】
前記生成工程における複数の部分塩基配列は、前記ターゲットの塩基配列の相補的な塩基配列から得られる部分塩基配列であることを特徴とする請求項１に記載のプローブ設計方法。
【請求項３】
前記生成工程における複数の部分塩基配列は、前記ターゲットの塩基配列から得られる部分塩基配列であり、
前記決定工程は、前記評価工程の評価結果に基づいて部分塩基配列を選択し、選択された部分塩基配列の相補的な塩基配列をプローブとして使用すべき部分塩基配列として決定することを特徴とする請求項１に記載のプローブ設計方法。
【請求項４】
前記生成工程は、前記ターゲットの塩基配列に基づいて得られる全ての部分塩基配列を分類するツリーを生成することを特徴とする請求項１乃至３のいずれか１項に記載のプローブ設計方法。
【請求項５】
前記評価工程は、前記ツリー上のノードの前記選択されたノードが表す部分塩基配列の中央付近に、ターゲットの塩基配列に対する特異性が変化する塩基配列が存在する場合に、プローブとして適切であるとする評価関数を導入することを特徴とする請求項１乃至３のいずれか１項に記載のプローブ設計方法。
【請求項６】
前記評価工程は、前記選択されたノードと接続する前記ツリー上の経路に存在する各ノードについて、前記特異性としてエントロピーを算出し、該算出されたエントロピーの減少状態に基づいて評価を行なうことを特徴とする請求項１乃至３のいずれか１項に記載のプローブ設計方法。
【請求項７】
前記評価工程は、前記選択されたノードに対応する部分塩基配列の中央付近における前記エントロピーの変化量に対して他の部分よりも大きい重み値を与える評価関数を導入することを特徴とする請求項６に記載のプローブ設計方法。
【請求項８】
前記決定工程は、前記評価工程において前記評価関数によって求めた値が所定値を超えるノードに対応する部分塩基配列をプローブとして決定することを特徴とする請求項５または７に記載のプローブ設計方法。
【請求項９】
前記決定工程は、前記エントロピーの変化量が所定値を超えるノードに対応する部分塩基配列をプローブとして決定することを特徴とする請求項６に記載のプローブ設計方法。
【請求項１０】
前記複数の部分塩基配列を前記ターゲットの塩基配列に対する特異性によってグループ分けするグループ化工程を更に備え、
前記決定工程は、前記評価工程の評価結果に基づいて各グループよりプローブとして用いる部分塩基配列を決定することを特徴とする請求項１乃至９のいずれか１項に記載のプローブ設計方法。
【請求項１１】
グループ化手段が、前記評価工程によりプローブとして利用可能であると評価された複数の部分塩基配列を前記ターゲットの塩基配列に対する特異性によってグループ分けするグループ化工程と、
選択手段が、前記グループ化工程で得られたグループのうち、解析対象のターゲットの塩基配列に基づいて使用すべきグループを選択する選択工程とを更に備え、
前記決定工程は前記選択工程で選択されたグループのそれぞれより、前記評価工程による評価結果に基づいてプローブとして用いる部分塩基配列を決定することを特徴とする請求項１乃至９のいずれか１項に記載のプローブ設計方法。
【請求項１２】
前記選択工程は、情報的に完全に独立で、ターゲットの塩基配列の検出に必要十分なグループのみを選択することを特徴とする請求項１１に記載のプローブ設計方法。
【請求項１３】
前記選択工程は、少なくとも解析対象の全てのターゲットの塩基配列に対して特異性のないグループを排除することを特徴とする請求項１１に記載のプローブ設計方法。
【請求項１４】
前記ターゲットの塩基配列は複数の塩基配列パターンを含み、
前記グループ化工程は、前記複数の塩基配列パターンのそれぞれとの反応の有無が同じ部分塩基配列を同一グループに所属させることを特徴とする請求項１０乃至１３のいずれか１項に記載のプローブ設計方法。
【請求項１５】
前記ツリーにおいて、ノードの接続によって表わされる部分塩基配列の塩基配列順序が、前記ターゲット中の塩基配列順序と一致していることを特徴とする請求項１乃至１４のいずれか１項に記載のプローブ設計方法。
【請求項１６】
前記ツリーにおいて、ノードの接続によって表わされる部分塩基配列の塩基配列順序が、前記ターゲット中の対応する部分塩基配列の塩基配列順序をその中央部が先頭となるように並べ替えられていることを特徴とする請求項１乃至１４のいずれか１項に記載のプローブ設計方法。
【請求項１７】
前記評価工程は、予め指定された範囲の長さの部分塩基配列のみを評価することを特徴とする請求項１乃至１６のいずれか１項に記載のプローブ設計方法。
【請求項１８】
前記評価工程は、予め指定された範囲の融解温度条件に適合する部分塩基配列のみを評価すること特徴とする請求項１乃至１６のいずれか１項に記載のプローブ設計方法。
【請求項１９】
前記決定工程は、予め指定された範囲の長さの部分塩基配列の中から、前記評価工程による評価結果に基づいてプローブとしての部分塩基配列を決定することを特徴とする請求項１乃至１６のいずれか１項に記載のプローブ設計方法。
【請求項２０】
前記決定工程は、予め指定された範囲の融解温度条件に適合する部分塩基配列の中から、前記評価工程による評価結果に基づいてプローブとしての部分塩基配列を決定することを特徴とする請求項１乃至１６のいずれか１項に記載のプローブ設計方法。
【請求項２１】
前記特異性の算出は、エントロピーを算出する関数を用いるものであり、前記プローブとしての利用可能性を評価するための評価関数に、前記エントロピーの減少があると高いスコアになる評価関数を用いることを特徴とする請求項１に記載のプローブ設計方法。
【請求項２２】
請求項１乃至２１のいずれか１項に記載のプローブ設計方法を実現する情報処理装置。
【請求項２３】
請求項１乃至２１のいずれか１項に記載のプローブ設計方法をコンピュータに実現させるためのプログラム。
【請求項２４】
請求項１乃至２１のいずれか１項に記載のプローブ設計方法をコンピュータに実現させるためのプログラムを格納する記憶媒体。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は複数の核酸プローブを配置したマイクロアレイの設計支援に好適な情報処理装置及びプローブ設計方法に関する。
【０００２】
【従来の技術】
従来より、特開平10-272000号公報または特開平11-187900号公報にあるように、ＤＮＡマイクロアレイを用いて遺伝子の発現やシークエンスを決定するシステムが提案されている。このシステムは、あらかじめ試料とハイブリダイズさせる塩基配列断片（プローブ）を設計する必要がある。良いプローブセットを設計することができれば、非常に高い確率で試料の中に存在する塩基配列の断片に対する大量の情報が得られることが特徴である。
【０００３】
【発明が解決しようとする課題】
しかしながら、解析対象の試料（ターゲット）の中に存在する可能性のある塩基配列断片は、非常に多様であり、実験系が異なると全く異なる。この場合、全く異なるプローブを設計する必要がある。従来は、このプローブ設計を経験に基づいて人手で行なっていたが、大量の塩基配列が決定されつつある今では、人手でいちいちプローブ設計を行うことは、事実上不可能となりつつある。
【０００４】
【課題を解決するための手段】
本発明は、上記の課題に鑑みてなされたものであり、解析対象となるターゲットの塩基配列に応じて、解析に適切なプローブを自動的に選択することを可能とし、プローブ設計を効果的に支援することを目的とする。
【０００５】
上記の目的を達成するための本発明によるプローブ設計方法は、
未知の核酸断片とハイブリダイゼーション反応させて遺伝子解析を行うためのプローブに用いられる塩基配列を設計するプローブ設計装置を用いてプローブを設計するプローブ設計方法であって、
生成手段が、ターゲットの塩基配列に基づいて得られる複数の部分塩基配列をノードの経路上に配置したツリーを生成する生成工程と、前記ツリーにおいて、各階層は部分塩基配列における各塩基位置と１対１に対応し、各ノードは部分塩基配列における所定位置の塩基に対応し、
評価手段が、前記ツリーから選択されたノードに関して、該選択されたノードと接続する前記ツリー上の経路に存在する各ノードに配置された部分塩基配列の数に基づいて、該選択されたノードが表す部分塩基配列の各位置における特異性を算出し、算出された複数の特異性とそれらの位置とに基づいてプローブとしての利用可能性を評価する評価工程と、
決定手段が、前記評価工程による評価結果に基づいてプローブとして使用すべき部分塩基配列を決定する決定工程とを備える。
【０００６】
また、上記の目的を達成するための本発明の他の態様によるプローブ設計方法は、
未知の核酸断片とハイブリダイゼーション反応させて遺伝子解析を行うためのプローブに用いられる塩基配列を設計するプローブ設計装置を用いてプローブを設計するプローブ設計方法であって、
生成手段が、ターゲットとなる塩基配列データより得られる複数の部分塩基配列のリストを分類する判別ツリーを生成する生成工程と、前記ツリーにおいて、各階層は部分塩基配列における各塩基位置と１対１に対応し、各ノードは部分塩基配列における所定位置の塩基に対応し、
評価手段が、前記判別ツリーに存在するプローブ候補に対して、前記プローブ候補に対応する前記ツリー上の経路に存在する各ノードに配置された部分塩基配列の数に基づいて、該プローブ候補が有する部分塩基配列の各位置における特異性を算出し、算出された複数の特異性とそれらの位置とに基づいてプローブとしての利用可能性を評価する評価工程と、
選択手段が、前記評価工程の評価結果に基づいて使用すべきプローブを選択する選択工程とを備える。
【０００７】
【発明の実施の形態】
以下、添付の図面に基づいて本発明の好適な実施形態を説明する。
【０００８】
＜第１実施形態＞
［装置及び動作の概要］
以下の実施形態では、いわゆるＤＮＡマイクロアレイを用いた核酸配列解析システムに用いる最適なオリゴヌクレオチドプローブを計算機を用いて設計するプローブ設計方法とそれを実行する装置を説明する。
【０００９】
図２は、本実施形態によるプローブ設計方法が適用される情報処理装置の構成を示すブロック図である。本実施形態では、外部記憶装置２０１、中央処理装置（ＣＰＵ）２０２、メモリ２０３、入出力装置２０４から構成される情報処理装置に、以下に説明するプローブ設計方法を実行させる。このような構成の情報処理装置としては、例えばパーソナルコンピュータが挙げられる。
【００１０】
外部記憶装置２０１は、本実施形態によるプローブ設計方法を実現するプログラムや、ターゲットの塩基配列データ及びパラメータを保持する。また、本実施形態によって導かれたプローブ配列を保持するのにも使用される。中央処理装置（ＣＰＵ）２０２は、プローブ設計のプログラムを実行したり、すべての装置の制御を行なったりする。
【００１１】
メモリ２０３は中央処理装置（ＣＰＵ）２０２が実行するプログラム、及びサブルーチンやデータを一時的に記録する。外部記憶装置２０１に格納されたプローブ設計方法の制御プログラムは、メモリ２０３にロードされ、中央処理装置２０２により実行される。入出力装置２０４は、ユーザとのインタラクションを行う。多くの場合、以下に説明するプローブ設計方法を実現するプログラム実行のトリガはこの入出力装置を介してユーザーが出すことになる。また、ユーザーが結果を見たり、プログラムのパラメータ制御をこの入出力装置を介して行う。
【００１２】
図１は、第１の実施形態によるプローブ設計方法の処理手順を説明するフローチャートである。本実施形態のプローブ設計方法処理概要を説明すると、まず、ステップＳ１０１において、外部記憶装置２０１からターゲットの塩基配列データを取得する。次にステップＳ１０２において、ステップＳ１０１で取得した塩基配列データに基づいて判別ツリーを作成する。判別ツリー及びその作成については後述する。そして、ステップＳ１０３において、ステップＳ１０２で作成した判別ツリーを外部記憶装置２０１に格納する。
【００１３】
ステップＳ１０４では、ステップＳ１０３で格納した判別ツリーを用いて、ツリー上に表現されている各プローブ（部分塩基配列）の評価を行なう。すなわち、判別ツリーの中に存在するプローブ候補の中からプローブとしてふさわしいものを選択するための評価を行なう。ステップＳ１０５では、上記評価の結果をもとに、最終結果としての最適プローブセットを選択する。ステップＳ１０６では最適なプローブを出力する。出力先としては、外部記憶装置２０１、入出力装置２０４に含まれる表示装置或いはプリンタ等が挙げられる。
【００１４】
以下、各処理について詳細に説明する。
【００１５】
［判別ツリーの生成について（ステップＳ１０１〜Ｓ１０３）］
図３は、ステップＳ１０１で取得する塩基配列データの例を示す図である。図３にある塩基配列情報はＨＬＡまたはＭＨＣと呼ばれる遺伝子群の一部である。この部分は、人間の抗体生成に関連する遺伝子をコーディングしており個体差がある。図３に類似しているが完全に同一ではない塩基配列があるのは、このことを意味している。このＨＬＡに関して現在人類で確認されている遺伝子タイプは１００種類以上あり、この種類が異なると生体間移植の際に致命的な拒絶反応がおこる可能性が高い。現在では、患者と臓器提供者のＤＮＡをシークエンスし、完全に塩基配列を決定することによって、このＨＬＡのタイピング（分類）を行っている。この作業を適切なＤＮＡマイクロアレイを用いて行うことによって、簡単かつ高速に行うことができるようになる。
【００１６】
図４はＤＮＡマイクロアレイ上のハイブリダイゼーションの様子を示した図である。生体内でほとんどの場合、ＤＮＡ塩基配列は２重らせん構造をしていて、その２本鎖の間の結合は塩基間の水素結合で実現されている。一方、ＲＮＡ塩基配列は１本で存在する場合が多い。塩基の種類はＤＮＡの場合はＡＣＧＴの４種類、ＲＮＡの場合はＡＣＧＵの４種類であり、それぞれ水素結合ができる塩基対はＡ−Ｔ（Ｕ）、Ｇ−Ｃのペアとなっている。以下では、説明のため、ＤＮＡ塩基配列に関するプローブの作成を例に挙げる。
【００１７】
ハイブリダイゼーションとは、１本鎖状態の塩基配列分子同士がある部分の相補的な塩基配列を介して部分的に結合して２本鎖状態を形成することをいう。本実施形態で想定している反応は、図４の上側の基板にくっついた塩基配列（プローブ配列）の方が下側の試料中にある塩基配列分子より短い。よって、試料中に存在する塩基配列分子がプローブ配列を含む場合は、このハイブリダイゼーション反応はうまくいき、試料中のターゲット塩基配列分子はトラップされることとなる。
【００１８】
但し、プローブ配列の全ての領域が相補的な場合のみ、ハイブリダイゼーションが可能なわけではなく、一部にペアを形成しない部分があっても、ターゲット塩基配列分子はトラップされることがある。特に、プローブの末端だけの塩基が結合できないような場合は、ハイブリダイゼーションを起こしてしまう可能性が高く、ターゲット塩基配列の正確な選別が出来ない。従って、末端だけが異なるプローブ配列セットは、ハイブリダイゼーションを使った実験には不適格といえる。すなわち、塩基配列中の中央部にターゲットの塩基配列と特異となる配列が存在するものが好ましい。また、ハイブリダイゼーション反応の強度はプローブの塩基配列の長さが長いほど高くなる。よって、理想的にはＤＮＡマイクロアレイ上に配置されるプローブ配列はハイブリダイゼーション強度が似ているもの、すなわち含まれる塩基数が近い塩基配列を選ぶほうが良い。
【００１９】
上述のように、最終出力であるプローブ塩基配列は、試料中の部分配列と相補的な塩基配列である必要がある。そこで、本実施形態では、ステップＳ１０２の判別ツリー作成処理において、予め全ての部分塩基配列を用意する。この様子を図にしたのが図５である。
【００２０】
図５はターゲットとなる塩基配列から得られる部分塩基配列（プローブ）を示す図である。図５で列挙している塩基配列群は、図３のDRB1*0101というタイプの配列の部分配列群の一部である。図５のように、末端を１塩基ずつシフトすることによって、全く異なる部分塩基配列が得られ、この全てがプローブとハイブリダイゼーション反応を行う配列候補となる。あるタイプの対象配列数がｎベースだったとすると、その部分塩基配列は全部でｎ−１個存在することとなる。本実施形態では全ての部分塩基配列を分類する判別ツリーを生成するが、予め所定数の深さを設定し、所定数以上の塩基を有する部分塩基配列について判別ツリーを生成するようにしてもよい。
【００２１】
次に、この全ての部分配列を判別／分類する判別ツリーを作成する。図６Ａは判別ツリーを説明する図である。ツリーを構成する全てのノードは、４つの子ノード（Ａ，Ｔ，Ｇ，Ｃ）を持つ。それぞれ４つの子ノードは、注目している位置の塩基配列がＡＴＧＣのいずれかで分類される。まず、一番簡単な例として、注目している位置を１つずつずらす方法を説明する。
【００２２】
図６Ｂは判別ツリーにおける各ノードへの部分塩基配列の分類例を示す図である。ルートノードで、各部分塩基配列の左端の塩基配列を見る。そして、その塩基がどの種類かで部分塩基配列を分類する。例えば、図５に示したＨＬＡのDRB1*0101由来の塩基配列ならば、図５の２番の部分塩基配列はＡのノード６１１に、１，３，５番の部分塩基配列はＧのノード６１２に、４番の部分塩基配列はＣのノード６１３にそれぞれ分類される。そして、１，３，５番の部分配列は、更に深いノードで細分されることとなる。具体的には、１番目の部分塩基配列の次のノードはＡであり、３番目の部分塩基配列の次のノードはＣであり、５番目の部分塩基配列の次のノードはＧである。よって、それぞれ、ノード６１４、ノード６１５、ノード６１６に格納されることになる。
【００２３】
以上の様な処理により、全てのターゲットの全ての部分塩基配列について各ノードへの登録を行なって、判別ツリーを作成する。この判別ツリー作成処理では、図５に示したような部分配列に展開したデータを全て対象として判別ツリーを作成するので、作成された判別ツリーの上のどこかに任意の部分塩基配列が登録されることとなる。上記例のように位置を左から一つずつずらす方法で判別ツリーを作成した場合には、ある部分塩基配列の登録されたノードをルートノードへ向けて逆に辿ると、登録された塩基配列の左側の部分がわかることとなる。
【００２４】
図７及び図８は判別ツリー上のあるノードに登録された部分塩基配列の様子を示す図である。例として、図３に示した４種類の５０ベースからなる塩基配列が試料に入っていると仮定する。図６Ａで示した判別ツリーの中で、ＣＡＧの順でノードを辿っていくと、図６Ａのノード６０１には、図７の表の右側に示した情報が格納されている。すなわち、ノード６０１には、DRB1*0101…「NULL」、DRB1*04011…「21」、DRB1*07011…「10,37」、DRB1*15011…「21,37」が格納されている。この情報によれば、DRB1*0101にはＣＡＧの部分塩基配列は存在せず（NULL）、DRB1*04011には２１番目の場所にＣＡＧが存在し、DRB1*07011には１０番目と３７番目、DRB1*15011には２１番目にそれぞれＣＡＧが存在することが示されている。
【００２５】
図７で示したノードに接続される子ノードの中でＴとなるノード、すなわち図６Ａのノード６０２の内容を示したのが図８である。このノード６０２に対応する部分塩基配列はＣＡＧＴであり、これにマッチングする部分配列の位置は、DRB1*07011の１０番目の部分配列のみである。
【００２６】
以上示してきた例は、塩基の種類を左から１つずつずらした位置で見てきた例であり、ノードの順序と塩基配列とが一致している。しかしながら、この位置はアルゴリズムによって変わってもよい。ハイブリダイゼーション反応の場合、中心の塩基配列の相違が重要である。よって、判別ツリーを作成する際に、注目する塩基配列の位置を中心から周辺へ移していくように並べ替えることが一つの有力な方法としてあげられる。
【００２７】
例えば、図７の左からＣ→Ａ→Ｇとう並びの部分塩基配列に相当するようなノードを、「真中のＡ」→「右のＧ」→「左のＣ」というふうに見ていくことがその方法にあたる。この場合、左から１つずつ見るのではなく、真中、右端、左端の順番で見ていくことになる。従って、左から見る方法だと図７のノード６０１は｛ルートノード→Ｃのノード→Ａのノード→Ｇのノード｝という位置に存在するが、中心から周辺へ見る方法だと、｛ルートノード→Ａのノード→Ｇのノード→Ｃのノード｝という位置に存在することになる。このことは、注目する塩基配列（上記例ではＣＡＧ）の位置を変えると図６Ａにおいてノード６０１がノード６０３（Ａ→Ｇ→Ｃ）に変わることを意味する。この場合、ノード６０１だけが変わるのではなく、判別ツリー全体の構造が変化する。
【００２８】
以上のような並べ替えを一般化して示せば、ｎ個の塩基配列について、ｎが奇数の場合は(n+1)/2を、ｎが偶数の場合はn/2をａとした場合に、「a, a+1, a-1, a+2, a-2, … 」というように並べ替えて判別ツリーを構成することになる。もちろん、中心、左端、右端の順番で見ていっても問題ない。
【００２９】
以上説明してきた通り、判別ツリー作成処理で判別ツリーが作成されるが、ここで注意するのは、ノードの数は指数関数的に増えていくので、ノードの節約を行わないと必要となるコンピュータリソースが不足することである。具体的には、想定する試料（ターゲット)中にノードに対応する塩基配列が存在しなくなれば、その子ノードの伸張はストップさせるという処理が必要となる。この処理によって、ある以上の深さを超えると、増加するノードは、定常的になる。例えば、判別すべき塩基配列のタイプの数をＴ、それぞれの全部分配列をＮとすると、最悪のケースでもＴ×Ｎ個のノードのみが増えることとなる。
【００３０】
［プローブの評価方法について（ステップＳ１０４）］
次に図１のプローブ評価処理（Ｓ１０４）を説明する。前述した判別ツリーの全てのノードは試料中に存在しうる部分塩基配列のどれかをコーディングしているので、判別ツリーの全てのノードはプローブの候補となりうる。しかし、ハイブリダイゼーション反応の条件よりプローブとして有効に働く部分塩基配列は限られてくる。このプローブとしての評価をするのが以下に説明するプローブ評価処理である。
【００３１】
まず、チップ上に載るプローブ群はハイブリダイゼーション強度が似ているほうがよい。よって、あらかじめプローブの塩基配列の数を決定する方法が考えられる。この場合、判別ツリーのある深さのノードのみがプローブ候補となる。もちろん、この塩基配列の数に幅を持たせてプローブ候補を絞ってもよい。いうまでもなく、この場合はある深さからある深さの範囲のノードがプローブ候補となる。
【００３２】
また、厳密に言えば、塩基配列の数が同じでもハイブリダイゼーション強度は異なる。これを計算する方法として、例えば、特開平８−３１７７９０「オリゴリボヌクレオチド／オリゴデオキシリボヌクレオチド相補塩基対構造の安定性解析方法」が提案されている。この方法によれば、塩基配列が決まればそのハイブリダイゼーション強度は推定できる。従って、各ノードに対応する塩基配列から各ノードのハイブリダイゼーション強度を推定し、この推定値を元に、判別ツリー上のノード群に絞込を行なう方法もある。但し、この場合、例えば文献“Predicting DNA duplex stability from the base sequence”(Proc. Natl. Acad. Sci. USA Vol. 83, pp. 3746-3750, June 1986 Biochmistry)に記載の方法で、判別ツリー上のノードの塩基配列より、塩基配列の融解温度を求めることになる。
【００３３】
本実施形態の非常に優れた点の１つは、判別ツリーが任意の長さの部分塩基配列をコーディングするので、ある範囲の長さやある範囲の融解温度を持ったプローブ候補を選ぶ処理を非常に高速に行なえる点にある。
【００３４】
さて、上述の絞込により、例えばある特定の長さ（例えば２０ベース）の部分配列を全て含むプローブを設計しても良いのであるが、全てのターゲットに含まれる部分塩基配列には特異性が全くなく、プローブとして意味がない。すなわち、このプローブを使って各ターゲットとハイブリダイゼーション反応を行なっても、何の情報も得ることはできない。よって、判別ツリーのノードによって特定される部分塩基配列の特異性を測る必要がある。この特異性の尺度として通常用いられるのがエントロピーである。
【００３５】
例えばターゲットとなる塩基配列タイプが全部でｔタイプあり、あるノードに含まれるタイプ毎の部分塩基配列の数がＮ１，Ｎ２，Ｎ３，．．．，Ｎｔであるとすると、そのノードのエントロピーは次の式で計算されることとなる。
【００３６】
【数１】
0004109873.tif 000002
【００３７】
もちろん、ここで説明したエントロピー以外の方法で塩基配列の特異性を測ってもよい。
【００３８】
さて、上記の式によって算出されたエントロピーを用いて最適なプローブを探すことになるが、その方法を以下に説明する。図９は、プローブとして典型的に最適なエントロピーの推移をグラフに示した図である。グラフの横軸はノードの深さを示し、この例では１５ベースの長さの塩基配列を調べた結果となっている。判別ツリーのノードを辿っていくと、通常、エントロピーは減少していく。例えば、図６Ａのノード６０１よりノード６０２の方が特異性が高く、エントロピーは減少する。図６Ａのノード６０１、６０２が図７、図８のような部分塩基配列のセットを保持するとすると、ノード６０２はDRB1*07011に固有の部分配列に相当し、エントロピーは０となる。
【００３９】
図９のグラフで８番目の塩基と１０番目の塩基の場所でエントロピーの減少が見られる。つまり、このノードのコーディングする部分塩基配列の８番目と１０番目があるタイプセットに特異的な塩基配列であることがグラフからわかる。図９にあるように、中央付近で急激なエントロピー減少があるノードがプローブとして最適となる。
【００４０】
本実施形態の非常に優れた点の１つは、このエントロピー減少の様子がノードを深さ方向に辿ることによって、自然と把握することができ、この解析によって、確実に中心位置での塩基特異性を持ったプローブを選べることにある。
【００４１】
以上の様なプローブの選択処理を実現するための具体的なプローブ評価関数としては、中心位置にエントロピー減少があると高いスコアになるようなものであればよい。例えば次に示す評価関数を使用することが可能である。
【００４２】
【数２】
評価スコア=(エントロピーの減少量)*(1/(中心位置からの距離+1))n
但し、この方法では、予めプローブの長さを想定しなければならい。なぜならば、中心位置はプローブの全長を決定して初めて決定されるからである。
【００４３】
予めプローブの長さを制限せずに評価を行なう方法として、判別ツリー作成ステップでも述べた中心から周辺へ塩基のバリエーションを見ていく方法が挙げられる。この方法で判別ツリーを作成した場合、判別ツリーの最初の方で、急激なエントロピー減少があった方がよく、それに対応した評価関数を用意することとなる。
【００４４】
例えば、次の評価関数を使用することが可能である。
【００４５】
【数３】
評価スコア=(エントロピーの減少量)*(1/(ルートノードからの距離))n
この方法だと任意の長さのプローブ設計に対応できることとなる。
【００４６】
以上のような評価関数を導入してプローブの評価を行ない、その評価値を各ノードに格納しておく。
【００４７】
図１０は本実施形態によるプローブの評価処理（ステップＳ１０４）とプローブセットの選択処理（ステップＳ１０５）の詳細を示すフローチャートである。ステップＳ２０１では、上述したように、塩基数（ノード深さ）や塩基配列に基づく融解温度により、判別ツリー上よりプローブ候補の絞込を行なう。そして、ステップＳ２０２において、上記絞込によって得られたプローブ候補のツリー上のパスにある各ノードについて、＜数１＞によりエントロピーを計算し、各ノードに格納する。そして、ステップＳ２０３において、プローブ候補のノードについて、＜数２＞或いは＜数３＞に示した評価関数を用いて評価値を算出し、これを当該ノードに格納する。
【００４８】
なお、図１０において、上記評価関数を用いたプローブの評価は、所定長の或いは所定長範囲の部分塩基配列を抽出して得られたプローブ候補、或いは所定範囲の融解温度を持ったプローブ候補について行うようにしているがこれに限らない。評価を判別ツリー上の全てのプローブについて行い、以下に説明するプローブセットの選択処理時に所定長範囲或いは所定融解温度範囲で絞込を行うようにしてもよい。
【００４９】
［プローブセットの選択処理について（ステップＳ１０５）］
次に図１のプローブセット選択処理（ステップＳ１０５）を説明する。図１０のステップＳ２０４〜Ｓ２０６はプローブセット選択処理の手順を示している。なお、必要があれば、プローブ候補を抽出して不適当なプローブを排除する。前述のプローブ評価ステップにおいて、ある範囲の長さ、または、ある範囲の融解温度のプローブのみを選んで得られたプローブ候補を評価したが、例えば、このプローブセット選択ステップにおいて、その判断を行うようにしてもよい。その場合、判別ツリー上の全てのプローブについて上述したステップＳ２０２とＳ２０３のプローブ評価を行なうことになる。
【００５０】
ステップＳ２０４において、プローブ候補となっている塩基配列の集合を複数のサブセットに分ける。これは冗長なプローブセットを選ぶことを防ぐためである。プローブの特異性を示す情報は、それぞれのターゲット塩基配列タイプに含まれるか含まれないかのビット列で表現できる。例えば、塩基配列タイプがｔ個存在するとすれば、ｔビットの情報でプローブの特異性が表現できるわけである。塩基配列タイプが含む場合ビットを１に設定し、含まない場合ビットを０に設定したとすると、例えば全てのタイプの塩基配列に含まれるプローブは全てが１であるｔビットの値で表現されることとなる。このｔビット列の異なるプローブを異なるグループとするとプローブの特異性を反映したサブセットが出来上がる。この処理によって同じサブセットに含まれる全てのプローブは、同タイプの特異性を持ち、情報的にはプローブとして同じ機能を持つこととなる。従って、各タイプのサブセットから上述したプローブ評価における評価値の高いプローブを選択してこれをマイクロアレイに用いるようにしてもよいが、本実施形態ではこのプローブセットを単位として更なる絞込を行なう。
【００５１】
上述のようにプローブ候補の塩基配列群を複数のサブセットに分けると、この複数のサブセットの中から必要でないサブセットを除くことが可能である。この処理を行うと、必要十分な独立プローブセットが出来上がる。以下に、具体的に例を示してこの処理を説明する。
【００５２】
ターゲットとなる塩基配列タイプがＡ、Ｂ、Ｃ、Ｄの４種類だったとする。いま、プローブ候補の塩基配列がその特異性に関して５つのサブグループに分けられるとすると、それぞれのサブグループの特異性が前述の４ビットの値で表せることとなる。ここで、例えば５つのサブグループが、1010,1100,0001,1110,0111というようなビット列でその特異性が表せたとする。すると実は、情報的には最初の２つのサブグループ（1010,1100）で４つのタイプを決定することは必要十分にできるのである。
【００５３】
この場合、最初の２つのサブグループの中から選んだプローブをα、βとすると、もし試料がＡタイプの塩基配列ならば、αβともにハイブリダイゼーション反応する。同様に、試料がＢタイプならαにはハイブリダイゼーション反応は起こらずβに対してのみ反応が起こる。試料がＣタイプならαにのみハイブリダイゼーション反応が起き、βには起きない。試料がＤタイプならα、βともにハイブリダイゼーション反応は起こらない。
【００５４】
もし実験が完全に制御できて、上記のハイブリダイゼーション反応の発生、非発生情報が非常に高い確率で得られたとすると、ターゲットとなるＡ、Ｂ、Ｃ、Ｄの４種類のタイプを判別する塩基プローブは、１番目と２番目のサブグループに含まれるプローブで必要十分となる。このような状況を情報的な観点でいうと、１番目と２番目のサブグループが独立で、その生成する情報空間に３番目、４番目、５番目のサブグループが従属しているという。よって、３番目、４番目、５番目のサブグループに属するプローブ候補は厳密に言えば必要なくなる。
【００５５】
なお、この必要十分な独立サブセットの組を選ぶ処理は非常に時間のかかる処理となる。また、実際は、最終的に選択するプローブセットは、ある程度冗長に選ぶ必要がある場合が多い。というのも、コンピュータの扱うデータと違って、生体物質の反応の結果として得られるデータは、実験の途中で多くのノイズが入る可能性が高く、必要最低限のプローブセットしか用意していないと、このノイズの多い場合に的確な実験結果を再現できないこととなるためである。よって、実際は上で述べたような必要十分な独立サブセットのみを選ぶ必要はない場合が多い。
【００５６】
ステップＳ２０５では以上のようにしてプローブセットを選択し、最終的に、対象となるプローブ候補のサブセットを得る。そして、選択されたそれぞれのサブセットの中で、プローブ評価ステップで得られた評価値の高いプローブを最終的なプローブセットとして選択する。なお、本実施形態では、上記手順で得られる塩基配列はターゲットの塩基配列から抽出された部分塩基配列そのものの配列順序を有するので、ステップＳ１０６で、選択された部分塩基配列を相補的な塩基配列に変換して出力することになる。例えば、ステップＳ１０５で部分塩基配列「ＧＡＧＣＧ」が選択されていれば、ステップＳ１０６では、これに対応するプローブとして「ＣＴＣＧＣ」が出力されることになる。もちろん、ステップＳ１０２において判別ツリーを作成する時点で、ターゲットの塩基配列と相補的な塩基配列を用いてもよい。この場合は、ステップＳ１０５で選択されるプローブの塩基配列がそのままステップＳ１０６で出力されることになる。
【００５７】
以上のようにして選択されたプローブを用いることにより、ＤＮＡマイクロアレイシステムに最適なオリゴヌクレオチドプローブの設計が実現できる。そして、このことにより、より確かな遺伝子発現情報や、個体識別情報が得られるという効果がある。
なお、上記実施形態において説明した評価については、種々の変形が考えられる。例えば、評価関数によって求めた値が所定値を超える場合に、即座にそのノードに対応する部分塩基配列をプローブとして決定するように処理を簡略化することも可能である。或いは＜数１＞で算出したエントロピーの変化量が所定値を超える部分を有する部分塩基配列をプローブとして決定するように処理を簡略化することも可能である。
【００５８】
＜第２の実施形態＞
第１の実施形態ではプローブ選択のために判別ツリーを生成し、これに基づいてプローブ評価を行なった。第２の実施形態では、判別ツリーの代わりにハッシュ表を用いる。以下、第２の実施形態の処理を説明する。なお、装置構成は第１の実施形態（図２）と同様であるので説明を省略する。
【００５９】
図１１は第２の実施形態によるプローブ設計方法の処理手順を説明するフローチャートである。ステップＳ１００１においてターゲットの塩基配列データを取得し、ステップＳ１００２において、取得された塩基配列データに含まれるある長さの範囲の部分塩基配列を分類する部分配列ハッシュ表を作成する。そして、ステップＳ１００３において、ステップＳ１００２で作成された部分塩基配列ハッシュ表を外部記憶装置２０１に格納する。
【００６０】
ステップＳ１００４は部分塩基配列ハッシュ表に存在するプローブ候補について、プローブとしてふさわしいかどうかを評価するプローブ評価処理を実行する。ステップＳ１００５では、更にプローブセットの選択を行ない、上記プローブ評価の結果を元に最終結果の最適プローブを選択する。ステップＳ１００６は、最終結果として得られた最適プローブである。上記処理において、入力はターゲットとなる塩基配列データであり、出力はマイクロアレイに搭載するのに適したプローブである。
【００６１】
［部分塩基配列ハッシュ表の生成について（Ｓ１００２、Ｓ１００３）］
第２の実施形態では、第１の実施形態で説明した判別ツリーを作成する替わりに、部分塩基配列ハッシュ表を作成する。これは判別ツリーのある深さのノードを作成するのに対応する。ハッシュ表とはソフトウェア技術として非常によく使われる方法で、大量のバリエーションのありうる情報列をコンパクトに保存する方法である。
【００６２】
例えば、２０ベースの長さの部分塩基配列がターゲットとなる塩基配列の中にどのように含まれるかを見る場合を考えてみる。２０ベースの塩基配列の取りうるバリエーションは、理論的には４20（＝４０ビット）となる。これだけの大きな空間を用意することは事実上不可能である。一方、試料中に存在しうる塩基配列のバリエーションは、この４０ビットの空間と比べると、極めて小さい。よって、２０ベースの部分塩基配列を４０ビットの値にコーディングしておいて、それを例えば十分大きな素数でわることによって、事実上重複のないコンパクトな空間に収めることができる。この方法をハッシュ表を用いたデータ格納整理方法と呼んでいる。なお、仮に実際は異なる塩基配列にもかかわらず、重複が生じた場合は、２重にハッシュ表を用意するとか、ハッシュ表の次の番地に格納するというような回避方法がある。
【００６３】
ハッシュ表を用いて、ターゲットの塩基配列の中の特定の長さの全部分配列が格納整理できると、それが含まれる塩基配列の種類と位置も自動的に特定できる。これは、図７や図８の左のカラムがハッシュ表のそれぞれのエントリーに登録されている状態と同じである。このことによって、実施形態１の判別ツリーのある深さのノード情報が全て得られたことになる。
【００６４】
この操作を複数の長さの部分塩基配列に対して実施することによって、第１の実施形態である深さからある深さまでの全てのノード展開ができたことと同じになる。第１の実施形態で述べたとおり、ある範囲の長さのプローブやある範囲の融解温度のプローブを設計する場合は、これでプローブ候補が全て得られたこととなる。
【００６５】
［プローブ候補の評価について（Ｓ１００４）］
次に、この全てのプローブ候補から、プローブとして適切か否かを判定するためのプローブ評価について説明する。第２の実施形態が第１の実施形態と決定的に異なるところは、このプローブ評価処理である。第２の実施形態は第１の実施形態より若干性能が劣る。つまり、第１の実施形態の方がより高性能なプローブを設計できる場合が多い。しかしながら、第２の実施形態は、処理の簡易さ、迅速さにおいて第１の実施形態を上回るものであり、その利用価値は高いものである。
【００６６】
第１の実施形態でも述べたとおり、プローブの塩基配列の特異性は中央部分に存在するほうがプローブとしては質が高い。よって、その中央部の特異性を調べたいのであるが、第１の実施形態と違って、塩基配列中の部分的な情報が全くないので非常に困難である。そこで、第２の実施形態では、ターゲットの塩基配列中に存在する特異点（特異塩基位置という）を予め求める。この様子を示したのが図１２である。図１２の表の一番下に表示したのが特異塩基位置である。＊のマークのある場所が、ターゲットの塩基配列の全てが共通とはなっていない箇所、すなわちあるタイプ間に特異性がある場所（特異塩基位置）である。そして、この位置を中央に含むプローブ候補が良質なプローブとなる。この特異塩基位置のプローブ候補中での場所を評価する関数を使って、プローブ候補の質を評価する。
【００６７】
具体的には、例えば、
【数４】
評価スコア=Σ{1/(中心位置からの特異塩基位置の距離+1)}n}
というような評価関数を利用して判定することができる。
【００６８】
また、上述の例では、特異塩基位置の判定において、特異性があるかないかの２値で評価していたが、実施形態１に示したようなエントロピーのような尺度を使用して、連続値で評価しても問題ない。例えば、図１２の＊のマークのあるところは、特異塩基配列と呼ばれるが、場所によってはほとんどがＧで、一種類のタイプだけＡである場所（エントロピーの減少量が低い）と、半分のタイプがＧで、残りの半分のタイプがＡである場所（エントロピーの減少量が高い）が混じっている。上述したように、エントロピーの減少量が多い場所が中央にあったほうが判別に有利なので、＜数４＞に対してそのような重み付けを行なうこと、例えば、特異塩基配列の位置におけるエントロピーの減少量を＜数４＞に積算するというような重み付けをすることで、より正確なプローブ候補の質評価が可能となる。
【００６９】
［プローブセット選択処理（Ｓ１００５）］
ステップＳ１００５で実行されるプローブセット選択処理は、第１の実施形態（ステップＳ１０６）と同様の処理を行うものであり、説明を省略する。
【００７０】
以上のように第２の実施形態によれば、ハッシュ表を用いてプローブ候補を選択するので、処理の簡素化、高速化が達成される。
【００７１】
以上説明したように、第１或いは第２の実施形態のプローブ設計処理に従って得られた塩基プローブを基板上に配置することで、ターゲットのタイピングに好適な塩基プローブを含むＤＮＡマイクロアレイが得られることになる。なお、塩基プローブが決定されれば、周知のＤＮＡマイクロアレイの製造方法を用いて、所定のターゲットに対して好適なＤＮＡマイクロアレイを製造することができる。また、このようなＤＮＡマイクロアレイを用いる遺伝子検査装置を提供することで、ターゲットのタイピングを良好に実現することができる。なお、ＤＮＡマイクロアレイの製造方法と構造、およびこれを用いた遺伝子検査装置そのものは周知であるので、その詳細な説明を省略する。
【００７２】
＜第３実施形態＞
上述の各実施形態においては、ＰＣＲなどの手法を用いて、調べたいＤＮＡ、あるいは、ＲＮＡ（ｃＤＮＡ）の塩基配列を増幅して検体溶液を作成するのが、最も一般的であった。この場合、検体中に存在する塩基配列は、調べたい塩基配列をその部分配列として包含している可能性が非常に高かった。即ち、第１及び第２実施形態は、ＰＣＲ等の実験手法を用いた検体調製をした場合に、特に好適である。
【００７３】
ところが、例えば細胞中に存在するＲＮＡをまとめて抽出し、その中の一部のＲＮＡの存在量を測定するといった目的の場合、作成したＤＮＡマイクロアレイと反応する可能性のある塩基配列は、細胞中に存在する全てのＲＮＡとなる。そして、これらのＲＮＡ、あるいはそれから作成されたｃＤＮＡの塩基配列の中に調べたい塩基配列が存在する確率は極めて低い。
【００７４】
これを図を使って説明すると、図１３のようになる。第１及び第２実施形態では、「検体中に存在する可能性のある塩基配列、あるいはその部分塩基配列」と、「調べたい塩基配列、あるいはその部分塩基配列」は、ほとんど同じ集合であった。これを第１及び第２実施形態においては、ターゲット塩基配列と呼んでいた。例えば、ＭＨＣなどの人体特異性を示す遺伝子ＤＮＡの解析などは、この場合に当てはまる。
【００７５】
これに対して、上述のようなＲＮＡの発現量を測定する実験の場合、検体中に存在する塩基配列は、調べたい塩基配列に対して、圧倒的に種類が多いのが一般的である。よって、この場合は、検体中に存在する可能性のある全ての塩基配列が特異性の対象となるターゲット塩基配列となる。そして、調べたい塩基配列は、そのターゲット塩基配列の一部となる。
【００７６】
図１３で説明すると、第１実施形態の判別ツリー、または、第２実施形態のハッシュテーブルを作成するに際して、使用される塩基配列は、図１３の「検体中に存在する可能性のある塩基配列の全部分塩基配列」（１３０１）となる。そして、調べたい塩基配列の部分配列の特異性は、全体であるターゲット塩基配列との特異性に基づいて導かれ、プローブ候補が決定されることとなる。
【００７７】
以下、図１４、図１５を用いて、第３実施形態によるプローブ設計方法を具体的に説明する。
【００７８】
まず、第１実施形態と同様の判別ツリーを用いたプローブ設計手順について図１４を用いて説明する。ステップＳ１４０１において、検体中に存在する可能性のある全塩基配列であるターゲット塩基配列１４１１を取得する。この中には、調べたい塩基配列データ１４１０が含まれる。ターゲット塩基配列１４１１は、図１３に示されるように、調べたい塩基配列１３０２を含む塩基配列集合１３０１である。
【００７９】
ステップＳ１４０２では、ステップＳ１４０１で得られたターゲット塩基配列データ１４１１より得られる複数の部分塩基配列を分類する判別ツリー１４１２を生成する。判別ツリーの生成については第１実施形態（ステップＳ１０２、Ｓ１０３）で説明したとおりである。
【００８０】
次に、ステップＳ１４０３において、判別ツリー１４１２の中に存在するプローブ候補についてプローブとしてどの程度ふさわしいかを評価する。評価の方法は第１の実施形態（ステップＳ１０４）と同様である。但し、プローブとして選ばれる部分塩基配列は、調べたい塩基配列データ１４１０の部分配列の中から選ばれ、その特異性は、ターゲット塩基配列データ１４１１をもとに求められることとなる。
【００８１】
続いて、ステップＳ１４０４において、ステップＳ１４０３におけるプローブの評価結果に基づいて最終結果の最適プローブセットを選択する。この処理は第１の実施形態のステップＳ１０５と同様である。そして、ステップＳ１４０５において、ステップＳ１４０４で選択されたプローブが最適プローブとして出力される。本実施形態の入力は、調べたい塩基配列データであり、出力は調べたい塩基配列に含まれる最適プローブセット、またはその相補鎖となる。
【００８２】
次に、第２実施形態と同様のハッシュテーブルを用いたプローブ設計手順について図１５を参照して説明する。
【００８３】
ステップＳ１５０１において、検体中に存在する可能性のある全塩基配列であるターゲット塩基配列１５１１を取得する。この中には、調べたい塩基配列データ１５１０が含まれる。ターゲット塩基配列１５１１は、図１３に示されるように、調べたい塩基配列１３０２を含む塩基配列集合１３０１である。
【００８４】
ステップＳ１５０２では、ステップＳ１５０１で得られたターゲット塩基配列データ１５１１より得られる複数の部分塩基配列を分類する部分塩基配列ハッシュ表１５１２を生成する。部分塩基配列ハッシュ表１５１２の生成については第２実施形態（ステップＳ１００２、Ｓ１００３）で説明したとおりである。
【００８５】
次に、ステップＳ１５０３において、部分塩基配列ハッシュ表１５１２の中に存在するプローブ候補についてプローブとしてどの程度ふさわしいかを評価する。評価の方法は第２の実施形態（ステップＳ１００４）と同様である。但し、プローブとして選ばれる部分塩基配列は、調べたい塩基配列データ１５１０の部分配列の中から選ばれ、その特異性は、ターゲット塩基配列データ１５１１をもとに求められることとなる。
【００８６】
続いて、ステップＳ１５０４において、ステップＳ１５０３におけるプローブの評価結果に基づいて最終結果の最適プローブセットを選択する。この処理は第２の実施形態のステップＳ１００５と同様である。そして、ステップＳ１５０５において、ステップＳ１５０４で選択されたプローブが最適プローブとして出力される。本実施形態の入力は、調べたい塩基配列データであり、出力は調べたい塩基配列に含まれる最適プローブセット、またはその相補鎖となる。
【００８７】
以上のように、本実施形態によれば、例えば細胞中に存在するＲＮＡをまとめて抽出し、その中の一部のＲＮＡの存在量を測定するといった用途に対して、適切なマイクロアレイを設計することが可能となる。
【００８８】
また、本発明の目的は、前述した実施形態の機能を実現するソフトウェアのプログラムコードを記録した記憶媒体を、システムあるいは装置に供給し、そのシステムあるいは装置のコンピュータ（またはＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読出し実行することによっても、達成されることは言うまでもない。
【００８９】
この場合、記憶媒体から読出されたプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコードを記憶した記憶媒体は本発明を構成することになる。
【００９０】
プログラムコードを供給するための記憶媒体としては、例えば、フロッピディスク，ハードディスク，光ディスク，光磁気ディスク，ＣＤ−ＲＯＭ，ＣＤ−Ｒ，磁気テープ，不揮発性のメモリカード，ＲＯＭなどを用いることができる。
【００９１】
また、コンピュータが読出したプログラムコードを実行することにより、前述した実施形態の機能が実現されるだけでなく、そのプログラムコードの指示に基づき、コンピュータ上で稼働しているＯＳ（オペレーティングシステム）などが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【００９２】
さらに、記憶媒体から読出されたプログラムコードが、コンピュータに挿入された機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに書込まれた後、そのプログラムコードの指示に基づき、その機能拡張ボードや機能拡張ユニットに備わるＣＰＵなどが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【００９３】
【発明の効果】
以上説明したように、本発明によれば、解析対象となるターゲットの塩基配列に応じて、解析に適切なプローブを自動的に選択することが可能となり、プローブ設計の効果的な支援が実現される。
【００９４】
【配列表】
0004109873.tif 000003
0004109873.tif 000004
0004109873.tif 000005
【図面の簡単な説明】
【図１】第１の実施形態によるプローブ設計方法の処理手順を説明するフローチャートである。
【図２】本実施形態によるプローブ設計方法が適用される情報処理装置の構成を示すブロック図である。
【図３】ステップＳ１０１で取得する塩基配列データの例を示す図である。
【図４】ＤＮＡマイクロアレイ上のハイブリダイゼーションの様子を示した図である。
【図５】ターゲットとなる塩基配列から得られる部分塩基配列（プローブ）を示す図である。
【図６Ａ】第１の実施形態における判別ツリーを説明する図である。
【図６Ｂ】判別ツリーにおける各ノードの内容例を示す図である。
【図７】判別ツリー上における、あるノードに登録された部分塩基配列の様子を示す図である。
【図８】判別ツリー上における、あるノードに登録された部分塩基配列の様子を示す図である。
【図９】プローブとして典型的に最適なエントロピーの推移をグラフに示した図である。
【図１０】第１の実施形態によるプローブの評価処理（ステップＳ１０４）とプローブセットの選択処理（ステップＳ１０５）の詳細を示すフローチャートである。
【図１１】第２の実施形態によるプローブ設計方法の処理手順を説明するフローチャートである。
【図１２】第２の実施形態による、ターゲット中の特異塩基位置を説明する図である。
【図１３】第３の実施形態によるプローブ設計方法が適用される塩基配列データの状態を説明する図である。
【図１４】第３の実施形態による、判別ツリーを用いたプローブ設計方法を説明するフローチャートである。
【図１５】第３の実施形態による、ハッシュ表を用いたプローブ設計方法を説明するフローチャートである。
【図１】
0004109873.tif 000006
【図２】
0004109873.tif 000007
【図３】
0004109873.tif 000008
【図４】
0004109873.tif 000009
【図５】
0004109873.tif 000010
【図６Ａ】
0004109873.tif 000011
【図６Ｂ】
0004109873.tif 000012
【図７】
0004109873.tif 000013
【図８】
0004109873.tif 000014
【図９】
0004109873.tif 000015
【図１０】
0004109873.tif 000016
【図１１】
0004109873.tif 000017
【図１２】
0004109873.tif 000018
【図１３】
0004109873.tif 000019
【図１４】
0004109873.tif 000020
【図１５】
0004109873.tif 000021