(57)【要約】【課題】コンテクストに依存して顕在化する遺伝子の発現関係を抽出すること。【解決手段】ネットワーク抽出部１１０が発現量行列から相関ネットワークを作成するとともに同期遺伝子集団を抽出し、コンテクスト生成部１２０が、特定の同期遺伝子集団に属する遺伝子の発現量に基づいてコンテクストを生成し、ネットワーク抽出部１１０がコンテクスト生成部１２０により生成されたコンテクストに対応する発現量行列から相関ネットワークを作成して同期遺伝子集団を抽出する。そして、コンテクスト生成部１２０によるコンテクストの生成とネットワーク抽出部１１０による相関ネットワークの作成および同期遺伝子集団の抽出を繰り返す。【選択図】        図２
2006260073.tif 000002
【特許請求の範囲】
【請求項１】
複数の遺伝子の発現量に関するデータである発現データに基づいて遺伝子間の関係を抽出する遺伝子関係抽出プログラムであって、
前記複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて遺伝子の発現環境であるコンテクストを生成するコンテクスト生成手順と、
前記コンテクスト生成手順により生成されたコンテクストにおける遺伝子間の関係を特定する関係特定手順と、
をコンピュータに実行させることを特徴とする遺伝子関係抽出プログラム。
【請求項２】
前記コンテクスト生成手順は、前記関係特定手順により特定された関係が所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成し、
前記関係特定手順は、前記コンテクスト生成手順により生成された新たなコンテクストにおける遺伝子間の関係を特定し、
前記コンテクスト生成手順と前記関係特定手順を所定の条件が満たされるまで繰り返すことを特徴とする請求項１に記載の遺伝子関係抽出プログラム。
【請求項３】
前記関係特定手順により複数のコンテクストにおいてそれぞれ特定された複数の遺伝子間の関係を比較する関係比較手順をさらにコンピュータに実行させることを特徴とする請求項１または２に記載の遺伝子関係抽出プログラム。
【請求項４】
前記関係特定手順により特定される遺伝子間の関係は二つの遺伝子の発現量間の相関関係であり、前記コンテクスト生成手順は、前記関係特定手順により特定される相関係数の値が所定の閾値より大きい遺伝子間を接続して構成される相関ネットワークのうちの部分ネットワークに含まれる複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成することを特徴とする請求項２または３に記載の遺伝子関係抽出プログラム。
【請求項５】
前記コンテクスト生成手順は、前記関係特定手順により特定された関係が所定の関係を満たす複数の遺伝子の発現量の平均値について前記複数の発現データを対象としてヒストグラムを作成し、該作成したヒストグラムを用いて新たなコンテクストを生成することを特徴とする請求項２〜４のいずれか一つに記載の遺伝子関係抽出プログラム。
【請求項６】
前記コンテクスト生成手順は、前記関係特定手順により特定された関係が所定の関係を満たす複数の遺伝子の発現量の平均値を用いて前記複数の発現データに対してクラスタリングを行って新たなコンテクストを生成することを特徴とする請求項２〜４のいずれか一つに記載の遺伝子関係抽出プログラム。
【請求項７】
複数の遺伝子の発現量に関するデータである発現データに基づいて遺伝子間の関係を抽出する遺伝子関係抽出方法であって、
前記複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて遺伝子の発現環境であるコンテクストを生成するコンテクスト生成工程と、
前記コンテクスト生成工程により生成されたコンテクストにおける遺伝子間の関係を特定する関係特定工程と、
を含んだことを特徴とする遺伝子関係抽出方法。
【請求項８】
前記コンテクスト生成工程は、前記関係特定工程により特定された関係が所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成し、
前記関係特定工程は、前記コンテクスト生成工程により生成された新たなコンテクストにおける遺伝子間の関係を特定し、
前記コンテクスト生成工程と前記関係特定工程を所定の条件が満たされるまで繰り返すことを特徴とする請求項７に記載の遺伝子関係抽出方法。
【請求項９】
複数の遺伝子の発現量に関するデータである発現データに基づいて遺伝子間の関係を抽出する遺伝子関係抽出装置であって、
前記複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて遺伝子の発現環境であるコンテクストを生成するコンテクスト生成手段と、
前記コンテクスト生成手段により生成されたコンテクストにおける遺伝子間の関係を特定する関係特定手段と、
を備えたことを特徴とする遺伝子関係抽出装置。
【請求項１０】
前記コンテクスト生成手段は、前記関係特定手段により特定された関係が所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成し、
前記関係特定手段は、前記コンテクスト生成手段により生成された新たなコンテクストにおける遺伝子間の関係を特定し、
前記コンテクスト生成手段によるコンテクストの生成と前記関係特定手段による遺伝子間の関係の特定とを所定の条件が満たされるまで繰り返すことを特徴とする請求項９に記載の遺伝子関係抽出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、複数の遺伝子の発現量に関するデータである発現データに基づいて遺伝子間の関係を抽出する遺伝子関係抽出プログラム、遺伝子関係抽出方法および遺伝子関係抽出装置に関し、特に、複数の発現データから遺伝子の発現環境としてのコンテクストを生成し、特定のコンテクストにおいて顕在化する関係を抽出することができる遺伝子関係抽出プログラム、遺伝子関係抽出方法および遺伝子関係抽出装置に関するものである。
【背景技術】
【０００２】
近年、遺伝子解析技術の進歩により、数千〜数万種類におよぶ遺伝子の発現状態を一度に知ることができるようになった。そこで、多数の遺伝子の発現状態を利用して遺伝子間の発現関係を抽出しようとする技術の開発が進められている（例えば、特許文献１および非特許文献１参照。）。
【０００３】
抽出する遺伝子間の発現関係としては、例えば、ある遺伝子Ａの発現量が多くなると他の遺伝子Ｂの発現量も多くなる、あるいは、ある遺伝子Ａの発現量が多くなると他の遺伝子Ｂの発現量が少なくなる、といった遺伝子間の促進・抑制などの発現関係があり、かかる発現関係を特定することによって、病気の原因の究明や病気の治療に役立てることができる。
【０００４】
【特許文献１】再公表ＷＯ２００２／０４８９１５号公報
【非特許文献１】Homin K. Lee, Amy K. Hsu, Jon Sajdak, Jie Qin and Paul Pavlidis, "Coexpression Analysis of Human Genes Across Many Microarray Data Sets," Genome Research 14:1085-1094, 2004.
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、遺伝子間の発現関係は、特定のコンテクスト（遺伝子の発現環境）において顕在化すると思われ、多様なコンテクストの下で得られた発現データをむやみに解析するだけでは、発現関係の抽出が困難であるという問題がある。ここで、コンテクストとしては、例えば、組織や細胞内の部位などの空間的コンテクストや、発達時期や細胞周期などの時間的コンテクストがあるが、遺伝子の発現環境としてのコンテクストは、多数の要因が相互に影響し合う複雑なものと考えられる。
【０００６】
この発明は、上述した従来技術による問題点を解消するためになされたものであり、複数の発現データから遺伝子の発現環境としてのコンテクストを生成し、特定のコンテクストにおいて顕在化する関係を抽出することができる遺伝子関係抽出プログラム、遺伝子関係抽出方法および遺伝子関係抽出装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上述した課題を解決し、目的を達成するため、請求項１の発明に係る遺伝子関係抽出プログラムは、複数の遺伝子の発現量に関するデータである発現データに基づいて遺伝子間の関係を抽出する遺伝子関係抽出プログラムであって、前記複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて遺伝子の発現環境であるコンテクストを生成するコンテクスト生成手順と、前記コンテクスト生成手順により生成されたコンテクストにおける遺伝子間の関係を特定する関係特定手順と、をコンピュータに実行させることを特徴とする。
【０００８】
この請求項１の発明によれば、複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて遺伝子の発現環境であるコンテクストを生成し、生成したコンテクストにおける遺伝子間の関係を特定するよう構成したので、コンテクストに依存する遺伝子間の関係を特定することができる。
【０００９】
また、請求項２の発明に係る遺伝子関係抽出プログラムは、請求項１の発明において、前記コンテクスト生成手順は、前記関係特定手順により特定された関係が所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成し、前記関係特定手順は、前記コンテクスト生成手順により生成された新たなコンテクストにおける遺伝子間の関係を特定し、前記コンテクスト生成手順と前記関係特定手順を所定の条件が満たされるまで繰り返すことを特徴とする。
【００１０】
この請求項２の発明によれば、特定した関係が所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成し、生成した新たなコンテクストにおける遺伝子間の関係を特定することを所定の条件が満たされるまで繰り返すよう構成したので、様々なコンテクストにおける遺伝子間の関係を特定することができる。
【００１１】
また、請求項３の発明に係る遺伝子関係抽出プログラムは、請求項１または２の発明において、前記関係特定手順により複数のコンテクストにおいてそれぞれ特定された複数の遺伝子間の関係を比較する関係比較手順をさらにコンピュータに実行させることを特徴とする。
【００１２】
この請求項３の発明によれば、複数のコンテクストにおいてそれぞれ特定された複数の遺伝子間の関係を比較するよう構成したので、コンテクストに依存する遺伝子間の関係の特定を容易にすることができる。
【００１３】
また、請求項４の発明に係る遺伝子関係抽出プログラムは、請求項２または３の発明において、前記関係特定手順により特定される遺伝子間の関係は二つの遺伝子の発現量間の相関関係であり、前記コンテクスト生成手順は、前記関係特定手順により特定される相関係数の値が所定の閾値より大きい遺伝子間を接続して構成される相関ネットワークのうちの部分ネットワークに含まれる複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成することを特徴とする。
【００１４】
この請求項４の発明によれば、特定する遺伝子間の関係は二つの遺伝子の発現量間の相関関係であり、特定した相関係数の値が所定の閾値より大きい遺伝子間を接続して構成される相関ネットワークのうちの部分ネットワークに含まれる複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成するよう構成したので、遺伝子の発現量間の相関関係に基づいてコンテクストを生成することができる。
【００１５】
また、請求項５の発明に係る遺伝子関係抽出プログラムは、請求項２〜４の発明において、前記コンテクスト生成手順は、前記関係特定手順により特定された関係が所定の関係を満たす複数の遺伝子の発現量の平均値について前記複数の発現データを対象としてヒストグラムを作成し、該作成したヒストグラムを用いて新たなコンテクストを生成することを特徴とする。
【００１６】
この請求項５の発明によれば、特定した関係が所定の関係を満たす複数の遺伝子の発現量の平均値について複数の発現データを対象としてヒストグラムを作成し、作成したヒストグラムを用いて新たなコンテクストを生成するよう構成したので、所定の関係を満たす複数の遺伝子の発現量に基づいてコンテクストを生成することができる。
【００１７】
また、請求項６の発明に係る遺伝子関係抽出プログラムは、請求項２〜４の発明において、前記コンテクスト生成手順は、前記関係特定手順により特定された関係が所定の関係を満たす複数の遺伝子の発現量の平均値を用いて前記複数の発現データに対してクラスタリングを行って新たなコンテクストを生成することを特徴とする。
【００１８】
この請求項６の発明によれば、特定した関係が所定の関係を満たす複数の遺伝子の発現量の平均値を用いて複数の発現データに対してクラスタリングを行って新たなコンテクストを生成するよう構成したので、所定の関係を満たす複数の遺伝子の発現量に基づいてコンテクストを生成することができる。
【００１９】
また、請求項７の発明に係る遺伝子関係抽出方法は、複数の遺伝子の発現量に関するデータである発現データに基づいて遺伝子間の関係を抽出する遺伝子関係抽出方法であって、前記複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて遺伝子の発現環境であるコンテクストを生成するコンテクスト生成工程と、前記コンテクスト生成工程により生成されたコンテクストにおける遺伝子間の関係を特定する関係特定工程と、を含んだことを特徴とする。
【００２０】
この請求項７の発明によれば、複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて遺伝子の発現環境であるコンテクストを生成し、生成したコンテクストにおける遺伝子間の関係を特定するよう構成したので、コンテクストに依存する遺伝子間の関係を特定することができる。
【００２１】
また、請求項８の発明に係る遺伝子関係抽出方法は、請求項７の発明において、前記コンテクスト生成工程は、前記関係特定工程により特定された関係が所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成し、前記関係特定工程は、前記コンテクスト生成工程により生成された新たなコンテクストにおける遺伝子間の関係を特定し、前記コンテクスト生成工程と前記関係特定工程を所定の条件が満たされるまで繰り返すことを特徴とする。
【００２２】
この請求項８の発明によれば、特定した関係が所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成し、生成した新たなコンテクストにおける遺伝子間の関係を特定することを所定の条件が満たされるまで繰り返すよう構成したので、様々なコンテクストにおける遺伝子間の関係を特定することができる。
【００２３】
また、請求項９の発明に係る遺伝子関係抽出装置は、複数の遺伝子の発現量に関するデータである発現データに基づいて遺伝子間の関係を抽出する遺伝子関係抽出装置であって、前記複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて遺伝子の発現環境であるコンテクストを生成するコンテクスト生成手段と、前記コンテクスト生成手段により生成されたコンテクストにおける遺伝子間の関係を特定する関係特定手段と、を備えたことを特徴とする。
【００２４】
この請求項９の発明によれば、複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて遺伝子の発現環境であるコンテクストを生成し、生成したコンテクストにおける遺伝子間の関係を特定するよう構成したので、コンテクストに依存する遺伝子間の関係を特定することができる。
【００２５】
また、請求項１０の発明に係る遺伝子関係抽出装置は、請求項９の発明において、前記コンテクスト生成手段は、前記関係特定手段により特定された関係が所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成し、前記関係特定手段は、前記コンテクスト生成手段により生成された新たなコンテクストにおける遺伝子間の関係を特定し、前記コンテクスト生成手段によるコンテクストの生成と前記関係特定手段による遺伝子間の関係の特定とを所定の条件が満たされるまで繰り返すことを特徴とする。
【００２６】
この請求項１０の発明によれば、特定した関係が所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成し、生成した新たなコンテクストにおける遺伝子間の関係を特定することを所定の条件が満たされるまで繰り返すよう構成したので、様々なコンテクストにおける遺伝子間の関係を特定することができる。
【発明の効果】
【００２７】
請求項１、７および９の発明によれば、コンテクストに依存する遺伝子間の関係を特定するので、特定のコンテクストにおいて顕在化する関係を抽出することができるという効果を奏する。
【００２８】
また、請求項２、８および１０の発明によれば、様々なコンテクストにおける遺伝子間の関係を特定するので、様々なコンテクストにおいて顕在化する関係を抽出することができるという効果を奏する。
【００２９】
また、請求項３の発明によれば、コンテクストに依存する遺伝子間の関係の特定を容易にするので、特定のコンテクストにおいて顕在化する関係の抽出を容易にすることができるという効果を奏する。
【００３０】
また、請求項４の発明によれば、遺伝子の発現量間の相関関係に基づいてコンテクストを生成するので、遺伝子の発現量間の相関関係に基づいて顕在化する関係を抽出することができるという効果を奏する。
【００３１】
また、請求項５および６の発明によれば、所定の関係を満たす複数の遺伝子の発現量に基づいてコンテクストを生成するので、適切なコンテクストを生成することができるという効果を奏する。
【発明を実施するための最良の形態】
【００３２】
以下に添付図面を参照して、この発明に係る遺伝子関係抽出プログラム、遺伝子関係抽出方法および遺伝子関係抽出装置の好適な実施例を詳細に説明する。
【実施例】
【００３３】
まず、本実施例に係る遺伝子関係抽出装置によるコンテクストを用いた発現関係抽出の概念について説明する。図１は、本実施例に係る遺伝子関係抽出装置によるコンテクストを用いた発現関係抽出の概念を説明するための説明図である。
【００３４】
同図（ａ）は、特定のコンテクストで顕在化する発現関係の例を示し、遺伝子Ｇｃが休止しているというコンテクストでは、遺伝子Ｇａと遺伝子Ｇｂとの間に正の相関関係、すなわちＧａの発現量が多い場合にはＧｂの発現量も多いという関係が出現し、遺伝子Ｇｃが活性化しているというコンテクストでは、遺伝子Ｇａと遺伝子Ｇｂとの間に負の相関関係、すなわちＧａの発現量が多い場合にはＧｂの発現量が少ないという関係が出現することを示している。
【００３５】
また、同図（ｂ）は、コンテクストの組合せによる発現関係の限定の例を示し、コンテクストＡでＧａ、ＧｂおよびＧｃの間に正の相関関係が観測され、コンテクストＢでＧａとＧｂの間に正の相関関係が観測された場合には、コンテクストＡにおける発現関係は、ＧａとＧｂを含むものである可能性が増大することを示している。
【００３６】
このように、本実施例に係る遺伝子関係抽出装置は、コンテクストを用いることによって、特定のコンテクストにおいて顕在化する遺伝子間の発現関係を抽出する。ここで、コンテクストを用いて発現関係を抽出するためには、コンテクストをどのように選定するかが重要となる。
【００３７】
そこで、本実施例に係る遺伝子関係抽出装置では、多くの遺伝子が同期して発現される場合、すなわち、部分遺伝子集団の同期的な発現に基づいてコンテクストを選定することとしている。なお、部分遺伝子集団の同期的な発現の代わりに同期的な抑制に基づいてコンテクストを選定することもできる。
【００３８】
次に、本実施例に係る遺伝子関係抽出装置の構成について説明する。図２は、本実施例に係る遺伝子関係抽出装置の構成を示す機能ブロック図である。同図に示すように、この遺伝子関係抽出装置１００は、ネットワーク抽出部１１０と、コンテクスト生成部１２０と、ネットワーク比較部１３０とを有する。
【００３９】
ネットワーク抽出部１１０は、複数のサンプル（発現データ）から構成される発現量行列を入力して遺伝子の発現量間の相関関係に基づく相関ネットワークを生成し、生成した相関ネットワークから部分ネットワークを抽出する処理部である。
【００４０】
図３は、ネットワーク抽出部１１０が入力する発現量行列を示す図である。同図に示すように、発現量行列は、遺伝子ｍ（１≦ｍ≦Ｍ）のサンプルｎ（１≦ｎ≦Ｎ）における発現量をｘmnとする行列であり、ここでは、Ｍ＝２２２８３、Ｎ＝１５００である場合を示している。
【００４１】
また、コンテクストは、複数のサンプルから構成され、発現量行列において、複数の列から構成される部分行列に対応する。なお、図３では、コンテクストが列の連続するサンプルから構成される例を示しているが、コンテクストを構成する列の列番号が連続している必要はない。
【００４２】
図４は、ネットワーク抽出部１１０により部分ネットワークが抽出される相関ネットワークの一例を示す図である。同図において、「ＧＭＦＧ」、「ＣＯＲＯ１Ａ」、「ＴＪＰ１」、「ＮＣＫＡＰ１」などは遺伝子の名前であり、二つの遺伝子間の線は、この例では、二つの遺伝子の発現量間の相関係数の絶対値が所定の値以上であることを示している。
【００４３】
また、相関ネットワークでは、二つの遺伝子の発現量間の相関係数の値が大きいほど、すなわち二つの遺伝子の相関が強いほど遺伝子が近くに配置される。
【００４４】
例えば、「ＴＪＰ１」と「ＮＣＫＡＰ１」の相関係数は、０．８７０９であり、「ＴＪＰ１」と「ＮＣＫＡＰ１」は近くに配置されている。なお、遺伝子ｉにおけるサンプルｊの発現量をｘijとすると、遺伝子ペア（α、β）の相関係数ｒαβは式（１）のようになる。
【数１】
2006260073.tif 000003
【００４５】
ネットワーク抽出部１１０は、全体の遺伝子から相関係数の値が所定値（例えば０．８）以上である遺伝子の集まりである同期遺伝子集団を抽出する。ここで、二つの遺伝子の発現量間の相関係数の値が所定値以上である場合に二つの遺伝子間を線で結ぶ相関ネットワークを作成すると、同期遺伝子集団は、相関ネットワークの部分ネットワークに対応し、ネットワーク抽出部１１０は、相関ネットワークから複数の部分ネットワークを抽出することとなる。
【００４６】
なお、ここでは、二つの遺伝子の発現量間の相関係数の値が所定値以上である遺伝子の集まり、すなわち正の相関関係が強い遺伝子の集まりを同期遺伝子集団としたが、負の相関関係が強い遺伝子の集まりを同期遺伝子集団とすることもできる。
【００４７】
また、ここでは、ネットワーク抽出部１１０は、遺伝子ペアごとに相関係数を計算し、その値が閾値以上のペアをリンクとした相関ネットワークを作ることで同期遺伝子集団を抽出しているが、相関係数についてクラスタリングを行うことによって同期遺伝子集団を抽出することもできる。
【００４８】
コンテクスト生成部１２０は、ネットワーク抽出部１１０により抽出された部分ネットワークに含まれる遺伝子、すなわち、同期遺伝子集団に属する遺伝子についての代表的発現量を算出し、算出した代表的発現量に基づいてコンテクストを生成する処理部である。ここで、コンテクストは、サンプル全体の集合を複数に分割して得られる部分集合である部分サンプルに対応する。
【００４９】
図５は、コンテクスト生成部１２０によるコンテクスト生成を説明するための説明図である。同図は、コンテクスト生成部１２０が各サンプルに対して同期遺伝子集団に属する遺伝子の発現量の平均値を代表的発現量として算出し、算出した平均値に関して全サンプルを対象として作成したヒストグラムを示す。コンテクスト生成部１２０は、このヒストグラムに基づいてサンプルを分割することによってコンテクストを生成する。
【００５０】
例えば、同図では、同期遺伝子集団ａに属する遺伝子の発現量の平均値に基づいて作成したヒストグラムを二つの山に分割することが可能であるため、コンテクスト生成部１２０は、それぞれの山に対応するコンテクストを生成する。
【００５１】
このように、このコンテクスト生成部１２０が同期遺伝子集団に属する遺伝子の発現量の平均値に基づいてサンプルを分割してコンテクストを生成することによって、特定のコンテクストにおける遺伝子の発現関係を抽出することができる。
【００５２】
そして、コンテクスト生成部１２０によって生成した各コンテクストに対応する部分サンプルに対してネットワーク抽出部１１０により同期遺伝子集団を抽出し、ネットワーク抽出部１１０によって抽出した同期遺伝子集団に属する遺伝子の発現量の平均値に基づいてコンテクストを生成する過程を繰り返すことによって、様々なコンテクストを生成するとともに、様々なコンテクストにおける遺伝子の発現関係を抽出することができる。
【００５３】
なお、同期遺伝子集団の抽出とコンテクストの生成の繰り返しは、生成されたコンテクストに属するサンプルの数が一定の値以下になった場合など、所定の条件をみたした場合に終了することができる。あるいは、ユーザの指示によって終了することもできる。
【００５４】
また、ここでは、コンテクスト生成部１２０は、同期遺伝子集団に属する遺伝子の代表的発現量として遺伝子の発現量の平均値を用いることとしたが、代表的発現量としては、特異値分解で得られる第一主成分などの他の値を用いることもできる。
【００５５】
また、ここでは、コンテクスト生成部１２０は、発現量の平均値についてヒストグラムを作成してサンプルを分割することとしたが、代表的発現量に対して、クラスタリングや二値化処理などを適用することで二つまたは三つ以上のコンテクストを得ることもできる。
【００５６】
ネットワーク比較部１３０は、ネットワーク抽出部１１０により抽出された様々な部分ネットワークを比較して表示する処理部である。すなわち、このネットワーク比較部１３０は、様々なコンテクストにおける相関ネットワークを比較することによってコンテクストの比較を行う。
【００５７】
図６は、ネットワーク比較部１３０によるコンテクスト比較の一例を示す図である。同図に示すように、ネットワーク比較部１３０は、例えば、コンテクストＡ、コンテクストＢおよびコンテクストＡ＋Ｂにおける相関ネットワークを比較して表示する。
【００５８】
次に、ネットワーク抽出部１１０の処理手順について説明する。図７は、ネットワーク抽出部１１０の処理手順を示す図である。同図に示すように、このネットワーク抽出部１１０は、全体の発現量行列から特定のコンテクストに対応する部分行列を抽出する（ステップＳ１０１）。なお、全サンプルから部分ネットワークを抽出する場合には、全体の発現量行列が抽出された部分行列となる。
【００５９】
そして、抽出した部分行列を用いて二つの遺伝子の発現量間の相関係数を計算するペア相関計算を行い（ステップＳ１０２）、遺伝子間相関行列を作成する。ここで、遺伝子間相関行列とは、遺伝子ペア（α、β）の相関係数ｒαβを要素とする行列である。
【００６０】
そして、作成した遺伝子間相関行列に基づいて全体の遺伝子集団から同期遺伝子集団を抽出する（ステップＳ１０３）。すなわち、遺伝子間相関行列に基づいて相関ネットワークを作成し、作成した相関ネットワークから部分ネットワークを抽出する。
【００６１】
このように、ネットワーク抽出部１１０が特定のコンテクストに対して相関ネットワークを作成することによって、特定のコンテクストにおいて顕在化する発現関係を抽出することができる。また、作成した相関ネットワークから部分ネットワークを抽出することによって、新たなコンテクストの生成に用いる同期遺伝子集団を抽出することができる。
【００６２】
次に、コンテクスト生成部１２０の処理手順について説明する。図８は、コンテクスト生成部１２０の処理手順を示す図である。同図に示すように、このコンテクスト生成部１２０は、全体の発現量行列、コンテクストを生成する元のコンテクストおよび特定の同期遺伝子集団を用いて新たなコンテクストを生成する（ステップＳ２０１）。
【００６３】
また、コンテクストを生成するとともに、生成した各コンテクストに対する評価値も算出する。ここで、評価値としては、コンテクストに含まれるサンプル数、他のコンテクストからの分離性、元のコンテクストからの変化量などを用いることができる。
【００６４】
そして、生成したコンテクストを評価値に基づいてランキングし（ステップＳ２０２）、評価値とともにユーザに提示する（ステップＳ２０３）。そして、ユーザにコンテクストを選択させ（ステップＳ２０４）、ユーザにより選択されたコンテクストを新たなコンテクストとしてネットワーク抽出部１１０に渡す。
【００６５】
このように、コンテクスト生成部１２０が全体の発現量行列、コンテクストを生成する元のコンテクストおよび特定の同期遺伝子集団を用いて新たなコンテクストを生成することによって、特定のコンテクストにおいて顕在化する発現関係を抽出することができる。
【００６６】
なお、ここでは、生成したコンテクストおよび評価値をユーザに提示して選択させることとしたが、所定の値以上の評価値を有するコンテクストを自動的に選択してネットワーク抽出部１１０に渡すこともできる。
【００６７】
次に、コンテクスト生成部１２０によるコンテクスト生成処理の処理手順について説明する。図９は、コンテクスト生成部１２０によるコンテクスト生成処理の処理手順を示す図である。なお、このコンテクスト生成処理は、図８のステップＳ２０１の処理に対応する。
【００６８】
図９に示すように、このコンテクスト生成処理では、全体の発現量行列から元のコンテクストに含まれるサンプルおよび同期遺伝子集団に含まれる遺伝子の発現量で構成される部分行列を抽出し（ステップＳ３０１）、抽出した部分行列の各サンプルに対して代表発現量、すなわち発現量の平均値を計算する（ステップＳ３０２）。
【００６９】
そして、各サンプルに対して計算した平均値に基づいてヒストグラムを作成し、作成したヒストグラムからサンプルを分割してコンテクストを生成する（ステップＳ３０３）。また、生成した各コンテクストに対して評価値を計算し（ステップＳ３０４）、コンテクストとともに蓄積する（ステップＳ３０５）。
【００７０】
このように、コンテクスト生成部１２０は、同期遺伝子集団に含まれる遺伝子の発現量の平均値を各サンプルに対して計算し、計算した平均値に基づいてサンプルを分割してコンテクストを生成する。
【００７１】
上述してきたように、本実施例では、ネットワーク抽出部１１０が発現量行列から相関ネットワークを作成するとともに同期遺伝子集団を抽出し、コンテクスト生成部１２０が、特定の同期遺伝子集団に属する遺伝子の発現量に基づいてコンテクストを生成し、ネットワーク抽出部１１０がコンテクスト生成部１２０により生成されたコンテクストに対応する発現量行列から相関ネットワークを作成することによって、特定のコンテクストにおいて顕在化する発現関係を抽出することができる。
【００７２】
また、本実施例では、コンテクスト生成部１２０によるコンテクストの生成とネットワーク抽出部１１０による同期遺伝子集団の抽出を繰り返すことによって、様々なコンテクストを生成し、様々なコンテクストにおいて顕在化する発現関係を抽出することができる。
【００７３】
なお、本実施例では、遺伝子関係抽出装置について説明したが、遺伝子関係抽出装置が有する構成をソフトウェアによって実現することで、同様の機能を有する遺伝子関係抽出プログラムを得ることができる。そこで、この遺伝子関係抽出プログラムを実行するコンピュータについて説明する。
【００７４】
図１０は、本実施例に係る遺伝子関係抽出プログラムを実行するコンピュータの構成を示す機能ブロック図である。同図に示すように、このコンピュータ２００は、ＲＡＭ２１０と、ＣＰＵ２２０と、ＨＤＤ２３０と、ＬＡＮインタフェース２４０と、入出力インタフェース２５０と、ＤＶＤドライブ２６０とを有する。
【００７５】
ＲＡＭ２１０は、プログラムやプログラムの実行途中結果などを記憶するメモリであり、ＣＰＵ２２０は、ＲＡＭ２１０からプログラムを読み出して実行する中央処理装置である。
【００７６】
ＨＤＤ２３０は、プログラムやデータを格納するディスク装置であり、ＬＡＮインタフェース２４０は、コンピュータ２００をＬＡＮ経由で他のコンピュータに接続するためのインタフェースである。
【００７７】
入出力インタフェース２５０は、マウスやキーボードなどの入力装置および表示装置を接続するためのインタフェースであり、ＤＶＤドライブ２６０は、ＤＶＤの読み書きを行う装置である。
【００７８】
そして、このコンピュータ２００において実行される遺伝子関係抽出プログラム２１１は、ＤＶＤに記憶され、ＤＶＤドライブ２６０によってＤＶＤから読み出されてコンピュータ２００にインストールされる。
【００７９】
あるいは、この遺伝子関係抽出プログラム２１１は、ＬＡＮインタフェース２４０を介して接続された他のコンピュータシステムのデータベースなどに記憶され、これらのデータベースから読み出されてコンピュータ２００にインストールされる。
【００８０】
そして、インストールされた遺伝子関係抽出プログラム２１１は、ＨＤＤ２３０に記憶され、ＲＡＭ２１０に読み出されてＣＰＵ２２０によって遺伝子関係抽出プロセス２２１として実行される。
【００８１】
（付記１）複数の遺伝子の発現量に関するデータである発現データに基づいて遺伝子間の関係を抽出する遺伝子関係抽出プログラムであって、
前記複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて遺伝子の発現環境であるコンテクストを生成するコンテクスト生成手順と、
前記コンテクスト生成手順により生成されたコンテクストにおける遺伝子間の関係を特定する関係特定手順と、
をコンピュータに実行させることを特徴とする遺伝子関係抽出プログラム。
【００８２】
（付記２）前記コンテクスト生成手順は、前記関係特定手順により特定された関係が所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成し、
前記関係特定手順は、前記コンテクスト生成手順により生成された新たなコンテクストにおける遺伝子間の関係を特定し、
前記コンテクスト生成手順と前記関係特定手順を所定の条件が満たされるまで繰り返すことを特徴とする付記１に記載の遺伝子関係抽出プログラム。
【００８３】
（付記３）前記関係特定手順により複数のコンテクストにおいてそれぞれ特定された複数の遺伝子間の関係を比較する関係比較手順をさらにコンピュータに実行させることを特徴とする付記１または２に記載の遺伝子関係抽出プログラム。
【００８４】
（付記４）前記関係特定手順により特定される遺伝子間の関係は二つの遺伝子の発現量間の相関関係であり、前記コンテクスト生成手順は、前記関係特定手順により特定される相関係数の値が所定の閾値より大きい遺伝子間を接続して構成される相関ネットワークのうちの部分ネットワークに含まれる複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成することを特徴とする付記２または３に記載の遺伝子関係抽出プログラム。
【００８５】
（付記５）前記コンテクスト生成手順は、前記関係特定手順により特定された関係が所定の関係を満たす複数の遺伝子の発現量の平均値について前記複数の発現データを対象としてヒストグラムを作成し、該作成したヒストグラムを用いて新たなコンテクストを生成することを特徴とする付記２〜４のいずれか一つに記載の遺伝子関係抽出プログラム。
【００８６】
（付記６）前記コンテクスト生成手順は、前記関係特定手順により特定された関係が所定の関係を満たす複数の遺伝子の発現量の平均値を用いて前記複数の発現データに対してクラスタリングを行って新たなコンテクストを生成することを特徴とする付記２〜４のいずれか一つに記載の遺伝子関係抽出プログラム。
【００８７】
（付記７）前記コンテクスト生成手順は、前記関係特定手順により特定された関係が所定の関係を満たす複数の遺伝子の発現量の平均値を用いて前記複数の発現データの二値化を行って新たなコンテクストを生成することを特徴とする付記２〜４のいずれか一つに記載の遺伝子関係抽出プログラム。
【００８８】
（付記８）前記コンテクスト生成手順は、前記関係特定手順により特定された関係が所定の関係を満たす複数の遺伝子の発現量の特異値分解で得られる第一主成分について複数の発現データを対象としてヒストグラムを作成し、該作成したヒストグラムを用いて新たなコンテクストを生成することを特徴とする付記２〜４のいずれか一つに記載の遺伝子関係抽出プログラム。
【００８９】
（付記９）前記コンテクスト生成手順は、前記複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて複数のコンテクスト候補とともに各コンテクスト候補の評価値を生成してユーザに提示し、提示した複数のコンテクスト候補および評価値に基づいてユーザに選択させることによってコンテクストを生成することを特徴とする付記１〜８のいずれか一つに記載の遺伝子関係抽出プログラム。
【００９０】
（付記１０）前記コンテクスト生成手順は、前記複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて複数のコンテクスト候補とともに各コンテクスト候補の評価値を生成し、該生成した評価値が所定の値以上であるコンテクスト候補を選択してコンテクストを生成することを特徴とする付記１〜８のいずれか一つに記載の遺伝子関係抽出プログラム。
【００９１】
（付記１１）前記コンテクスト生成手順が生成する評価値は、元のコンテクストからの変化量であることを特徴とする付記９または１０に記載の遺伝子関係抽出プログラム。
【００９２】
（付記１２）複数の遺伝子の発現量に関するデータである発現データに基づいて遺伝子間の関係を抽出する遺伝子関係抽出方法であって、
前記複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて遺伝子の発現環境であるコンテクストを生成するコンテクスト生成工程と、
前記コンテクスト生成工程により生成されたコンテクストにおける遺伝子間の関係を特定する関係特定工程と、
を含んだことを特徴とする遺伝子関係抽出方法。
【００９３】
（付記１３）前記コンテクスト生成工程は、前記関係特定工程により特定された関係が所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成し、
前記関係特定工程は、前記コンテクスト生成工程により生成された新たなコンテクストにおける遺伝子間の関係を特定し、
前記コンテクスト生成工程と前記関係特定工程を所定の条件が満たされるまで繰り返すことを特徴とする付記１２に記載の遺伝子関係抽出方法。
【００９４】
（付記１４）前記関係特定工程により特定される遺伝子間の関係は二つの遺伝子の発現量間の相関関係であり、前記コンテクスト生成工程は、前記関係特定工程により特定される相関係数の値が所定の閾値より大きい遺伝子間を接続して構成される相関ネットワークのうちの部分ネットワークに含まれる複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成することを特徴とする付記１３に記載の遺伝子関係抽出方法。
【００９５】
（付記１５）前記コンテクスト生成工程は、前記関係特定工程により特定された関係が所定の関係を満たす複数の遺伝子の発現量の平均値について前記複数の発現データを対象としてヒストグラムを作成し、該作成したヒストグラムを用いて新たなコンテクストを生成することを特徴とする付記１３または１４に記載の遺伝子関係抽出方法。
【００９６】
（付記１６）複数の遺伝子の発現量に関するデータである発現データに基づいて遺伝子間の関係を抽出する遺伝子関係抽出装置であって、
前記複数の遺伝子のうち所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて遺伝子の発現環境であるコンテクストを生成するコンテクスト生成手段と、
前記コンテクスト生成手段により生成されたコンテクストにおける遺伝子間の関係を特定する関係特定手段と、
を備えたことを特徴とする遺伝子関係抽出装置。
【００９７】
（付記１７）前記コンテクスト生成手段は、前記関係特定手段により特定された関係が所定の関係を満たす複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成し、
前記関係特定手段は、前記コンテクスト生成手段により生成された新たなコンテクストにおける遺伝子間の関係を特定し、
前記コンテクスト生成手段によるコンテクストの生成と前記関係特定手段による遺伝子間の関係の特定とを所定の条件が満たされるまで繰り返すことを特徴とする付記１６に記載の遺伝子関係抽出装置。
【００９８】
（付記１８）前記関係特定手段により特定される遺伝子間の関係は二つの遺伝子の発現量間の相関関係であり、前記コンテクスト生成手段は、前記関係特定手段により特定される相関係数の値が所定の閾値より大きい遺伝子間を接続して構成される相関ネットワークのうちの部分ネットワークに含まれる複数の遺伝子についての複数の発現データに基づいて新たなコンテクストを生成することを特徴とする付記１７に記載の遺伝子関係抽出装置。
【００９９】
（付記１９）前記コンテクスト生成手段は、前記関係特定手段により特定された関係が所定の関係を満たす複数の遺伝子の発現量の平均値について前記複数の発現データを対象としてヒストグラムを作成し、該作成したヒストグラムを用いて新たなコンテクストを生成することを特徴とする付記１７または１８に記載の遺伝子関係抽出装置。
【産業上の利用可能性】
【０１００】
以上のように、本発明に係る遺伝子関係抽出プログラム、遺伝子関係抽出方法および遺伝子関係抽出装置は、遺伝子解析に有用であり、特に、遺伝子間の発現関係の解析に適している。
【図面の簡単な説明】
【０１０１】
【図１】本実施例に係る遺伝子関係抽出装置によるコンテクストを用いた発現関係抽出の概念を説明するための説明図である。
【図２】本実施例に係る遺伝子関係抽出装置の構成を示す機能ブロック図である。
【図３】ネットワーク抽出部が入力する発現量行列を示す図である。
【図４】ネットワーク抽出部により部分ネットワークが抽出される相関ネットワークの一例を示す図である。
【図５】コンテクスト生成部によるコンテクスト生成を説明するための説明図である。
【図６】ネットワーク比較部によるコンテクスト比較の一例を示す図である。
【図７】ネットワーク抽出部の処理手順を示す図である。
【図８】コンテクスト生成部の処理手順を示す図である。
【図９】コンテクスト生成部によるコンテクスト生成処理の処理手順を示す図である。
【図１０】本実施例に係る遺伝子関係抽出プログラムを実行するコンピュータの構成を示す機能ブロック図である。
【符号の説明】
【０１０２】
１００    遺伝子関係抽出装置
１１０    ネットワーク抽出部
１２０    コンテクスト生成部
１３０    ネットワーク比較部
２００    コンピュータ
２１０    ＲＡＭ
２１１    遺伝子関係抽出プログラム
２２０    ＣＰＵ
２２１    遺伝子関係抽出プロセス
２３０    ＨＤＤ
２４０    ＬＡＮインタフェース
２５０    入出力インタフェース
２６０    ＤＶＤドライブ
【図１】
2006260073.tif 000004
【図２】
2006260073.tif 000005
【図３】
2006260073.tif 000006
【図４】
2006260073.tif 000007
【図５】
2006260073.tif 000008
【図６】
2006260073.tif 000009
【図７】
2006260073.tif 000010
【図８】
2006260073.tif 000011
【図９】
2006260073.tif 000012
【図１０】
2006260073.tif 000013