"生物ゲノム中の標的遺伝子を含む遺伝子クラスタ及び/または該遺伝子クラスタ中の標的遺伝子を探索する方法であって、生物細胞の生理状態変化を生じる条件とコントロール条件下において生じたゲノム遺伝子の発現量変動比を、ゲノムDNA上に配列する複数の遺伝子により構成される仮想の遺伝子クラスタ単位の発現量変動比として合算することにより、仮想の遺伝子クラスタ単位毎にスコアリングし、得られたスコアに基づき、上記生理状態変化の原因遺伝子である標的遺伝子を含む遺伝子クラスタ及び/または該遺伝子クラスタ中の標的遺伝子を探索する方法。生物細胞の生理状態変化を生じる条件とコントロール条件下とを一の対比条件セットとして、該対比条件セットが一種以上設定されていることを特徴とする請求項1に記載の方法。生理状態変化を生じる条件とコントロール条件が、少なくとも代謝産物の産生誘導条件下と非誘導条件下あるいは代謝産物の産生抑制条件下と非抑制条件下との対比条件セットを含むことを特徴とする請求項1または2に記載の方法。代謝物産生に関与する遺伝子が2次代謝物産生に関与する遺伝子であることを特徴とする請求項3に記載の方法。仮想の各遺伝子クラスタは、ゲノムDNA上に連続する遺伝子を2個から遺伝子数を一つずつ増やして、想定される遺伝子クラスタに含まれる最大限のゲノム遺伝子数になるまで抽出し、かつ該抽出において、抽出する遺伝子の各個数毎に、直鎖状DNAからなるゲノムの場合には該DNAのいずれかの末端から、環状DNAからなるゲノムの場合には任意の遺伝子を起点として順にゲノムDNA上に配列する遺伝子を一つずつずらしながら抽出された各遺伝子群からなることを特徴とする、上記請求項1~4のいずれかに記載の方法。スコアリングされる仮想の各遺伝子クラスタの集合体が、ゲノムDNA上に連続する遺伝子を2個から遺伝子数を一つずつ増やし想定される遺伝子クラスタに含まれる最大限のゲノム遺伝子数になるまで抽出し、かつ該抽出において、抽出する遺伝子の各個数毎に、直鎖状DNAからなるゲノムの場合には該DNAのいずれかの末端から、環状DNAからなるゲノムの場合には任意の遺伝子を起点として順にゲノムDNA上に配列する遺伝子を一つずつずらしながら抽出された各遺伝子群からなる仮想の各遺伝子クラスタの集合からなり、ゲノム上に存在する遺伝子クラスタの全てが仮想の遺伝子クラスタの集合体中に含まれるように構成されていることを特徴とする請求項1~5のいずれかに記載の方法。仮想の各遺伝子クラスタのスコアリングが以下の計算式aによりなされることを特徴とする、請求項1~6のいずれかに記載の方法。計算式aゲノムDNA上に配列する遺伝子が、標的とする遺伝子機能を有すると推定される場合、あるいは標的とする遺伝子機能を有する可能性が低いか若しくはその可能性がないと推定される場合において、当該ゲノムDNA上に配列する遺伝子については、以下の重み付け計算が適用されることを特徴とする、請求項7の方法。ゲノムDNA上に配列する遺伝子が、標的とする遺伝子機能を有すると推定される場合において、標的とする遺伝子機能を有すると推定された遺伝子を含む仮想の遺伝子クラスタを選出し、選出された仮想の遺伝子クラスタについて、スコアリングすることを特徴とする、請求項7に記載の方法。仮想の遺伝子クラスタが、ゲノムにおいて近傍に存在することを条件として、以下の1~3の内の1以上の遺伝子のみから、あるいは該遺伝子を少なくとも含む1以上の遺伝子から構築されることを特徴とする、上記4に記載の方法。12次代謝物産生に関与していると想定される酵素種に属する酵素遺伝子。2トランスポーター遺伝子3転写因子をコードする遺伝子仮想の各遺伝子クラスタのスコアリングが以下の計算式aによりなされることを特徴とする、上記10に記載の方法。計算式a仮想の遺伝子クラスタ全体のスコアの分布から乖離して存在するスコアを有する仮想の遺伝子クラスタを、標的の遺伝子クラスタ候補として選定することを特徴とする、請求項1~11のいずれかに記載の方法。仮想の遺伝子クラスタ全体のスコアの分布からの乖離の程度を示す判定値Iχを、以下の計算式bにより算出し、算出された該判定値Iχに基づき仮想の遺伝子クラスタを標的の遺伝子クラスタ候補として選定することを特徴とする、請求項12に記載の方法。計算式b仮想の遺伝子クラスタ全体のスコアの分布からの乖離の程度を示す判定値IIυを、以下の計算式cにより算出し、算出された判定値IIを求め、該判定値に基づき、予め、標的とする遺伝子クラスタがゲノム中に存在するか否かあるいは、標的クラスタが存在する場合のその遺伝子サイズを予測することを特徴とする、上記方法。計算式a計算式e遺伝子数がk個のときのε値εkと、その前後数のときのε値εk-1、εk1が、以下の関係にあるとき、標的とする遺伝子クラスタがゲノム中に存在すると判定し、標的遺伝子クラスタに含まれる遺伝子数をk個と予想することを特徴とする、請求項16に記載の方法。生物ゲノム中の標的遺伝子を含む遺伝子クラスタ及び/または該遺伝子クラスタ中の標的遺伝子を探索する装置であって、a生物細胞の生理状態変化を生じる条件とコントロール条件下におけるゲノムDNA上に配列する各遺伝子の発現量データに基づき算出された上記2つの条件下における上記各遺伝子の発現量変動比を記憶する手段、bゲノムDNA上に配列する複数の遺伝子を組み合わせて仮想の遺伝子クラスタを構築する手段、c該算出され、記憶されたゲノムDNA上に配列する各遺伝子の発現量変動比を複数の遺伝子により構築された上記仮想の遺伝子クラスタ単位の発現量変動比として合算し、仮想の遺伝子クラスタ単位毎にスコアリングし、仮想の各遺伝子クラスタのスコアを記憶する手段、及びd得られたスコアに基づき上記生理状態変化の原因遺伝子である標的遺伝子を含む遺伝子クラスタを選定する手段を有するか、あるいはさらにe選定された遺伝子クラスタ中に含まれる遺伝子を表示する手段を有することを特徴とする、上記装置。発現量データが、遺伝子発現量測定用DNAマイクロアレイによる蛍光強度情報であることを特徴とする請求項18に記載の装置。蛍光強度情報が、蛍光強度を読み取り、数値化する手段を有する蛍光強度読み取り装置により出力される数値データであることを特徴とする、請求項19に記載の装置。生物細胞の生理状態変化を生じる条件とコントロール条件とを1の対比条件セットとして1以上設定されている場合において、各対比条件セットに含まれる条件毎に各遺伝子の発現量データが入力され、各対比条件セットにおける同一遺伝子の発現量変動比が算出されることを特徴とする、請求項18~20のいずれかに記載の装置。標的遺伝子が代謝物産生に関与する遺伝子であることを特徴とする、請求項18~21のいずれかに記載の装置。代謝物産生に関与する遺伝子が2次代謝物産生に関与する遺伝子であることを特徴とする、請求項22に記載の装置。設定される対比条件セットが、少なくとも代謝産物の産生誘導条件下と非誘導条件下あるいは代謝産物の産生抑制条件下と非抑制条件下との対比条件セットを含むことを特徴とする請求項22に記載の装置。代謝産物が2次代謝産物であることを特徴とする、請求項24に記載の装置。仮想の各遺伝子クラスタの構築手段が、ゲノムDNA上に連続する遺伝子を2個から遺伝子数を1ずつ増やして、想定される遺伝子クラスタに含まれる最大限のゲノム遺伝子数になるまで抽出し、かつ該抽出において、抽出する遺伝子の各個数毎に、直鎖状DNAからなるゲノムの場合には該DNAのいずれかの末端から、環状DNAからなるゲノムの場合には任意の遺伝子を起点として順にゲノムDNA上に配列する遺伝子を一つずつずらしながら抽出した各遺伝子群により構築する手段であることを特徴とする、請求項18~25のいずれかに記載の装置。仮想の各遺伝子クラスタのスコアリングが以下の計算式aによりなされることを特徴とする、請求項18~26のいずれかに記載の装置。計算式aゲノムDNA上に配列する各遺伝子中の特定の遺伝子を選定するためのアノテーション付与手段を有し、上記遺伝子クラスタのスコアリングにおいて、付与されたアノテーションに基づき選定された遺伝子についての発現量変動比計算を以下の重み付け計算式により行うことを特徴とする、請求項27に記載の装置。アノテーション付与手段が、それぞれ遺伝子機能の種類毎に異なるアノテーションを付与する手段であることを特徴とする請求項28に記載の装置。アノテーションに基づき選定される遺伝子が、1~3のうちの1以上の遺伝子であることを特徴とする、請求項29に記載の装置12次代謝物産生に関与していると想定される酵素種に属する酵素遺伝子。2トランスポーター遺伝子3転写因子をコードする遺伝子上記請求項28~30のいずれかに記載のアノテーション付与手段と、構築された仮想の遺伝子クラスタから、アノテーションに基づき選出された遺伝子を含む仮想の遺伝子クラスタを選出する手段を有し、選出された仮想の遺伝子クラスタについてスコアリングすることを特徴とする、請求項27に記載の装置。ゲノムDNA上に配列する各遺伝子中の特定遺伝子を選定するためのアノテーション付与手段を有し、ゲノムDNA上において近傍に位置することを条件として、アノテーションに基づき選定された遺伝子のみから、あるいは該遺伝子を少なくとも含む1以上の遺伝子から仮想の遺伝子クラスタを構築する手段を有することを特徴とする、請求項18~25のいずれかに記載の装置。請求項32に記載のアノテーション付与手段が、それぞれ遺伝子機能の種類に応じたアノテーションを付与する手段であることを特徴とする請求項32に記載の装置。アノテーション付与に基づき選定される遺伝子が、1~3のうちの1以上の遺伝子であることを特徴とする、請求項33に記載の装置12次代謝物産生に関与していると想定される酵素種に属する酵素遺伝子。2トランスポーター遺伝子3転写因子をコードする遺伝子仮想の各遺伝子クラスタのスコアリングが以下の計算式aによりなされることを特徴とする、請求項32~34のいずれかに記載の装置。計算式a仮想の遺伝子クラスタ全体のスコアの分布から乖離して存在するスコアを有する仮想の遺伝子クラスタを、標的の遺伝子クラスタ候補として選定する手段を有することを特徴とする、請求項18~35のいずれかに記載の装置。標的の遺伝子クラスタ候補として選定する手段として、仮想の遺伝子クラスタ全体のスコアの分布からの乖離の程度を示す判定値Iχを、以下の計算式bにより算出するプログラムが格納されていることを特徴とする、請求項36に記載の装置。計算式b標的の遺伝子クラスタ候補として選定する手段として、遺伝子クラスタ全体のスコアの分布からの乖離の程度を示す判定値IIυを、以下の計算式cにより算出するプログラムが格納されていることを特徴とする、請求項37に記載の装置。計算式cさらに、以下の計算式dの算出結果に基づき、bが100未満の仮想のクラスタを少なくとも除外し、標的の遺伝子クラスタ候補をさらに絞り込むプログラムが格納されていることを特徴とする、請求項37または38に記載の装置。計算式da生物細胞の生理状態変化を生じる条件とコントロール条件下において生じたゲノムDNA上に配列する各遺伝子の発現量を入力する手段、b入力された上記2つの条件下における同一遺伝子の発現量の比を算出する発現量変動比算出手段、c該算出されたゲノムDNA上に配列する各遺伝子の発現量変動比を複数の遺伝子により構築された上記仮想の遺伝子クラスタ単位の発現量変動比として合算し、仮想の遺伝子クラスタ単位毎にスコアリングする手段、及びd得られた仮想の遺伝子クラスタのスコアから遺伝子クラスタに含まれる遺伝子数単位毎の遺伝子クラスタ分布判定値εを算出する手段を有し、該遺伝子クラスタ分布判定値εから、標的とする遺伝子クラスタがゲノム中に存在する否かあるいは、標的遺伝子クラスタが存在する場合の遺伝子サイズを予測する装置であって、仮想の遺伝子クラスタの構築手段が、ゲノムDNA上に連続する遺伝子を2個から遺伝子数を一つずつ増やし想定される遺伝子クラスタに含まれる最大限のゲノム遺伝子数になるまで抽出し、かつ該抽出において抽出する遺伝子の各個数毎に、直鎖状DNAからなるゲノムの場合には該DNAのいずれかの末端から、あるいは環状DNAからなるゲノムの場合には任意の遺伝子を起点として順にゲノムDNA上に配列する遺伝子を一つずつずらしながら抽出された各遺伝子群を仮想の各遺伝子クラスタとする手段であり、上記仮想の遺伝子クラスタ単位のスコアリング手段は以下の計算式aによる演算手段からなるとともに、上記遺伝子クラスタ分布判定値εの算出手段が、以下の計算式eによるものであることを特徴とする、上記装置。計算式a計算式e遺伝子数がk個のときの遺伝子クラスタ分布判定値ε値εkと、その前後数のときの同ε値εk-1、εk1が、以下の関係にあるとき、標的とする遺伝子クラスタがゲノム中に存在すると判定し、標的遺伝子クラスタに含まれる遺伝子数をk個とする予想値を出力することを特徴とする、請求項40に記載の装置。請求項26に記載の仮想の遺伝子クラスタの構築手段を実行するプログラムであって、ゲノム遺伝子の位置情報に基づき、以下の1または2の手段を実行することを特徴とする、仮想の遺伝子クラスタ構築プログラム。1ゲノム遺伝子が直鎖状ゲノムの場合、a.ゲノムDNAの一方の末端に位置する遺伝子を起点として、他方の末端方向に、順次、ゲノムDNA上に連続する遺伝子を同一方向に2個から一つずつ増やして想定される遺伝子クラスタに含まれる遺伝子数の最大限になるまで組み合わせ、起点とした遺伝子を含み、かつ遺伝子の個数の異なる複数の遺伝子群を構成する手段。b.起点を、順次、他方の末端方向に一遺伝子ずつずらせながら、上記a.と同様の処理を行い、新たな起点遺伝子を含みかつ遺伝子の個数が異なる複数の遺伝子群を構成し、a.の遺伝子群と併せて、複数の遺伝子を組み合わせた遺伝子群からなる仮想の遺伝子クラスタを構築する手段。2ゲノム遺伝子が環状の場合、ゲノムDNA上の任意の遺伝子を起点として、上記1a.及びb.と同様の処理を順次行い、最初に起点とした遺伝子が起点となる時点で処理を終了する手段。請求項42のプログラムにより構築された仮想の遺伝子クラスタについて、以下の計算式aによるスコアリングを実行することを特徴とする、仮想の遺伝子クラスタのスコアリングプログラム。計算式a上記遺伝子クラスタのスコアリングにおいて、付与されたアノテーションに基づきゲノム遺伝子を選定し、選定された遺伝子についての発現量変動比計算を以下の重み付け計算式により行うことを特徴とする、請求項43に記載のスコアリングプログラム。上記遺伝子クラスタのスコアリングにおいて、付与されたアノテーションに基づきゲノム遺伝子を選定し、構築された遺伝子クラスタの中から、該選定されたゲノム遺伝子を含む仮想の遺伝子クラスタを選出し、選出された仮想の遺伝子クラスタについてスコアリングを実行することを特徴とする、請求項43に記載のスコアリングプログラム。上記請求項32に記載の仮想の遺伝子クラスタの構築手段を実行するプログラムであって、ゲノムDNA上において近傍に位置することを条件として、アノテーションに基づき選定された遺伝子のみから、あるいは該遺伝子を少なくとも含む1以上の遺伝子から仮想の遺伝子クラスタを構築することを特徴とする、仮想の遺伝子クラタの構築プログラム。請求項46のプログラムにより構築された仮想の遺伝子クラスタについて、以下の計算式aによるスコアリングを実行することを特徴とする、仮想の遺伝子クラスタのスコアリングプログラム。計算式a請求項43~45又は47のいずれかに記載のスコアリングプログラムにより算出された各仮想の遺伝子クラスタのスコアについて、仮想の遺伝子クラスタ全体のスコアの分布からの乖離の程度を算出するプログラムであって、以下の計算式bにより、判定値Iχを算出することを特徴とする、上記プログラム。計算式b請求項43~45又は47のいずれかに記載のスコアリングプログラムにより算出された各仮想の遺伝子クラスタのスコアについて、仮想の遺伝子クラスタ全体のスコア分布からの乖離の程度を算出するプログラムであって、以下の計算式cにより判定値IIυの算出を実行する、上記プログラム。計算式c生物細胞の生理状態変化を生じる条件とコントロール条件下とにおけるゲノムDNA上に配列する各遺伝子の発現量変動比を複数の遺伝子により構築された上記仮想の遺伝子クラスタ単位の発現量変動比として合算し、仮想の遺伝子クラスタ単位毎にスコアリングする手段、及び得られた仮想の遺伝子クラスタのスコアから遺伝子クラスタに含まれる遺伝子数単位毎の遺伝子クラスタ分布判定値εを算出し、該遺伝子クラスタ分布判定値εから、標的とする遺伝子クラスタがゲノム中に存在する否かあるいは、標的遺伝子クラスタが存在する場合の遺伝子サイズを予測する手段に用いるプログラムであって、少なくとも以下A~Cの手段を実行するプログラム。Aゲノム遺伝子の位置情報に基づき、以下の1または2の手段により仮想の遺伝子クラスタを構築する手段、1ゲノム遺伝子が直鎖状の場合、a.ゲノムDNAの一方の末端に位置する遺伝子を起点として、他方の末端方向に、順次、ゲノムDNA上に連続する遺伝子を同一方向に2個から一つずつ増やして想定される遺伝子クラスタに含まれる遺伝子数の最大限になるまで組み合わせ、起点とした遺伝子を含み、かつ遺伝子の個数の異なる複数の遺伝子群を構成する手段。b.起点を、順次、他方の末端方向に一遺伝子ずつずらせながら、上記a.と同様の処理を行い、新たな起点遺伝子を含みかつ遺伝子の個数が異なる複数の遺伝子群を構成し、a.の遺伝子群と併せて、複数の遺伝子の組み合わせた遺伝子群からなる仮想の遺伝子クラスタを構築する手段。2ゲノム遺伝子が環状の場合、ゲノムDNA上の任意の遺伝子を起点として、上記1a.及びb.と同様の処理を順次行い、最初に起点とした遺伝子が起点となる時点で処理を終了する手段。B上記Aの手段により構築された仮想の遺伝子クラスタについて、以下の計算式aにより仮想の遺伝子クラスタ単位毎にスコアリングする手段。計算式aC上記Bの手段により得られた仮想の遺伝子クラスタのスコアから、以下の計算式eにより仮想の遺伝子クラスタに含まれる遺伝子数単位毎の遺伝子クラスタ分布判定値εを算出する手段。計算式e遺伝子数がk個のときの遺伝子クラスタ分布判定値ε値εkと、その前後数のときの同ε値εk-1、εk1が、以下の関係にあるとき、標的とする遺伝子クラスタがゲノム中に存在すると判定し、標的遺伝子クラスタに含まれる遺伝子数をk個とする予想値を出力することを特徴とする、請求項50に記載のプログラム。本発明は、上記従来技術における、代謝産物の産生に関与する遺伝子等の有用遺伝子の探索において、上記従来技術にみられる研究者の知識や経験等に大きく依存することなく、また、遺伝子破壊実験を逐次行わずとも、論理的、システマチックに極めて短時間で効率的に、有用遺伝子の探索、同定を行う方法、およびそのための装置を提供することにあり、これにより、今後増大するゲノム情報を活用して、新たな有用遺伝子の探索を加速させ、ゲノム中の遺伝子配列と有用遺伝子の対応関係について、詳細かつ膨大な情報を収集してデータベースなどを構築することを可能にし、多数の有用な遺伝子産物の発見に資することを課題とするものである。本発明者は、上記課題を解決するため鋭意研究の結果、従来におけるマイクロアレイのゲノム遺伝子の発現誘導あるいは破壊実験等による有用遺伝子探索にみられるような、ゲノム中の個々の遺伝子についての発現変動情報から直接標的遺伝子を絞りこむのではなく、マイクロアレイ等によるゲノム上の各遺伝子の発現変動情報を、複数の遺伝子により構成される仮想の遺伝子クラスタ単位の発現変動情報として合算して、仮想の各遺伝子クラスタをスコアリングし、この仮想の遺伝子クラスタの中から、有用遺伝子を含む遺伝子クラスタ及び該クラスタに含まれる有用遺伝子を見いだすことにより、上記従来の有用遺伝子探索法に比べ、遙かに正確かつ効率的に有用遺伝子の探索、同定が可能となることを見いだし、本発明を完成するに至った。すなわち、本発明は、以下のとおりである。1本発明は、以下の有用遺伝子の探索、同定を行う方法を提供する。1生物ゲノム中の標的遺伝子を含む遺伝子クラスタ及び/または該遺伝子クラスタ中の標的遺伝子を探索する方法であって、生物細胞の生理状態変化を生じる条件とコントロール条件下において生じたゲノム遺伝子の発現量変動比を、ゲノムDNA上に配列する複数の遺伝子により構成される仮想の遺伝子クラスタ単位の発現量変動比として合算することにより、仮想の遺伝子クラスタ単位毎にスコアリングし、得られたスコアに基づき、上記生理状態変化の原因遺伝子である標的遺伝子を含む遺伝子クラスタ及び/または該遺伝子クラスタ中の標的遺伝子を探索する方法。2生物細胞の生理状態変化を生じる条件とコントロール条件下とを一の対比条件セットとして、該対比条件セットが一種以上設定されていることを特徴とする上記1に記載の方法。3生理状態変化を生じる条件とコントロール条件が、少なくとも代謝産物の産生誘導条件下と非誘導条件下あるいは代謝産物の産生抑制条件下と非抑制条件下との対比条件セットを含むことを特徴とする上記1または2に記載の方法。4代謝物産生に関与する遺伝子が2次代謝物産生に関与する遺伝子であることを特徴とする上記3に記載の方法。5仮想の各遺伝子クラスタは、ゲノムDNA上に連続する遺伝子を2個から遺伝子数を一つずつ増やして、想定される遺伝子クラスタに含まれる最大限のゲノム遺伝子数になるまで抽出し、かつ該抽出において、抽出する遺伝子の各個数毎に、直鎖状DNAからなるゲノムの場合には該DNAのいずれかの末端から、環状DNAからなるゲノムの場合には任意の遺伝子を起点として順にゲノムDNA上に配列する遺伝子を一つずつずらしながら抽出された各遺伝子群からなることを特徴とする、上記1~4のいずれかに記載の方法。6スコアリングされる仮想の各遺伝子クラスタの集合体が、ゲノムDNA上に連続する遺伝子を2個から遺伝子数を一つずつ増やし想定される遺伝子クラスタに含まれる最大限のゲノム遺伝子数になるまで抽出し、かつ該抽出において、抽出する遺伝子の各個数毎に、直鎖状DNAからなるゲノムの場合には該DNAのいずれかの末端から、環状DNAからなるゲノムの場合には任意の遺伝子を起点として順にゲノムDNA上に配列する遺伝子を一つずつらしながら抽出された各遺伝子群からなる仮想の各遺伝子クラスタの集合からなり、ゲノム上に存在する遺伝子クラスタの全てが仮想の遺伝子クラスタの集合体中に含まれるように構成されていることを特徴とする上記1~5のいずれかに記載の方法。7仮想の各遺伝子クラスタのスコアリングが以下の計算式aによりなされることを特徴とする、上記1~6のいずれかに記載の方法。計算式a8ゲノムDNA上に配列する遺伝子が、標的とする遺伝子機能を有すると推定される場合、あるいは標的とする遺伝子機能を有する可能性が低いか若しくはその可能性がないと推定される場合において、当該ゲノムDNA上に配列する遺伝子については、以下の重み付け計算が適用されることを特徴とする、上記7の方法。9ゲノムDNA上に配列する遺伝子が、標的とする遺伝子機能を有すると推定される場合において、標的とする遺伝子機能を有すると推定された遺伝子を含む仮想の遺伝子クラスタを選出し、選出された仮想の遺伝子クラスタについて、スコアリングすることを特徴とする、上記7に記載の方法。10仮想の遺伝子クラスタが、ゲノムにおいて、近傍に存在することを条件として、以下の1~3の内の1以上の遺伝子のみから、あるいは該遺伝子を少なくとも含む1以上の遺伝子から構築されることを特徴とする、上記4に記載の方法。12次代謝物産生に関与していると想定される酵素種に属する酵素遺伝子。2トランスポーター遺伝子3転写因子をコードする遺伝子11仮想の各遺伝子クラスタのスコアリングが以下の計算式aによりなされることを特徴とする、上記10に記載の方法。計算式a12仮想の遺伝子クラスタ全体のスコアの分布から乖離して存在するスコアを有する仮想の遺伝子クラスタを、標的の遺伝子クラスタ候補として選定することを特徴とする、上記1~11のいずれかに記載の方法。13仮想の遺伝子クラスタ全体のスコアの分布からの乖離の程度を示す判定値Iχを、以下の計算式bにより算出し、算出された該判定値Iχに基づき仮想の遺伝子クラスタを標的の遺伝子クラスタ候補として選定することを特徴とする、上記12に記載の方法。計算式b14仮想の遺伝子クラスタ全体のスコアの分布からの乖離の程度を示す判定値IIυを、以下の計算式cにより算出し、算出された判定値IIを求め、該判定値に基づき、予め、標的とする遺伝子クラスタがゲノム中に存在するか否かあるいは、標的クラスタが存在する場合のその遺伝子サイズを予測することを特徴とする、上記方法。計算式a計算式e17遺伝子数がk個のときのε値εkと、その前後数のときのε値εk-1、εk1が、以下の関係にあるとき、標的とする遺伝子クラスタがゲノム中に存在すると判定し、標的遺伝子クラスタに含まれる遺伝子数をk個と予想することを特徴とする、上記16に記載の方法。2本発明は、また、以下の有用遺伝子の探索、同定を行うための装置、およびそのためのプログラムを提供する。18生物ゲノム中の標的遺伝子を含む遺伝子クラスタ及び/または該遺伝子クラスタ中の標的遺伝子を探索する装置であって、a生物細胞の生理状態変化を生じる条件とコントロール条件下におけるゲノムDNA上に配列する各遺伝子の発現量データに基づき算出された上記2つの条件下における上記各遺伝子の発現量変動比を記憶する手段、bゲノムDNA上に配列する複数の遺伝子を組み合わせて仮想の遺伝子クラスタを構築する手段、c該算出され、記憶されたゲノムDNA上に配列する各遺伝子の発現量変動比を複数の遺伝子により構築された上記仮想の遺伝子クラスタ単位の発現量変動比として合算し、仮想の遺伝子クラスタ単位毎にスコアリングし、仮想の各遺伝子クラスタのスコアを記憶する手段、及びd得られたスコアに基づき上記生理状態変化の原因遺伝子である標的遺伝子を含む遺伝子クラスタを選定する手段を有するか、あるいはさらにe選定された遺伝子クラスタ中に含まれる遺伝子を表示する手段を有することを特徴とする、上記装置。19発現量データが、遺伝子発現量測定用DNAマイクロアレイによる蛍光強度情報であることを特徴とする上記18に記載の装置。20蛍光強度情報が、蛍光強度を読み取り、数値化する手段を有する蛍光強度読み取り装置により出力される数値データであることを特徴とする、上記19に記載の装置。21生物細胞の生理状態変化を生じる条件とコントロール条件とを1の対比条件セットとして1以上設定されている場合において、各対比条件セットに含まれる条件毎に各遺伝子の発現量データが入力され、各対比条件セットにおける同一遺伝子の発現量変動比が算出されることを特徴とする、上記18~20のいずれかに記載の装置。22標的遺伝子が代謝物産生に関与する遺伝子であることを特徴とする、上記18~21のいずれかに記載の装置。23代謝物産生に関与する遺伝子が2次代謝物産生に関与する遺伝子であることを特徴とする、上記22に記載の装置。24設定される対比条件セットが、少なくとも代謝産物の産生誘導条件下と非誘導条件下あるいは代謝産物の産生抑制条件下と非抑制条件下との対比条件セットを含むことを特徴とする上記22に記載の装置。25代謝産物が2次代謝産物であることを特徴とする、上記24に記載の装置。26仮想の各遺伝子クラスタの構築手段が、ゲノムDNA上に連続する遺伝子を2個から遺伝子数を一つずつ増やして、想定される遺伝子クラスタに含まれる最大限のゲノム遺伝子数になるまで抽出し、かつ該抽出において、抽出する遺伝子の各個数毎に、直鎖状DNAからなるゲノムの場合には該DNAのいずれかの末端から、環状DNAからなるゲノムの場合には任意の遺伝子を起点として順にゲノムDNA上に配列する遺伝子を一つずつずらしながら抽出した各遺伝子群により構築する手段であることを特徴とする、上記18~25のいずれかに記載の装置。27仮想の各遺伝子クラスタのスコアリングが以下の計算式aによりなされることを特徴とする、上記18~26のいずれかに記載の装置。計算式a28ゲノムDNA上に配列する各遺伝子中の特定の遺伝子を選定するためのアノテーション付与手段を有し、上記遺伝子クラスタのスコアリングにおいて、付与されたアノテーションに基づき選定された遺伝子についての発現量変動比計算を以下の重み付け計算式により行うことを特徴とする、上記27に記載の装置。29アノテーション付与手段が、それぞれ遺伝子機能の種類毎に異なるアノテーションを付与する手段であることを特徴とする上記28に記載の装置。30アノテーションに基づき選定される遺伝子が、1~3のうちの1以上の遺伝子であることを特徴とする、上記29に記載の装置12次代謝物産生に関与していると想定される酵素種に属する酵素遺伝子。2トランスポーター遺伝子3転写因子をコードする遺伝子31上記28~30のいずれかに記載のアノテーション付与手段と、構築された仮想の遺伝子クラスタから、アノテーションに基づき選出された遺伝子を含む仮想の遺伝子クラスタを選出する手段を有し、選出された仮想の遺伝子クラスタについてスコアリングすることを特徴とする、上記27に記載の装置。32ゲノムDNA上に配列する各遺伝子中の特定遺伝子を選定するためのアノテーション付与手段を有し、ゲノムDNA上において近傍に位置することを条件として、アノテーションに基づき選定された遺伝子により、あるいは該遺伝子を少なくとも含む1以上の遺伝子から仮想の遺伝子クラスタを構築する手段を有することを特徴とする、上記18~25に記載の装置。33上記32に記載のアノテーション付与手段が、それぞれ遺伝子機能の種類に応じたアノテーションを付与する手段であることを特徴とする上記32に記載の装置。34アノテーション付与に基づき選定される遺伝子が、1~3のうちの1以上の遺伝子であることを特徴とする、上記33に記載の装置12次代謝物産生に関与していると想定される酵素種に属する酵素遺伝子。2トランスポーター遺伝子3転写因子をコードする遺伝子35仮想の各遺伝子クラスタのスコアリングが以下の計算式aによりなされることを特徴とする、上記32~34のいずれかに記載の装置。計算式a36仮想の遺伝子クラスタ全体のスコアの分布から乖離して存在するスコアを有する仮想の遺伝子クラスタを、標的の遺伝子クラスタ候補として選定する手段を有することを特徴とする、上記18~35のいずれかに記載の装置。37標的の遺伝子クラスタ候補として選定する手段として、仮想の遺伝子クラスタ全体のスコアの分布からの乖離の程度を示す判定値Iχを、以下の計算式bにより算出するプログラムが格納されていることを特徴とする、上記36に記載の装置。計算式b38標的の遺伝子クラスタ候補として選定する手段として、遺伝子クラスタ全体のスコアの分布からの乖離の程度を示す判定値IIυを、以下の計算式cにより算出するプログラムが格納されていることを特徴とする、上記36に記載の装置。計算式c39さらに、以下の計算式dの算出結果に基づき、bが100未満の仮想のクラスタを少なくとも除外し、標的の遺伝子クラスタ候補をさらに絞り込むプログラムが格納されていることを特徴とする、上記37または38に記載の装置。計算式d40a生物細胞の生理状態変化を生じる条件とコントロール条件下において生じたゲノムDNA上に配列する各遺伝子の発現量を入力する手段、b入力された上記2つの条件下における同一遺伝子の発現量の比を算出する発現量変動比算出手段、c該算出されたゲノムDNA上に配列する各遺伝子の発現量変動比を複数の遺伝子により構築された上記仮想の遺伝子クラスタ単位の発現量変動比として合算し、仮想の遺伝子クラスタ単位毎にスコアリングする手段、及びd得られた仮想の遺伝子クラスタのスコアから遺伝子クラスタに含まれる遺伝子数単位毎の遺伝子クラスタ分布判定値εを算出する手段を有し、該遺伝子クラスタ分布判定値εから、標的とする遺伝子クラスタがゲノム中に存在する否かあるいは、標的遺伝子クラスタが存在する場合の遺伝子サイズを予測する装置であって、仮想の遺伝子クラスタの構築手段が、ゲノムDNA上に連続する遺伝子を2個から遺伝子数を一つずつ増やし想定される遺伝子クラスタに含まれる最大限のゲノム遺伝子数になるまで抽出し、かつ該抽出において抽出する遺伝子の各個数毎に、直鎖状DNAからなるゲノムの場合には該DNAのいずれかの末端から、あるいは環状DNAからなるゲノムの場合には任意の遺伝子を起点として順にゲノムDNA上に配列する遺伝子を一つずつずらしながら抽出された各遺伝子群を仮想の各遺伝子クラスタとする手段であり、上記仮想の遺伝子クラスタ単位のスコアリング手段は以下の計算式aによる演算手段からなるとともに、上記遺伝子クラスタ分布判定値εの算出手段が、以下の計算式eによるものであることを特徴とする、上記装置。計算式a計算式e41遺伝子数がk個のときの遺伝子クラスタ分布判定値ε値εkと、その前後数のときの同ε値εk-1、εk1が、以下の関係にあるとき、標的とする遺伝子クラスタがゲノム中に存在すると判定し、標的遺伝子クラスタに含まれる遺伝子数をk個とする予想値を出力することを特徴とする、上記40に記載の装置。42上記26に記載の仮想の遺伝子クラスタの構築手段を実行するプログラムであって、ゲノム遺伝子の位置情報に基づき、以下の1または2の手段を実行することを特徴とする、仮想の遺伝子クラスタ構築プログラム。1ゲノム遺伝子が直鎖状ゲノムの場合、a.ゲノムDNAの一方の末端に位置する遺伝子を起点として、他方の末端方向に、順次、ゲノムDNA上に連続する遺伝子を同一方向に2個から一つずつ増やして想定される遺伝子クラスタに含まれる遺伝子数の最大限になるまで組み合わせ、起点とした遺伝子を含み、かつ遺伝子の個数の異なる複数の遺伝子群を構成する手段。b.起点を、順次、他方の末端方向に一遺伝子ずつずらせながら、上記a.と同様の処理を行い、新たな起点遺伝子を含みかつ遺伝子の個数が異なる複数の遺伝子群を構成し、a.の遺伝子群と併せて、複数の遺伝子を組み合わせた遺伝子群からなる仮想の遺伝子クラスタを構築する手段。2ゲノム遺伝子が環状の場合、ゲノムDNA上の任意の遺伝子を起点として、上記1a.及びb.と同様の処理を順次行い、最初に起点とした遺伝子が起点となる時点で処理を終了する手段。43上記42のプログラムにより構築された仮想の遺伝子クラスタについて、以下の計算式aによるスコアリングを実行することを特徴とする、仮想の遺伝子クラスタのスコアリングプログラム。計算式a44上記遺伝子クラスタのスコアリングにおいて、付与されたアノテーションに基づきゲノム遺伝子を選定し、選定された遺伝子についての発現量変動比計算を以下の重み付け計算式により行うことを特徴とする、上記43に記載のプログラム。45上記遺伝子クラスタのスコアリングにおいて、付与されたアノテーションに基づきゲノム遺伝子を選定し、構築された遺伝子クラスタの中から、該選定されたゲノム遺伝子を含む仮想の遺伝子クラスタを選出し、選出された仮想の遺伝子クラスタについてスコアリングを実行することを特徴とする、上記43に記載のスコアプログラム。46上記32に記載の仮想の遺伝子クラスタの構築手段を実行するプログラムであって、ゲノムDNA上において近傍に位置することを条件として、アノテーションに基づき選定された遺伝子により、あるいは該遺伝子を少なくとも含む1以上の遺伝子から仮想の遺伝子クラスタを構築することを特徴とする、仮想の遺伝子クラタの構築プログラム。47上記46のプログラムにより構築された仮想の遺伝子クラスタについて、以下の計算式aによるスコアリングを実行することを特徴とする、仮想の遺伝子クラスタのスコアリングプログラム。計算式a48上記43~45又は47のいずれかに記載のスコアリングプログラムにより算出された各仮想の遺伝子クラスタのスコアについて、仮想の遺伝子クラスタ全体のスコアの分布からの乖離の程度を算出するプログラムであって、以下の計算式bにより、判定値Iχを算出することを特徴とする、上記プログラム。計算式b49上記43~45又は47のいずれかに記載のスコアリングプログラムにより算出された各仮想の遺伝子クラスタのスコアについて、仮想の遺伝子クラスタ全体のスコア分布からの乖離の程度を算出するプログラムであって、以下の計算式cにより判定値IIυの算出を実行する、上記プログラム。計算式c50生物細胞の生理状態変化を生じる条件とコントロール条件下とにおけるゲノムDNA上に配列する各遺伝子の発現量変動比を複数の遺伝子により構築された上記仮想の遺伝子クラスタ単位の発現量変動比として合算し、仮想の遺伝子クラスタ単位毎にスコアリングする手段、及び得られた仮想の遺伝子クラスタのスコアから遺伝子クラスタに含まれる遺伝子数単位毎の遺伝子クラスタ分布判定値εを算出し、該遺伝子クラスタ分布判定値εから、標的とする遺伝子クラスタがゲノム中に存在する否かあるいは、標的遺伝子クラスタが存在する場合の遺伝子サイズを予測する手段に用いるプログラムであって、少なくとも以下A~Cの手段を実行するプログラム。Aゲノム遺伝子の位置情報に基づき、以下の1または2の手段により仮想の遺伝子クラスタを構築する手段、1ゲノム遺伝子が直鎖状の場合、a.ゲノムDNAの一方の末端に位置する遺伝子を起点として、他方の末端方向に、順次、ゲノムDNA上に連続する遺伝子を同一方向に2個から一つずつ増やして想定される遺伝子クラスタに含まれる遺伝子数の最大限になるまで組み合わせ、起点とした遺伝子を含み、かつ遺伝子の個数の異なる複数の遺伝子群を構成する手段。b.起点を、順次、他方の末端方向に遺伝子一つずつずらせながら、上記a.と同様の処理を行い、新たな起点遺伝子を含みかつ遺伝子の個数が異なる複数の遺伝子群を構成し、a.の遺伝子群と併せて、複数の遺伝子の組み合わせた遺伝子群からなる仮想の遺伝子クラスタを構築する手段。2ゲノム遺伝子が環状の場合、ゲノムDNA上の任意の遺伝子を起点として、上記1a.及びb.と同様の処理を順次行い、最初に起点とした遺伝子が起点となる時点で処理を終了する手段。B上記Aの手段により構築された仮想の遺伝子クラスタについて、以下の計算式aにより仮想の遺伝子クラスタ単位毎にスコアリングする手段。計算式aC上記Bの手段により得られた仮想の遺伝子クラスタのスコアから、以下の計算式eにより仮想の遺伝子クラスタに含まれる遺伝子数単位毎の遺伝子クラスタ分布判定値εを算出する手段。計算式e51遺伝子数がk個のときの遺伝子クラスタ分布判定値ε値εkと、その前後数のときの同ε値εk-1、εk1が、以下の関係にあるとき、標的とする遺伝子クラスタがゲノム中に存在すると判定し、標的遺伝子クラスタに含まれる遺伝子数をk個とする予想値を出力することを特徴とする、上記50に記載のプログラム。従来技術では主にDNAマイクロアレイなどを用いることにより、例えば、代謝産物の産生に関与する遺伝子を探索する場合、着目する化合物が産生される、あるいは着目する活性が観察される条件における発現の誘導や強い発現強度を有することなどを指標として標的となる遺伝子の同定を行ってきた。しかし、生物情報特有のデータの曖昧性、誤り、複雑性などによって正しい遺伝子を高い精度で予測することは極めて困難であった。これに対して、本発明の遺伝子探索方法および装置は、隣接あるいは近傍に位置する複数の遺伝子から、仮想の遺伝子クラスタを構成し、この仮想の遺伝子クラスタをまず探索対象として、有用遺伝子の探索を行うものであり、その手法自体極めて論理的、機械的であり、従来のDNAマイクロアレイによる解析にみられるような研究者の知識や経験等に大きく依存することなく、コンピューターを使用して、迅速、正確に有用遺伝子を特定することが可能となるとともに、同時に該遺伝子を含む遺伝子クラスタも特定できる。一方、本発明の遺伝子探索法においては、探索条件に誤りがある場合には、取得されたデータ自体から把握でき、この場合には探索条件を再設定し、探索をやり直すことができる。これに対し、従来法においては、解析結果が誤りであるか否かの判断には、遺伝子破壊実験等の検証実験を必要とし、膨大な費用、手間をかけざるを得ない。したがって、本発明の遺伝子探索法および探索装置の有利性は明らかである。また、本発明の遺伝子探索法および装置は、代謝産物産生遺伝子、とりわけ従来困難であった2次代謝物産生遺伝子の探索に極めて適している。これは、2次代謝物質の産生に関与する遺伝子が遺伝子クラスタを構成していることが多いからである。さらに、このようにして探索され、特定された2次代謝物質産生遺伝子等の有用遺伝子等の配列情報を利用すれば、新たな類似遺伝子の取得も可能となる。しかし、本発明の遺伝子探索法および装置によれば、このような代謝産物産生遺伝子の探索のみではなく、広く普遍性を有し、生物の様々な生理状態変化をもたらす原因遺伝子、それのみならず、同時に生理状態変化に関与する遺伝子クラスタも探索可能であり、これにより該原因遺伝子と協働する他の遺伝子も特定することが可能となる。したがって、本発明は、例えば、代謝産物、特に2次代謝物質産生遺伝子、様々な疾病の原因遺伝子、あるいはこれらと共働する遺伝子等の探索に極めて有効であり、新たな有用な化合物の取得、その大量生産あるいは医薬品開発等において、その技術を飛躍的に向上させることができる。本発明は、生物細胞の生理状態変化を生じる条件とコントロール条件下において生じたゲノムDNA上に配列する遺伝子の発現量変動比を、ゲノムDNA上に配列する複数の遺伝子により構成される仮想の遺伝子クラスタ単位の発現量変動比として合算することにより、仮想の遺伝子クラスタ単位毎にスコアリングし、得られたスコアに基づき、まず、上記生理状態変化の原因遺伝子である標的遺伝子を含む遺伝子クラスタを特定し、さらに該クラスタから標的遺伝子を特定する方法である。本発明は、また、上記方法を基本原理とし、生物ゲノム中の標的遺伝子を含む遺伝子クラスタ及び/または該遺伝子クラスタ中の標的遺伝子を探索する装置以下、単に、本発明の遺伝子探索装置という場合がある。に関するものであり、さらに該装置の一部を応用した遺伝子クラスタの有無及びそのサイズを予測する装置に関する。本発明の探索法および探索装置においては、真核生物、原核生物を問わず、あらゆる生物種について、ゲノム中の有用遺伝子を含有する遺伝子クラスタを探索対象とすることができる。また、本発明によれば、ゲノムの配列が明らかになっているものであれば、遺伝子クラスタの境界が明らかになっていない場合であっても本発明の手法および装置を適用でき、遺伝子クラスタ及び該クラスタ中の有用遺伝子を探索することができる。本発明における生理状態変化とは、例えば、生物の代謝物産生量の変化、分泌物質の種類と量の変化、増殖速度など増殖相グロースフェーズの違い、静止期間期などの細胞の分裂状態の違い、細胞の形態や機能の違い菌糸、分生子など分化状態の違いを含む等をいい、本発明においては、これら生理状態変化を生じる条件とコントロール条件とを一の対比条件セットとして、該対比条件セットを1種あるいは2種以上設定し、それぞれの対比条件セットの各条件下におけるゲノム遺伝子の発現量を測定し、その比発現量変動量比を求める。生理状態変化を生じる条件とは、例えば、薬剤の使用、温度、栄養源、培地、培養時間等の調整により、人為的に生理状態変化を誘導する場合の他、特にこのような誘導をせず、経時的に生理状態変化が生じる場合の時間条件も含める。コントロール条件とは、生理状態変化を生じないかあるいは生じても変化が少なく、生理状態変化を生じる条件下での生理状態変化と対比しうるものをいう。例えば、2次代謝物の産生に関与する遺伝子クラスタあるいは遺伝子を探索する場合、2次代謝物産生誘導条件下あるいは抑制条件下とコントロール条件としての2次代謝物産生非誘導条件下あるいは産生条件下におけるゲノム遺伝子の発現量を測定する。比較する上記2次代謝物産生誘導条件と2次代謝物産生非誘導条件、あるいは2次代謝物産生抑制条件と2次代謝物産生条件とは、代謝物産生速度、量等に差が生じる条件であればよく、例えば、薬剤の使用、温度、栄養源、培地等の調製の有無等、あるいは特にこのような誘導をせず、経時的に2次代謝物産生量に生じる場合の時間条件も含まれる。本発明における遺伝子クラスタ及び遺伝子探索法の全体的な流れを、図1に示した。このうちグレーで示した大きな四角の内部白い二つの四角を含むが、本発明の特徴部分である。本発明のプロセスにおいては、ゲノムDNA上に配列する各遺伝子の発現量の測定は、例えばマイクロアレイ等により行うが、その他のプロセスは、ゲノムDNA上に配列する遺伝子の発現量データに基づき、数学的データ処理により行うことができ、実験を必要とせず、また、上記発現量測定対象とするゲノム遺伝子の選定等も機械的に、あるいは研究者の特別な知識あるいは勘にほとんど左右されることがなく行うことができる。したがって、本発明の探索法は、コンピューター利用に極めて適しており、本発明によれば、迅速、効率的に有用遺伝子が探索可能となり、従来困難であって、代謝物、とりわけ2次代謝物産生に関与する遺伝子及び該遺伝子が含まれる遺伝子クラスタの探索に特に効力を発揮する。以下、本発明のプロセスについて、さらに具体的に説明する。本発明における上記仮想の遺伝子クラスタの構成手法としては、例えば、AゲノムDNA上に配列する複数の遺伝子を配列順に組み合わせ、各サイズの異なる仮想の遺伝子クラスタを構成する手法、B近傍に位置し機能的に遺伝子クラスタを構成する可能性のある複数の遺伝子から構成する手法が挙げられる。この2つの手法は、発現量を測定する遺伝子の対象範囲が異なるため、使用する発現量変動比データ、仮想の遺伝子クラスタの構成ゲノム遺伝子が異なるが、構成された仮想の遺伝子クラスタスコアリング等その他の数学的処理プロセス自体は共通している。以下に、本発明のプロセスについて、順次具体的に説明する図1参照。1上記Aの手法による場合の発現量の測定及び発現量変動比データの取得、Aの手法による場合、原則、ゲノムDNA上に配列する各遺伝子全てについて、生理状態変化を生じる条件とコントロール条件下とにおいて、それぞれ発現量を測定し、両条件下における発現量の比を求め、発現量変動比生理状態変化条件下での発現量を分子、コントロール条件下での発現量を分母として算出した値とする。発現量の測定は、例えば、ゲノムDNA上に配列する各遺伝子に特異的なプローブを有するマイクロアレイを用いてそれ自体周知の方法で行うことができる。例えば、代謝産物、特に2次代謝物の産生に関与する有用遺伝子を標的とする場合、1以上の2次代謝物産生誘導条件下あるいは抑制条件下で細胞を培養し、細胞からゲノムRNAを抽出し、ゲノムDNA上の各遺伝子に特異的なプローブを有するマイクロアレイでゲノムDNA上の各遺伝子の発現量を測定する。一方、コントロール条件として、上記2次代謝物の産生非誘導条件下の場合における発現量を測定し、両条件下における発現量の比をとり、これを発現量変動比とする。各遺伝子発現量の測定は、例えば、上記培養細胞からmRNAを抽出して、色素等でラベリングし、各遺伝子クラスタにおける上記各遺伝子中のDNA配列の一部を有するオリゴDNAをプローブとして基板に固定化したアレイを用い、上記該ラベリングしたmRNAを各オリゴDNAにハイブリダイズさせ、洗浄した後、発光強度等を測定することにより行う。2上記Aの手法による仮想の遺伝子クラスタの構築仮想の各遺伝子クラスタは、ゲノムDNA上に連続する遺伝子を2個から遺伝子数を1ずつ増やして、想定される遺伝子クラスタに含まれる最大限の遺伝子数になるまで抽出し、かつ該抽出において、抽出する遺伝子の各個数毎に、直鎖状DNAからなるゲノムの場合には該DNAのいずれかの末端から、環状DNAからなるゲノムの場合には任意の遺伝子を起点として順にゲノムDNA上に配列する遺伝子を一つずつずらしながら抽出された各遺伝子群から構成される。この仮想の遺伝子クラスタの構築手法をより具体的に示すと例えば以下の手法が挙げられる。1ゲノム遺伝子が直鎖状ゲノムの場合、aゲノムDNAの一方の末端に位置する遺伝子を起点として、他方の末端方向に、順次、ゲノムDNA上に連続する遺伝子を同一方向に2個から一つずつ増やしてN1、想定される遺伝子クラスタに含まれる遺伝子数の最大限nclになるまで組み合わせ、起点とした遺伝子を含み、かつ遺伝子の個数の異なる複数の遺伝子群を構成する。b起点を、順次、他方の末端方向に一遺伝子づつずらしながら、上記aと同様の処理を行い、新たな起点遺伝子を含みかつ遺伝子の個数が異なる複数の遺伝子群を構成し、aの遺伝子群と併せて、複数の遺伝子の組み合わせた遺伝子群からなる仮想の遺伝子クラスタを構築する。2ゲノム遺伝子が環状の場合、ゲノムDNA上の任意の遺伝子を起点として、上記1a及びbと同様の処理を順次行い、最初に起点とした遺伝子が起点となる時点で処理を終了する最初に起点として遺伝子に基づく仮想の遺伝子クラスタの構築は再度行わない。。上記仮想の遺伝子クラスタの構築においては、仮想の遺伝子クラスタが複数の遺伝子から構成される点で、遺伝子2個から一つずつ増やす手法が採用されるが、本発明は、遺伝子1個から一つずつ増やす手法を排除するものではない。すなわち、この場合、遺伝子一個の場合が構築される仮想の遺伝子クラスタに混入することになるが、本発明においては、この混入した遺伝子を含む遺伝子2個以上の組み合わせからなる仮想の遺伝子遺伝子クラスタが必ず構築され、また、仮想の遺伝子クラスタのスコアは、組み合わせた各遺伝子の発現量変動比の合算であるから、ゲノム中に標的遺伝子が存在した場合、この標的遺伝子単独のスコアに比べ、これを含む仮想の遺伝子クラスタのスコアは少なくとも同等以上になり、上記混入は実質的な問題ではない。したがって、仮想の遺伝子構築において、遺伝子2個から一つずつ増やす手法を含む限り、遺伝子1個から一つずつ増やしたとしても、本発明に包含される。例えば、以下のように、ゲノムDNA上に配列する遺伝子がA~Jの10個である場合、構築される仮想の遺伝子クラスタは、表1の示される各遺伝子群から構成される。すなわち、上記抽出によって構築される仮想の各遺伝子クラスタは、以下の各遺伝子群からなる。遺伝子の2個の各仮想の遺伝子クラスタ9個;ABBCCDDEEFFGGHHIIJ同3個の仮想の各遺伝子クラスタ8個;ABCBCDCDE,DEFEFGFGHGHIIJK同4個の仮想の各遺伝子クラスタ7個;ABCDBCDECDEFDEFGEFGHFGHIGHIJ同5個の仮想の各遺伝子クラスタ6個;ABCDEBCDEFCDEFGDEFGHEFGHIFGHIJ同6個の仮想の各遺伝子クラスタ5個;ABCDEFBCDEFGCDEFGHDEFGHI,EFGHIJ同7個の仮想の各遺伝子クラスタ4個;ABCDEFGBCDEFGHCDEFGHIDEFGHIJ同8個の仮想の各遺伝子クラスタ3個;ABCDEFGHBCDEFGHICDEFGHIJ同9個の仮想の各遺伝子クラスタ2個;ABCDEFGHI、BCDEFGHIJ同10個の仮想の各遺伝子クラスタ1個;ABCDEFGHIJしたがって、この場合、仮想の各遺伝子クラスタの構築数は45個であるが、これらの各遺伝子クラスタはデータ上において構築されるだけであって、実験によって実際に構築されるものではない。なお、実際のゲノムDNA上の遺伝子数は、麹菌の場合、外部データベースDOGANに登録されているもので12084個であり、これよりも遺伝子の定義を緩めてDNAマイクロアレイのプラットフォーム作成に使用されたものの場合14032個である。このうち連続していることが判明しているゲノム上の領域から、仮想の遺伝子クラスタを構築する。抽出する遺伝子の数の最大限は論理上ゲノム中の遺伝子の数とすることができるが、想定される遺伝子クラスタサイズの最大限の遺伝子数でよく、実際問題として、遺伝子クラスタを構成する遺伝子の数は、最大でも30個程度であり、これを超える必要は通常ない。1’上記Bの手法による場合の発現量の測定及び発現量変動比データの取得このB手法は、上記Aの手法に比べ簡便であり、2次代謝産物の産生に関与する遺伝子クラスタ及び該クラスタ中の2次代謝産物産生遺伝子の探索に特に適している。この手法は、ゲノムDNAの配列中、12次代謝に関与していると想定される酵素種に属する酵素遺伝子、2トランスポーター遺伝子、3転写因子をコードする遺伝子の内、1種以上、好ましくは2種以上が近傍に位置する場合、これらの遺伝子から、あるいはこれら遺伝子が含まれるようにゲノム遺伝子を組み合わせて仮想の遺伝子クラスタとするものであり、この場合において、近傍に位置する具体的な条件は、ゲノム上に配列する遺伝子数でいえば、上限30程度以内に存在すればよい。上記遺伝子の発現量の測定は、上記Aの手法と同様に、例えば2次代謝物産生誘導条件下あるいは抑制条件下で細胞を培養し、細胞からゲノムRNAを抽出し、ゲノムDNA上の各遺伝子に特異的なプローブを有するマイクロアレイを使用して、ゲノムの各遺伝子の発現量を測定し、上記2次代謝物産生非誘導条件下あるいは産生条件下の場合と比較して、発現量変動比を求める。この手法の場合、マイクロアレイにおける発現量測定においては、ゲノムDNA上の全遺伝子について行うが発現変動量の抽出対象遺伝子が絞られているため、これら遺伝子に対応する配列を有するプローブを用いたマイクロアレイのみを用いても良い。比較する上記2次代謝物産生誘導条件と2次代謝物産生非誘導条件、あるいは2次代謝物産生抑制条件と2次代謝物産生条件とは、代謝物産生速度、量等に差が生じる条件であればよく、例えば、薬剤の使用、温度、栄養源、培地等の調製の有無等、あるいは特にこのような誘導をせず、経時的に2次代謝物産生量に生じる場合の時間条件も含まれる。なお、この手法においても、上記Aの手法と同様に発現変動量の測定の他は格別の実験を必要とせず、数学的データ処理によりなされる。一方、ゲノム配列中の12次代謝に関与していると想定される酵素種に属する酵素遺伝子、2トランスポーター遺伝子、3転写因子をコードする遺伝子の判別は、既知の同じ酵素種の遺伝子との相同性あるいはモチーフ等により判別すればよく、例えば、各仮想の遺伝子クラスタ中の遺伝子配列中にこれらの遺伝子が存在するか否かは、上記酵素種に属する酵素、トランスポーター、転写因子の各アミノ酸配列に特有なモチーフと共通するアミノ酸配列をコードする塩基配列が、上記遺伝子クラスタ中に存在するか否かにより識別できる。これらについては市販ソフトウェアを使用することができる。すなわち、上記機能遺伝子の選定、及び以下に示す仮想の遺伝子クラスタのスコアリングにおいて重み付けする遺伝子の選定においては、アノテーション機能注釈付与を行い、これに基づき対象遺伝子を選定することが有効である。このようなアノテーション付与は、探索対象ゲノム上の各遺伝子の塩基配列情報等を基に、記憶部に記憶された検索対象のゲノム上の各遺伝子の位置情報中の遺伝子について行うもので、コンピューターにより自動で行うことができる。このようなアノテーション付与は、装置使用者が、検索探索対象ゲノム上の遺伝子について予め相同検索あるいはモチーフ検索等の結果に基づき、上記記憶されたゲノム上の各遺伝子の位置情報中の遺伝子を逐一指定し、この指定された遺伝子にアノテーションが付与されるように構成しても良いが、ゲノム上の遺伝子の数は極めて多数であり、上記モチーフ検索を行う市販のソフトウェアを付属のモチーフ情報とともにコンピューターに格納するか本発明装置に格納するか、あるいは該ソフトウェアをモチーフ情報とともに格納した外部コンピューターを利用することが好ましい。これにより、探索対象ゲノム上の各遺伝子の塩基配列情報を上記コンピューターあるいは外部コンピューターに入力することにより、期待される機能に対応したモチーフについて検索を行い、アノテーション付与する遺伝子を自動で選定することができる。また、他のアノテーション付与手段として、上記モチーフ検索により探索対象のゲノム上の全遺伝子にアノテーションを付与した後に、付与されたアノテーションの種類遺伝子機能から、期待される機能と一致する遺伝子を選定してもよい。このようにすれば、アノテーション付与を研究者の手を煩わすことなく自動で行うことができる。アノテーション付与は、機能が同様なゲノム遺伝子に付与しても良いし、機能の種類が異なる複数種の遺伝子に付与しても良い。機能が異なる複数種のゲノム遺伝子にアノテーションを付与する場合には、ゲノム遺伝子の機能毎に識別可能なように付与する。アノテーションによる選定の対象となる遺伝子は、例えば、2次代謝物産生に関与する遺伝子クラスタあるいはその中の遺伝子を標的とする場合、ゲノムDNAの配列中、12次代謝に関与していると想定される酵素種に属する酵素遺伝子、2トランスポーター遺伝子、3転写因子をコードする遺伝子を選定可能である。上記1の酵素遺伝子の判別において、酵素種は、2次代謝物の化学構造、前駆体、関与しうる補酵素、化学的物理的性質、既知の酵素反応の事例、生産効率速度等からその産生反応を推測し、関与する酵素種を想定するが、この酵素種の想定においては、実際にその反応に関与したであろう特定酵素のレベルまで想定しなければならないというわけではなく、該反応に関与することがより確実なレベルの酵素種でよい。例えば、オキシゲナーゼに属する酵素であることは分かるが、その下位概念の酵素種まで特定できないときは、オキシゲナーゼのレベルを酵素種として選定して、ゲノム上の各遺伝子の配列を探索し、その範疇に属する全てのゲノム遺伝子のそれぞれを、各仮想の遺伝子クラスタの構成遺伝子とすればよい。ただし、下位概念の酵素種を選定できた場合には、探索対象とする仮想の遺伝子クラスタの範囲が狭まる可能性があり、その分探索が効率的となる。また、2次代謝物産生反応において複数の酵素が関与していると想定できる場合には、その複数の酵素種を選定することも可能である。トランスポーター遺伝子及び転写因子遺伝子においても同様で、標的とする2次代謝物産生に直接関与している遺伝子を特定しなければならないというわけではない。2’上記Bの手法による仮想の遺伝子クラスタの構築上記Bの手法による場合、近傍に位置する2次代謝に関与していると想定される酵素種に属する酵素遺伝子、2トランスポーター遺伝子、3転写因子をコードする遺伝子のうち少なくとも1種以上、好ましくは2種以上の遺伝子を抽出し、これらを組み合わせることにより、あるいはこれら遺伝子が含まれるようにゲノムDNA上の遺伝子を抽出して仮想の遺伝子クラスタとする。例えば、以下のように、ゲノムDNA上に配列する遺伝子がA~Jの10個である場合、は該当する酵素種をコードする遺伝子、はトランスポーター遺伝子前者の場合、仮想の遺伝子クラスタは、AC及びGJとにより構成される。一方、後者の場合は、ABC及びGHIJで構成してもよく、さらにABCDEあるいはFGHIJのように、各仮想の遺伝子クラスタが一定数の遺伝子により構成されるようにゲノムを分割して、各仮想の遺伝子クラスタを構成しても良い。3仮想の遺伝子クラスタのスコアリング上記1のプロセスにより取得されたゲノムDNA上に配列する各遺伝子の発現量変動比は、各対比条件セット毎に、正規化され、上記2のプロセスにより構築された仮想の各遺伝子クラスタ単位で、以下の計算式aにより合算され、算出された値を、仮想の各遺伝子クラスタのスコアとする。計算式aなお、上記全ての仮想の遺伝子クラスタに含まれる全遺伝子とは、全ての仮想の遺伝子クラスタを構成するために抽出されたゲノムDNA上の全ての遺伝子をいう。一方、1’のプロセスにより取得された各遺伝子の発現量変動比も、同様に各対比条件セット毎に、正規化され、上記2のプロセスにより構築された仮想の各遺伝子クラスタ単位で合算されるが、この手法はアノテーション付与により選定された特定の遺伝子のみの発現量変動比を用いるため、計算式aの定義が異なる。すなわち上記式中、Mは各仮想の遺伝子クラスタのスコア、mはスコアリングされる仮想の各遺伝子クラスタに含まれるアノテーション付与に基づき選定された各遺伝子の発現量変動比、m-は全ての仮想の遺伝子クラスタに含まれるアノテーション付与に基づき選定された全遺伝子の発現量変動比m値の平均、smは全ての仮想の遺伝子クラスタに含まれるアノテーション付与に基づき選定された全遺伝子の発現量変動比m値の標準偏差を表す。本発明によれば、このようにして得られた一群の仮想の遺伝子クラスタのスコアに対する出現頻度分布をみる場合、全体としては大凡正規分布となるが、このような全体のスコア分布から離れて存在する仮想の遺伝子クラスタが存在すれば、少なくとも標的の遺伝子クラスタと対応していると判定できる。すなわち、この仮想の遺伝子クラスタは、該クラスタ中の少なくとも2つの遺伝子が、代謝物産生誘導条件下協働した結果、発現変動量の総量であるスコアが増大したものであり、標的の遺伝子クラスタとみなすことができ、この仮想の遺伝子クラスタ中の遺伝子は少なくとも実際の遺伝子クラスタ中に存在する代謝物産生に関与する遺伝子であると同定することができる。さらに、仮想の遺伝子クラスタ中の遺伝子及び必要に応じ代謝産物の産生機構を検討すれば、直接代謝産物の産生に関与する標的遺伝子のみではなく、未知の機能を有する遺伝子の発見も期待でき、さらに代謝物産生機構の全体像も明らかにすることができる。一方、Aの手法において、ゲノムDNA上に配列する遺伝子が、標的とする遺伝子機能を有すると推定される場合、あるいは標的とする遺伝子機能を有する可能性が低いか若しくはその可能性がないと推定できる場合においては、当該遺伝子については、以下の計算式により重み付けをしておくことができる。重みwの設定は、上記標的とする遺伝子機能を有すると推定される場合、1を超えるように設定し、標的とする遺伝子機能を有する可能性が低いか若しくは可能性がないと推定できる場合は、0以上1未満になるように設定する。標的とする遺伝子機能を有するか、あるいはその可能性が低いかの推定は上記と同様に既知の遺伝子との相同性あるいはモチーフ等により判別でき、上記したアノテーション付与手段を利用することができる。また、ゲノムDNA上に配列する遺伝子が、標的とする遺伝子機能を有すると推定される場合においは、Aの手法により構築された仮想の遺伝子クラスタの中から、標的とする遺伝子機能を有すると推定された遺伝子を含む仮想の遺伝子クラスタを選出し、選出された仮想の遺伝子クラスタのみについて、スコアリングすることも可能である。標的とする遺伝子機能を有するか否かについての推定においては、上記したアノテーション付与手段の全てを利用することができる。この手法によれば、スコアリングする対象となる仮想の遺伝子クラスタの数を低減することができる。また、この手法により選定された仮想の遺伝子クラスタは、結果として上記手法Bにより構築された仮想の遺伝子クラスタと同様になる場合があるが、この手法による場合、一度Aの手法による網羅的な仮想の遺伝子クラスタ群を構築しておけば、自由に標的とする遺伝子あるいはこれを含む遺伝子クラスタの機能を変更でき、機能選択的な遺伝子解析が容易に行える点で有利である。また、該当するアノテーションが付与されなかった遺伝子のスコアを考慮に含めることができるため、機能未知の遺伝子の影響が大きい場合などに柔軟に対応できる。本発明は、ゲノムDNA上の複数の遺伝子を組み合わせて仮想の遺伝子クラスタを構成し、これら複数の遺伝子の生理状態変化条件下の発現量変動比を合算して各仮想の遺伝子クラスタをスコアリングし、これに基づき、まず、標的の遺伝子クラスタを探索する方法である。スコアリングされてスコアが高いものが得られた場合は、仮想の遺伝子クラスタに含まれる複数の遺伝子が協働した結果であり、各遺伝子単独の発現量変動比スコアをみるよりも、全体のスコア分布に対する特異性がより鮮明となる。一方、従来のように一つ一つの遺伝子の発現変動量のみから有用遺伝子を検出する場合には、正解の遺伝子であっても、全体のスコア分布中に吸収されてしまい、高い順位の遺伝子であっても、目的の遺伝子であるか否かの遺伝子破壊実験等の検証を必要とする。これに加え、上記したように重み付けをした遺伝子についての発現量変動比は、Aの手法において構築される仮想の各遺伝子クラスタのスコアリングにおいて、他の遺伝子の発現量変動比と合算され、上記標的とする遺伝子機能を有すると推定される遺伝子を含む仮想の各遺伝子クラスタのスコアはより高くなり、反対に標的とする遺伝子機能を有する可能性が低いか若しくは可能性がないと推定される遺伝子を含む仮想の遺伝子クラスタのスコアはより低くなり、全体のスコア分布との乖離が明瞭となる。したがって、これにより、標的とする遺伝子機能を有する遺伝子あるいはこれを含む遺伝子クラスタの探索がより効率的になる。4全体のスコア分布からの乖離の程度の算出仮想の遺伝子クラスタ全体のスコアの分布からの乖離の程度を示す判定値は、上記3のプロセスにより算出されたスコアに基づき、例えば、以下の計算式bあるいはcから算出される。計算式b上記計算式b中のスコアMの出現頻度は、仮想の遺伝子クラスタの全てを含む集団における各スコアの出現頻度Pの累計を1としたときの値であるため、1を超えることはなく、したがってlogPは正になることはない。また、出現頻度が低いものほどlogPは-∞に近づくため、頻度の低いスコア値を持つ遺伝子クラスタほどlogPの絶対値は大きくなる。したがって、上記計算式bにおいては、logPと仮想の各遺伝子クラスタのスコアを掛け合わせて-1を乗算することにより、頻度が低くかつスコアの高いものが、より大きな判定値Iχを持つこととなる。上記計算式bによれば、判定値Iχが0を超え、高い値を示す仮想の遺伝子クラスタは、仮想の各遺伝子クラスタのスコアに対する出現頻度分布から離れており、高い判定値Iを示した仮想の遺伝子クラスタを標的の遺伝子クラスタあるいは標的の遺伝子クラスタに対応する候補として選定することができる。候補の選定は、例えば判定値Iが高い順に仮想の遺伝子クラスタを一定数選定するかあるいは判定値Iが一定値以上を示した仮想のクラスタを選定するか等により行う。計算式cこの判定値IIυは、仮想の各遺伝子クラスタのスコアについて、仮想の遺伝子クラスタ全体の平均スコアからのずれを、上記標準偏差の実数倍で割ったものを次元数d’乗したもので、正規分布様のスコアに対する出現頻度分布から乖離するスコアを有する仮想の遺伝子クラスタにおいて大きな値となる。上記式中d’は任意に設定できる正の整数たる次元数であり、値が大きくなるほど平均スコアからの隔たりが強調されることになる。あまり大きくしすぎると、平均スコアから大きく外れたものの値が強調されて相対的にそれ以外の値が小さくなるため、通常2または4に設定する。外れたものをより鋭敏に検出したい場合は、6以上の偶数とする。また式中のaは外れ度を表す係数で、この値を調節することにより、上記正規分布様分布からどの程度乖離したものをとるかを調節することができる。1を超えて大きく設定するほど、平均スコアから大きく外れたもの以外のυ値はゼロに近付くため、このa値は通常1~2に設定する。逆に1未満の場合、より外れ方の小さなものも拾うことができる。この計算式cによる場合も上記判定値Iと同様に、υが0を超え、高い値を示す仮想の遺伝子クラスタを標的の遺伝子クラスタあるいは標的の遺伝子クラスタに対応する候補として選定することができる。候補の選定は、例えば判定値IIが高い順に仮想の遺伝子クラスタを一定数選定するかあるいは判定値IIが一定値以上を示した仮想のクラスタを選定するか等により行う。5遺伝子クラスタ候補の絞り込み上記計算式b、cにより算出された判定値χあるいはυにより、標的の遺伝子クラスタ候補となった仮想の遺伝子クラスタの数が多く、さらに候補を絞り込みたい場合は、以下の計算式dの算出結果に基づき、bが100未満の仮想のクラスタを少なくとも除外することにより、標的の遺伝子クラスタ候補をさらに絞り込むことが可能である。計算式d上記計算式d中、bはどの程度の遺伝子クラスタ候補を絞り込むかを決定するための閾値であり、bを大きくとるほど候補の絞り込み効果がより高くなる。また小さくとるほど多くの候補遺伝子クラスタを選択することができる。bの値の設定は対象とする生物種や培養条件に依存する。すなわち、候補遺伝子クラスタが強くかつ多く発現している系であれば値を高くする必要があるが、逆に発現強度が弱くかつ数が少なければ値を低くしなければ候補遺伝子が出現しない。前者の場合、例えば5000~10000あるいは10000~30000の範囲内の任意の数値に設定し、後者の場合、通常100以上、例えば1000~2000、あるいは2000~5000の範囲内の任意の数値に設定する。6標的とする遺伝子クラスタの有無及び標的とする遺伝子クラスタが存在する場合のサイズの推定本発明においては、予めゲノム中に標的の遺伝子クラスタが存在するか否か及び標的遺伝子クラスタが存在する場合の遺伝子サイズクラスタを構成する遺伝子数;nclを推定することができる。この手法は、まず、生物細胞の生理状態変化を生じる条件とコントロール条件下において生じたゲノムDNA上に配列する遺伝子の発現量変動比を合算し、仮想の遺伝子クラスタのスコアとするが、上記発現量の測定、発現量変動比データの取得、仮想の遺伝子クラスタの構築、及び仮想の各遺伝子クラスタをスコアリングの各プロセスは、上記Aの手法中1~3のプロセスと同様なプロセスである。すなわち、この手法においては、生物細胞の生理状態変化を生じる条件とコントロール条件下において生じたゲノムDNA上の各遺伝子の発現量変動比を、ゲノムDNA上の複数の遺伝子により構成される仮想の遺伝子クラスタ単位の発現量変動比として合算することにより、仮想の遺伝子クラスタ単位毎にスコアリングするが、この仮想の各遺伝子クラスタは、ゲノムDNA上に連続する遺伝子を2個から遺伝子数を一つずつ増やし想定される遺伝子クラスタに含まれる最大限のゲノム遺伝子数になるまで抽出し、かつ該抽出において抽出する遺伝子の各個数毎に、直鎖状DNAからなるゲノムの場合には該DNAのいずれかの末端から、あるいは環状DNAからなるゲノムの場合には任意の遺伝子を起点として順にゲノムDNA上に配列する遺伝子を一つずつずらしながら抽出された各遺伝子群から構成する。このように構成された各遺伝子クラスタは、上記Aの手法中3のプロセスと同様に以下の計算式aにより、そのスコアが算出される。計算式aついで、このスコアを仮想の各遺伝子クラスタに含まれる遺伝子数毎に分け、以下の計算式eにより、各遺伝子数単位毎に遺伝子クラスタスコア分布判定値を求める。計算式eこの計算式eによれば、仮想の遺伝子クラスタが、実際のゲノムDNAにおいてクラスタを形成していなければ、仮想の遺伝子クラスタ中に含まれる、標的の生理状態変化に関与せず、発現量変動のない遺伝子の影響を受けるので、仮想の遺伝子クラスタのスコアMはサイズ遺伝子数;nclが大きくなるほど平均化され、すなわちスコアの平均値に近づくためサイズの増大に伴いε値は単調減少する図2上から1および3番目の曲線参照。しかし、仮想の遺伝子クラスタのあるサイズのものがクラスタを形成している場合、そのサイズにおいて、分布の偏りεは大きくなり、上記単調減少曲線とはならず、ε値はそのサイズにおいて特異点を形成する図2矢印で示した点参照。したがって、ε値が特異点を形成するか否か、及び特異点を形成した遺伝子クラスタのサイズから、遺伝子クラスタの存在及びそのサイズを推定することができる。具体的には、仮想の遺伝子クラスタの該クラスタに含まれる遺伝子数毎の集計において、ある遺伝子数k個のときのε値εkと、その前後数のときのε値εk-1、εk1が以下の関係にあれば、標的とする遺伝子クラスタがゲノム中に存在すると判定し、標的遺伝子クラスタに含まれる遺伝子数をk個と予想することができる。この手法は、本発明による標的の遺伝子クラスタの探索法、特にBの手法を行うに際し、予め行う手法として有効である。すなわち、遺伝子クラスタが存在し、そのサイズを予想できれば、予想されるサイズ内に、標的とする酵素種に属する酵素遺伝子、2トランスポーター遺伝子、3転写因子をコードする遺伝子が存在するゲノム配列のみを対象として、上記仮想の遺伝子クラスタとして探索すればよい。また、この手法は、ある条件下で細胞が何らかの生理的状態変化を起こす場合においては、どのような生理的状態変化であっても変化を対比する条件が設定できれば、その原因遺伝子はもちろんその変化を生じる機構そのものが全く不明である場合においても、その変化原因が遺伝子クラスタ中の遺伝子の連関にあるのか否か、遺伝子クラスタ中の遺伝子の連関による場合該クラスタの遺伝子サイズも容易に予測できる。すなわちこの手法は、生物の生理的変化が、極めて探索の難しい複数の遺伝子の連関によって生じている場合において、その原因が遺伝子クラスタ中の遺伝子の共働によるものであることを明らかにでき、かつそのサイズも予測できる点で極めて有用である。7本発明の手法により解が得られなかった場合一方、本発明の手法を行った結果、仮に仮想の遺伝子クラスタ全体のスコア分布から乖離したスコアの遺伝子クラスタが見いだされなかった場合、設定する生理状態変化条件、重み付けするゲノムDNA上の遺伝子の選定、あるいは上記Bの手法による仮想の遺伝子クラスタ構築のためのゲノムDNA上の遺伝子の選定等の探索条件設定に問題点がある。したがって、このような場合には、探索条件を再設定して、バックグランドの分布から離れたスコアの遺伝子クラスタが見つかるまで、上記した遺伝子クラスタの探索法を繰り返し行えばよい。すなわち本発明においては、得られたデータのみから、探索条件設定の問題点を把握できる。これに対して、上記したような従来法の場合には、もともと正解の遺伝子であっても、遺伝子全体の発現量についての分布中に埋もれてしまうので、得られたデータからでは正解か否かは不明であり、結果的に無意味かもしれない検証実験を繰り返さなければならない。次に、本発明のプロセスを実施するために用いられる、本発明の遺伝子探索装置について、説明する。本発明の遺伝子探索装置は、ゲノムDNA上に配列する遺伝子の発現量データに基づき、数学的データ処理を行うもので、研究者の特別な知識あるいは勘にほとんど左右されることがなく、迅速、効率的に有用遺伝子が探索可能となり、従来困難であった代謝物、とりわけ2次代謝物産生に関与する遺伝子及び該遺伝子が含まれる遺伝子クラスタの探索に特に効力を発揮する。本発明の遺伝子探索装置は少なくとも以下のa~fの手段により構成される。a生物細胞の生理状態変化を生じる条件とコントロール条件下におけるゲノムDNA上に配列する各遺伝子の発現量データを入力する手段。b入力された上記2つの条件下における各遺伝子の発現量の比を算出する手段。cゲノムDNA上に配列する複数の遺伝子を組み合わせて仮想の遺伝子クラスタを構築する手段。d該算出されたゲノムDNA上に配列する各遺伝子の発現量変動比を複数の遺伝子により構築された上記仮想の遺伝子クラスタ単位の発現量変動比として合算し、仮想の遺伝子クラスタ単位毎にスコアリングする手段。e得られたスコアに基づき上記生理状態変化の原因遺伝子である標的遺伝子を含む遺伝子クラスタを選定する手段。あるいはさらにf選定された遺伝子クラスタ中に含まれる遺伝子を表示する手段。このような手段を伴う本発明の装置の概要は図2に示される。なお、図2中、点線部は、本発明装置において、さらに記憶されるのが好ましいデータ及び該データに関連する処理部を示す。本発明装置は、データの入出力部キーボード、マウス、ディスプレイ等、該入出力部の制御を行う入出力制御インターフェース、記憶部ハードディスク、主記憶部メモリ、制御演算部CPU、外部ネットワークと接続する通信制御インターフェースを含む。本装置の記憶部には、各遺伝子の発現量データ、該発現量変動比データ、ゲノム上の各遺伝子位置データ、及び仮想の遺伝子クラスタのスコアデータが記憶され、さらに必要に応じて、塩基配列に対する遺伝子機能の対応データ、各遺伝子のアノテーションデータ、仮想の遺伝子クラスタのスコア乖離度データが順次格納される。また、制御演算部には、ゲノム各遺伝子の発現量変動比の算出部、ゲノム上の遺伝子の位置情報に基づき仮想の遺伝子クラスタの構築を行う仮想の遺伝子クラスタの構築部、及び上記算出された発現量変動比を合算し、仮想の遺伝子クラスタのスコアリングを行う、仮想の遺伝子クラスタのスコアリング部を少なくとも設ける。また、さらに必要に応じて、各遺伝子へのアノテーション付与部、アノテーションに応じて、仮想の遺伝子遺伝子のスコアリングにおいて重み付けを行う重み付け付与部、仮想の遺伝子クラスタ構築を選定された機能遺伝子に限定して行うための機能遺伝子選択部、仮想の遺伝子クラスタの全体分布からの乖離度を算出する仮想の遺伝子クラスタの乖離度算出部を設け、さらに算出された乖離度では、遺伝子クラスタ候補の選定が充分できない場合に遺伝子クラスタ候補の絞り込みを行う遺伝子クラスタ候補の絞り込み部を設けてもよい。一方、本発明の遺伝子探索装置においては、装置構成はそのままで、標的とする遺伝子クラスタの存在の有無、存在する場合そのサイズを予測する機能をさらに保有させることが可能であり、この場合には、仮想の遺伝子クラスタのサイズ毎にスコアリングするサイズスコアリング部、及び仮想の遺伝子クラスタ分布判定値ε算出部を設ける。本装置は、特別なコンピューターを必要とせず、一般的な、制御演算処理装置CPU、主記憶装置メモリ、記憶装置ハードディスク、入出力装置キーボード、マウス、ディスプレイからなるもので構成可能である。オペレーティングシステムは、Linux、Windows、Macのいずれも使用可能であるが、メモリ空間を考慮すると、64bitのものがより望ましい。メモリは、本装置が生物のゲノム全体を対象とすることを考慮して、できれば2GB以上のものが望ましいが、1GB程度であっても、微生物であれば可能である。なお、ゲノム上の各遺伝子の位置情報、および機能に対応した塩基配列のデータベースは、NCBIhttpwwwncbinlmnihgovやInterproScanhttpwwwebiacukToolsInterProScanなどの、外部データベースを利用することができる。以下に、本発明の装置について、その処理プロセスに従い、具体的に説明する。A遺伝子探索装置1ゲノムDNA上に配列する各遺伝子の発現量データ入力及び発現量変動比算出本発明装置の場合、原則、ゲノムDNA上に配列する各遺伝子全てについて、生理状態変化条件下とコントロール条件下における発現量を測定し、得られた各遺伝子の発現量データを本発明装置の入力手段に入力し、入力された各遺伝子の発現量データに基づき、発現量変動比が算出される。発現量の測定は、例えば、ゲノムDNA上に配列する各遺伝子に特異的なプローブを有するマイクロアレイを用いてそれ自体周知の手段により行うことができる。例えば、代謝産物、特に2次代謝物の産生に関与する有用遺伝子を標的とする場合、1以上の2次代謝物産生誘導条件下あるいは抑制条件下で細胞を培養し、細胞からゲノムRNAを抽出し、ゲノムDNA上の各遺伝子に特異的なプローブを有するマイクロアレイでゲノムDNA上の各遺伝子の発現量を測定する。一方、コントロール条件として、上記2次代謝物の産生非誘導条件下の場合における発現量を測定し、両条件下における発現量の比をとり、これを発現量変動比とする。各遺伝子発現量の測定は、例えば、上記培養細胞からmRNAを抽出して、色素等でラベリングし、上記各遺伝子中のDNA配列の一部を有するオリゴDNAをプローブとして基板に固定化したアレイを用い、上記該ラベリングしたmRNAを各オリゴDNAにハイブリダイズさせ、洗浄した後、発光強度等を測定することにより行う。マイクロアレイ中の各遺伝子の発光強度は、例えば、マイクロアレイ読み取り装置における走査手段を伴う画像読み取り手段により読み取り、読み取った発光強度を数値化して、上記aの入力手段により本発明の装置に入力する。このような画像読み取り装置は、市販されている装置が使用できるが、このような読み取り装置の手段全部あるいは例えば数値化手段等の一部手段を本発明の装置に組み込むか、あるいは該読み取り装置が出力する数値データを介して本発明装置の入力手段に自動で入力可能なように設計しても良い。本発明装置に入力された、上記両条件下に遺伝子の発光強度についての数値化データは、それぞれ本発明装置の記憶部に記憶しておき、この記憶された各条件下における数値化データを各遺伝子について呼び出し、発現量変動比生理状態変化条件下での発現量を分子、コントロール条件下での発現量を分母として算出した値算出プログラムを有する発現量変動比算出手段により、各遺伝子同一遺伝子毎にその発現量変動比を算出する。この算出には、必要に応じて、各遺伝子の発現強度による歪みの補正も含まれる。すなわち、遺伝子の発現量変動比は、発現の強度に依存して、ノイズの影響で値が強調されることがあるため、発現量変動比の分布がどの発現の強度に対してもほぼ一定となるようなバックグラウンド補正を行う。これらの発現量変動比算出プロセスには、フリーのソフトウェアであるR内のRowessアルゴリズム等が利用できる。算出された各遺伝子の発現量変動比は、本発明装置の記憶部に記憶される。一方、この各遺伝子の発現量変動比は、上記両条件下の発現量データから予め発現量変動比を求めておき、この発現量変動量を本装置に入力し、本装置の記憶装置に記憶させてもよい。2仮想の遺伝子クラスタの構築a本発明の遺伝子探索装置においては、この仮想の遺伝子クラスタの構築手段として、ゲノム上での遺伝子の連続情報及び/又は位置番号を含むゲノム上の各遺伝子の位置情報、及び仮想の遺伝子クラスタの構築を実行する仮想の遺伝子構築プログラムが格納される。仮想の各遺伝子クラスタは、上記ゲノム上の各遺伝子の位置情報に基づき、上記仮想の遺伝子クラスタ構築プログラムを実行することにより構築される。すなわち、仮想の遺伝子クラスタは、ゲノムDNA上に連続する遺伝子を同一方向に2個から遺伝子数を一つずつ増やして、想定される遺伝子クラスタに含まれる最大限の遺伝子数になるまで抽出され、かつ該抽出において、抽出する遺伝子の各個数毎に、直鎖状DNAからなるゲノムの場合には該DNAのいずれかの末端から、環状DNAからなるゲノムの場合には任意の遺伝子を起点として順にゲノムDNA上に配列する遺伝子を一つずつずらしながら抽出された各遺伝子群からなるが、このような仮想の遺伝子クラスタを構築するため、上記仮想の遺伝子クラスタの構築プログラムは、本発明装置の記憶装置に記憶されたゲノムDNA上の各遺伝子の位置情報に基づき、以下の処理手段を実行する。その手順を図3に示す。なお、図3中、Nは、仮想の遺伝子クラスタを構成する遺伝子数を表す。1ゲノム遺伝子が直鎖状ゲノムの場合、aゲノムDNAの一方の末端に位置する遺伝子を起点として、他方の末端方向に、順次、ゲノムDNA上に連続する遺伝子を同一方向に2個から一つずつ増やしてN1、想定される遺伝子クラスタに含まれる遺伝子数の最大限nclになるまで組み合わせ、起点とした遺伝子を含み、かつ遺伝子の個数の異なる複数の遺伝子群を構成する。b起点を、順次、他方の末端方向に一遺伝子ずつずらせながら、上記aと同様の処理を行い、新たな起点遺伝子を含みかつ遺伝子の個数が異なる複数の遺伝子群を構成し、aの遺伝子群と併せて、複数の遺伝子の組み合わせた遺伝子群からなる仮想の遺伝子クラスタを構築する。2ゲノム遺伝子が環状の場合、ゲノムDNA上の任意の遺伝子を起点として、上記1a及びbと同様の処理を順次行い、最初に起点とした遺伝子が起点となる時点で処理を終了する最初に起点として遺伝子に基づく仮想の遺伝子クラスタの構築は再度行わない。。上記仮想の遺伝子クラスタの構築においては、仮想の遺伝子クラスタが複数の遺伝子から構成される点で、遺伝子2個から一つずつ増やす手法が採用されるが、本発明は、遺伝子1個から一つずつ増やす手法を排除するものではない。すなわち、この場合、遺伝子一個の場合が構築される仮想の遺伝子クラスタに混入することになるが、本発明においては、この混入した遺伝子を含む遺伝子2個以上の組み合わせからなる仮想の遺伝子遺伝子クラスタが必ず構築され、また、仮想の遺伝子クラスタのスコアは、組み合わせた各遺伝子の発現量変動比の合算であるから、ゲノム中に標的遺伝子が存在した場合、この標的遺伝子単独のスコアに比べ、これを含む仮想の遺伝子クラスタのスコアは少なくとも同等以上になり、上記混入は実質的な問題ではない。したがって、仮想の遺伝子構築において、遺伝子2個から一つずつ増やす手法を含む限り、遺伝子1個から一つずつ増やしたとしても、本発明に包含される。なお、上記ゲノム上の各遺伝子の位置情報は、マイクロアレイによる発現量データにも同様な位置情報を付与することにより、以下の仮想の遺伝子クラスタのスコアリングの際の遺伝子照合に用いられる他、特定の遺伝子の重み付けあるいは特定の遺伝子で仮想の遺伝子クラスタを選定する際の識別手段ともなる。一方、上記のようにゲノム上の各遺伝子の各位置情報を格納しなくとも、例えば予めマイクロアレイ上の各DNAをゲノムDNA上の配列順に整列させておくことにより、ゲノムDNA上の遺伝子の配列順に従いそのまま入力して、入力された遺伝子の順序を遺伝子位置番号として記憶し、該位置番号を用いて仮想遺伝子クラスタを構築することもできる。この仮想の遺伝子クラスタ構築プログラムは、コマンドに基づき組み合わせる遺伝子数の上限を設定可能なようにしてもよい。上限は検索対象の遺伝子クラスタにもよるが、ほとんどの場合、最大30個で充分である。このようにして構築された仮想の遺伝子クラスタは記憶部に記憶される。構築される仮想の遺伝子クラスタは、例えば、次のように、ゲノムDNA上に配列する遺伝子がA~Jの10個である場合、以下の遺伝子群からなる表1。したがって、この場合、仮想の各遺伝子クラスタの構築数は45個であるが、これらの各遺伝子クラスタは本発明の装置内でデータ処理に基づき構築されるだけであって、実験によって実際に構築されるものではない。なお、実際のゲノムDNA上の遺伝子数は、麹菌の場合、外部データベースDOGANに登録されているもので12084個であり、これよりも遺伝子の定義を緩めてDNAマイクロアレイのプラットフォーム作成に使用されたものの場合14032個である。このうち連続していることが判明しているゲノム上の領域から、仮想の遺伝子クラスタを構築する。抽出する遺伝子の数の最大限は論理上ゲノム中の遺伝子の数とすることができるが、想定される遺伝子クラスタサイズの最大限の遺伝子数でよく、実際問題として、遺伝子クラスタを構成する遺伝子の数は、最大でも30個程度であり,これを超えて遺伝子クラスタを構築する必要は通常ない。3仮想の遺伝子クラスタのスコアリング上記のように構築された仮想の各遺伝子クラスタは、本発明装置のスコアリング手段により、スコアリングされる。該スコアリング手段は、本装置の処理演算部に格納されているスコアリングプログラムにより実行される図4。該プログラムは、記憶部に記憶されている、ゲノムDNA上の各遺伝子の発現量変動比データと上記構築された仮想の遺伝子クラスタ情報を呼び出して、各仮想の遺伝子クラスタを構成する遺伝子と各発現量変動比データの遺伝子を照合して、以下の計算式aを使用して各遺伝子の発現量変動比を合算して、各仮想の遺伝子クラスタのスコアを算出する手段を実行する。得られた各仮想の遺伝子クラスタのスコアは、出力されるか及び/又は記憶部に記憶される。計算式aなお、上記式の定義中、全ての仮想の遺伝子クラスタに含まれる全遺伝子とは、全ての仮想の遺伝子クラスタを構成するために抽出されたゲノムDNA上の全ての遺伝子をいう。本発明によれば、このようにして得られた一群の仮想の遺伝子クラスタのスコアに対する出現頻度分布をみる場合、全体としては大凡正規分布となるが、このような全体のスコア分布から離れて存在する仮想の遺伝子クラスタが存在すれば、少なくとも標的の遺伝子クラスタと対応していると判定できる。すなわち、この仮想の遺伝子クラスタは、該クラスタ中の少なくとも2つの遺伝子が、代謝物産生誘導などの生理状態変化条件下で協働した結果、発現変動量の総量であるスコアが増大したものであり、標的の遺伝子クラスタとみなすことができ、この仮想の遺伝子クラスタ中の遺伝子は少なくとも実際の遺伝子クラスタ中に存在する代謝物産生などの生理状態変化に関与する遺伝子であると同定することができる。さらに、例えば、仮想の遺伝子クラスタ中の遺伝子及び必要に応じ代謝産物の産生機構を検討すれば、直接代謝産物の産生に関与する標的遺伝子のみではなく、未知の機能を有する遺伝子の発見も期待でき、さらに代謝物産生機構の全体像も明らかにすることができる。4アノテーション付与本発明の遺伝子探索装置においては、入力されたゲノム上の各遺伝子にアノテーションを付与する手段を設けることができる。アノテーション付与は、ゲノム上の遺伝子が、標的とする遺伝子機能を有すると推定される場合、あるいは標的とする遺伝子機能を有する可能性が低いか若しくはその可能性がないと推定できる場合において行う。このようなアノテーション付与は、探索対象ゲノム上の各遺伝子の塩基配列情報等を基に、記憶部に記憶されたゲノム上の各遺伝子の位置情報中の遺伝子について行う。このアノテーション付与手段は、装置使用者が、検索探索対象ゲノム上の遺伝子について予め相同検索あるいはモチーフ検索等の結果に基づき、上記記憶されたゲノム上の各遺伝子の位置情報中の遺伝子を逐一指定し、この指定された遺伝子にアノテーションが付与されるように構成しても良いが、ゲノム上の遺伝子の数は極めて多数であり、上記モチーフ検索を行う市販のソフトウェアを付属のモチーフ情報とともに本発明装置に格納するか、あるいは該ソフトウェアをモチーフ情報とともに格納した外部コンピューターと接続可能にすることが好ましい。これにより、探索対象ゲノム上の各遺伝子の塩基配列情報を本発明装置の入力手段に入力するかあるいは外部コンピューターに入力することにより、期待される機能に対応したモチーフについて検索を行い、アノテーション付与する遺伝子を自動で選定することができる。また、他のアノテーション付与手段として、上記モチーフ検索により探索対象のゲノム上の全遺伝子にてアノテーションを付与した後に、付与されたアノテーションの種類遺伝子機能から、期待される機能と一致する遺伝子を選定してもよい。選定された遺伝子は、本発明装置の記憶部に記憶されたゲノム上の遺伝子の位置情報における各遺伝子と照合される。このようなシステムによれば、アノテーション付与を研究者の手を煩わすことなく自動で行うことができる。アノテーション付与は、機能が同様なゲノム遺伝子に付与しても良いし、機能の種類が異なる複数種の遺伝子に付与しても良い。機能が異なる複数種のゲノム遺伝子にアノテーションを付与する場合には、ゲノム遺伝子の機能毎に識別可能なように付与する。アノテーションによる選定の対象となる遺伝子は、例えば、2次代謝物産生に関与する遺伝子クラスタあるいはその中の遺伝子を標的とする場合、ゲノムDNAの配列中、12次代謝に関与していると想定される酵素種に属する酵素遺伝子、2トランスポーター遺伝子、3転写因子をコードする遺伝子を選定可能である。5アノテーションを付与した場合の仮想の遺伝子クラスタスコアリング11本発明の遺伝子探索装置においては、各仮想の各遺伝子クラスタのスコアリングにおいて、各仮想の遺伝子クラスタ中に該当する機能に関するアノテーションが付与された遺伝子がある場合、その遺伝子についての発現量変動比に重み付けを実行する重み付けスコアリングプログラム図5を格納することができる。これにより、アノテーションに基づき選定されたゲノム遺伝子についての発現量変動比は重み付けがなされ、仮想の各遺伝子クラスタのスコアリングがなされる。この重み付けスコアリングプログラムは、各仮想の遺伝子クラスタのスコアリングにおいて、アノテーションに基づき選定された遺伝子について、以下の計算式による重み付け計算手段を実行する他は、上記3のスコアリングプログラムと同様の手段を実行する。1重みwの設定は、上記標的とする遺伝子機能を有すると推定される場合、1を超えるように設定し、標的とする遺伝子機能を有する可能性が低いか若しくは可能性がないと推定できる場合は、0以上1未満になるように設定する。標的とする遺伝子機能を有するか、あるいはその可能性が低いかの推定は上記と同様に既知の遺伝子との相同性あるいはモチーフ等により判別すればよい。2一方、本発明の遺伝子探索装置においては、上記重み付けの代わりに、構築された仮想の遺伝子クラスタの中から、アノテーションに基づき選定された遺伝子を含む仮想の遺伝子クラスタを選出し、この選出された仮想の遺伝子クラスタについてのスコアリングを実行するプログラムを格納してもよい。このような手段は、上記標的とする遺伝子機能を有すると推定される場合に有効であり、例えば上記した2次代謝物産生に関与する機能遺伝子の探索等においては、特に有効である。これにより、スコアリングする仮想の遺伝子クラスタの数を削減できスコアリング時間を短縮することができる。例えば、上記表1において、アノテーション付与された遺伝子がAとCである場合、遺伝子AとCを含む仮想の遺伝子クラスタのスコアリングは、合計8個ですむ。また、この手法により選定された仮想の遺伝子クラスタは、結果として、後記する5アノテーションに基づき遺伝子を選定した場合の仮想の遺伝子クラスタスコアリング2において示される、選定された機能遺伝子により構築された仮想の遺伝子クラスタと同様なものになる場合があるが、この手法による場合、後記する一度Aの手法による網羅的な仮想の遺伝子クラスタ群を構築しておけば、自由に標的とする遺伝子あるいはこれを含む遺伝子クラスタの機能を変更でき、種々の機能選択的な遺伝子解析が容易に行える点で有利である。また、該当するアノテーションが付与されなかった遺伝子のスコアを考慮することもできるため、機能未知の遺伝子の影響が大きい場合などに柔軟に対応できる。本発明は、ゲノムDNA上の複数の遺伝子を組み合わせて仮想の遺伝子クラスタを構成し、これら複数の遺伝子の生理状態変化条件下の発現量変動比を合算して各仮想の遺伝子クラスタをスコアリングし、これに基づき、まず、標的の遺伝子クラスタを探索する装置に関する。スコアリングされてスコアが高いものが得られた場合は、仮想の遺伝子クラスタに含まれる複数の遺伝子が協働した結果であり、各遺伝子単独の発現量変動比スコアをみるよりも、全体のスコア分布に対する特異性がより鮮明となる。一方、従来のように一つ一つの遺伝子の発現変動量のみから有用遺伝子を検出する場合には、正解の遺伝子であっても、全体のスコア分布中に吸収されてしまい、高い順位の遺伝子であっても、目的の遺伝子であるか否かの遺伝子破壊実験等の検証を必要とする。これに加え、上記したように重み付けをした遺伝子についての発現量変動比は、仮想の各遺伝子クラスタのスコアリングにおいて、他の遺伝子の発現量変動比と合算され、上記標的とする遺伝子機能を有すると推定される遺伝子を含む仮想の各遺伝子クラスタのスコアはより高くなり、反対に標的とする遺伝子機能を有する可能性が低いか若しくは可能性がないと推定される遺伝子を含む仮想の遺伝子クラスタのスコアはより低くなり、全体のスコア分布との乖離が明瞭となる。したがって、これにより、標的とする遺伝子機能を有する遺伝子あるいはこれを含む遺伝子クラスタの探索がより効率的になる。5アノテーション付与により遺伝子を選定した場合の仮想の遺伝子クラスタスコアリング2一方、本発明の遺伝子探索装置においては、ゲノム上近傍に存在する遺伝子について、アノテーションの種類毎に、機能遺伝子を一種以上、好ましくは2種以上抽出するか、あるいはこれら遺伝子が含まれるようにゲノムDNA上の遺伝子を抽出して仮想の遺伝子クラスタとする、仮想の遺伝子をクラスタの構築手段を設けることができる。これによればスコアリングの対象となる遺伝子クラスタの数を大幅に減らすことができ、処理データ量が少なく簡便であり、2次代謝産物の産生に関与する遺伝子クラスタ及び該クラスタ中の2次代謝産物産生遺伝子の探索に特に適している。このような処理を実行するプログラム図6は、記憶部に記憶されたゲノム上の遺伝子の位置情報に基づき、アノテーションにより選定された遺伝子について、ゲノムDNA上において近傍に位置することを条件として、上記選定された遺伝子を1種以上、好ましくは2種以上抽出し、これら抽出した遺伝子により、仮想の遺伝子クラスタを構築するかあるいは少なくともこれら選定された遺伝子が含まれるようにゲノム遺伝子を抽出して仮想の遺伝子クラスタとする手段を実行する。例えば、これら仮想の遺伝子クラスタの構築において機能遺伝子のみを組み合わせる場合、ゲノム上に配列する遺伝子数で上限30程度の範囲にある遺伝子であり、本発明装置においては、組み合わせる機能遺伝子の範囲の入力、設定手段を設けるとともに、上記プログラムはこれに基づき組み合わせる機能遺伝子を選択する。該プログラムは、上記遺伝子に付与されたアノテーションの種類と上記記憶部に記憶された上記ゲノム上の各遺伝子の位置情報中の位置番号により組みあわせる遺伝子を選択する。二次代謝産物の産生に関与する遺伝子クラスタ及び該クラスタ中の二次代謝産物産生遺伝子を探索する場合、アノテーションによる選定は、例えば、ゲノムDNAの配列中、12次代謝に関与していると想定される酵素種に属する酵素遺伝子、2トランスポーター遺伝子、3転写因子をコードする遺伝子を対象に行う。例えば、以下のように、ゲノムDNA上に配列する遺伝子がA~jの10個である場合、は該当する酵素種をコードする遺伝子、はトランスポーター遺伝子仮想の遺伝子クラスタは、AC及びGJとにより構成してもよく、また、これら遺伝子が含まれるように、ABC及びGHIJで構成してもよく、さらにABCDEあるいはFGHIJのように、各仮想の遺伝子クラスタが一定数の遺伝子により構成されるようにゲノムを分割して、各仮想の遺伝子クラスタを構成しても良い。ゲノム配列中の1二次代謝に関与していると想定される酵素種に属する酵素遺伝子、2トランスポーター遺伝子、3転写因子をコードする遺伝子の判別は、既知の同じ酵素種の遺伝子との相同性あるいはモチーフ等により判別すればよく、例えば、各仮想の遺伝子クラスタ中の遺伝子配列中にこれらの遺伝子が存在するか否かは、上記酵素種に属する酵素、トランスポーター、転写因子の各アミノ酸配列に特有なモチーフと共通するアミノ酸配列をコードする塩基配列が存在するか否かにより識別でき、それぞれ種類の異なるアノテーションが付与されるが、このような識別、アノテーション付与は、上記4アノテーション付与において述べた手法を用いればよい。上記1の酵素遺伝子の判別において、酵素種は、二次代謝物の化学構造、前駆体、関与しうる補酵素、化学的物理的性質、既知の酵素反応の事例、生産効率速度等からその産生反応を推測し、関与する酵素種を想定するが、この酵素種の想定においては、実際にその反応に関与したであろう特定酵素のレベルまで想定しなければならないというわけではなく、該反応に関与することがより確実なレベルの酵素種でよい。例えば、オキシゲナーゼに属する酵素であることは分かるが、その下位概念の酵素種まで特定できないときは、オキシゲナーゼのレベルを酵素種として選定して、ゲノム上の各遺伝子の配列を探索し、その範疇に属する全てのゲノム遺伝子のそれぞれを、各仮想の遺伝子クラスタの構成遺伝子とすればよい。ただし、下位概念の酵素種を選定できた場合には、探索対象とする仮想の遺伝子クラスタの範囲が狭まる可能性があり、その分探索が効率的となる。また、二次代謝物産生反応において複数の酵素が関与していると想定できる場合には、その複数の酵素種を選定することも可能である。また、このような機能遺伝子を組み合わせた仮想の各遺伝子クラスタのスコアリングも、上記1aの計算式による計算において、選定された機能遺伝子の発現量変動比のみを用いて行うだけでよく、このような設定を行うだけで、上記3仮想の遺伝子クラスタのスコアリングにおいて説明したスコアリングプログラムを使用できる。すなわち、この場合、計算式1aの定義は、上記式中、Mは各仮想の遺伝子クラスタのスコア、mはスコアリングされる仮想の各遺伝子クラスタに含まれるアノテーション付与に基づき選定された各遺伝子の発現量変動比、m-は全ての仮想の遺伝子クラスタに含まれるアノテーション付与に基づき選定された全遺伝子の発現量変動比m値の平均、sは全ての仮想の遺伝子クラスタに含まれるアノテーション付与に基づき選定された全遺伝子の発現量変動比m値の標準偏差を表す。になる。6スコアリング結果の表示本発明の遺伝子探索装置においては、上記したように仮想の遺伝子クラスタンスコアリングにより算出されたスコアあるいはこれを加工した形態で、画面表示及び/または紙等の表示媒体に出力する手段を設けることができる。表示手段としては、例えば、スコアの高い順に仮想の各遺伝子クラスタを表示したり、あるいは仮想の遺伝子クラスタのスコアの分布状態を表すグラフ等があげられ、さらに仮想の遺伝子クラスタに含まれる遺伝子を表示する手段を設けることもでき、これらに基づき、仮想の遺伝子クラスタを選定することができる。一方、スコアが高く全体分布と乖離している仮想の遺伝子クラスタは、実際に存在する標的の遺伝子クラスタに一致ないし対応する仮想の遺伝子クラスタの可能性が高い。以下に示す7あるいは8の手段は、仮想の各遺伝子クラスタのスコアの全体のスコアからの乖離の程度をみることにより、標的の遺伝子クラスタ候補を選定するか、あるいはさらに候補の絞り込みを行うための手段であり、本発明装置にこれら7あるいはさらに8の手段を設けて、乖離度を示す判定値Iχ、判定値IIυあるいは絞り込み結果b値を上記選定された仮想の遺伝子クラスタ及びその中に含まれる遺伝子とともに表示することができる。これらにより、標的の遺伝子クラスタ及び該遺伝子クラスタに含まれる標的の遺伝子を特定できる。7全体のスコア分布からの乖離の程度の算出上記スコアリング結果の表示から、標的とする遺伝子クラスタあるいはその中の標的遺伝子を見いだすことは十分可能と考えられるが、より客観性及び効率性を高めるため、本発明装置においては、さらに仮想の遺伝子クラスタ全体のスコアの分布から乖離して存在するスコアを有する仮想の遺伝子クラスタを、標的の遺伝子クラスタ候補として選定する手段を設けることができる。本発明の装置における、このような、仮想の遺伝子クラスターのスコアの、全体分布からの乖離度を判定する手順について、図7に示す。この候補選定手段には、仮想の遺伝子クラスタ全体のスコアの分布からの乖離の程度を示す判定値を算出する、乖離度判定プログラムが格納されている。この乖離度判定プログラムは2種あり、上記仮想の遺伝子クラスタのスコアリングプロセスにより算出されたスコアに基づき、例えば、以下の計算式bあるいはcに基づき、それぞれ判定値Iχあるいは判定値IIυを算出し、判定値Iχあるいは判定値IIυが、例えば、予め設定した一定値以上を示した仮想のクラスタを標的の遺伝子クラスタの候補として選定する手段を実行する図7。選定結果は判定値とともに出力されるが、併せて乖離度の平均値等も出力するようにしてもよい。これら2種のプログラムは、本発明装置にともに格納しても良いが、そのうち1種のみを格納しても良い。計算式b上記計算式b中のスコアMの出現頻度は、仮想の遺伝子クラスタの全てを含む集団における各スコアの出現頻度Pの累計を1としたときの値であるため、1を超えることはなく、したがってlogPは正になることはない。また、出現頻度が低いものほどlogPは-∞に近づくため、頻度の低いスコア値を持つ遺伝子クラスタほどlogPの絶対値は大きくなる。したがって、上記計算式bにおいては、logPと仮想の各遺伝子クラスタのスコアを掛け合わせて-1を乗算することにより、頻度が低くかつスコアの高いものが、より大きな判定値Iχを持つこととなる。逆に、頻度が低くかつスコアの低いものは、より小さな負の判定値Iχを持つ。上記計算式bによれば、判定値Iχが0を超え、その絶対値が高い値を示す仮想の遺伝子クラスタは、仮想の各遺伝子クラスタのスコアに対する出現頻度分布から離れており、その絶対値が高い判定値Iを示した仮想の遺伝子クラスタを標的の遺伝子クラスタあるいは標的の遺伝子クラスタに対応する候補として選定することができる。計算式cこの判定値IIυは、仮想の各遺伝子クラスタのスコアについて、仮想の遺伝子クラスタ全体の平均スコアからのずれを、上記標準偏差の実数倍で割ったものを次元数d’乗したもので、正規分布様のスコアに対する出現頻度分布から乖離するスコアを有する仮想の遺伝子クラスタにおいて大きな値となる。上記式中d’は任意に設定できる正の偶数たる次元数であり、値が大きくなるほど平均スコアからの隔たりが強調されることになる。あまり大きくしすぎると、平均スコアから大きく外れたものの値が強調されて相対的にそれ以外の値が小さくなるため、通常2または4に設定する。外れたものをより鋭敏に検出したい場合は、6以上の偶数とする。また式中のaは外れ度を表す係数で、この値を調節することにより、上記正規分布様分布からどの程度乖離したものをとるかを調節することができる。1を超えて大きく設定するほど、平均スコアから大きく外れたもの以外のυ値はゼロに近付くため、このa値は通常1~2に設定する。逆に1未満の場合、より外れ方の小さなものも拾うことができる。この計算式cによる場合も上記判定値Iと同様に、υが0を超え、高い値を示す仮想の遺伝子クラスタを標的の遺伝子クラスタあるいは標的の遺伝子クラスタに対応する候補として選定することができる。8遺伝子クラスタ候補の絞り込み上記計算式b、cにより算出された判定値χあるいはυにより、標的の遺伝子クラスタ候補となった仮想の遺伝子クラスタの数が多く、さらに候補を絞り込みたい場合に備えて、本発明装置においては、遺伝子クラスタ候補絞り込み手段として以下の計算式dによる計算を行う、候補絞り込みプログラムを格納することができる図8。すなわち、各仮想の遺伝子クラスタについて、判定値IおよびIIの積を取った値について、bが100未満の仮想のクラスタを少なくとも除外することにより、標的の遺伝子クラスタ候補をさらに絞り込むことが可能である。計算式d上記計算式d中、bはどの程度の遺伝子クラスタ候補を絞り込むかを決定するための閾値であり、bを大きくとるほど候補の絞り込み効果がより高くなる。また小さくとるほど多くの候補遺伝子クラスタを選択することができる。bの値の設定は対象とする生物種や培養条件に依存する。すなわち、候補遺伝子クラスタが強くかつ多く発現している系であれば値を高くする必要があるが、逆に発現強度が弱くかつ数が少なければ値を低くしなければ候補遺伝子が出現しない。前者の場合、例えば5000~10000あるいは10000~30000の範囲内の任意の数値に設定し、後者の場合、通常100以上、例えば1000~2000、あるいは2000~5000の範囲内の任意の数値に設定する。9本発明の装置を用いて正解が得られなかった場合一方、本発明の手法を行った結果、仮に仮想の遺伝子クラスタ全体のスコア分布から乖離したスコアの遺伝子クラスタが見いだされなかった場合、設定する生理状態変化条件、重み付けするゲノムDNA上の遺伝子の選定、あるいは上記Bの手法による仮想の遺伝子クラスタ構築のためのゲノムDNA上の遺伝子の選定等の探索条件設定に問題点がある。したがって、このような場合には、探索条件を再設定して、バックグランドの分布から離れたスコアの遺伝子クラスタが見つかるまで、上記した遺伝子クラスタの探索法を繰り返し行えばよい。すなわち本発明においては、得られたデータのみから、探索条件設定の問題点を把握できる。これに対して、上記したような従来法の場合には、もともと正解の遺伝子であっても、遺伝子全体の発現量についての分布中に埋もれてしまうので、得られたデータからでは正解か否かは不明であり、結果的に無意味かもしれない検証実験を繰り返さなければならない。B遺伝子クラスタ予測装置一方、本発明における上記仮想の遺伝子クラスタの構築手段及びそのスコアリング手段を用いた他の態様として、標的とする遺伝子クラスタの有無及び標的とする遺伝子クラスタが存在する場合のサイズクラスタを構成する遺伝子数;nclを推定する装置以下、遺伝子クラスタ予測装置という。を挙げることができる。本発明の装置における、この遺伝子クラスタ予測装置の概要を図9に示す。この遺伝子クラスタ予測装置においては、まず、生物細胞の生理状態変化を生じる条件とコントロール条件下において生じたゲノムDNA上に配列する遺伝子の発現量変動比を合算し、仮想の遺伝子クラスタのスコアとするが、ゲノムDNA上に配列する各遺伝子の発現量データの入力、発現量変動比の計算、仮想の遺伝子クラスタの構築、及び仮想の各遺伝子クラスタのスコアリングの各手段は、上記1~3に記載した手段と同様である。すなわち、この装置は、本発明の上記遺伝子探索装置における、a生物細胞の生理状態変化を生じる条件とコントロール条件下において生じたゲノムDNA上に配列する各遺伝子の発現量を入力する手段、b入力された上記2つの条件下における同一遺伝子の発現量の比を算出する発現量変動比算出手段、cゲノムDNA上に配列する各遺伝子の発現量変動比を複数の遺伝子により構築された上記仮想の遺伝子クラスタ単位の発現量変動比として合算し、仮想の遺伝子クラスタ単位毎にスコアリングする手段を有するものであって、仮想の遺伝子クラスタの構築手段が、ゲノムDNA上に連続する遺伝子を2個から遺伝子数を一つずつ増やし想定される遺伝子クラスタに含まれる最大限のゲノム遺伝子数になるまで抽出し、かつ該抽出において抽出する遺伝子の各個数毎に、直鎖状DNAからなるゲノムの場合には該DNAのいずれかの末端から、あるいは環状DNAからなるゲノムの場合には任意の遺伝子を起点として順にゲノムDNA上に配列する遺伝子を一つずつずらしながら抽出された各遺伝子群を仮想の各遺伝子クラスタとする手段であり、スコアリング手段として、以下の計算式aによる計算を行うプログラムを格納している点では、本発明の遺伝子探索装置と共通する。この装置の特徴点は、上記1~3の手段のプロセスを行い、出力された仮想の各遺伝子クラスタのスコアに基づき、d仮想の遺伝子クラスタに含まれる遺伝子数単位毎の遺伝子クラスタ分布判定値εを算出する手段にあり、この手段を実行するプログラムとして遺伝子クラスタ分布判定値ε値算出プログラムを格納した点にある図9。計算式aこの遺伝子クラスタ分布判定値は、以下の計算式eにより求められる。計算式eこの計算式eによれば、仮想の遺伝子クラスタが、実際のゲノムDNAにおいてクラスタを形成していなければ、仮想の遺伝子クラスタ中に含まれる、標的の生理状態変化に関与せず、発現量変動のない遺伝子の影響を受けるので、仮想の遺伝子クラスタのスコアMはサイズ遺伝子数;nclが大きくなるほど平均化され、すなわちスコアの平均値に近づくためサイズの増大に伴いε値は単調減少する図10上から1および3番目の曲線参照。しかし、仮想の遺伝子クラスタのあるサイズのものがクラスタを形成している場合、そのサイズにおいて、分布の偏りεは大きくなり、上記単調減少曲線とはならず、ε値はそのサイズにおいて特異点を形成する図10矢印で示した点参照。したがって、ε値が特異点を形成するか否か、及び特異点を形成した遺伝子クラスタのサイズから、遺伝子クラスタの存在及びそのサイズを推定することができる。具体的には、仮想の遺伝子クラスタの該クラスタに含まれる遺伝子数毎の集計において、ある遺伝子数k個のときのε値εと、その前後数のときのε値ε、εが以下の関係にあれば、標的とする遺伝子クラスタがゲノム中に存在すると判定し、標的遺伝子クラスタに含まれる遺伝子数をk個と予想することができる。本発明の遺伝子クラスタ予測装置は、上記a~dの手段を具備した独立の装置としても構成しても良いが、上記a~cの手段が、本発明の遺伝子探索装置と共通するため遺伝子数単位毎の遺伝子クラスタ分布判定値εを算出する手段を、本発明の遺伝子探索装置にさらに設けて、本発明の遺伝子探索装置に標的の遺伝子クラスタの有無および遺伝子クラスタのサイズ予測機能を付加しても良い。このような予測機能は、本発明の遺伝子探索装置を用いて、選定した複数種の機能遺伝子を組み合わせて仮想の遺伝子クラスタを構築し、そのスコアリングを行う際、予め行う手法として有効である。すなわち、遺伝子クラスタが存在し、そのサイズを予想できれば、予想されるサイズ内に、標的とする酵素種に属する酵素遺伝子、2トランスポーター遺伝子、3転写因子をコードする遺伝子が存在するゲノム配列のみを対象として、上記仮想の遺伝子クラスタとして探索することが可能となる。また、この遺伝子クラスタ予測装置によれば、ある条件下で細胞が何らかの生理的状態変化を起こす場合においては、どのような生理的状態変化であっても変化を対比する条件が設定できれば、その原因遺伝子はもちろんその変化を生じる機構そのものが全く不明である場合においても、その変化原因が遺伝子クラスタ中の遺伝子の連関にあるのか否か、遺伝子クラスタ中の遺伝子の連関による場合該クラスタの遺伝子サイズも容易に予測できる。すなわちこの手法は、生物の生理的変化が、極めて探索の難しい複数の遺伝子の連関によって生じている場合において、その原因が遺伝子クラスタ中の遺伝子の共働によるものであることを明らかにでき、かつそのサイズも予測できる点で極めて有用である。参考例1コウジ酸の産生に必須の遺伝子の同定本参考例は、本発明による遺伝子の探索、同定の有利性を明らかにするため、まず従来法によるアスペルギルスオリゼのコウジ酸産生遺伝子の探索、同定手法について示すものである。アスペルギルスオリゼAspergillusoryzaeの菌株RIB40以下、単にアスペルギルスオリゼと書いた場合にはこの菌株を指すは、以下の組成の液体培地中で、30°C、150回転/毎分の条件下で生育させた場合、コウジ酸を培地中に産生する。500mLのこぶつき三角フラスコ中に250mLの培地を入れ、アスペルギルスオリゼの胞子懸濁液を105-107/mLになるように接種する。培地組成:以下コウジ酸産生培地と呼ぶ10WVグルコース025WVイーストエクストラクトYeastExtract01WVK2HPO4005WVMgSO47H2OpHを60に調整後、オートクレーブにより滅菌する。アスペルギルスオリゼが上記の培養によりコウジ酸を産生することは、コウジ酸と塩化第二鉄とのキレート化合物の生成による赤色の発色により検出することが可能である。また培養の上清などを適宜希釈した試料に、最終濃度10mM程度になるように高濃度の塩化第二鉄溶液を添加した液を作成し、波長500nmの吸光度を測定することで、コウジ酸量の定量測定が可能である。この波長500nmの吸光度は、01~10程度の範囲でコウジ酸濃度に比例する。このような検出法によれば、接種後3または4日目には産生を検出することが可能であり、少なくとも7日目には十分な速度をもってコウジ酸の産生が行われている。またコウジ酸の産生は、上記の産生培地に01W/V以上の硝酸ナトリウムを加えることで阻害される。この硝酸ナトリウムによる阻害は可逆的である。硝酸ナトリウムの添加によって阻害された菌糸を、培地成分の洗浄後、新たに用意した産生条件を満たす培地に移すことによって、菌はコウジ酸の産生を開始する。アスペルギルスオリゼのコウジ酸の産生量が異なる下記に記載の条件からなるC1~C3の3つの系で、ゲノム中にコードされたほとんどの遺伝子の網羅的な発現の解析を、DNAマイクロアレイを用いた実験により比較した。C1上記コウジ酸産生培地で、4日間および2日間生育させた菌体の遺伝子の発現を比較した4日目/2日目。C2上記コウジ酸産生培地で、7日間および4日間生育させた菌体の遺伝子の発現を比較した7日目/4日目。C3上記コウジ酸産生培地に03W/Vの硝酸ナトリウムを添加してコウジ酸産生を阻害した菌体と、上記コウジ酸産生培地で生育させた菌体を比較した。どちらも4日間、30°C、150回転/分の条件で生育させたNO3なしあり。上記各系での菌体の遺伝子発現を、DNAマイクロアレイを用いて解析した結果、系C1~C3のそれぞれにおいて、比較する条件下で培養した2つの菌体での各遺伝子の発現量の比、および発現の強度に相当する値が得られた。各々比較する条件間で、コウジ酸の産生がより顕著な条件で、発現がより顕著になっている遺伝子を抽出するために、以下の手続きにより、候補を抽出した。発現量の比、および発現の強度に相当する値は、それぞれで正規分布に近い分布をするが、値の絶対値は大きな違いがある。この両者を統合して候補を抽出するために、発現量の比、発現の強度は、それぞれで値の正規化を実施した後で、比較した。それぞれ正規化した発現量の比、および発現の強度に相当する値の積を作成した。その積が高いほど、コウジ酸の産生に関係する可能性が高いと考え、それぞれの実験で高い積の値をもつもの上位5つを選び出した表2。アスペルギルスオリゼのDNAマイクロアレイにおける発現比と発現量の積からなるスコア上位遺伝子表2に示した遺伝子は、それぞれの比較する二つの条件下で、コウジ酸の産生条件において顕著に発現が高くなっている遺伝子である。すなわちコウジ酸の産生に必須の遺伝子である可能性が高い遺伝子である。これらの遺伝子について、上位のものから遺伝子欠失破壊実験を行った。ここで上記C1~C3の3つの系は、いずれもコウジ酸の産生量が有意に異なる2つの条件を比較したものである。したがって理想的には、どの系においてもコウジ酸産生に必須の遺伝子が上位に現れると予想した。しかし現実には、3つの系全てにおいて上位にくる遺伝子はなかった。したがっていずれの系においても、上位に来るものはコウジ酸の産生に必須であるか、または各条件に特異的に誘導される遺伝子である可能性の両者を含む。これらの中からコウジ酸の産生に必須の遺伝子を選び出すために、各系において上位に来ている候補遺伝子のいずれかを破壊して変異体を作製し、当該変異体のコウジ酸産生能を解析した。その結果、AO090113000136およびAO090113000138の両遺伝子は、破壊によりコウジ酸の産生が著しく低下することが判明した。上記2つの遺伝子は、他の生物種のゲノム中にある機能が既知の遺伝子とのオーソロガスな関係を持たないため、ゲノム情報から両遺伝子の機能を知ることはできなかった。ただし、該遺伝子のアミノ酸配列には散在する既知の配列モチーフが存在し、機能の概略を予測することが可能であった。AO090113000136の遺伝子は、FAD依存性の酸化還元酵素のモチーフを持っている。これはグルコースからのコウジ酸への変換を考えたとき、変換の過程で複数の酸化還元反応が関係していると予想されていることから、この遺伝子がコウジ酸の生合成における酵素であることを強く示唆している。一方で、AO090113000138は膜輸送に関わる配列モチーフを持っており、Majorfacilitatorsuperfamilyと分類される。コウジ酸の生合成に伴って産生されたコウジ酸が培地中に分泌されることは明確であり、この遺伝子がコウジ酸の産生に必須であることを示唆している。この両遺伝子はゲノム上で近傍に位置する。間には1遺伝子しか存在せず、そのAO090113000137遺伝子のアミノ酸配列は転写因子のモチーフを持つことが判明した。この遺伝子の破壊によってもコウジ酸の産生が著しく低下することが判明した。以上の解析により、AO090113000136、AO090113000137、AO090113000138の3つの遺伝子がコウジ酸の産生に必須の遺伝子であると同定された。本同定過程には、培養条件の検討などを除いて、およそ1年の時間を要した。このようにして同定されたコウジ酸の産生に必須の3つの遺伝子について、系C1~C3におけるDNAマイクロアレイの結果において、その発現量変動比m値が全遺伝子中どの位置に来るかを表3にまとめた。アスペルギルスオリザにおけるコウジ酸産生に必須の3つの遺伝子のスコアm値とその順位また系C1~C3における各分布を図11から図13に示した。表3にあるように、系C2では該当する3つの遺伝子は1位、6位、71位と上位におり、この系のアレイであれば必須遺伝子の同定は比較的容易である。一方系C3では、コウジ酸の産生が顕著に見られるにも関わらず、必須の遺伝子の値は最高でも2658位であるなど、上位には見られない。このアレイをもとにした場合、従来法で遺伝子を特定することは事実上不可能である。その上、必須の遺伝子が分かっていない状況では、どのアレイが正解を与えうるものであるかの判断ですら困難である。ここに示した3つのアレイデータだけを用いて3つの遺伝子がコウジ酸の産生に必須であることを同定することは、上記に示した方法で可能ではあったが、偶然的な幸運の要素が大きく一般性は低い。機能の注釈を元に推定する場合でも、100以上の遺伝子を破壊してみなければ分からない可能性があり、この場合、検証に通常3年程度以上はかかる。実施例1アスペルギルスオリゼにおける遺伝子クラスタスコアリングによるコウジ酸合成遺伝子同定本特許で出願する該当遺伝子の同定手法に従い、本発明装置を使用して、アスペルギルスオリゼのコウジ酸産生関連遺伝子からなる遺伝子クラスタを同定した。この実験に使用した装置は、データの入出力装置、入出力インターフェース、記憶装置、制御演算装置CPUから構成され、上記制御演算装置は、発現量変動比算出部、仮想の遺伝子クラスタ構築部、仮想の遺伝子クラスタのスコアリング部、仮想の遺伝子クラスタの乖離度判定値算出部、遺伝子クラスタ候補絞り込み部、及び遺伝子クラスタ予測部を有し、これら各部には、それぞれ順に、発現量変動比算出プログラム、仮想の遺伝子クラスタ構築プログラム、仮想の遺伝子クラスタスコアリングプログラム、乖離度判定値χおよびυ算出プログラム、候補絞り込みプログラム並びに遺伝子クラスタ分布判定値ε算出プログラムが格納されている。また、これら各部での計算は、Linuxオペレーティングシステム上で、フリーソフトウェアR、およびプログラム言語Perlを用いて行った。DNAマイクロアレイのデータは、参考例1と同様のものを使用した。すなわち以下の、コウジ酸を産生する培養条件を分子に、コントロールとなる培養条件を分母にして測定した、以下のC1~C3の系における二色法データである。C14日目/2日目C27日目/4日目C3NO3なしありこれらは各々、産生条件と非産生条件の下で生育させた菌体からmRNAを取り出し、それぞれ色素でラベリングした後でアレイ上のオリゴDNAにハイブリダイズすることによりデータを得て、そこから各遺伝子の発現量変動比にしたがって遺伝子クラスタの判定値χを算出した図16。具体的には、本発明の装置に記憶されている系C2における各仮想の遺伝子クラスタのスコアに、仮想の遺伝子クラスタの乖離度算出部に格納された仮想の遺伝子乖離度判定プログラムのうちχ値算出プログラムを適用し、計算式bにしたがって、各仮想の遺伝子クラスタについての判定値χを算出した。なお、図16中の各折れ線は、仮想の遺伝子クラスタ構築において起点となる遺伝子が共通する遺伝子サイズ1~30の各仮想の遺伝子クラスタの判定値を結んだものである図17、18、21、23、30~32、35~37も同様。。ここで、ncl=1のときの値がncl=2のときの値より大きい仮想の遺伝子クラスタは、本手法におけるクラスタスコアリングによってスコアを上げているわけではないので、該当しない。またncl=1のときの値が負の仮想の遺伝子クラスタは、本手法におけるクラスタスコアリングにおいて、そのスコアの上昇に寄与しないため、該当しない。そこで図13においては、これらのものを除外してある。図16にあるように、1つの仮想遺伝子クラスタがncl=3で極大かつ最大値を取ることが分かる。この仮想の遺伝子クラスタは、コウジ酸産生に必須の3つの遺伝子、AO090113000136、AO090113000137、AO090113000138のみを含んだものであった。この結果は参考例の結果と一致し、上記遺伝子クラスタ分布判定値ε算出による予測結果が正しいことが分かる。また、判定値χによって、標的とする遺伝子クラスタ及び該クラスタに含まれる遺伝子も同定可能なことが明らかとなった。続いてもう一つの遺伝子クラスタ判定値であるυを、上記と同様の各仮想の遺伝子クラスタのスコアに、各仮想の遺伝子クラスタについて、仮想の遺伝子クラスタの乖離度算出部に格納された遺伝子乖離度判定プログラムのうちυ値算出プログラムを適用し、計算式cにしたがって算出した図17。ここでも、χ値と同様、ncl=1のときの値がncl=2のときの値より大きい仮想の遺伝子クラスタは除外してある。ここで次元数d’は2、係数aは1を採用した。すると図のように、ncl=3のとき極大かつ最大値をとる仮想遺伝子クラスタが一つある。これはχ値のときと同様、コウジ酸産生に必須の3つの遺伝子のみを含む遺伝子クラスタである。したがって、判定値υによっても標的とする遺伝子クラスタ及び該クラスタに含まれる遺伝子が同定可能であることが明らかとなった。こうして得られたχ値およびε値に対し、遺伝子クラスタ絞り込み部に格納された候補絞り込みプログラムを適用して、計算式dに従って二つの値の積から遺伝子クラスタ評価値を算出した図18。図18をみると明らかなように、5000以上の値で最大値を取る仮想遺伝子クラスタが極めて明確に一つ存在しており、ncl=3で極大値をとっている。これがコウジ酸産生に必須の3つの遺伝子、AO090113000136、AO090113000137、AO090113000138のみを含んだものである。このように本発明の手法および装置を用いることで、標的とする生合成遺伝子を同定することができた。また計算式dにおける閾値bを例えば2000とすれば、該当する遺伝子クラスタは4個しかなく、実験系による検証する場合においても容易に行える数値である。χ値図16とυ値図17を乗算することで、それぞれにおいて存在していた多くのピークがキャンセルされ、探索対象に該当するもののみ高い値を示している。以上から、本発明の手法および装置は、DNAマイクロアレイデータのみを用いて、ゲノム上に集合して機能を果たす生合成遺伝子の探索、同定を可能とする実効的な手段であることが示された。実施例2アスペルギルスオリゼにおけるアノテーション機能注釈による重み付けを行った場合の仮想の遺伝子クラスタスコアリングによるコウジ酸合成遺伝子の探索アスペルギルスオリゼのコウジ酸産生関連遺伝子からなる遺伝子クラスタを同定することを目的として、予測される機能に関連した注釈のついた遺伝子のm値に重み付けをした後、該当遺伝子の同定を行った。この実験に使用した装置は、上記実施例1に記載した装置と基本的に同様であるが、アノテーションによる遺伝子選定部、選定された遺伝子に対する発現量変動比の重み付与部を有している点で異なる。コウジ酸産生に必要な機能は以下の3つの機能を選出した。膜輸送体:transporterまたはmajorfacilitator転写制御因子:transcription酸化還元酵素:oxidoreductaseまたはdehydrogenaseなお、上記英単語は、アノテーションによる遺伝子選定に用いたキーワードである。これらを選出した理由として、コウジ酸の生合成がグルコースから酸化により変換されていると推定されること、産生されたコウジ酸の膜輸送による培地中への分泌に膜輸送体が必要と推定されること、および関与する遺伝子の転写制御に転写因子が必要と推定されることが挙げられる。Aアノテーション機能注釈による重み付けとクラスタスコアリング一般に利用されているアノテーション推定ソフトウェアシステムの一つであるInterproscanhttpwwwebiacukToolsInterProScanを用いて、アスペルギルスオリゼのゲノムDNA上の各遺伝子についてアノテーションを付与し、その結果付与されたアノテーションに基づき、上記3つの機能に該当する遺伝子を選出した。具体的には、該各遺伝子についてのアノテーションデータを本装置の入力装置に入力し、記憶装置に記憶した。記憶したアノテーションデータのデータを呼び出し、上記3種の機能を有する遺伝子を機能遺伝子選択部の選択プログラムを適用して選定した。なお、選定は各遺伝子について付与されたアノテーション中に上記3つの機能群に対応する英単語が含まれるかどうかで行い、その結果、該当した遺伝子は5179個のうち709個であった。続いてこれらの該当する注釈のついた遺伝子に対して、実施例1に記載した3つのアレイ測定の系C1~3のそれぞれについて、その発現量変動比m値に正規化後重みw=20を積算した後、計算式aにしたがってncl=1~30でクラスタスコアリングを行い、各仮想の遺伝子クラスタのM値を得た。具体的には、このように選定された遺伝子についての発現量変動比は、重み付与部により重み付け[数2]参照がなされ、この重み付けされた発現量変動比を用いて、各仮想の遺伝子クラスタのスコアが算出される。アノテーションに基づき選定された遺伝子の発現量変動比に重み付けする以外は、仮想の遺伝子クラスタの構築プログラム、スコアリングプログラム自体、実施例1と相違しない。この実験においては、発現量変動比m値に正規化後、アノテーションにより選定された遺伝子の発現量変動比には、重みw=20を積算し、仮想の遺伝子クラスタのスコアリング部に格納されたスコアリングプログラムを適用し、計算式aにしたがってncl=1~30でクラスタスコアリングを行って各仮想の遺伝子クラスタのスコアM値を得た。なお、算出された各仮想の遺伝子クラスタのスコアは、本発明装置の記憶装置に記憶した。図19は、上記算出された各仮想の遺伝子クラスタスコアのヒストグラムである。左の拡大図を図14と比較すると、重み付けによってより高いスコアが出現したために、相対的にゼロを中心とした山型の分布がより尖ってみえ、かつ山の中心がより左側にずれていることが分かる。Bデータ判定続いて計算式eにしたがって、系C1~3におけるスコア分布評価値εを算出した図20。具体的には、上記Aにより算出され、記憶された各仮想の遺伝子クラスタのスコアを呼び出し、遺伝子クラスタ予測部に格納されている遺伝子クラスタ分布判定値ε算出プログラムを適用し、計算を行った。ここで実施例1と同様に、仮想遺伝子クラスタ数nは5179、次元数dは6を採用した。図20にあるように、系C1およびC3ではε値は基本的に単調減少しているのに対し、系C2ではncl=3のときε値は大きく増加し極大値を示す。その値は実施例1におけるもの図15の10倍以上であった。すなわち機能の注釈による重み付けによって、より機能限定的に高精度で該当遺伝子クラスタの存在及びその遺伝子数が予測可能であることを示している。この実験によって、系C2のマイクロアレイデータ中に標的とする遺伝子クラスタによるものと推定されるデータが存在することが強く示唆されたため、続いてC2のDNAマイクロアレイデータを用いて以下の検証、同定実験を行った。C遺伝子クラスタ判定上記Aで得られた系C2についての注釈重み付け後の仮想に遺伝子クラスタのスコアM値から、計算式bにしたがって遺伝子クラスタ判定値χを算出した図21。具体的には、記憶されている系C2における各仮想の遺伝子クラスタのスコアを呼び出し、仮想の遺伝子クラスタの乖離度算出部に格納された仮想の遺伝子乖離度判定プログラムのうちχ値算出プログラムを適用し、計算式bにしたがって、各仮想の遺伝子クラスタについての判定値χを算出した。なお、図21においても、実施例1の図16と同様に、ncl=1のときの値がncl=2のときの値より大きい仮想の遺伝子クラスタ、及びncl=1のときの値が負の仮想の遺伝子クラスタは除外している。結果は、図21に示されるように、実施例1と同様、1つの仮想遺伝子クラスタがncl=3で極大かつ最大値を示した。これがコウジ酸産生に必須の3つの遺伝子、AO090113000136、AO090113000137、AO090113000138のみを含んだものである点は、実施例1図16と同様だが、ここでは上位のχ値が重み付けによって図16におけるものより2倍程度高くなり、他との差が広がっている。すなわち注釈による重み付けによって、機能に即した該当遺伝子クラスタの検出精度が向上しているといえる。続いて遺伝子クラスタ判定値υを、各仮想の遺伝子クラスタについて計算式cにしたがって算出した。具体的には、仮想の遺伝子クラスタの乖離度算出部に格納されたυ値算出プログラムを適用し、各仮想の遺伝子クラスタについて計算式cにしたがって判定値υを算出した。次元数d’および係数aは、実施例1と同様に、それぞれ2および1を採用した。結果は図22に示される。なお、図22においても、実施例1の図17と同様に、ncl=1のときの値がncl=2のときの値より大きい仮想の遺伝子クラスタは除外している。結果を図22に示す。実施例1と同様、ncl=3のとき極大かつ最大値をとる仮想の遺伝子クラスタが一つあり、これがコウジ酸産生に必須の3つの遺伝子のみを含む遺伝子クラスタである。その他にncl=2に小さなピークを持つものが1つ見受けられるが、これはコウジ酸産生関連遺伝子の3つのうちの2つAO090113000137、AO090113000138からなるものである。図22を図17と比較すると分かるように、アノテーション付与により選定された機能を有する遺伝子の発現量変動比に重み付けを行うことで、標的とする遺伝子クラスタのスコアが嵩上げされて浮き彫りになり、探索対象の遺伝子クラスタをより高精度に検出可能となっている。続いて計算式dに従って、χ値およびυ値の積から遺伝子クラスタ評価値を算出した図23。具体的には、こうして得られたχ値およびυ値に対し、遺伝子クラスタ絞り込み部に格納された候補絞り込みプログラムを適用して、計算を行った。図23をみると明らかなように、結果二つの仮想の遺伝子クラスタが、10000以上の突出して大きな値をとっており、それと比較してその他のものは相対的に非常に小さなピークしか示さない。このうちncl=3で極大および最大値をとる仮想の遺伝子クラスタが、実施例1と同様、コウジ酸産生に必須の3つの遺伝子、AO090113000136、AO090113000137、AO090113000138のみを含んだものである。もう一つ顕著なものとして、ncl=2で極大値を持つものがあるが、これはコウジ酸産生に必須の3つの遺伝子のうちの2つAO090113000137およびAO090113000138からなるものである。それ以外の仮想の遺伝子クラスタは、相対的にほぼゼロと見なせる。この結果から明らかなように、アノテーションにより選定された遺伝子の発現量変動比に重み付けを行うことで標的とする遺伝子クラスタに対応する仮想の遺伝子クラスタのスコアをより顕著に高めていることは、重み付けを行わないもの図18と比較すると、明らかである。以上より、該当アノテーションにより選択されたる遺伝子の発現量変動比に重み付けを行うことによって、より高精度に機能に即した形で該当遺伝子クラスタを検出、同定できることが示された。実施例3アスペルギルスオリゼにおける特定機能を有するゲノム遺伝子により、仮想の遺伝子クラスタを構築し、スコアリングした場合の、コウジ酸生合成遺伝子の探索本実施例は、アスペルギルスオリゼのゲノム遺伝子中、特定の機能を持った遺伝子によって仮想の遺伝子クラスタを構築し、仮想の遺伝子クラスタのスコアを解析することにより、コウジ酸産生に必須の遺伝子を探索しうることを検証するための実験である。本実施例では、仮想の遺伝子クラスタのサイズnclを5として、アスペルギルスオリゼのゲノム配列より14032個の仮想遺伝子クラスタを作成した。実施例1と同様、途中抜けているものやゲノム断片の末端に位置する仮想の遺伝子クラスタは、ncl個より少ない遺伝子よりなるものとして構成した。この実験においては、実施例2の装置を用いた。ただし、実施例2における実験系C1からC3の3種のアレイデータについては、足し合わせて一つの発現量変動比m値にまとめたものを使用した。また、重み付けの代わりに、仮想の遺伝子クラスタのサイズを遺伝子数で5と設定し、構築された仮想の遺伝子クラスタ中から、アノテーションにより選定された複数種の機能遺伝子が含まれていることを条件として、仮想の遺伝子クラスタを選出し、該選出された仮想の遺伝子クラスタを、スコアリングする対象の仮想の遺伝子クラスタとするように、システムを変更した。その他は実施例2と同様である。すなわち、ゲノム上近傍に位置する条件として、仮想の遺伝子クラスタのサイズnclを5と設定して、記憶装置に記憶されたアスペルギルスオリゼのゲノム上の位置情報に基づき、14032個の仮想遺伝子クラスタを作成した。この場合において実施例1と同様、途中抜けているものやゲノム断片の末端に位置する仮想の遺伝子クラスタは、ncl個より少ない遺伝子で構成した。これらの仮想の遺伝子クラスタのうち、実施例2と同様にして、特定の機能の遺伝子を含むものを、機能に即したモチーフとの配列相同性を検索することで選び出した。特定の機能とはすなわち、以下の3つである。膜輸送体:transporterまたはmajorfacilitator転写制御因子:transcription酸化還元酵素:oxidoreductaseまたはdehydrogenase続いて、総数で14032個ある仮想の遺伝子クラスタの中から、該当する機能の注釈を持つ遺伝子を含むものを選び出した。その数のベン図を図24に示す。上記3つの因子膜輸送体、転写制御因子、酸化還元酵素の全てをもっている仮想遺伝子クラスタは、14032個のうち176個であった。また、上記3つの内、酸化還元酵素を除いて2つの要素膜輸送体、転写制御因子をもつものとすると、該仮想遺伝子クラスタは636個であった。上記手順は、具体的には、実施例2と同様にして、記憶装置に記憶されたアノテーションデータの中から以下の3種の機能を有する遺伝子を機能遺伝子選択部の選択プログラムを適用して選定し、さらに、構築された総数で14032個ある仮想の遺伝子クラスタの中から、選定された機能遺伝子を含むものを選出することにより行った。続いて、選び出した各仮想遺伝子クラスタに対して、クラスタスコアリングを行った。なお、アレイのデータは、参考例1および実施例1~2で述べたものと同様、系C1~C3におけるものを二色法によって測定したものであり、産生条件と非産生条件の下で生育させた菌体からmRNAを取り出し、それぞれ色素でラベリングした後でアレイ上のオリゴDNAにハイブリダイズすることによりデータを得て、そこから各遺伝子の発現量変動比より、GSE15435のIDで登録されたものの一部を用いた参考文献1。すなわちこのデータを、遺伝子発現量入力部を通して記憶部に保存した。このアレイデータは実施例1~4と異なり、一色法で測定されている。そこで各ゲノム遺伝子の発現量変動比m値を得るため、以下のように二次代謝産物をより多く産生すると考えられる条件とそうでない条件を比較し、前者を分子に後者を分母にした値をm値として算出した。検討した系は全部で2つである。C1培養開始後96時間/同18時間C2培養中、育成温度28°C/同37°C以降これら二つの系をそれぞれ系C1、C2とする。遺伝子は、2つの系それぞれにおいて12955個含まれている。Aクラスタスコアリング系C1、C2のそれぞれについて、実施例1と同様に計算式aに従って、仮想の遺伝子クラスタのサイズncl=1~30でクラスタスコアリングを行い、各仮想の遺伝子クラスタのスコアM値を得た。図28右は、各仮想の遺伝子クラスタの各サイズ毎にスコアの分布状態を示したヒストグラムである。図28左のグラフはヒストグラムの一部拡大図である。これをみると分かるように、ゼロを中心とした山型の正規分布様集団から外れて高いスコアM値を持つ仮想の遺伝子クラスタがあると、全体を表したヒストグラムにおいて山の中心が左側にずれるが、一見して、系C2において、nclの増加とともに、山の中心が左にずれていくことが分かる。Bデータ判定実施例1と同様にして、計算式eにしたがって、系C1およびC2におけるスコア分布評価値εを算出した図29。ここで仮想の遺伝子クラスタ数nは各々のクラスタサイズについて12955、次元数dは6を採用した。なお実施例1および2同様、ncl=1のときの値に応じて該当しない仮想の遺伝子クラスタは除外してある。図29示されるように、系C1ではε値はほぼゼロであるのに対し、系C2ではncl=18において極大かつ最大値を示す。これは、アフラトキシン産生至適温度25°Cであるため、系C2における温度に関する生理状態変化条件の設定が適切であり、他方系C1の変化条件設定は適切ではなかったことを反映している。すなわち系C2に基づく発現量変動比データを用いたクラスタスコアリングによって、ε値を増加させる同定対象の遺伝子クラスタが存在すること、及び、そのクラスタサイズは20前後であることが推定することができた。なおアスペルギルスフラバスが最も強く産生する二次代謝産物の一つは前述のアフラトキシンであり、その生合成遺伝子は29の遺伝子AFLA139100AFLA139440からなる遺伝子クラスタを形成していることが知られている参考文献2。ただしこの全てが同時に発現しているわけではなく、環境等によってその発現強度は変化する。本結果においてncl=20程度の大きなクラスタサイズの位置にピークが存在することは、アフラトキシンの生合成遺伝子クラスタの発現と対応していると考えられる。また本図ではε値が10の4乗のオーダーの値を示しているが、二次代謝産物の発現の弱い種であるアスペルギルスオリゼのε値は、図15にあるように10の3乗オーダーである。これは、アスペルギルスフラバスが同オリゼと比較して二次代謝産物を非常に強く発現する事実と一致する。以上より、系C2の発現量変動比データを用いた仮想の遺伝子クラスタ中のスコアから、構築された仮想の遺伝子クラスタ中に、標的とする遺伝子クラスタが含まれると予測できたため、系C2のDNAマイクロアレイデータセットを用いて、以下の実験を行った。C遺伝子クラスタ判定系C2に基づく仮想の各遺伝子クラスタのスコアM値から、実施例1と同様にして、各仮想の遺伝子クラスタについて、計算式bにしたがって遺伝子クラスタ判定値χを算出した。なお、この算出においては、実施例1および2同様、ncl=1のときの値に応じて該当しない仮想の遺伝子クラスタは除外してある。結果を図30に示すが、この折れ線グラフは、実施例1Cに記載したように、各仮想の遺伝子クラスタの構築において起点とした遺伝子が共通する、遺伝子サイズ1~30の各仮想の遺伝子クラスタの判定値を結んだものである。図30の結果から明らかなように、起点を同じくする仮想遺伝子クラスタの各折れ線グラフは、あるサイズでχ値の極大値をとっている。この起点を同じくする仮想の遺伝子クラスタの各折れ線グラフにおいて、χ値の極大値が高く、150程度を示すものは、大凡4種のサイズに分けられるが、このうち、極大値が高いピークを最も多く含むものは、サイズnclが20付近のものである。アスペルギルスフラバスのアフラトキシン合成に関与する遺伝子クラスタは既知であり、このサイズ20付近のピークの高いもの上位10個の各仮想の遺伝子クラスタ中の遺伝子について機能注釈を参照した結果、いずれもアフラトキシン合成に関与する遺伝子を含むことが明らかとなった。この結果は上記bの予測結果と一致し、このχ値の算出により、アスペルギルスフラバスにおけるアフラトキシン生合成に関与する遺伝子クラスタとその中に含まれるアフラトキシン生合成遺伝子を、ある程度特定可能であることを示す。一方、図30では他にも大きな値をとる仮想の遺伝子クラスタが複数存在するが、これは、その推定機能の注釈をみても、未知の二次代謝産物合成に関与する遺伝子クラスタである可能性が高い。次に、系C2の仮想の各遺伝子クラスタについて、実施例1と同様にして、遺伝子クラスタ判定値υを計算式cにしたがって算出した。ここで次元数d’は2、係数aは1を採用した。また実施例1および2同様、ncl=1のときの値に応じて該当しない仮想の遺伝子クラスタは除外してある。結果を図31に示すが、図31も図30と同様に、各仮想の遺伝子クラスタの構築において起点とした遺伝子が共通する、遺伝子サイズ1~30の各仮想の遺伝子クラスタの判定値を結んだものである。図31に示されるように、多くの仮想の遺伝子クラスタが極大値を示している。このうちυ値が200前後を示すものは、4つのサイズに分けられるが、このうち極大値が高いピークを最も多く含むものは、χ値の場合と同様にサイズが20付近ものであり、この各ピークの上位10個の仮想の遺伝子クラスタは、いずれも上述のアフラトキシン合成遺伝子を含んでいた。中には上述のアフラトキシン合成遺伝子クラスタを含むものが含まれる。すなわち本評価値υによっても、アフラトキシン生合成に関与する遺伝子クラスタ及びその中に含まれるアフラトキシン生合成遺伝子をある程度特定することができた。こうして得られたχ値およびυ値を元に、さらに遺伝子クラスタ候補を絞り込むため、実施例1と同様にして計算式dに従って二つの値の積から遺伝子クラスタ判定評価値を算出した。図32は、この算出結果に基づき、仮想の遺伝子クラスタサイズとχ×υ値の関係をグラフ化したものである。図32から明らかなように、多くの仮想遺伝子クラスタが特定のnclにおいて極大値を示すことが分かる。このうち、ncl=18で最大値を取る仮想の遺伝子クラスタは、AFLA139150AFLA139220、AFLA139240AFLA139280、AFLA139300AFLA139320からなるもので、そのすべてが既知であるアフラトキシン生合成遺伝子クラスタに含まれる遺伝子である。またその他に25000以上の値を示す仮想の遺伝子クラスタの機能注釈を見てみると、代表的な二次代謝産物関連遺伝子機能であるNRPSやP450といったものが挙げられており、これらも未知の二次代謝産物合成遺伝子クラスタである可能性が高い。さらに値の大きさを実施例1のアスペルギルスオリゼ図18と比較してみると、3倍近く同フラバスが高いことが分かる。これは、アスペルギルスフラバスが二次代謝産物を非常に活発に産生する種であるという事実と対応している。以上より、本発明がDNAマイクロアレイデータからゲノム上に集合して機能を果たす生合成遺伝子を同定する有効な手段であることが示された。参考文献1BeyondaflatoxinfourdistinctexpressionpatternsandfunctionalrolesassociatedwithAspergillusflavussecondarymetabolismgeneclustersDRYANGEORGIANNAら、MOLECULARPLANTPATHOLOGY11213226参考文献2GeneticregulationofaflatoxinbiosynthesisfromgenetogenomeDRYANGEORGIANNAら、FungalGeneticsandBiology46113125実施例6アスペルギルスニガーにおける遺伝子クラスタスコアリングによる生合成遺伝子推定本発明の同定手法に従って、アスペルギルスニガーの二次代謝産物を合成する遺伝子クラスタを推定した。この実験の使用装置は、実施例1の装置と同様である。DNAマイクロアレイのデータは、遺伝子発現解析データの公共データベースであるNCBBIのGEOより、GSE17329のIDで登録されたものの一部を用いた。すなわち、このデータを、遺伝子発現量データ入力部を通じて記憶部にゲノム遺伝子発現量データとして保存した。このアレイデータは実施例1~4で用いたアスペルギルスオリゼのものとは異なり、一色法で測定されている。そこで、ゲノム遺伝子発現量変動比m値を得るため、遺伝子発現量変動比算出部において、以下のように生理状態が変化する条件として以下の条件を設定し、前者を分子に後者を分母にした値をm値として算出した。検討した系は以下の2つである。なお、これらの系は、炭素源欠乏条件下で何らかの二次代謝関連遺伝子クラスタが関与していることを期待したもので、例えば上記したコウジ酸あるいはアフラトキシン産生等の特定の機能を標的としたものではない。C1培養中、炭素源枯渇後5555時間/同5時間C2培養中、炭素源枯渇後24時間/炭素源枯渇前35時間以下、上記2つの生理状態が変化する条件をそれぞれ系C1、C2とする。なお、発現量変動比は、2つの系それぞれにおいて14509個の遺伝子について算出した。Aクラスタスコアリング系C1~2のそれぞれについて、実施例1と同様に計算式aに従ってncl=1~30でクラスタスコアリングを行い、各仮想遺伝子クラスタのM値を得た。図33右は、各仮想の遺伝子クラスタの各サイズ毎にスコアの分布状態を示したヒストグラムである。図33左のグラフはヒストグラムの一部拡大図である。図33左の拡大図をみると分かるように、ゼロ付近を中心とした山型の正規分布様集団から外れて高いM値を持つ仮想遺伝子クラスタがあると、全体を表したヒストグラムにおいて山の中心が左側にずれるが、系C2において、ncl=5付近で山の中心が左にずれていることが分かる。Bデータ判定実施例1、2、5と同様にして、計算式eにしたがって、系C1~2におけるスコア分布評価値εを算出した図34。ここで仮想の遺伝子クラスタ数nは14509、次元数dは6を採用した。図34に示されるように、系C1はncl=8、系C2はncl=5においてそれぞれ極大値を示している。したがって二つの系においてともに、クラスタスコアリングによる平均化の方向に反して値を増加させる仮想の遺伝子クラスタが存在することを意味する。すなわちこの二つの系C1、C2における発現量変動比データを用いた、仮想の遺伝子クラスタのスコアリングによって、ε値を増大させる遺伝子クラスタが存在すること、及びその遺伝子クラスタのサイズが8前後あるいは5前後と推定された。ただし、この実験においては、上記したように、生理状態変化状件として炭素源欠乏条件を設定したものであり、該条件は特定の遺伝子クラスタを標的としたものではないので、極めて多数の遺伝子クラスタが関与していることが予想され、ε値によるサイズの予想は確定的なものではない。この点をふまえ、さらに以下の実験を行った。C遺伝子クラスタ判定系C1およびC2のそれぞれについて、系C1およびC2のDNAマイクロアレイデータから、実施例1と同様に計算式bにしたがって、各仮想遺伝子クラスタについて遺伝子クラスタ判定値χ算出した図35a;C1、同b;C2。なお実施例1、2、5と同様に、ncl=1のときの値に応じて該当しない仮想の遺伝子クラスタは除外してある。図35の結果から明らかなように、系C1、C2の双方において、多くの仮想遺伝子クラスタが極大値を示している。これより、アスペルギルスニガーには、系C2、C3の生理状態変化条件で変動する遺伝子クラスタが存在すると考えられ、これは既存の事実と一致する参考文献3。次に、系C1およびC2の各仮想遺伝子クラスタについて、実施例1と同様にして遺伝子クラスタ評価値υを計算式cにしたがって算出した図36a;C1、同b;C2。ここで次元数d’は2、係数aは1を採用した。ここでも実施例1、2、5と同様に、ncl=1のときの値に応じて該当しない仮想の遺伝子クラスタは除外してある。図36の結果に示されるように、系C1およびC2の双方において、複数の仮想の遺伝子クラスタが極大値を示す。ただしυ値の上位と下位の差はχ値図35に比べて増大しており、本実験においてはυ値の方が、少数の仮想の遺伝子クラスタを抽出するためにはより有利である。例えば、系C1においてυ値が100以上のものとした場合、該当する仮想遺伝子クラスタは1つのみである。系C2の場合、υ値が60以上の仮想遺伝子クラスタは3つのみである。こうして得られたχ値およびυ値を元に、実施例1と同様にして、計算式dに従って二つの値の積から遺伝子クラスタ判定評価値を算出した図37a;C1、同b;C2。図37の結果から明らかなように、系C1において、ncl=3で極大かつ最大値をとる仮想遺伝子クラスタが一つ存在することが分かる。系C2においてもいくつかの顕著なピークが見られ、例えば値4000以上とした場合、該当する仮想遺伝子クラスタは4個である。これらの仮想の遺伝子クラスタを構成する遺伝子について、その配列に基づきモチーフ検索による推定機能の注釈を見たところ、その多くが機能未知であり、該当する機能遺伝子は見いだせなかった。しかし本評価値において高い値を示すことを考えると、未知の該当遺伝子クラスタである可能性が高い。参考文献3ReviewofsecondarymetabolitesandmycotoxinsfromtheAspergillusnigergroupKFOGNIELSENら、AnalyticalandBioanalyticalChemistry39512251242実施例7アノテーション機能注釈に基づき選定された遺伝子を一以上含むことを条件として、仮想の遺伝子クラスタを構築した場合における、コウジ酸合成遺伝子の探索アスペルギルスオリゼのコウジ酸産生関連遺伝子からなる遺伝子クラスタを同定することを目的として、予測される機能に関連した注釈のついた遺伝子を選定したのち、該遺伝子を1以上含むように仮想の遺伝子クラスタを構築し、構築された各仮想の遺伝子クラスタをスコアリングして、該当遺伝子の同定を行った。この実験に使用した手法は、実施例1と基本的には同様であるが、実施例1においては、仮想の遺伝子クラスタの構築をする際、仮想の遺伝子クラスタのサイズを1~30と設定し、ゲノム遺伝子の配列順に全ての遺伝子が含まれるように仮想の遺伝子クラスタを構築したが、本実施例においては、ゲノム位置情報配列情報中に、アノテーション付与に基づき選定された機能遺伝子が出現したとき、その機能遺伝子を起点とする仮想の遺伝子クラスタを構築するように変更した点、構築された仮想の遺伝子クラスタのスコアリングにおいて、選定された機能遺伝子以外の遺伝子についての発現量変動比m値は無視し、選定された遺伝子の発現量変動比のみを用いるように変更した点で異なる。なお、遺伝子サイズについては、実施例1と同様にゲノムの遺伝子配列順に1~30個と設定した。具体的には、この実験に使用した装置は、実施例1に記載した装置と基本的に同様であるが、仮想の遺伝子クラスタ構築プログラムにおいて、ゲノム位置情報配列情報中に、アノテーション付与に基づく遺伝子選定部において選定された機能遺伝子が出現したとき、その機能遺伝子を起点とする仮想の遺伝子クラスタを構築するように変更した点、構築された仮想の遺伝子クラスタのスコアリングにおいて、選定された機能遺伝子以外の遺伝子についての発現量変動比m値は無視し、選定された遺伝子の発現量変動比のみを用いるように変更した点で異なる。なお、遺伝子サイズについては、実施例1と同様にゲノムの遺伝子配列順に1~30個と設定した。なお、本実施例では、実施例1のデータ判定結果において該当遺伝子クラスターが含まれていると予測された、系C27日目/4日目のアレイデータのみを用いて実験を行った。また、実施例2と同様、コウジ酸産生に必要な機能として以下の3つの機能を選出した。膜輸送体:transporterまたはmajorfacilitator転写制御因子:transcription酸化還元酵素:oxidoreductaseまたはdehydrogenaseこれらを選出した理由として、コウジ酸の生合成がグルコースから酸化により変換されていると推定されること、産生されたコウジ酸の膜輸送による培地中への分泌に膜輸送体が必要と推定されること、および関与する遺伝子の転写制御に転写因子が必要と推定されることが挙げられる。なお、上記の英単語は、アノテーションによる遺伝子選定に用いたキーワードである。アノテーション機能注釈による遺伝子の選定及び仮想の遺伝子クラスタの構築スコアリング一般に利用されているアノテーション推定プログラムの一つであるInterproscanhttpwwwebiacukToolsInterProScanを用いて、アスペルギルスオリゼのゲノムDNA上の各遺伝子についてアノテーションを付与し、上記3種の機能を有する遺伝子を選定した。具体的には、各遺伝子についてのアノテーションデータを本装置の入力装置に入力し、記憶装置に記憶した。記憶したアノテーションデータのデータを呼び出し、上記3種の機能を有する遺伝子を機能遺伝子選択部の選択プログラムを適用して選定した。なお、選定は各遺伝子について付与されたアノテーション中に上記3つの機能群に対応するキーワードが含まれるかどうかで行い、その結果、選定された遺伝子は、系C2において有効に遺伝子発現データを取得できた5595個の遺伝子のうち、796個であった。仮想の遺伝子クラスタの構築においては、上記変更プログラムを適用し、ゲノム遺伝子の位置情報に基づき、ゲノムの遺伝子配列中にこのように選定された機能遺伝子が出現した場合、この選定された遺伝子を起点遺伝子として、ゲノム遺伝子の配列順にクラスタサイズを1から30まで変化させて、仮想の遺伝子クラスタを構築した。すなわち、これにより、構築される仮想の遺伝子サイズ中には、必ず付与されたアノテーションに基づき選定された遺伝子が一つ以上含まれ、選定された機能遺伝子が含まれない仮想の遺伝子クラスタは構築されないが、構築された遺伝子クラスタにおいては、上記選定された機能遺伝子以外の遺伝子も含まれる。このようにした理由は、実施例1の装置に格納された仮想の遺伝子構築プログラムの変更を出来るだけ最小限にしたためである。しかし、構築された仮想の遺伝子クラスタのスコアリングにおいては、上記選定された機能遺伝子以外の遺伝子についての発現量変動比については無視し、選定された機能遺伝子の発現量変動比のみを用いて、計算式aによる計算を行った。これによれば、仮想の遺伝子クラスタのスコアは、仮想の遺伝子クラスタを、上記選定された機能遺伝子のみから構築した場合のスコアと全く同じになる。このようにして、得られた各仮想の遺伝子クラスタのスコアは、本発明装置の記憶部に記憶した。なお、本実施例では、構築された仮想の遺伝子クラスタに1つの遺伝子のみしか含まれない場合が含まれ、また、本実施例においては、実施例1~4と同様に、ゲノム上の末端側に位置する遺伝子については、組み合わせうる最大個数の遺伝子で仮想の遺伝子クラスタを構築したが、クラスタスコアリングの性質上、これらによる遺伝子クラスタの探索についての影響はない。このようにして構築された仮想の遺伝子クラスタは、各クラスタサイズについてそれぞれ796個である。つづいて、構築した各々の仮想の遺伝子クラスタについて、計算式aにしたがって、ncl=1~30でクラスタスコアリングを行って各仮想の遺伝子クラスタのスコアM値を得た。遺伝子クラスタ判定算出した各仮想の遺伝子クラスタのスコアM値に基づいて、計算式bにしたがって、各仮想の遺伝子クラスタについての判定値χを算出した。具体的には、記憶されている各仮想の遺伝子クラスタのスコアを呼び出し、仮想の遺伝子クラスタの乖離度算出部に格納された仮想の遺伝子乖離度判定プログラムのうちχ値算出プログラムを適用し、計算式bにしたがって、各仮想の遺伝子クラスタについての判定値χを算出した。図38は、各仮想の遺伝子クラスタについて、起点遺伝子を同じくする仮想の遺伝子クラスタの判定値χを、横軸をクラスタサイズとして結んで描いたものである。ここで、クラスタスコアリングによって絶対値を増加させない仮想の遺伝子クラスタは該遺伝子クラスタではないため、ncl=1のときの絶対値がncl=2のときの絶対値より大きい仮想の遺伝子クラスタは除外している。図をみると、多くの仮想の遺伝子クラスタの判定値χがゼロ付近に位置するのに対し、起点を同じくする仮想の遺伝子クラスタの3組が、大きな値を取っていることが分かる。その中で最大のものは、ncl=4のときに極大かつ最大値を取っている。これは、コウジ酸産生に必須の3つの遺伝子、AO090113000136、AO090113000137、AO090113000138に加えて、その隣に位置し、本実施例において遺伝子選定の対象であるアノテーションとしてmajorfacilitator膜輸送体を持つAO090113000139を含むものである。すなわち、本実施例では、選定対象のアノテーションがついた遺伝子の発現量変動比のみを用いて仮想の遺伝子クラスタのスコアリングを行っているため、スコアリングの対象となる要素が極度にそぎ落とされており、その結果、該当遺伝子クラスタの近傍にアノテーションによって選定された遺伝子がある場合、それを含んだ遺伝子クラスタが高い値をとりうる。しかし、この最大値を示す仮想の遺伝子クラスタが、コウジ酸の産生に必須の3つの遺伝子を含むことから、本手法は遺伝子クラスタの探索手法として有効である。実際、図38において最大値を示す起点遺伝子を同じくする仮想の遺伝子クラスタの組において、コウジ酸産生に必須の3つの遺伝子のみからなるncl=3のものと、隣接したAO090113000139を含むncl=4のものは、値がそれほど大きく変わらない。なお、その他にゼロから外れて大きな値を示す2組の仮想の遺伝子クラスタが存在するが、これらはコウジ酸産生に必須の3つの遺伝子のうち、AO090113000136を含まないものである。また、同様に、各仮想の遺伝子クラスタについて、計算式cにしたがって判定値υを算出した。具体的には、仮想の遺伝子クラスタの乖離度算出部に格納されたυ値算出プログラムを適用し、各仮想の遺伝子クラスタについて計算式cにしたがって判定値υを算出した。次元数d’および係数aは、実施例1と同様に、それぞれ2および1を採用した。起点遺伝子を同じくする仮想の遺伝子クラスタの判定値υをncl=1~30で結んだ結果を図39に示す。なお図39においても、ncl=1のときの値がncl=2のときの値より大きい仮想の遺伝子クラスタは除外している。判定値χのときと同様、ncl=4のとき極大かつ最大値をとる仮想の遺伝子クラスタが一つあり、これがコウジ酸産生に必須の3つの遺伝子に加えてもう一つの遺伝子AO090113000139を含む遺伝子クラスタである。図17と比較すると分かるように、アノテーションによる遺伝子の選定によって候補となる仮想の遺伝子クラスタの数が大きく減り、該当遺伝子クラスタが存在する場合、他のゼロ付近の値を持つものとの差がより明確になる。こうして得られたχ値およびυ値に対し、遺伝子クラスタ絞り込み部に格納された候補絞り込みプログラムを適用して、計算式dに従って二つの値の積から遺伝子クラスタ評価値を算出した図40。図をみると明らかなように、一つの仮想の遺伝子クラスタが、ncl=4のとき6000以上の大きな値をとっており、これがコウジ酸産生に必須の3つの遺伝子を含むものである。その他にゼロ付近の集団から離れて比較的大きな値を示す2つの仮想の遺伝子クラスタが存在するが、これはコウジ酸産生に必須の3つの遺伝子のうちAO090113000136を含まないものである。図40を図38および図39と比較すると明らかなように、二つの判定値を掛け合わせることで、該当遺伝子クラスタがより明確に大きな値をとり、該当遺伝子クラスタの予測精度が向上することが分かる。図41は、上記遺伝子クラスタ評価値を、横軸を遺伝子クラスタ番号として、各クラスタサイズについてプロットしたものである。クラスタサイズに対応した各図は、縦軸のスケールを合わせてある。この図でncl=4のときに最大値をとり、いずれのクラスタサイズにおいても突出して高い値を示しているのが、コウジ酸産生に必須の3つの遺伝子の3つあるいは2つを含む遺伝子クラスタである。このように本手法によって、本実施例の予測の対象であるコウジ酸産生遺伝子クラスタは鋭敏に検出されていることが分かる。以上の実験結果により、アノテーションにより選定されたる遺伝子を含む仮想の遺伝子クラスタ構築し、該選定された遺伝子の発現量変動比を用いてクラスタスコアリングを行うことで、高感度に、目的とする遺伝子クラスタ及びその中に含まれる遺伝子を探索できることが示された。また、この実験結果からいえば、アノテーションにより選定された1以上の遺伝子を組み合わせて仮想の遺伝子クラスタを構築し、スコアリングしても同様な結果が得られることは明らかである。本手法は強いフィルタリング操作を伴うものであり、該当するアノテーションを持つ遺伝子のm値を過度に反映する場合もある。しかし逆に、遺伝子間の発現変動比が比較的小さい場合などには、目的の遺伝子クラスタを鋭敏に予測できる手法である。実施例8フザリウムバーティシリオイデスにおける遺伝子クラスタスコアリングによる二次代謝産物生合成遺伝子の予測と検証本発明の同定手法に従って、菌類であるフザリウム属の一種、フザリウムバーティシリオイデスの二次代謝産物を合成する遺伝子クラスタを予測した。フザリウム属は、実施例1~6で用いた真菌類アスペルギルス属とは、進化系統樹的に遠い菌類である参考文献4。またフモニシンを始めとするマイコトキシンを産生することで知られており、その他多くの二次代謝産物生合成遺伝子クラスタを有すると考えられる参考文献5。DNAマイクロアレイのデータは、米国国立生物工学情報センターNCBIが提供する遺伝子発現解析データの公共データベースGEOより、GSE16900のIDで登録されたものの一部を用いた。このアレイデータは、フモニシン産生培地における培養時間が24487296時間である培養条件のそれぞれについて、一色法にて遺伝子の発現量を測定したものである。そこで、発現量変動比m値を得るため、以下のように二次代謝産物をより多く産生すると考えられる条件とそうでない条件を比較し、前者を分子に後者を分母にした値をm値として算出した。検討した系は2つである。C1培養時間72時間同24時間C2培養時間96時間同48時間以降これらの系をそれぞれC1C2とする。本発現情報には、遺伝子クラスタを構成するのに用いられる遺伝子が12230個含まれている。また、元のアレイデータでは各培養時間について3つのデータがとられているため、各遺伝子について3つのデータ間で発現量を平均化した後、以下の手順に進んだ。クラスタスコアリング系C1C2のそれぞれについて、計算式aに従ってncl=1~30でクラスタスコアリングを行い、各仮想遺伝子クラスタのM値を得た。図42はそのヒストグラムである。左の拡大図をみると分かるように、ゼロを中心とした山型の正規分布様集団から外れて高いM値を持つ仮想遺伝子クラスタがあると、各nclにおけるM値のヒストグラムにおいて、山の中心が左側にずれる。図右側の各nclにおけるヒストグラムを上から見ていくと、nclの増加とともに山の中心が左にずれていくことが分かる。なおクラスタスコアリングを行う際に必要となるゲノム上の遺伝子の連続情報は、アメリカ合衆国の研究機関であるBroadInstituteがウェブ上で公開しているデータベースFusariumComparativeSequencingProjectBroadInstituteofHarvardandMIT中の、fusariumverticillioides3genomesummarypergenetxtを用いた。データ判定計算式eにしたがって、系C1C2におけるスコア分布評価値eを算出した図43。ここで、仮想の遺伝子クラスタ数nは各nclにおいて12230であり、次元数dは6を採用した。図にあるように、系C1はncl=14、系C2はncl=5において極大値を示す。これは、二つの系においてともに、クラスタスコアリングによる平均化の方向に反して値を増加させる仮想の遺伝子クラスタが存在することを意味する。すなわちこの二つの系C1、C2において、クラスタスコアリングによって値を増加させる、本提案の同定対象に該当する遺伝子クラスタが存在すると判定できる。以上より、系C1、C2双方の遺伝子発現情報を用いて、以下の同定過程へ進む。遺伝子クラスタ判定系C1およびC2のDNAマイクロアレイデータから、各仮想遺伝子クラスタについて、計算式bにしたがって遺伝子クラスタ判定値cを算出した図44。ここで、遺伝子クラスタを検出するという目的に鑑み、ncl=1のときの絶対値がncl=2のときのものより大きな仮想の遺伝子クラスタについては除外してある。すると図にあるように、系C1、C2の双方において、ncl=1以外で複数の仮想遺伝子クラスタが極大値、極小値を示している。これより、フザリウムバーティシリオイデスには、複数の二次代謝関連遺伝子クラスタが存在すると考えられる。これは既存の事実と一致する参考文献5。次に、系C1およびC2の各仮想遺伝子クラスタについて、遺伝子クラスタ評価値uを計算式cにしたがって算出した図45。ここで次元数d’は2、係数aは1を採用した。ここでcと同様、ncl=1のときの値がncl=2のときのものより大きな仮想の遺伝子クラスタについては除外してある。すると図のように、系C1およびC2の双方において、複数の仮想の遺伝子クラスタが極大値を示す。ただし極大値におけるu値の上位と下位の差はc値図44に比べて増大しており、u値によって上位の少数の仮想の遺伝子クラスタを抽出することはc値よりさらに容易である。例えば系C1においてu値が100以上のものとした場合、該当する仮想遺伝子クラスタは1つのみである。系C2の場合、u値が150以上の仮想遺伝子クラスタは3つのみである。このように、評価値uを用いることで、該当遺伝子クラスタを高く評価することが出来ると考えられる。こうして得られたc値およびu値を元に、計算式dに従って二つの値の積から遺伝子クラスタ判定評価値を算出した。図46は、アレイデータに含まれるフザリウムバーティシリオイデスの遺伝子12230個のそれぞれを起点とし、ncl=1~30で構成した30個の仮想の遺伝子クラスタのうち、最大の評価値を取るものを、横軸をゲノム上の起点となる遺伝子としてプロットした図である。図において、C1とC2の縦軸の縮尺は合わせてある。系C1において、突出して高い値をとる3つの仮想の遺伝子クラスタがある。これらはそれぞれ、遺伝子FVEG00316、FVEG08708、FVEG12519を起点とし、クラスタサイズが14516のものである。これらの仮想の遺伝子クラスタを構成する遺伝子について、遺伝子配列相同性検索blastを行った結果を表3に示す。データベースは、NCBIが提供するもので、微生物を含む多くの生物種の遺伝子配列を格納するNRNonRedundanthttpwwwncbinlmnihgovstafftaoURLAPIblastdbhtmlを使用した。表では、相同性の高さを評価する値Evalueが10の100乗以下のもののうち、ベストヒットのものを抜粋してある。既知の二次代謝物質であるフモニシンの生合成遺伝子は、15個からなるクラスタであることが、フザリウムバーティシリオイデスの完全世代であるジベレリンモニリフォルミスにおいて報告されている参考文献678。また、フザリウムバーティシリオイデスにおいては、これまでFUM1FUM6FUM7FUM8FUM9の5つがフモニシン生合成遺伝子と同定されている参考文献:9。表をみると、Aとラベルした遺伝子クラスタの構成遺伝子14個が、フモニシン生合成遺伝子15個のうちの14個であることが分かるラベルにFumと併記した。すなわち、本発明手法によって、二次代謝物質フモニシンの生合成遺伝子クラスタをほぼ正確に予測できることが示された。なお今回の結果に含まれなかった残り一つのフモニシン生合成遺伝子は、遺伝子相同性検索の結果によると、FVEG00316の一つ前のFVEG00315である。仮想の遺伝子クラスタの本判定評価値をみると、クラスタサイズ15のときのFVEG00315を起点とするものの値は9242であり、一方クラスタサイズ14のときのFVEG00316を起点とするものの値は9763である。図をみるとわかるようにこれら2点は近接しており図46のC1、遺伝子クラスタAのピークにおける2点、僅差でFVEG00315がもれたことが分かる。従って、同一の遺伝子が含まれる予測されたクラスタにおいては、判定評価値の最大値が近接している場合、仮想クラスタサイズが最大となる物を選択することによって、最も正確な予測結果が得られると考えられる。系C2においても複数の顕著なピークがみられるが図46、なかでもFVEG03696を起点とするクラスタサイズ4の仮想の遺伝子クラスタが、値10000を超える最大の正のピークを示している。このピークは、培養開始72時間後を24時間後と比較した系C1では見られなかったもので、培養開始96時間後に初めて発現してくる遺伝子クラスタの存在を示唆している。また系C2では、FVEG08709を起点とするクラスタサイズ4の仮想の遺伝子クラスタが、大きな負の値をとっている。これは、系C1においては正の値を示していたものと、起点は一つずつずれているものの同等であり、培養開始72時間後には発現していたものが、96時間後には発現を止めた遺伝子クラスタであると推測される。このように、本手法を用いる際に、比較する系を目的に応じて選ぶことで、一つの生物種であっても異なる遺伝子クラスタを検出することが可能である。これらの該当遺伝子クラスタ候補の機能については、blast検索の結果をみても判然としないが表4、仮想の遺伝子クラスタCに含まれるFVEG12523は二次代謝物質生合成遺伝子の一つであるpolyketidesynthaseと高い配列相同性を示しており、これまでに知られていない新規な二次代謝物質生合成遺伝子が検出されたと期待される。以上より本提案方法が、アスペルギルス属からは進化系統樹的に遠い菌類であるフザリウムバーティシリオイデスにおいても、アスペルギルス属のものと同様、全遺伝子の発現情報からゲノム上に集合して機能を果たす生合成遺伝子を同定する実効的な手段であることが示された。参考文献4EvolutionoftheFot1transposonsinthegenusFusariumdiscontinuousdistributionandepigeneticinactivationMJDaboussiら、MolecularBiologyandEvolution19510520参考文献5BiochemistryandgeneticsofFusariumtoxinsAEDesjardinsら、FusariumPaulENelsonSymposiumAPSPress参考文献6LinkageamonggenesresponsibleforfumonisinbiosynthesisinGibberellafujikuroimatingpopulationADesjardinsら、AppliedandEnvironmentalMicrobiology6225712576参考文献7ApolyketidesynthasegenerequiredforbiosynthesisoffumonisinmycotoxinsinGibberellafujikuroimatingpopulationARHProctorら、FungalGeneticsandBiology27100112参考文献8Coexpressionof15contiguousgenesdelineatesafumonisinbiosyntheticgeneclusterinGibberellamoniliformisRHProctorら、FungalGeneticsandBiology38237249参考文献9CharacterizationoffourclusteredandcoregulatedgenesassociatedwithFumonisinbiosynthesisinFusariumverticillioidesJASeoら、FungalGeneticsandBiology34155165実施例9大腸菌における遺伝子クラスタスコアリングによるラクトースオペロンの検出と検証本発明の同定手法に従って、大腸菌のラクトースオペロンを検出した。大腸菌は原核生物であり、実施例1~8までで本発明手法の検証に用いた真核生物とは生物の分類上大きく異なる。大腸菌はオペロンの存在が実証された最初の生物である。オペロンとは、ゲノム上に集合して機能を果たす一つの制御単位であり、複数の遺伝子がゲノム上にまとまって存在し高発現して機能するという性質上、本発明の同定対象に該当する。ここで、本実施例で実証するラクトースオペロンについて説明する。ラクトースオペロンは、リプレッサータンパク質をコードするlacIに続いて、プロモーター配列lacP、オペレーター配列lacO、そしてラクトースを代謝する3つの遺伝子lacZlacYlacAlacZYAから構成される。lacIは常時発現しており、lacO領域と強く結合するため、通常はその下流lacZYAは翻訳されない。ところが、lacIに翻訳されるリプレッサータンパク質は、ラクトースが異性化したような誘導物質が存在すると、高次構造を変化させて、lacO領域から遊離する。これによってラクトース代謝系であるlacZYAが翻訳され、ラクトースの代謝が可能となる参考文献10。DNAマイクロアレイのデータは、米国国立生物工学情報センターNCBIが提供する遺伝子発現解析データの公共データベースGEOより、GSE7265のIDで登録されたものを用いた参考文献1112。このアレイデータは、大腸菌MG1655株およびその変異株を用いて、グルコースとラクトースの2つの栄養源を含む培地上で培養したときの遺伝子発現の変化を、分刻みで追ったものである。この2つの栄養源を含む培地上では、大腸菌はまずグルコースを代謝し、グルコースがなくなった後ラクトースを代謝する。すなわち、グルコースからラクトースへと栄養源を切り替える際に、最初に実証されたオペロンであるラクトースオペロンが発現する。ここではこのデータセットのうち、野生株のデータを用いる。野生株のデータには、培養開始後17段階におけるデータセットが含まれており、それぞれ、培養開始7808308618698788888989089199299399699991035104910701089分後に取られたものである。各データは、対数増殖初期780分後のデータを分母とした発現誘導比の形で記載されているため、そのまま本手法に適応可能である。ただし、各測定段階について3~4つのデータが取られているため、各遺伝子について3~4つのデータ間で値を平均化した後、以下の手順に進んだ。なおデータに含まれる遺伝子数は4102である。クラスタスコアリング17の各測定段階における系のそれぞれについて、計算式aに従ってncl=1~30でクラスタスコアリングを行い、各仮想遺伝子クラスタのM値を得た。なおクラスタスコアリングを行う際に必要となるゲノム上の遺伝子の連続情報は、公共の学術データベースであるNCBIに登録されている、大腸菌MG1655株のゲノム情報IDNC000913httpwwwncbinlmnihgovnuccoreNC000913を用いた。ここで大腸菌は環状ゲノムであるため、起点を上記ゲノム情報の中でb0001と名付けられた遺伝子とし、全ての遺伝子は連続しているものとして扱った。また、ラクトースオペロンを構成するlacIlacZlacYlacAの4つの遺伝子は、本ゲノム情報では向きが逆になっており、lacAlacYlacZlacIの順に並んでいる。これらの遺伝子IDはそれぞれ、b0342b0343b0344b0345である。図47は、各仮想の遺伝子クラスタのM値のヒストグラムの一部である。左の拡大図をみると分かるように、ゼロを中心とした山型の正規分布様集団から外れて高いM値を持つ仮想遺伝子クラスタがあると、各nclにおけるM値のヒストグラムにおいて、山の中心が左側にずれる。図右側の各nclにおけるヒストグラムを上から見ていくと、nclの増加とともに山の中心が左または右にずれていくことが分かる。これは、正規分布から外れて高い低い値を持つ仮想の遺伝子クラスタの存在を示している。データ判定計算式eにしたがって、17つの各系におけるスコア分布評価値eを算出した図48。ここで、仮想の遺伝子クラスタ数nは各nclにおいて4102であり、次元数dは6を採用した。図にあるように、e値は、培養開始後878888898分後、および104910701089分後の6つの系において、大きな極大値を示す。ここでこの結果を、大腸菌の成長速度と照らし合わせる。図49は、本アレイデータに関する文献参考文献11に記載されている、培養開始後の大腸菌の増殖を表す濁度の時系列変化である。前培養をどこからとるかでアレイデータとの時間のラベルがずれているが、図49での開始点が、アレイデータの780分にあたり、以降の各点がアレイデータの各17段階に順次相当する。すると、図をみると分かるように、スコア分布評価値eが大きな極大値を示す878888898分後789点目および104910701089分後151617点目のデータはすべて、図49において濁度の上昇が留まっている箇所、すなわち増殖の停滞期にあたる。このうち最初の停滞期は、グルコースを全て消費して栄養源をラクトースに切り替えようとしている段階である。この段階では増殖を一時停止するため、ゲノム上でまとまって存在する増殖に必須のリボソーム遺伝子群が強く抑制される一方、ラクトースを消費するためにラクトースオペロンを発現する。したがって、この段階でe値が大きな極大値を示していることは、リボソーム遺伝子群の抑制およびラクトースオペロンの発現という現象と一致する。二つ目の停滞期はラクトースも枯渇した段階であり、成長そのものが停滞するために、ここでも増殖に必須のリボソーム遺伝子が強く抑制される参考文献13。この段階でのe値の大きな極大値は、このリボソーム遺伝子の抑制を検出していると考えられる。以上より、e値によって、ゲノム上でまとまって発現または抑制して機能する遺伝子群の存在を感度よく判定できることが示された。本実施例では、すでに同定されているラクトースオペロンを本手法によって検出できることを示すことが目的であるため、引き続き17段階すべてのデータを用いて以下の手順に進む。遺伝子クラスタ判定大腸菌MG1655株の培養開始後17段階におけるDNAマイクロアレイデータから、各仮想遺伝子クラスタについて、計算式bにしたがって遺伝子クラスタ判定値cを算出した図50。図において、一つの仮想の遺伝子クラスタにつき一本の線がグレーで描かれており、このうち黒の太線で描かれた線が、ラクトースオペロンを構成する4つの遺伝子のゲノム情報上の最初の遺伝子lacAb0342を起点とする遺伝子クラスタである。この遺伝子クラスタは、869分の系から徐々に上昇を始め、908919分の系で、極大値を示す仮想の遺伝子クラスタの中で最大値を示している。極大値を示す点はクラスタサイズ3のときであり、lacZYAからなるものである。一方この遺伝子クラスタにlacIは含まれないが、これは、lacIがラクトースオペロンの発現に関係なく常時発現しているという事実と一致する。次に同様に、17の各系における各仮想遺伝子クラスタについて、遺伝子クラスタ評価値uを計算式cにしたがって算出した図51。ここで次元数d’は2、係数aは1を採用した。すると図のように、c値のときと同様、黒い太線で示したlacAb0342から始まる遺伝子クラスタは、869分の系から徐々に値の上昇を始め、908919分の系で、全ての仮想の遺伝子クラスタの中で最大の極大値をクラスタサイズ3で示す。これは、グルコースが枯渇してラクトース代謝系が動き出すときに、ラクトース代謝遺伝子群lacZYAが発現するという事実と一致する。またこの段階において、黒の太線とその他のグレーで示した線の差は、図50におけるよりも拡大しており、評価値uを用いることで、該当遺伝子クラスタをc値よりもさらに高く評価することが出来ることが示された。こうして得られたc値およびu値を元に、計算式dに従って二つの値の積から遺伝子クラスタ判定評価値c ́uを算出した図52。cおよびu値のみの場合図50、図51よりもより鋭敏に、908分の系で、lacZYAである黒の太線が検出されていることが分かる。図53は、評価値c ́uについて、各々の起点遺伝子について、ncl=1~30で構成した30個の仮想の遺伝子クラスタのうちの最大値を、横軸を起点遺伝子IDとして描いた図である。図において、17の系全ての縦軸の縮尺は合わせてある。黒矢印で示したラクトースオペロンが、908分の系で最大の値を示している。また、白抜きの矢印で示したリボソーム遺伝子群が、増殖停滞期である878888898分、そして104910701089分の系で、強く負の値を示している。これらより、本評価値が、ゲノム上で集合して機能する遺伝子群を、細胞の状態に応じて正確に検出できることが示された。以上より本提案方法が、真核生物だけでなく原核生物においても、ゲノム上に集合して機能を果たす遺伝子群を検出する実効的な手段であることが示された。参考文献10ThelactoserepressorsystemparadigmsforregulationallostericbehaviorandproteinfoldingCJWilsonら、CellularandMolecularLifeSciences64316参考文献11GeneexpressionprofilingofEscherichiacoligrowthtransitionsanexpandedstringentresponsemodelDongEunChangら、MolecularMicrobiology45289306参考文献12Guanosine3’5’bispyrophosphatecoordinatesglobalgeneexpressionduringglucoselactosediauxieinEscherichiacoliMatthewFTraxlerら、ProceedingsoftheNationalAcademyofSciences10323742379参考文献13ControlofproteinsynthesisinEscherichiacoliIITranslationanddegradationoflactoseoperonmessengerribonucleicacidafterenergysourceshiftdownKCWestoverら、JournalofBiologicalChemistry24962806287"