"ダイナミックパスウェイマップDPMの生成方法であって、複数のパスウェイ要素を格納するパスウェイ要素データベースへのアクセスを提供するステップであって、各パスウェイ要素が少なくとも1つのパスウェイへのその関与により特徴付けられるステップと、前記パスウェイ要素データベースに結合された操作エンジンへのアクセスを提供するステップと、前記操作エンジンを使用して第1のパスウェイ要素を少なくとも1つの先験的に既知の属性と関連付けるステップと、前記操作エンジンを使用して第2のパスウェイ要素を少なくとも1つの仮の属性と関連付けるステップと、前記操作エンジンを使用して、前記第1のパスウェイ要素及び前記第2のパスウェイ要素を相互相関させ、少なくとも1つのパスウェイについてのその影響レベルを、それぞれ前記既知の属性及び前記仮の属性を使用して割り当てることにより、確率的パスウェイモデルを形成するステップと、分析エンジンを介して前記確率的パスウェイモデルを用いることにより、患者試料の複数の要素についての複数の計測属性から、特定のパスウェイについての参照パスウェイ活性情報を有するDPMを導き出すステップと、を含む方法。前記パスウェイが調節パスウェイネットワーク内にある、請求項1に記載の方法。前記調節パスウェイネットワークが、加齢パスウェイネットワーク、アポトーシスパスウェイネットワーク、恒常性パスウェイネットワーク、代謝パスウェイネットワーク、複製パスウェイネットワーク、及び免疫応答パスウェイネットワークからなる群から選択される、請求項2に記載の方法。前記パスウェイが、シグナル伝達パスウェイネットワーク内にある、及び個別的パスウェイネットワークのネットワーク内にある、からなる群から選択される、請求項1に記載の方法。前記シグナル伝達パスウェイネットワークが、カルシウム/カルモジュリン依存性シグナル伝達パスウェイネットワーク、サイトカイン媒介性シグナル伝達パスウェイネットワーク、ケモカイン媒介性シグナル伝達パスウェイネットワーク、成長因子シグナル伝達パスウェイネットワーク、ホルモンシグナル伝達パスウェイネットワーク、MAPキナーゼシグナル伝達パスウェイネットワーク、ホスファターゼ媒介性シグナル伝達パスウェイネットワーク、Rasスーパーファミリー媒介性シグナル伝達パスウェイネットワーク、及び転写因子媒介性シグナル伝達パスウェイネットワークからなる群から選択される、請求項4に記載の方法。前記パスウェイ要素がタンパク質である、請求項1に記載の方法。前記タンパク質が、受容体、ホルモン結合タンパク質、キナーゼ、転写因子、メチラーゼ、ヒストンアセチラーゼ、及びヒストンデアセチラーゼからなる群から選択される、請求項6に記載の方法。前記パスウェイ要素が核酸である、請求項1に記載の方法。前記核酸が、タンパク質コード配列、ゲノム調節配列、調節RNA、及びトランス活性化配列からなる群から選択される、請求項8に記載の方法。前記参照パスウェイ活性情報が、正常組織、罹患組織、加齢組織、又は回復組織に関して特定のものである、請求項1に記載の方法。前記既知の属性が、化合物属性、クラス属性、遺伝子コピー数、転写レベル、翻訳レベル、及びタンパク質活性からなる群から選択される、請求項1に記載の方法。前記仮の属性が、化合物属性、クラス属性、遺伝子コピー数、転写レベル、翻訳レベル、及びタンパク質活性からなる群から選択される、請求項1に記載の方法。前記計測属性が、突然変異、示差的遺伝子配列オブジェクト、遺伝子コピー数、転写レベル、翻訳レベル、タンパク質活性、及びタンパク質相互作用からなる群から選択される、請求項1に記載の方法。ダイナミックパスウェイマップDPMの生成方法であって、複数のパスウェイ要素を含む確率的パスウェイモデルを格納するモデルデータベースへのアクセスを提供するステップであって、前記複数のパスウェイ要素の第1の数値が相互相関され、且つ既知の属性に基づき少なくとも1つのパスウェイについての影響レベルを割り当てられ、前記複数のパスウェイ要素の第2の数値が相互相関され、且つ仮の属性に基づき少なくとも1つのパスウェイについての影響レベルを割り当てられる、ステップと、分析エンジンを介して、患者試料の複数の要素についての複数の計測属性を用いることにより、前記確率的パスウェイモデルを操作して前記DPMを得るステップであって、前記DPMは特定のパスウェイについての参照パスウェイ活性情報を有する、ステップと、を含む方法。前記パスウェイが、調節パスウェイネットワーク、シグナル伝達パスウェイネットワーク、又は個別的パスウェイネットワークのネットワーク内にある、請求項14に記載の方法。前記パスウェイ要素が、受容体、ホルモン結合タンパク質、キナーゼ、転写因子、メチラーゼ、ヒストンアセチラーゼ、及びヒストンデアセチラーゼからなる群から選択されるタンパク質であり、又は核酸が、ゲノム調節配列、調節RNA、及びトランス活性化配列からなる群から選択される、請求項14に記載の方法。前記参照パスウェイ活性情報が、正常組織、罹患組織、加齢組織、又は回復組織に関して特定のものである、請求項14に記載の方法。前記既知の属性が、化合物属性、クラス属性、遺伝子コピー数、転写レベル、翻訳レベル、及びタンパク質活性からなる群から選択される、請求項14に記載の方法。前記仮の属性が、化合物属性、クラス属性、遺伝子コピー数、転写レベル、翻訳レベル、及びタンパク質活性からなる群から選択される、請求項14に記載の方法。前記計測属性が、突然変異、示差的遺伝子配列オブジェクト、遺伝子コピー数、転写レベル、翻訳レベル、タンパク質活性、及びタンパク質相互作用からなる群から選択される、請求項14に記載の方法。生物学的に関連性のある情報を分析する方法であって、ダイナミックパスウェイマップDPMを格納するモデルデータベースへのアクセスを提供するステップであって、前記DPMが、第1の細胞又は患者試料の複数の要素についての複数の計測属性によって確率的パスウェイモデルを操作することにより生成される、ステップと、第2の細胞又は患者試料の複数の要素についての複数の計測属性を得るステップと、分析エンジンを介して、前記DPM及び前記第2の細胞又は患者試料の前記複数の要素についての前記複数の計測属性を用いることにより、前記第2の細胞又は患者試料についての予測パスウェイ活性情報を決定するステップと、を含む方法。前記第1の細胞又は患者試料の前記複数の要素についての前記計測属性が、健常な細胞又は組織、細胞又は組織の特定の年齢、細胞又は組織の特定の疾患、罹患した細胞又は組織の特定の病期、特定の性別、特定の人種、特定の職業集団、及び特定の種に特徴的である、請求項21に記載の方法。前記第2の細胞又は患者試料の前記複数の要素についての前記計測属性が、突然変異、示差的遺伝子配列オブジェクト、遺伝子コピー数、転写レベル、翻訳レベル、タンパク質活性、及びタンパク質相互作用からなる群から選択される、請求項21に記載の方法。前記第1の試料及び前記第2の試料が同じ細胞又は患者から得られ、前記第2の細胞又は患者試料の前記複数の要素についての前記複数の計測属性を得る前に前記細胞又は患者に処置を提供するステップをさらに含む、請求項21に記載の方法。前記処置が、放射線照射、前記患者への医薬品の投与、及び前記細胞への候補分子の投与からなる群から選択される、請求項24に記載の方法。前記候補分子が、候補分子ライブラリのメンバーである、請求項25に記載の方法。前記予測パスウェイ活性情報が、ある要素を、少なくとも1つのパスウェイにおいて階層的に優位な要素として同定する、請求項21に記載の方法。前記予測パスウェイ活性情報が、ある要素を、疾患に関連する少なくとも1つのパスウェイにおいて疾患の決定因子である要素として同定する、請求項21に記載の方法。予測パスウェイ活性情報のグラフィック表現を生成するステップをさらに含む、請求項21に記載の方法。少なくとも一部において前記予測パスウェイ活性情報に基づく推奨処置を生成するステップをさらに含む、請求項21に記載の方法。前記予測パスウェイ活性情報を用いて、疾患についての診断、予後判定、又は処置オプションの選択肢、及び食事指導からなる群から選択される推奨を処方するステップをさらに含む、請求項21に記載の方法。前記予測パスウェイ活性情報を用いて、後成的要因、ストレス適応、生物体の状態、及び修復又は治癒の状態を特定するステップをさらに含む、請求項21に記載の方法。現在、疾患及び障害の特性決定、診断、予防、治療、及び転帰の判断に用い得る方法を提供することが必要とされている。一実施形態において、本発明はダイナミックパスウェイマップDPMの生成方法を提供し、この方法は、複数のパスウェイ要素を格納するパスウェイ要素データベースへのアクセスを提供するステップであって、各パスウェイ要素が少なくとも1つのパスウェイへのその関与により特徴付けられるステップと;パスウェイ要素データベースに結合された操作エンジンへのアクセスを提供するステップと;操作エンジンを使用して第1のパスウェイ要素を少なくとも1つの先験的に既知の属性と関連付けるステップと;操作エンジンを使用して第2のパスウェイ要素を少なくとも1つの仮の属性と関連付けるステップと;操作エンジンを使用して、第1のパスウェイ要素及び第2のパスウェイ要素を相互相関させ、少なくとも1つのパスウェイについてのその影響レベルを、それぞれ既知の属性及び仮の属性を使用して割り当てることにより、確率的パスウェイモデルを形成するステップと;分析エンジンを介して確率的パスウェイモデルを用いることにより、患者試料の複数の要素についての複数の計測属性から、特定のパスウェイについての参照パスウェイ活性情報を有するDPMを導き出すステップとを含む。好ましい一実施形態において、パスウェイ要素はタンパク質である。より好ましい実施形態において、タンパク質は、受容体、ホルモン結合タンパク質、キナーゼ、転写因子、メチラーゼ、ヒストンアセチラーゼ、及びヒストンデアセチラーゼからなる群から選択される。代替的な好ましい実施形態において、パスウェイ要素は核酸である。より好ましい実施形態において、核酸は、タンパク質コード配列、ゲノム調節配列、調節RNA、及びトランス活性化配列からなる群から選択される。別のより好ましい実施形態において、参照パスウェイ活性情報は、正常組織、罹患組織、加齢組織、又は回復組織に関して特定のものである。好ましい実施形態において、既知の属性は、化合物属性、クラス属性、遺伝子コピー数、転写レベル、翻訳レベル、及びタンパク質活性からなる群から選択される。別の好ましい実施形態において、仮の属性は、化合物属性、クラス属性、遺伝子コピー数、転写レベル、翻訳レベル、及びタンパク質活性からなる群から選択される。別の代替的実施形態において、計測属性は、突然変異、示差的遺伝子配列オブジェクト、遺伝子コピー数、転写レベル、翻訳レベル、タンパク質活性、及びタンパク質相互作用からなる群から選択される。好ましい実施形態において、パスウェイは調節パスウェイネットワーク内にある。より好ましい実施形態において、調節パスウェイネットワークは、加齢パスウェイネットワーク、アポトーシスパスウェイネットワーク、恒常性パスウェイネットワーク、代謝パスウェイネットワーク、複製パスウェイネットワーク、及び免疫応答パスウェイネットワークからなる群から選択される。さらにより好ましい実施形態において、パスウェイはシグナル伝達パスウェイネットワーク内にある。代替的なさらにより好ましい実施形態において、パスウェイは個別的パスウェイネットワークのネットワーク内にある。最も好ましい実施形態において、シグナル伝達パスウェイネットワークは、カルシウム/カルモジュリン依存性シグナル伝達パスウェイネットワーク、サイトカイン媒介性シグナル伝達パスウェイネットワーク、ケモカイン媒介性シグナル伝達パスウェイネットワーク、成長因子シグナル伝達パスウェイネットワーク、ホルモンシグナル伝達パスウェイネットワーク、MAPキナーゼシグナル伝達パスウェイネットワーク、ホスファターゼ媒介性シグナル伝達パスウェイネットワーク、Rasスーパーファミリー媒介性シグナル伝達パスウェイネットワーク、及び転写因子媒介性シグナル伝達パスウェイネットワークからなる群から選択される。本発明はまた、ダイナミックパスウェイマップDPMの生成方法も提供し、この方法は、複数のパスウェイ要素を含む確率的パスウェイモデルを格納するモデルデータベースへのアクセスを提供するステップであって;複数のパスウェイ要素の第1の数値が相互相関され、且つ既知の属性に基づき少なくとも1つのパスウェイについての影響レベルを割り当てられ;複数のパスウェイ要素の第2の数値が相互相関され、且つ仮の属性に基づき少なくとも1つのパスウェイについての影響レベルを割り当てられる、ステップと;分析エンジンを介して、患者試料の複数の要素についての複数の計測属性を用いることにより、確率的パスウェイモデルを操作してDPMを得るステップであって、DPMは特定のパスウェイについての参照パスウェイ活性情報を有する、ステップとを含む。好ましい一実施形態において、パスウェイは、調節パスウェイネットワーク、シグナル伝達パスウェイネットワーク、又は個別的パスウェイネットワークのなかあるネットワーク内にある。別の好ましい実施形態において、パスウェイ要素は、受容体、ホルモン結合タンパク質、キナーゼ、転写因子、メチラーゼ、ヒストンアセチラーゼ、及びヒストンデアセチラーゼからなる群から選択されるタンパク質であるか、又は核酸が、ゲノム調節配列、調節RNA、及びトランス活性化配列からなる群から選択される。なおさらに好ましい実施形態において、参照パスウェイ活性情報は、正常組織、罹患組織、加齢組織、又は回復組織に関して特定のものである。別の好ましい実施形態において、既知の属性は、化合物属性、クラス属性、遺伝子コピー数、転写レベル、翻訳レベル、及びタンパク質活性からなる群から選択される。別の好ましい実施形態において、仮の属性は、化合物属性、クラス属性、遺伝子コピー数、転写レベル、翻訳レベル、及びタンパク質活性からなる群から選択される。なおさらに好ましい実施形態において、計測属性は、突然変異、示差的遺伝子配列オブジェクト、遺伝子コピー数、転写レベル、翻訳レベル、タンパク質活性、及びタンパク質相互作用からなる群から選択される。本発明はさらに、生物学的に関連性のある情報を分析する方法を提供し、この方法は、ダイナミックパスウェイマップDPMを格納するモデルデータベースへのアクセスを提供するステップであって、DPMが、第1の細胞又は患者試料の複数の要素についての複数の計測属性によって確率的パスウェイモデルを操作することにより生成される、ステップと;第2の細胞又は患者試料の複数の要素についての複数の計測属性を得るステップと;分析エンジンを介して、DPM及び第2の細胞又は患者試料の複数の要素についての複数の計測属性を用いることにより、第2の細胞又は患者試料についての予測パスウェイ活性情報を決定するステップとを含む。好ましい一実施形態において、第1の細胞又は患者試料の複数の要素についての計測属性は、健常な細胞又は組織、細胞又は組織の特定の年齢、細胞又は組織の特定の疾患、罹患した細胞又は組織の特定の病期、特定の性別、特定の人種、特定の職業集団、及び特定の種に特徴的である。別の好ましい実施形態において、第2の細胞又は患者試料の複数の要素についての計測属性は、突然変異、示差的遺伝子配列オブジェクト、遺伝子コピー数、転写レベル、翻訳レベル、タンパク質活性、及びタンパク質相互作用からなる群から選択される。代替的な好ましい実施形態において、第1の試料及び第2の試料は同じ細胞又は患者から得られ、第2の細胞又は患者試料の複数の要素についての複数の計測属性を得る前に細胞又は患者に処置を提供するステップをさらに含む。より好ましい実施形態において、処置は、放射線照射、患者への医薬品の投与、及び細胞への候補分子の投与からなる群から選択される。別のより好ましい実施形態において、候補分子は、候補分子ライブラリのメンバーである。別の好ましい実施形態において、予測パスウェイ活性情報は、ある要素を、少なくとも1つのパスウェイにおいて階層的に優位な要素として同定する。より好ましい実施形態において、予測パスウェイ活性情報は、ある要素を、疾患に関連する少なくとも1つのパスウェイにおいて疾患の決定因子である要素として同定する。代替的実施形態において、この方法は、予測パスウェイ活性情報のグラフィック表現を生成するステップをさらに含む。代替的実施形態において、この方法は、少なくとも一部において予測パスウェイ活性情報に基づく推奨処置を生成するステップをさらに含む。代替的実施形態において、この方法は、予測パスウェイ活性情報を用いて、疾患についての診断、予後判定、又は処置オプションの選択肢、及び食事指導からなる群から選択される推奨を処方するステップをさらに含む。代替的実施形態において、この方法は、予測パスウェイ活性情報を用いて、後成的要因、ストレス適応、生物体の状態、及び修復又は治癒の状態を特定するステップをさらに含む。別の実施形態において、本発明は、必要がある個体の臨床転帰を予測する統合パスウェイ活性IPAの行列を作成するための変形方法であって、iキュレーションされたパスウェイの集合を提供するステップであって、パスウェイが複数のエンティティを含む、ステップと;ii各キュレーションされたパスウェイを個別的な確率的グラフィカルモデルPGMに変換するステップであって、PGMが各キュレーションされたパスウェイの因子グラフから導かれる、ステップと、iii個体からの生体試料を提供するステップであって、生体試料が、キュレーションされたパスウェイの一つに含まれる少なくとも1つの内因性エンティティを含む、ステップと;iv生体試料中の内因性エンティティのレベルを決定するステップと;v内因性エンティティのレベルを、別の個体からの既に決定済みの対照試料におけるエンティティのレベルと比較するステップと;vi対照エンティティレベルに対する内因性エンティティのレベルが活性化しているか、基準どおりであるか、又は不活性化しているかを決定するステップと;vii内因性エンティティに数値ステートを割り当てるステップであって、活性化を表すステートが1であり、基準活性を表すステートが0であり、及び不活性化を表すステートが-1である、ステップと;viii別の内因性エンティティについてステップii~viを繰り返すステップと;x各内因性エンティティの数値ステートを統合パスウェイ活性IPAの行列にコンパイルするステップとを含み、x統合パスウェイ活性の行列がAであり、式中、Aijは、生体試料jにおけるエンティティiの推論された活性を表す、方法を提供し;この方法により、個体の臨床転帰を予測する統合パスウェイ活性行列がもたらされる。一実施形態において、IPA行列を作成する方法は、必要がある個体に対して、臨床転帰を予測するステップ、診断を提供するステップ、治療を提供するステップ、治療を送達するステップ、治療を投与するステップ、治療を実施するステップ、治療を管理するステップ、又は治療を投薬するステップを含む。別の実施形態において、キュレーションされたパスウェイの集合は、ヒト生物学の分析によるものである。さらに別の代替的実施形態において、キュレーションされたパスウェイの集合は、非ヒト生物学の分析によるものである。別の実施形態において、対照エンティティレベルに対する内因性エンティティのレベルを決定するステップは、スチューデントのt検定を用いて実行される。代替的実施形態において、対照エンティティレベルに対する内因性エンティティのレベルを決定するステップは、ANOVAを用いて実行される。別の実施形態において、変形方法は、2つ以上の個体からの複数の統合パスウェイ活性行列が組み合わされるステップであって、組み合わせた複数の行列からクラスターが得られるステップと、得られたクラスターの個々の行列間の距離が決定されるステップとを含む。一実施形態では、決定された距離がK平均法クラスター分析を用いて分析される。別の代替的実施形態では、決定された距離がK2平均法クラスター分析を用いて分析される。さらに他の実施形態では、変形方法は、生体試料中の内因性エンティティのレベルを決定するステップを含み、このステップは、抗体によって内因性エンティティを検出するステップであって、それにより内因性エンティティのレベルを決定するステップを含む。代替的実施形態において、生体試料中の内因性エンティティのレベルを決定するステップは、核酸プローブによって内因性エンティティを検出するステップであって、それにより内因性エンティティのレベルを決定するステップを含む。別の代替的実施形態において、生体試料中の内因性エンティティのレベルを決定するステップは、有機試薬によって内因性エンティティを検出するステップであって、有機試薬が内因性エンティティに結合し、それにより検出可能なシグナルが得られ、それにより内因性エンティティのレベルを決定するステップを含む。さらに別の代替的実施形態において、生体試料中の内因性エンティティのレベルを決定するステップは、無機試薬によって内因性エンティティを検出するステップであって、無機試薬が内因性エンティティに結合し、それにより検出可能なシグナルが得られ、それにより内因性エンティティのレベルを決定するステップを含む。別の代替的実施形態において、生体試料中の内因性エンティティのレベルを決定するステップは、有機試薬によって内因性エンティティを検出するステップであって、有機試薬が内因性エンティティと反応し、それにより検出可能なシグナルが得られ、それにより内因性エンティティのレベルを決定するステップを含む。別の代替的実施形態において、生体試料中の内因性エンティティのレベルを決定するステップは、無機試薬によって内因性エンティティを検出するステップであって、無機試薬が内因性エンティティと反応し、それにより検出可能なシグナルが得られ、それにより内因性エンティティのレベルを決定するステップを含む。好ましい実施形態において、生体試料中の内因性エンティティのレベルを決定するステップは、内因性エンティティに最適な波長で内因性エンティティの吸光度を計測するステップであって、それにより内因性エンティティのレベルを決定するステップを含む。代替的な好ましい実施形態において、生体試料中の内因性エンティティのレベルを決定するステップは、内因性エンティティに最適な波長で内因性エンティティの蛍光を計測するステップであって、それにより内因性エンティティのレベルを決定するステップを含む。さらに別の代替的な好ましい実施形態において、生体試料中の内因性エンティティのレベルを決定するステップは、内因性エンティティを酵素と反応させるステップであって、酵素が内因性エンティティを選択的に消化することにより少なくとも1つの産物が産生されるステップと、その少なくとも1つの産物を検出するステップであって、それにより内因性エンティティのレベルを決定するステップとを含む。より好ましい実施形態において、内因性エンティティを酵素と反応させるステップにより、少なくとも2つの産物の産生が生じる。さらにより好ましい実施形態において、少なくとも2つの産物を生じる内因性エンティティを酵素と反応させるステップの次に、産物を別の酵素により処理するステップであって、酵素が、産物の少なくとも1つを選択的に消化することにより少なくとも第3の産物が産生されるステップと、その少なくとも第3の産物を検出するステップであって、それにより内因性エンティティのレベルを決定するステップとが続く。別の好ましい実施形態において、個体は、健常な個体、無症候の個体、及び症候を示す個体の群から選択される。より好ましい実施形態において、個体は、ある病態と診断される個体からなる群から選択され、その病態は疾患及び障害からなる群から選択される。好ましい実施形態において、病態は、後天性免疫不全症候群AIDS、アジソン病、成人呼吸窮迫症候群、アレルギー、強直性脊椎炎、アミロイドーシス、貧血症、喘息、アテローム性動脈硬化症、自己免疫性溶血性貧血、自己免疫性甲状腺炎、良性前立腺肥大症、気管支炎、チェディアック東症候群、胆嚢炎、クローン病、アトピー性皮膚炎、皮膚筋炎dermnatomyositis、真性糖尿病、気腫、胎児赤芽球症、結節性紅斑、萎縮性胃炎、糸球体腎炎、グッドパスチャー症候群、痛風、慢性肉芽腫性疾患、グレーブス病、橋本甲状腺炎、過好酸球増加症、過敏性腸症候群、多発性硬化症、重症筋無力症、心筋又は心膜炎症、変形性関節症、骨粗鬆症、膵炎、多嚢胞性卵巣症候群、多発筋炎、乾癬、ライター症候群、関節リウマチ、強皮症、重症複合型免疫不全症SCID、シェーグレン症候群、全身性アナフィラキシー、全身性エリテマトーデス、全身性硬化症、血小板減少性紫斑病、潰瘍性大腸炎、ブドウ膜炎、ウェルナー症候群、癌の、血液透析の、及び体外循環の合併症、ウイルス、細菌、真菌、寄生虫、原虫、及び蠕虫の感染症;及び腺癌、白血病、リンパ腫、黒色腫、骨髄腫、肉腫、テラトカルシノーマ、特に、副腎癌、膀胱癌、骨癌、骨髄癌、脳癌、乳癌、子宮頸癌、胆嚢癌、ガングリオンの癌、消化管癌、心臓癌、腎癌、肝癌、肺癌、筋肉癌、卵巣癌、膵癌、副甲状腺癌、陰茎癌、前立腺癌、唾液腺癌、皮膚癌、脾臓癌、精巣癌、胸腺癌、甲状腺癌、及び子宮癌、静坐不能akathesia、アルツハイマー病、健忘、筋萎縮性側索硬化症ALS、失調、双極性障害、緊張病、脳性麻痺、脳血管疾患クロイツフェルトヤコブ病、認知症、うつ病、ダウン症候群、遅発性ジスキネジー、ジストニー、癲癇、ハンチントン病、多発性硬化症、筋ジストロフィー、神経痛、神経線維腫症、ニューロパチー、パーキンソン病、ピック病、網膜色素変性症、統合失調症、季節性情動障害、老年痴呆、脳卒中、トゥレット症候群並びに腺癌、黒色腫、及びテラトカルシノーマを含む癌、特に脳の癌からなる群から選択される。代替的な好ましい実施形態において、病態は、癌、例えば、腺癌、白血病、リンパ腫、黒色腫、骨髄腫、肉腫、テラトカルシノーマ、特に、副腎癌、膀胱癌、骨癌、骨髄癌、脳癌、乳癌、子宮頸癌、胆嚢癌、ガングリオンの癌、消化管癌、心臓癌、腎癌、肝癌、肺癌、筋肉癌、卵巣癌、膵癌、副甲状腺癌、陰茎癌、前立腺癌、唾液腺癌、皮膚癌、脾臓癌、精巣癌、胸腺癌、甲状腺癌、及び子宮癌;免疫障害、例えば、後天性免疫不全症候群AIDS、アジソン病、成人呼吸窮迫症候群、アレルギー、強直性脊椎炎、アミロイドーシス、貧血症、喘息、アテローム性動脈硬化症、自己免疫性溶血性貧血、自己免疫性甲状腺炎、気管支炎、胆嚢炎、接触性皮膚炎、クローン病、アトピー性皮膚炎、皮膚筋炎、真性糖尿病、気腫、リンパ球毒素性偶発性リンパ球減少症、胎児赤芽球症、結節性紅斑、萎縮性胃炎、糸球体腎炎、グッドパスチャー症候群、痛風、グレーブス病、橋本甲状腺炎、過好酸球増加症、過敏性腸症候群、多発性硬化症、重症筋無力症、心筋又は心膜炎症、変形性関節症、骨粗鬆症、膵炎、多発筋炎、乾癬、ライター症候群、関節リウマチ、強皮症、シェーグレン症候群、全身性アナフィラキシー、全身性エリテマトーデス、全身性硬化症、血小板減少性紫斑病、潰瘍性大腸炎、ブドウ膜炎、ウェルナー症候群、癌の、血液透析の、及び体外循環の合併症、ウイルス、細菌、真菌、寄生虫、原虫、及び蠕虫の感染症、外傷、ブルトン型X連鎖性無γグロブリン血症X-linkedagammaglobinemiaofBruton、分類不能型免疫不全症CVI、ディジョージ症候群胸腺低形成、胸腺異形成、IgA単独欠損症、重症複合型免疫不全症SCID、血小板減少症と湿疹を伴う免疫不全症ヴィスコットオールドリッチ症候群、チェディアック東症候群、慢性肉芽腫性疾患、遺伝性血管神経性浮腫、及びクッシング病に関連する免疫不全症;並びに発育障害、例えば、尿細管性アシドーシス、貧血症、クッシング症候群、軟骨形成不全性小人症、デュシェンヌ型及びベッカー型筋ジストロフィー、癲癇、性腺形成不全、WAGR症候群ウィルムス腫瘍、無虹彩、泌尿生殖器奇形、及び精神遅滞、スミスマジェニス症候群、骨髄異形成症候群、遺伝性粘膜上皮異形成、遺伝性角皮症、遺伝性ニューロパチー、例えばシャルコーマリートゥース病及び神経線維腫症、甲状腺機能低下症、水頭症、発作性障害、例えばシデナム舞踏病Syndenham’schorea及び脳性麻痺、脊椎披裂、無脳症、頭蓋脊椎披裂、先天性緑内障、白内障、感音難聴、並びに対象の任意の組織、器官、又は系、例えば、脳、副腎、腎臓、骨格系又は生殖器系が関与する細胞成長及び分化、胚発生、及び形態形成に関連する任意の障害からなる群から選択される。別の好ましい実施形態において、病態は、内分泌障害、例えば、性腺機能低下症、シーハン症候群、尿崩症、カルマン病、ハンドシュラークリスチャン病、レットレルジーベ病、サルコイドーシス、エンプティセラ症候群、及び小人症を含む下垂体機能低下に関連する障害;先端巨大症、巨人症、及び抗利尿ホルモンADH不適合分泌症候群SIADHを含む下垂体機能亢進症;並びに、甲状腺腫、粘液水腫、細菌感染症に関連する急性甲状腺炎、ウイルス感染症に関連する亜急性甲状腺炎、自己免疫性甲状腺炎橋本病、及びクレチン病を含む甲状腺機能低下に関連する障害;甲状腺中毒症及びその様々な型、グレーブス病、前脛骨粘液水腫、中毒性多結節性甲状腺腫、甲状腺癌、及びプランマー病を含む甲状腺機能亢進に関連する障害;並びに、コン病慢性高カルシウム血症chronichypercalemiaを含む副甲状腺機能亢進に関連する障害;呼吸器疾患、例えば、アレルギー、喘息、急性及び慢性炎症性肺疾患、ARDS、気腫、肺うっ血及び肺浮腫、COPD、間質性肺疾患、及び肺癌;癌、例えば、腺癌、白血病、リンパ腫、黒色腫、骨髄腫、肉腫、テラトカルシノーマ、特に、副腎癌、膀胱癌、骨癌、骨髄癌、脳癌、乳癌、子宮頸癌、胆嚢癌、ガングリオンの癌、消化管癌、心臓癌、腎癌、肝癌、肺癌、筋肉癌、卵巣癌、膵癌、副甲状腺癌、陰茎癌、前立腺癌、唾液腺癌、皮膚癌、脾臓癌、精巣癌、胸腺癌、甲状腺癌、及び子宮癌;及び免疫不全症、例えば、後天性免疫不全症候群AIDS、アジソン病、成人呼吸窮迫症候群、アレルギー、強直性脊椎炎、アミロイドーシス、貧血症、喘息、アテローム性動脈硬化症、自己免疫性溶血性貧血、自己免疫性甲状腺炎、気管支炎、胆嚢炎、接触性皮膚炎、クローン病、アトピー性皮膚炎、皮膚筋炎、真性糖尿病、気腫、リンパ球毒素性偶発性リンパ球減少症、胎児赤芽球症、結節性紅斑、萎縮性胃炎、糸球体腎炎、グッドパスチャー症候群、痛風、グレーブス病、橋本甲状腺炎、過好酸球増加症、過敏性腸症候群、多発性硬化症、重症筋無力症、心筋又は心膜炎症、変形性関節症、骨粗鬆症、膵炎、多発筋炎、乾癬、ライター症候群、関節リウマチ、強皮症、シェーグレン症候群、全身性アナフィラキシー、全身性エリテマトーデス、全身性硬化症、血小板減少性紫斑病、潰瘍性大腸炎、ブドウ膜炎、ウェルナー症候群、癌の、血液透析の、及び体外循環の合併症、ウイルス、細菌、真菌、寄生虫、原虫、及び蠕虫の感染症、並びに外傷からなる群から選択される。本発明はまた、本明細書に開示されるとおりの変形方法も提供し、ここでは次に元の構成データ集合の代わりに行列Aを使用して、臨床転帰との関連性を特定することができる。より好ましい実施形態において、キュレーションされたパスウェイは、生化学的パスウェイ、遺伝的パスウェイ、代謝パスウェイ、遺伝子調節パスウェイ、遺伝子転写パスウェイ、遺伝子翻訳パスウェイからなる群から選択される。別のより好ましい実施形態において、エンティティは、核酸、ペプチド、タンパク質、ペプチド核酸、炭水化物、脂質、プロテオグリカン、因子、補因子、生化学的代謝産物、有機組成物、無機組成物、及び塩からなる群から選択される。さらに他の好ましい実施形態において、生体試料は、患者試料、対照試料、実験的に処理された動物試料、実験的に処理された組織培養物試料、実験的に処理された細胞培養物試料、及び実験的に処理されたインビトロ生化学組成物試料からなる群から選択される。より好ましい実施形態において、生体試料は患者試料である。本発明はまた、患者試料において変化する分子パスウェイを推論する出力を有する確率的グラフィカルモデルPGMフレームワークも提供し、このPGMは複数の因子グラフを含み、ここで因子グラフは統合された生物学的データ集合を表し、及び患者試料において変化する推論された分子パスウェイはデータから分かる分子パスウェイを含み、及び前記分子パスウェイは臨床的又は非臨床的条件をもたらし、推論された分子パスウェイは臨床レジメン又は治療によって調整されることが既知であり、及び出力は臨床レジメンを指示する。好ましい実施形態において、データは、実験データ、臨床データ、疫学データ、及び現象論的データから選択される。別の好ましい実施形態において、病態は、疾患及び障害からなる群から選択される。より好ましい実施形態において、病態は、後天性免疫不全症候群AIDS、アジソン病、成人呼吸窮迫症候群、アレルギー、強直性脊椎炎、アミロイドーシス、貧血症、喘息、アテローム性動脈硬化症、自己免疫性溶血性貧血、自己免疫性甲状腺炎、良性前立腺肥大症、気管支炎、チェディアック東症候群、胆嚢炎、クローン病、アトピー性皮膚炎、皮膚筋炎dermnatomyositis、真性糖尿病、気腫、胎児赤芽球症、結節性紅斑、萎縮性胃炎、糸球体腎炎、グッドパスチャー症候群、痛風、慢性肉芽腫性疾患、グレーブス病、橋本甲状腺炎、過好酸球増加症、過敏性腸症候群、多発性硬化症、重症筋無力症、心筋又は心膜炎症、変形性関節症、骨粗鬆症、膵炎、多嚢胞性卵巣症候群、多発筋炎、乾癬、ライター症候群、関節リウマチ、強皮症、重症複合型免疫不全症SCID、シェーグレン症候群、全身性アナフィラキシー、全身性エリテマトーデス、全身性硬化症、血小板減少性紫斑病、潰瘍性大腸炎、ブドウ膜炎、ウェルナー症候群、癌の、血液透析の、及び体外循環の合併症、ウイルス、細菌、真菌、寄生虫、原虫、及び蠕虫の感染症;及び腺癌、白血病、リンパ腫、黒色腫、骨髄腫、肉腫、テラトカルシノーマ、特に、副腎癌、膀胱癌、骨癌、骨髄癌、脳癌、乳癌、子宮頸癌、胆嚢癌、ガングリオンの癌、消化管癌、心臓癌、腎癌、肝癌、肺癌、筋肉癌、卵巣癌、膵癌、副甲状腺癌、陰茎癌、前立腺癌、唾液腺癌、皮膚癌、脾臓癌、精巣癌、胸腺癌、甲状腺癌、及び子宮癌、静坐不能akathesia、アルツハイマー病、健忘、筋萎縮性側索硬化症ALS、失調、双極性障害、緊張病、脳性麻痺、脳血管疾患クロイツフェルトヤコブ病、認知症、うつ病、ダウン症候群、遅発性ジスキネジー、ジストニー、癲癇、ハンチントン病、多発性硬化症、筋ジストロフィー、神経痛、神経線維腫症、ニューロパチー、パーキンソン病、ピック病、網膜色素変性症、統合失調症、季節性情動障害、老年痴呆、脳卒中、トゥレット症候群並びに腺癌、黒色腫、及びテラトカルシノーマを含む癌、特に脳の癌からなる群から選択される。代替的なより好ましい実施形態において、病態は、癌、例えば、腺癌、白血病、リンパ腫、黒色腫、骨髄腫、肉腫、テラトカルシノーマ、特に、副腎癌、膀胱癌、骨癌、骨髄癌、脳癌、乳癌、子宮頸癌、胆嚢癌、ガングリオンの癌、消化管癌、心臓癌、腎癌、肝癌、肺癌、筋肉癌、卵巣癌、膵癌、副甲状腺癌、陰茎癌、前立腺癌、唾液腺癌、皮膚癌、脾臓癌、精巣癌、胸腺癌、甲状腺癌、及び子宮癌;免疫障害、例えば、後天性免疫不全症候群AIDS、アジソン病、成人呼吸窮迫症候群、アレルギー、強直性脊椎炎、アミロイドーシス、貧血症、喘息、アテローム性動脈硬化症、自己免疫性溶血性貧血、自己免疫性甲状腺炎、気管支炎、胆嚢炎、接触性皮膚炎、クローン病、アトピー性皮膚炎、皮膚筋炎、真性糖尿病、気腫、リンパ球毒素性偶発性リンパ球減少症、胎児赤芽球症、結節性紅斑、萎縮性胃炎、糸球体腎炎、グッドパスチャー症候群、痛風、グレーブス病、橋本甲状腺炎、過好酸球増加症、過敏性腸症候群、多発性硬化症、重症筋無力症、心筋又は心膜炎症、変形性関節症、骨粗鬆症、膵炎、多発筋炎、乾癬、ライター症候群、関節リウマチ、強皮症、シェーグレン症候群、全身性アナフィラキシー、全身性エリテマトーデス、全身性硬化症、血小板減少性紫斑病、潰瘍性大腸炎、ブドウ膜炎、ウェルナー症候群、癌の、血液透析の、及び体外循環の合併症、ウイルス、細菌、真菌、寄生虫、原虫、及び蠕虫の感染症、外傷、ブルトン型X連鎖性無γグロブリン血症X-linkedagammaglobinemiaofBruton、分類不能型免疫不全症CVI、ディジョージ症候群胸腺低形成、胸腺異形成、IgA単独欠損症、重症複合型免疫不全症SCID、血小板減少症と湿疹を伴う免疫不全症ヴィスコットオールドリッチ症候群、チェディアック東症候群、慢性肉芽腫性疾患、遺伝性血管神経性浮腫、及びクッシング病に関連する免疫不全症;並びに発育障害、例えば、尿細管性アシドーシス、貧血症、クッシング症候群、軟骨形成不全性小人症、デュシェンヌ型及びベッカー型筋ジストロフィー、癲癇、性腺形成不全、WAGR症候群ウィルムス腫瘍、無虹彩、泌尿生殖器奇形、及び精神遅滞、スミスマジェニス症候群、骨髄異形成症候群、遺伝性粘膜上皮異形成、遺伝性角皮症、遺伝性ニューロパチー、例えばシャルコーマリートゥース病及び神経線維腫症、甲状腺機能低下症、水頭症、発作性障害、例えばシデナム舞踏病Syndenham’schorea及び脳性麻痺、脊椎披裂、無脳症、頭蓋脊椎披裂、先天性緑内障、白内障、感音難聴、並びに対象の任意の組織、器官、又は系、例えば、脳、副腎、腎臓、骨格系又は生殖器系が関与する細胞成長及び分化、胚発生、及び形態形成に関連する任意の障害からなる群から選択される。さらに他のより好ましい実施形態において、病態は、内分泌障害、例えば、性腺機能低下症、シーハン症候群、尿崩症、カルマン病、ハンドシュラークリスチャン病、レットレルジーベ病、サルコイドーシス、エンプティセラ症候群、及び小人症を含む下垂体機能低下に関連する障害;先端巨大症、巨人症、及び抗利尿ホルモンADH不適合分泌症候群SIADHを含む下垂体機能亢進症;並びに、甲状腺腫、粘液水腫、細菌感染症に関連する急性甲状腺炎、ウイルス感染症に関連する亜急性甲状腺炎、自己免疫性甲状腺炎橋本病、及びクレチン病を含む甲状腺機能低下に関連する障害;甲状腺中毒症及びその様々な型、グレーブス病、前脛骨粘液水腫、中毒性多結節性甲状腺腫、甲状腺癌、及びプランマー病を含む甲状腺機能亢進に関連する障害;並びに、コン病慢性高カルシウム血症chronichypercalemiaを含む副甲状腺機能亢進に関連する障害;呼吸器疾患、例えば、アレルギー、喘息、急性及び慢性炎症性肺疾患、ARDS、気腫、肺うっ血及び肺浮腫、COPD、間質性肺疾患、及び肺癌;癌、例えば、腺癌、白血病、リンパ腫、黒色腫、骨髄腫、肉腫、テラトカルシノーマ、特に、副腎癌、膀胱癌、骨癌、骨髄癌、脳癌、乳癌、子宮頸癌、胆嚢癌、ガングリオンの癌、消化管癌、心臓癌、腎癌、肝癌、肺癌、筋肉癌、卵巣癌、膵癌、副甲状腺癌、陰茎癌、前立腺癌、唾液腺癌、皮膚癌、脾臓癌、精巣癌、胸腺癌、甲状腺癌、及び子宮癌;及び免疫不全症、例えば、後天性免疫不全症候群AIDS、アジソン病、成人呼吸窮迫症候群、アレルギー、強直性脊椎炎、アミロイドーシス、貧血症、喘息、アテローム性動脈硬化症、自己免疫性溶血性貧血、自己免疫性甲状腺炎、気管支炎、胆嚢炎、接触性皮膚炎、クローン病、アトピー性皮膚炎、皮膚筋炎、真性糖尿病、気腫、リンパ球毒素性偶発性リンパ球減少症、胎児赤芽球症、結節性紅斑、萎縮性胃炎、糸球体腎炎、グッドパスチャー症候群、痛風、グレーブス病、橋本甲状腺炎、過好酸球増加症、過敏性腸症候群、多発性硬化症、重症筋無力症、心筋又は心膜炎症、変形性関節症、骨粗鬆症、膵炎、多発筋炎、乾癬、ライター症候群、関節リウマチ、強皮症、シェーグレン症候群、全身性アナフィラキシー、全身性エリテマトーデス、全身性硬化症、血小板減少性紫斑病、潰瘍性大腸炎、ブドウ膜炎、ウェルナー症候群、癌の、血液透析の、及び体外循環の合併症、ウイルス、細菌、真菌、寄生虫、原虫、及び蠕虫の感染症、並びに外傷からなる群から選択される。本書に開示する実施形態は事例的且つ例示的なものであり、本発明を限定することは意図されない。本発明の特許請求の範囲から逸脱することなく他の実施形態を利用することができ、及び構造の変更を行うことができる。本明細書で使用されるとき、及び添付の特許請求の範囲において、文脈上特に明確に指示されない限り、単数形a、an、及びtheは複数形の参照を含む。従って、例えば、anmiRNAへの参照は、複数のかかるmiRNAを含み、医薬担体apharmaceuticalcarrierへの参照は、1つ又は複数の医薬担体及びその均等物への参照である等となる。本明細書で使用されるとき、用語キュレーションされたは、科学的及び/又は臨床的原理に従い、当該技術分野において公知の方法、例えば、分子生物学的技法、生化学的技法、生理学的技法、解剖学的技法、ゲノミクス技法、トランスクリプトミクス技法、プロテオミクス技法、メタボロミクス技法、ADME、及びバイオインフォマティクス技法などを用いて、試験、分析、及び同定されている生体分子及び/又は非生体分子の集合における関係を意味する。この関係は、生化学的パスウェイ、遺伝的パスウェイ、代謝パスウェイ、遺伝子調節パスウェイ、遺伝子転写パスウェイ、遺伝子翻訳パスウェイ、miRNA調節パスウェイ、偽遺伝子調節パスウェイなどの生化学的なものであってよい。ハイスループットデータは、癌組織における分子変化の包括的な概観を提供している。新しい技法が、腫瘍試料及び癌細胞株のゲノムコピー数変動、遺伝子発現、DNAメチル化、及びエピジェネティクスの状態の同時ゲノムワイドアッセイを可能にする。癌ゲノムアトラスTheCancerGenomeAtlas:TCGA、スタンドアップトゥキャンサーStandUpToCancer:SU2Cなど、及びさらに多くの研究が、様々な腫瘍について近い将来に計画されている。現行のデータ集合の解析から、患者間での遺伝的変化は異なり得るが、多くの場合に共通のパスウェイが関与することが分かる。従って、癌の進行に関与する関連パスウェイを特定し、それらが異なる患者でどのように変化するかを検出することが、決定的に重要である。本発明者らは、遺伝子間のキュレーションされたパスウェイ相互作用を包含する患者特異的な遺伝的活性を推論するための新規方法を提示する。遺伝子は、因子グラフにより、遺伝子及びその産物の発現及び既知の活性を符号化する相互に接続し合う変数の集合としてモデル化され、種々のオミクスデータをエビデンスとして取り込むことが可能になる。本方法は、確率的推論を用いて、患者においてパスウェイの活性例えば、内部遺伝子の状態、相互作用、又は高レベル出力が変化する程度を予測する。SPIAと呼ばれる競合するパスウェイ活性推論手法と比較して、本発明者らの方法は、限定はされないが多形性膠芽腫GBM及び乳癌の双方のデータ集合において、少ない偽陽性で癌関連パスウェイの活性の変化を特定する。ゲノムモデルに関するデータ統合を用いたパスウェイ認識アルゴリズムPathwayRecognitionAlgorithmusingDataintegrationonGenomicModels:PARADIGMは、遺伝子が独立して考慮される場合に見過ごされるGBM患者のサブセットについて、一貫したパスウェイレベル活性を特定した。さらに、このアルゴリズムを用いてGBM患者をその有意なパスウェイ摂動に基づきグループ化することで、GBM患者が、有意に異なる生存転帰を有する臨床的に関連性のあるサブグループに分けられる。これらの知見は、患者群又は個人の共通して摂動を受ける1つ又は複数のパスウェイにおいて重要な点にある遺伝子を標的化することのできる治療法が選択され得ることを示唆している。本発明者らは、任意の数のゲノム及び機能的ゲノムデータ集合を統合して、患者試料で変化する分子パスウェイを推論することのできる因子グラフKschischang:2001、上記をベースとする確率的グラフィカルモデルPGMフレームワークについて記載する。本発明者らは、膠芽腫及び乳癌の双方のデータ集合のコピー数変動及び遺伝子発現データを使用してモデルを検証した。構造化したパスウェイモデルを用いて推論される活性は、成功裏に膠芽腫患者を臨床的に関連性のあるサブタイプに層別化する。この結果から、パスウェイ情報に基づく推論が、遺伝子レベルのデータを独立して使用するより有益な情報を提供することが示唆される。より優れた予後判定及び診断を提供することに加え、統合パスウェイの活性化は、疾患の進行を抑止するために用いられ得る有望な治療に関して重要な手掛かりをもたらす。本発明者らは、統合された患者データから遺伝的パスウェイの活性を推論するPARADIGMPAthwayRecognitionAlgorithmusingDataIntegrationonGenomicModelsゲノムモデルに関するデータ統合を用いたパスウェイ認識アルゴリズムと呼ばれる手法を開発した。図1は、本手法の概要を示す。単一の患者試料に関する複数のゲノムスケールの計測値を組み合わせて、単一のNationalCancerInstituteNCIパスウェイの遺伝子、産物、及び抽象的プロセス入出力の活性を推論する。PARADIGMは統合パスウェイ活性IPAの行列Aを生成し、ここでAijは、患者試料jにおけるエンティティiの推論された活性を表す。次に、行列Aを元の構成データ集合の代わりに使用して、臨床転帰との関連性を特定することができる。本発明者らは、初めに各NCIパスウェイを個別的な確率モデルに変換した。p53アポトーシスパスウェイのごく一部の小例を図2cに示す。NCIのパスウェイ図を、隠れ状態と観測状態との双方を含む因子グラフに変換した図2。因子グラフは、エンティティ間の既知の相互作用を表す構造により、遺伝子関連及び生物学的プロセス関連の状態情報に関する観測を統合する。因子グラフで生物学的パスウェイを表現するため、本発明者らは、特定のmRNA又は複合体などの、細胞中のエンティティの状態を変数を用いて記述し、及びそれらのエンティティ間の相互作用及び情報フローを、因子を用いて表現する。これらの変数は各エンティティの示差的状態を、直接的な分子エンティティ濃度ではなく、対照又は正常レベルとの比較で表現する。この表現により、多くの場合に遺伝子の示差的状態を直接計測するか、或いは直接的な計測値をマッチングした対照との相対的な計測値に変換するかのいずれかであるDNAマイクロアレイによって検出される遺伝子発現などの、多くのハイスループットデータ集合をモデル化することが可能になる。また、遺伝子間の多様な調節関係についても可能になる。例えば、p53のMDM2媒介性ユビキチン依存性分解を表す相互作用は、p53タンパク質レベルを阻害する活性化MDM2としてモデル化することができる。一実施形態では、本方法を用いて、様々な診断上及び治療上の適用、例えば、癌組織の検出、癌組織の病期分類、転移性組織の検出など;神経障害、例えば、限定はされないが、アルツハイマー病、筋萎縮性側索硬化症ALS、パーキンソン病、統合失調症、癲癇、及びそれらの合併症;ディジョージ症候群、自閉症などの発達障害、多発性硬化症などの自己免疫障害、糖尿病などの検出;感染症、例えば、限定はされないが、ウイルス感染症、細菌感染症、真菌感染症、リーシュマニア、住血吸虫症、マラリア、条虫類、象皮病、線虫、ネマチンnematineによる感染症などの治療に用いることのできる臨床情報を提供し得る。一実施形態では、本方法を用いて、変化した遺伝子発現、mRNAの発現過剰に対する非存在/存在を検出及び定量化し、又は治療介入中のmRNAレベルをモニタするための臨床情報を提供し得る。発現の変化に関連する病態、疾患又は障害としては、後天性免疫不全症候群AIDS、アジソン病、成人呼吸窮迫症候群、アレルギー、強直性脊椎炎、アミロイドーシス、貧血症、喘息、アテローム性動脈硬化症、自己免疫性溶血性貧血、自己免疫性甲状腺炎、良性前立腺肥大症、気管支炎、チェディアック東症候群、胆嚢炎、クローン病、アトピー性皮膚炎、皮膚筋炎dermnatomyositis、真性糖尿病、気腫、胎児赤芽球症、結節性紅斑、萎縮性胃炎、糸球体腎炎、グッドパスチャー症候群、痛風、慢性肉芽腫性疾患、グレーブス病、橋本甲状腺炎、過好酸球増加症、過敏性腸症候群、多発性硬化症、重症筋無力症、心筋又は心膜炎症、変形性関節症、骨粗鬆症、膵炎、多嚢胞性卵巣症候群、多発筋炎、乾癬、ライター症候群、関節リウマチ、強皮症、重症複合型免疫不全症SCID、シェーグレン症候群、全身性アナフィラキシー、全身性エリテマトーデス、全身性硬化症、血小板減少性紫斑病、潰瘍性大腸炎、ブドウ膜炎、ウェルナー症候群、癌の、血液透析の、及び体外循環の合併症、ウイルス、細菌、真菌、寄生虫、原虫、及び蠕虫の感染症;及び腺癌、白血病、リンパ腫、黒色腫、骨髄腫、肉腫、テラトカルシノーマ、特に、副腎癌、膀胱癌、骨癌、骨髄癌、脳癌、乳癌、子宮頸癌、胆嚢癌、ガングリオンの癌、消化管癌、心臓癌、腎癌、肝癌、肺癌、筋肉癌、卵巣癌、膵癌、副甲状腺癌、陰茎癌、前立腺癌、唾液腺癌、皮膚癌、脾臓癌、精巣癌、胸腺癌、甲状腺癌、及び子宮癌が挙げられる。この診断アッセイは、遺伝子発現の変化を検出するため、ハイブリダイゼーション又は増幅技法を用いて患者の生体試料における遺伝子発現を標準試料と比較し得る。この比較のための定性的又は定量的方法は、当該技術分野において公知である。一実施形態では、本方法を用いて、変化した遺伝子発現;mRNAの非存在、存在、又は過剰発現を検出及び定量化し;又は治療介入中のmRNAレベルをモニタする臨床情報を提供し得る。発現の変化に関連する障害としては、静坐不能akathesia、アルツハイマー病、健忘、筋萎縮性側索硬化症ALS、失調、双極性障害、緊張病、脳性麻痺、脳血管疾患クロイツフェルトヤコブ病、認知症、うつ病、ダウン症候群、遅発性ジスキネジー、ジストニー、癲癇、ハンチントン病、多発性硬化症、筋ジストロフィー、神経痛、神経線維腫症、ニューロパチー、パーキンソン病、ピック病、網膜色素変性症、統合失調症、季節性情動障害、老年痴呆、脳卒中、トゥレット症候群並びに腺癌、黒色腫、及びテラトカルシノーマを含む癌、特に脳の癌が挙げられる。一実施形態では、本方法を用いて、哺乳動物タンパク質の発現又は活性の変化に関連する病態についての臨床情報を提供し得る。かかる病態の例としては、限定はされないが、後天性免疫不全症候群AIDS、アジソン病、成人呼吸窮迫症候群、アレルギー、強直性脊椎炎、アミロイドーシス、貧血症、喘息、アテローム性動脈硬化症、自己免疫性溶血性貧血、自己免疫性甲状腺炎、良性前立腺肥大症、気管支炎、チェディアック東症候群、胆嚢炎、クローン病、アトピー性皮膚炎、皮膚筋炎、真性糖尿病、気腫、胎児赤芽球症、結節性紅斑、萎縮性胃炎、糸球体腎炎、グッドパスチャー症候群、痛風、慢性肉芽腫性疾患、グレーブス病、橋本甲状腺炎、過好酸球増加症、過敏性腸症候群、多発性硬化症、重症筋無力症、心筋又は心膜炎症、変形性関節症、骨粗鬆症、膵炎、多嚢胞性卵巣症候群、多発筋炎、乾癬、ライター症候群、関節リウマチ、強皮症、重症複合型免疫不全症SCID、シェーグレン症候群、全身性アナフィラキシー、全身性エリテマトーデス、全身性硬化症、血小板減少性紫斑病、潰瘍性大腸炎、ブドウ膜炎、ウェルナー症候群、癌の、血液透析の、及び体外循環の合併症、ウイルス、細菌、真菌、寄生虫、原虫、及び蠕虫の感染症;及び腺癌、白血病、リンパ腫、黒色腫、骨髄腫、肉腫、テラトカルシノーマ、特に、副腎癌、膀胱癌、骨癌、骨髄癌、脳癌、乳癌、子宮頸癌、胆嚢癌、ガングリオンの癌、消化管癌、心臓癌、腎癌、肝癌、肺癌、筋肉癌、卵巣癌、膵癌、副甲状腺癌、陰茎癌、前立腺癌、唾液腺癌、皮膚癌、脾臓癌、精巣癌、胸腺癌、甲状腺癌、及び子宮癌、静坐不能akathesia、アルツハイマー病、健忘、筋萎縮性側索硬化症、失調、双極性障害、緊張病、脳性麻痺、脳血管疾患クロイツフェルトヤコブ病、認知症、うつ病、ダウン症候群、遅発性ジスキネジー、ジストニー、癲癇、ハンチントン病、多発性硬化症、筋ジストロフィー、神経痛、神経線維腫症、ニューロパチー、パーキンソン病、ピック病、網膜色素変性症、統合失調症、季節性情動障害、老年痴呆、脳卒中、トゥレット症候群並びに腺癌、黒色腫、及びテラトカルシノーマを含む癌、特に脳の癌が挙げられる。一実施形態では、本明細書にerein開示される方法を用いて、核酸配列の発現又は活性の低下に関連する障害を検出、病期分類、診断、及び/又は治療し得る。かかる障害の例としては、限定はされないが、癌、例えば、腺癌、白血病、リンパ腫、黒色腫、骨髄腫、肉腫、テラトカルシノーマ、特に、副腎癌、膀胱癌、骨癌、骨髄癌、脳癌、乳癌、子宮頸癌、胆嚢癌、ガングリオンの癌、消化管癌、心臓癌、腎癌、肝癌、肺癌、筋肉癌、卵巣癌、膵癌、副甲状腺癌、陰茎癌、前立腺癌、唾液腺癌、皮膚癌、脾臓癌、精巣癌、胸腺癌、甲状腺癌、及び子宮癌;免疫障害、例えば、後天性免疫不全症候群AIDS、アジソン病、成人呼吸窮迫症候群、アレルギー、強直性脊椎炎、アミロイドーシス、貧血症、喘息、アテローム性動脈硬化症、自己免疫性溶血性貧血、自己免疫性甲状腺炎、気管支炎、胆嚢炎、接触性皮膚炎、クローン病、アトピー性皮膚炎、皮膚筋炎、真性糖尿病、気腫、リンパ球毒素性偶発性リンパ球減少症、胎児赤芽球症、結節性紅斑、萎縮性胃炎、糸球体腎炎、グッドパスチャー症候群、痛風、グレーブス病、橋本甲状腺炎、過好酸球増加症、過敏性腸症候群、多発性硬化症、重症筋無力症、心筋又は心膜炎症、変形性関節症、骨粗鬆症、膵炎、多発筋炎、乾癬、ライター症候群、関節リウマチ、強皮症、シェーグレン症候群、全身性アナフィラキシー、全身性エリテマトーデス、全身性硬化症、血小板減少性紫斑病、潰瘍性大腸炎、ブドウ膜炎、ウェルナー症候群、癌の、血液透析の、及び体外循環の合併症、ウイルス、細菌、真菌、寄生虫、原虫、及び蠕虫の感染症、外傷、ブルトン型X連鎖性無γグロブリン血症X-linkedagammaglobinemiaofBruton、分類不能型免疫不全症CVI、ディジョージ症候群胸腺低形成、胸腺異形成、IgA単独欠損症、重症複合型免疫不全症SCID、血小板減少症と湿疹を伴う免疫不全症ヴィスコットオールドリッチ症候群、チェディアック東症候群、慢性肉芽腫性疾患、遺伝性血管神経性浮腫、及びクッシング病に関連する免疫不全症;並びに発育障害、例えば、尿細管性アシドーシス、貧血症、クッシング症候群、軟骨形成不全性小人症、デュシェンヌ型及びベッカー型筋ジストロフィー、癲癇、性腺形成不全、WAGR症候群ウィルムス腫瘍、無虹彩、泌尿生殖器奇形、及び精神遅滞、スミスマジェニス症候群、骨髄異形成症候群、遺伝性粘膜上皮異形成、遺伝性角皮症、遺伝性ニューロパチー、例えばシャルコーマリートゥース病及び神経線維腫症、甲状腺機能低下症、水頭症、発作性障害、例えばシデナム舞踏病Syndenham’schorea及び脳性麻痺、脊椎披裂、無脳症、頭蓋脊椎披裂、先天性緑内障、白内障、感音難聴、並びに対象の任意の組織、器官、又は系、例えば、脳、副腎、腎臓、骨格系又は生殖器系が関与する細胞成長及び分化、胚発生、及び形態形成に関連する任意の障害が挙げられる。一実施形態では、本明細書にerein開示される方法を用いて、核酸配列の発現に関連する障害を検出、病期分類、診断、及び/又は治療し得る。かかる障害の例としては、限定はされないが、内分泌障害、例えば、性腺機能低下症、シーハン症候群、尿崩症、カルマン病、ハンドシュラークリスチャン病、レットレルジーベ病、サルコイドーシス、エンプティセラ症候群、及び小人症を含む下垂体機能低下に関連する障害;先端巨大症、巨人症、及び抗利尿ホルモンADH不適合分泌症候群SIADHを含む下垂体機能亢進症;並びに、甲状腺腫、粘液水腫、細菌感染症に関連する急性甲状腺炎、ウイルス感染症に関連する亜急性甲状腺炎、自己免疫性甲状腺炎橋本病、及びクレチン病を含む甲状腺機能低下に関連する障害;甲状腺中毒症及びその様々な型、グレーブス病、前脛骨粘液水腫、中毒性多結節性甲状腺腫、甲状腺癌、及びプランマー病を含む甲状腺機能亢進に関連する障害;並びに、コン病慢性高カルシウム血症chronichypercalemiaを含む副甲状腺機能亢進に関連する障害;呼吸器疾患、例えば、アレルギー、喘息、急性及び慢性炎症性肺疾患、ARDS、気腫、肺うっ血及び肺浮腫、COPD、間質性肺疾患、及び肺癌;癌、例えば、腺癌、白血病、リンパ腫、黒色腫、骨髄腫、肉腫、テラトカルシノーマ、特に、副腎癌、膀胱癌、骨癌、骨髄癌、脳癌、乳癌、子宮頸癌、胆嚢癌、ガングリオンの癌、消化管癌、心臓癌、腎癌、肝癌、肺癌、筋肉癌、卵巣癌、膵癌、副甲状腺癌、陰茎癌、前立腺癌、唾液腺癌、皮膚癌、脾臓癌、精巣癌、胸腺癌、甲状腺癌、及び子宮癌;及び免疫不全症、例えば、後天性免疫不全症候群AIDS、アジソン病、成人呼吸窮迫症候群、アレルギー、強直性脊椎炎、アミロイドーシス、貧血症、喘息、アテローム性動脈硬化症、自己免疫性溶血性貧血、自己免疫性甲状腺炎、気管支炎、胆嚢炎、接触性皮膚炎、クローン病、アトピー性皮膚炎、皮膚筋炎、真性糖尿病、気腫、リンパ球毒素性偶発性リンパ球減少症、胎児赤芽球症、結節性紅斑、萎縮性胃炎、糸球体腎炎、グッドパスチャー症候群、痛風、グレーブス病、橋本甲状腺炎、過好酸球増加症、過敏性腸症候群、多発性硬化症、重症筋無力症、心筋又は心膜炎症、変形性関節症、骨粗鬆症、膵炎、多発筋炎、乾癬、ライター症候群、関節リウマチ、強皮症、シェーグレン症候群、全身性アナフィラキシー、全身性エリテマトーデス、全身性硬化症、血小板減少性紫斑病、潰瘍性大腸炎、ブドウ膜炎、ウェルナー症候群、癌の、血液透析の、及び体外循環の合併症、ウイルス、細菌、真菌、寄生虫、原虫、及び蠕虫の感染症、並びに外傷が挙げられる。ポリヌクレオチド配列は、サザン解析若しくはノーザン解析、ドットブロット、又は他の膜ベースの技法;PCR法;ディップスティックアッセイ、ピンアッセイ、及びELISAアッセイ;及び患者の体液又は組織を利用して核酸配列発現の変化を検出するマイクロアレイにおいて用いられ得る。かかる定性的又は定量的方法は、当該技術分野において公知である。本発明の特徴付け及び最良の態様PARADIGM:PARADIGMを使用した多次元癌ゲノミクスデータからの患者特異的パスウェイ活性の推論パスウェイに基づく手法の一つの仮説は、パスウェイデータベースに存在する遺伝的相互作用が、癌において検出される遺伝子発現変化の間の相関性を解釈するための情報を有するというものである。例えば、癌関連パスウェイが転写アクチベータAから標的遺伝子Tへのリンクを含む場合、Aの発現がTの発現と正の相関E2E相関を有することが予想される。同様に、Aのコピー数とTの発現との間に正の相関C2E相関が予想される。さらに、Aにおける増幅は、必ずしもAがより高度に発現することを含意するものでなく、それには次にBの上方制御が必要なため、C2E相関はE2E相関より弱いことが予想される。このように、パスウェイの各リンクがデータに関する予想を提供する;多くの一貫したリンクを含むパスウェイは、さらなる考慮事項について関連性を有し得る。本発明者らはこれらの仮定を検証し、NCIパスウェイが、最近のTCGAGBMデータを予測する多くの相互作用を含むことを見出したTheTCGAresearchnetwork2008。本発明者らは、統合された患者データから遺伝的パスウェイの活性を推論するPARADIGMPAthwayRecognitionAlgorithmusingDataIntegrationonGenomicModelsゲノムモデルに関するデータ統合を用いたパスウェイ認識アルゴリズムと呼ばれる手法を開発した。PARADIGM方法は、既知のシグナル伝達パスウェイによって多様なハイスループットゲノミクス情報を統合し、遺伝子活性、複合体、及び細胞プロセスの状態に関する患者特異的ゲノミック推論を提供する。本方法の核心部は、様々なデータソースを組み合わせるために、因子グラフを用いて推論を利用する。元のハイスループットデータ集合の代わりに、又はそれと併せてかかる推論を使用することで、試料を臨床的に関連性のあるサブタイプに分類する能力が向上する。PARADIGMで統合した活性に基づくGBM患者のクラスタリングにより、異なる生存プロファイルと相関する患者サブタイプが明らかとなった。対照的に、発現データ又はコピー数データのいずれかを用いる試料のクラスタリングでは、データ集合において有意なクラスターは何ら明らかとならなかった。PARADIGMは、GBM及び乳癌の双方の腫瘍試料で有意に変化した遺伝子活性のパスウェイ推論を生成する。SPIAと呼ばれる競合するパスウェイ活性推論手法と比較して、本発明者らの方法は、癌関連パスウェイにおいて変化した活性をより少ない偽陽性で特定する。計算効率上、PARADIGMは現在、NCIパスウェイをそのまま使用する。PARADIGMはEMを用いて隠れ量を推論するが、NCIパスウェイにまだ存在しない新しい相互作用を推論しようと試みるものではない。尤度関数を増加させる新しい相互作用を導入するように手法を拡張することが想像され得る。この問題は概して厄介であるが、構造的EMFriedman1997上記などのヒューリスティクスを用いることにより、計算機による探索戦略を用いて相互作用を特定することができる。デノボで新規連結を探索するのでなく、タンパク質間相互作用マップ又は有意な数の発現データ集合において相関する遺伝子ペアから得られた相互作用を提案することにより、探索速度を有意に上げることができる。このパスウェイベースの手法の力は、観測される生存の差の根底にある考えられ得る機構に関して手掛かりを提供し得ることである。情報提供性を有するIPAは、治療標的を示唆したり、又は臨床試験に最適な患者を選択したりするのに有用であり得る。例えば、ErbB2増幅は、薬物トラスツズマブにより治療可能な特定の乳癌型の周知のマーカーである。しかしながら、ErbB2増幅を有する患者の一部は、治療抵抗性の腫瘍を有する。CircleMap表示を調べることにより、ErbB2増幅を有するが、PARADIGMにより推論されるときIPAが不活性又は不変のいずれかである患者を特定し得る。ErbB2増幅を有するが、予測される活性がない患者には、代替治療が検討され得る。将来、より多次元のデータ集合が利用可能となるに伴い、かかるパスウェイ推論がコホート間で一般化されるロバストなバイオマーカーを提供するかどうかを検証することは、有益であろう。乳癌における抗癌化合物に対するサブタイプ及びパスウェイ特異的応答現在、ヒト悪性腫瘍の治療に向けて800種を超える小分子阻害薬及び生物製剤の開発が進められているNewMedicinesDatabase|PHRMA.http://newmeds.phrma.org/2010。これらの薬剤の多くは、腫瘍を正常細胞と区別すると考えられる分子的特徴を標的とし、抗代謝産物及びDNA架橋剤、例えばトラスツズマブ及びラパチニブを含めた、癌のサブセットで調節が解除される分子イベント及び経路を選択的に標的化する広域特異性の従来の治療薬からの範囲に及ぶ例えば、Slamon,D.J.etal.UseofchemotherapyplusamonoclonalantibodyagainstHER2formetastaticbreastcancerthatoverexpressesHER2.NEnglJMed344,783-7922001;Vogel,C.L.etal.Efficacyandsafetyoftrastuzumabasasingleagentinfirst-linetreatmentofHER2-overexpressingmetastaticbreastcancer.JClinOncol20,719-7262002;Rusnak,D.W.etal.Theeffectsofthenovel,reversibleepidermalgrowthfactorreceptor/ErbB-2tyrosinekinaseinhibitor,GW2016,onthegrowthofhumannormalandtumor-derivedcelllinesinvitroandinvivo.MolCancerTher1,85-942001を参照のこと。Effectsofchemotherapyandhormonaltherapyforearlybreastcanceronrecurrenceand15-yearsurvival:anoverviewoftherandomisedtrials.Lancet365,1687-17172005。今日の薬剤開発における一般的な傾向は、従来の薬物と比べて高い効力及び低い毒性を示す標的化薬剤に向かいつつあるSawyers,C.Targetedcancertherapy.Nature432,294-2972004。ERBB2/EGFR阻害薬ラパチニブなど、高い標的特異性を示す薬物もあれば、SRC阻害薬ダサチニブなど、広範囲のキナーゼを阻害する薬物もあるKaraman,M.W.etal.Aquantitativeanalysisofkinaseinhibitorselectivity.NatBiotechnol26,127-1322008。臨床試験は応答の予測因子を含むべきで、且つ試験参加患者を層別化すべきであるという認識が高まっている。多くの分子的に標的化された治療剤は、患者を分類するための明らかな分子特徴を提供するが、ほとんどはそうでない。さらに、腫瘍間での分子的及び生物学的違い、標的化された経路の複雑なクロスカップリング及びフィードバック調節並びに不正確な標的特異性により、基本的機構の予測は多くの場合に複雑となる。分子マーカーに基づく臨床試験の過程では、応答性を有するサブセットが同定され得るが、この手法は手配が困難で、高価であり、且つ最も応答性が高いと思われる特定の部分集団において実験的化合物を最初に試験することができない。実際、現在開発中の薬物の大半は今後乳癌で試験されることはなく、従って乳癌患者の部分集団においてのみ極めて有効な化合物は見逃される確率が高い。有望な手法は、前臨床モデルから得られる応答の予測因子を用いて臨床試験の参加患者を層別化することであり、それにより開発費用が低減され、患者サブセットで特に有効であり得る薬物が特定され得る。細胞株パネルにおける前臨床試験は、初期臨床試験の指針として応答性を有する分子サブタイプを早期に且つ効率的に同定することを可能にする見込みがある。この手法の有用性の根拠は、細胞株パネルがaゲフィチニブに応答性を有するものとしてEGFR突然変異を伴う肺癌を予測しPaez,J.G.etal.EGFRmutationsinlungcancer:correlationwithclinicalresponsetogefitinibtherapy.Science304,1497-15002004、bトラスツズマブ及び/又はラパチニブに応答性を有するものとしてHER2/ERBB2増幅を伴う乳癌を予測しNeve,R.M.etal.Acollectionofbreastcancercelllinesforthestudyoffunctionallydistinctcancersubtypes.CancerCell10,515-5272006;Konecny,G.E.etal.ActivityofthedualkinaseinhibitorlapatinibGW572016againstHER-2-overexpressingandtrastuzumab-treatedbreastcancercells.CancerRes66,1630-16392006、及びcメシル酸イマチニブに応答性を有するものとして変異又は増幅したBCR-ABLを伴う腫瘍を予測するScappini,B.etal.ChangesassociatedwiththedevelopmentofresistancetoimatinibSTI571intwoleukemiacelllinesexpressingp210Bcr/Ablprotein.Cancer100,1459-14712004ことを示す試験に基づく。NCIのディスカバリーセラピューティックプログラムDiscoveryTherapeuticProgramは、この手法を大規模に達成しようとするもので、約60種の癌細胞株のコレクションにおいて分子的特徴と100,000種超の化合物に対する応答との間の関連性を特定しているWeinstein,J.N.Spotlightonmolecularprofiling:IntegromicanalysisoftheNCI-60cancercelllines.MolCancerTher5,2601-26052006;Bussey,K.J.etal.IntegratingdataonDNAcopynumberwithgeneexpressionlevelsanddrugsensitivitiesintheNCI-60celllinepanel.MolCancerTher5,853-8672006。多様な応答を伴う化合物の検出に有用であるものの、NCI60パネルはコレクション内での特定の癌サブタイプの代表性が比較的疎であるため、サブタイプ特異的な応答の検出力に限りがある点は議論の余地がある。例えば、このコレクションには乳癌細胞株は6種しか含まれず、これは既知の多様性を適切に表現するには不十分である。従って本発明者らは、乳癌におけるインビトロ治療化合物応答と分子サブタイプと活性化したシグナル伝達パスウェイとの間の関連性を統計的にロバストに同定するため、約50種の乳癌細胞株のコレクションの使用を図った。ここで本発明者らは、FDA承認済み薬物及び治験中の化合物の双方を含む77種の化合物について、定量的な成長阻害応答とサブタイプを定義する分子的特徴と活性化したパスウェイとの間の関連性の評価を報告する。約半数が異常又はサブタイプ特異性を示す。本発明者らはまた、遺伝子発現及びコピー数データの統合的分析から、観測されるサブタイプ関連応答の一部を特定のパスウェイ活性により説明し得ることも示す。統合された分子プロファイルはDcisにおけるインターロイキンシグナル伝達の歪み及び浸潤性乳癌における予後判定力の向上を明らかにする様々なレベルにおけるハイスループット腫瘍分子プロファイルの蓄積は、世界的に時間及び費用のかかるプロセスとなっている。様々なレベルで遺伝子調節を複合的に分析することにより、複数の上皮癌で調節が解除される特定の生物学的機能及び分子パスウェイが示され、並びに個々に合わせた治療及びモニタリングを行うための新規患者サブグループが明らかとなり得る。本発明者らは、約110人の乳癌患者からの、原発腫瘍、対応する血液由来、且つ既知の微小転移状態有りの新鮮凍結試料から得られるいくつかの分子レベルのハイスループットデータを収集した以下、MicMaデータ集合と称する。これらの患者は、播種性腫瘍細胞DTCの存在、再発についての長期経過観察及び全生存に関する情報を含む900件を超える乳癌症例コホートの一部である。MicMa集合は、全ゲノムmRNA発現Naume,B.etal.,2007,Presenceofbonemarrowmicrometastasisisassociatedwithdifferentrecurrenceriskwithinmolecularsubtypesofbreastcancer,1:160-17、arrayCGHRussnes,H.G.etal.,2010,Genomicarchitecturecharacterizestumorprogressionpathsandfateinbreastcancerpatients,2:38ra472、DNAメチル化Ronneberg,J.A.etal.,2011,Methylationprofilingwithapanelofcancerrelatedgenes:associationwithestrogenreceptor,TP53mutationstatusandexpressionsubtypesinsporadicbreastcancer,5:61-76、全ゲノムSNP及びSNP-CGHVan,LooP.etal.,2010,Allele-specificcopynumberanalysisoftumors,107:16910-169154、全ゲノムmiRNA発現解析EnerlyE,SteinfeldI,KleiviK,LeivonenS,AureMR,RussnesHG,RoennebergJA,JohnsenH,NavonR,RoedlandE,MaekelaeR,NaumeB,PeraelaeM,KallioniemiO,KristensenVN,YakhiniZ,Boerresen-DaleA.miRNA-mRNAintegratedanalysisrevealsrolesformiRNAsinprimarybreasttumors.PLoSONE2011;62:e16915。TP53突然変異状態依存性パスウェイ及びハイスループットペアエンドシーケンシングStephens,P.J.etal.,2009,Complexlandscapesofsomaticrearrangementinhumanbreastcancergenomes,462:1005-1010の並行するパイロット研究で使用されている。これは、単一の研究室によって同じ原発性乳房腫瘍集合に対して行われる包括的なハイスループット分子データコレクションである。以下にこれらの研究の知見を要約し、研究の各々は、mRNA発現をDNAコピー数、DNAメチル化における調節解除又はmiRNA発現のいずれかと統合することを試みている。過去に本発明者ら及び他の研究者らが複数の分子レベルで乳癌機構を調べているが、mRNA、CNA、miRNA、及びメチル化をパスウェイの文脈でモデル化することによりこれらの観点を統合する試みは、ごく僅かであった。本論において本発明者らは、乳癌からのかかるデータを、摂動を受けるパスウェイ及び個別的な表現型特性を有する分子サブタイプの双方の検出と合わせて分析した。ここで考察するMicMaデータ集合において、本発明者らはメチル化プロファイルに基づき3つの大クラスター及び1つの小クラスターを特定した;大クラスターの1つは主として筋上皮起源の腫瘍からなり、他の2つの大クラスターは主に内腔上皮起源の腫瘍を含んだ。クラスターは、TP53突然変異及びER、及びErbB2発現状態、並びに悪性度に関して異なった。パスウェイ解析により、EGF、NGFR及びTNFなどの遺伝子を含むカノニカルなキュレーションしたパスウェイ、樹状細胞成熟及びNF-κBシグナル伝達パスウェイとの有意な関連性が特定された。DCIS及び浸潤癌の試料に対する候補遺伝子のパイロシーケンシングから、ABCB1、FOXC1、PPP2R2B及びPTENが、DCISでメチル化される新規遺伝子として同定された。どの病変が侵襲性となるリスクを伴うかをより良く理解するためには、これらの後成的変化がどのように腫瘍進行の惹起に関与するかを理解することが重要である。本発明者らはまた、MicMaデータ集合におけるmiRNA及びmRNA発現間の関係についても、それらの互いの、及び臨床的特徴との相関性の観点で調査した。本発明者らは、増殖、細胞接着及び免疫応答などのいくつかの細胞プロセスが一定のmiRNAと強く関連することを示すことができた。分子固有サブタイプ間、及び増殖レベルが異なる試料間において、統計的に有意なmiRNA発現差が認められた。本発明者らは、細胞株に対するハイスループットライセートマイクロアレイを用いて増殖の調節におけるmiRNAの役割を検証し、このプロセスの潜在的なドライバーを指摘したEnerlyetal.2001上記。この乳癌患者コホートにおいて10e-6のp値カットオフレベルでTP53突然変異状態によれば濃縮差を示す40個超のKEGGパスウェイが同定された。パスウェイの濃縮差はまた、2つの異なるマイクロアレイプラットフォームに基づく、187個の乳癌試料からなるクロスプラットフォームのデータ集合においても認められた。示差的に濃縮されたパスウェイには、TP53シグナル伝達及び細胞周期などのいくつかの既知の癌パスウェイ、免疫応答及びサイトカイン活性化を含むシグナル伝達パスウェイ、並びに脂肪酸代謝を含む代謝パスウェイJoshietal,2011上記が含まれた。先述の研究の各々は、ハイスループット分子データからペアワイズ式にCNA/mRNA、miRNA/mRNA、DNAメチル化/mRNA、TP53/mRNA生物学的相互作用を導き出そうと試みている。本研究では、本発明者らは、調節解除されたパスウェイに焦点を合わせ、全ての分子レベルを同時に考慮した統合予後判定指標を開発しようと試みた。本発明者らは、ゲノムモデルに関するデータ統合を用いたパスウェイ認識アルゴリズムPARADIGMを適用することにより、様々な遺伝的パスウェイの相対活性を明らかにし、共同の潜在的予後判定能力を評価した。次に、クラスター及びPARADIGMにより特定される調節解除されたパスウェイを別のデータ集合Chin,S.F.etal.,2007,Usingarray-comparativegenomichybridizationtodefinemolecularportraitsofprimarybreastcancers,26:1959-1970で検証し、また、DCIS非浸潤性乳管癌などの、前癌性の新形成のデータ集合Muggerud,A.A.etal.,2010,MoleculardiversityinductalcarcinomainsituDCISandearlyinvasivebreastcancer,4:357-368で試験した。卵巣漿液性癌において高頻度で変化したパスウェイコピー数及び遺伝子発現の双方の統合分析によって有意に変化したパスウェイを同定するため、本発明者らは最近開発したパスウェイ活性推論法PARADIGMPMID:20529912を適用した。計算モデルはコピー数変化、遺伝子発現データ、及びパスウェイ構造を組み込み、パスウェイデータベースに存在する遺伝子、複合体、及び遺伝的プロセス毎に統合パスウェイ活性IPAを生成する。本発明者らは用語エンティティを、それが遺伝子か、複合体か、又は小分子かを問わず、パスウェイにおける任意の分子を指して使用する。エンティティのIPAは、最終活性のみを指す。遺伝子について、IPAはタンパク質の活性状態の推論された活性のみを指し、これは、パスウェイにおける他の遺伝子のコピー数、遺伝子発現、及びシグナル伝達から推論される。本発明者らはPARADIGMを卵巣試料に適用し、NationalCancerInstituteのパスウェイ相互作用データベースNCI-PIDに含まれるパスウェイに存在する多くの異なる遺伝子及びプロセスの変化を見出した。本発明者らは1000回のランダムシミュレーションを用いて推論された変化の有意性を評価し、それらのシミュレーションでは、同じ構造を有するパスウェイが使用されるが、パスウェイにおける種々の点に任意の遺伝子を割り当てた。換言すれば、所与のパスウェイのあるランダムシミュレーションにおいて相互作用の集合は固定されたままであり、従って任意の遺伝子集合がそれらのパスウェイ相互作用によって互いに連結された。全ての試料のIPAの有意性は同じヌル分布に対して評価して、各試料におけるエンティティ毎の有意水準を得た。標準偏差が少なくとも0.1のIPAを、図28にヒートマップとして表示する。表3は、PARADIGMにより見出された順列化試料に関して少なくとも3標準偏差だけ変化するパスウェイを示す。FOXM1転写因子ネットワークは、検証した全てのパスウェイのなかで最多数の試料において変化した-試料間で平均化したとき67のエンティティが活性の変化を有した。それに対して、卵巣コホートにおいて次に最も高いレベルの活性の変化を有したパスウェイとしては、PLK1シグナル伝達イベント27、オーロラBシグナル伝達24、及びトロンボキサンA2受容体シグナル伝達20が含まれた。従って、NCI-PIDのパスウェイのなかでは、卵巣試料に関して、FOXM1ネットワークが他のパスウェイと比べて有意により大きく変化する活性を含む。FOXM1転写因子ネットワークが、最も高い割合の患者試料において、正常対照と比較して腫瘍試料で示差的に変化することが分かった図29。FOXM1は、3つの既知の支配的なスプライシング型を有する多機能性転写因子であり、スプライシング型の各々は、細胞増殖及びDNA修復において様々な役割を有する個別的な遺伝子サブセットを調節する。FOXM1cアイソフォームは、AUKB、PLK1、CDC25、及びBIRC5を含めた、細胞増殖において既知の役割を有するいくつかの標的を直接調節するPMID:15671063。他方で、FOXM1bアイソフォームは、DNA修復遺伝子BRCA2及びXRCC1を含む完全に異なる遺伝子サブセットを調節する。ATMの間接的な制御下にあるCHEK2は、FOXM1発現レベルを直接調節するPMID:17101782。本発明者らは、FOXM1転写因子それ自体のIPAが、他の転写因子のIPAと比べてより高度に変化するかという問いを立てた。FOXM1の活性レベルをNCI-PIDにおける他の203個の転写因子の全てと比較した。NCI集合の他の転写因子と比較しても、FOXM1転写因子は有意に高いレベルの活性を有したp<0.0001;K-S検定ことから、それが重要なシグネチャーであり得ることがさらに示唆される図30。FOXM1はまた、上皮起源の多くの異なる正常組織においても発現するため、本発明者らは、PARADIGMにより同定されるシグネチャーが、他の組織では正常と見なされ得る上皮シグネチャーに起因するかという問いを立てた。これに答えるため、卵管上皮及び卵巣腫瘍組織がマイクロダイセクションされたGEOからの独立したデータ集合をダウンロードしGSE10971PMID:18593983、遺伝子発現をアッセイした。FOXM1のレベルが正常試料と比較して腫瘍試料において有意に高いことが見出されたことから、実際に癌性組織においてはFOXM1調節が、正常上皮組織に見られるものを超えて亢進することが示唆される図31。TCGA卵巣コホート全体が高悪性度漿液性腫瘍に由来する試料を含んだため、本発明者らは、FOXM1シグネチャーが高悪性度漿液性に特異的であるのかという問いを立てた。低悪性度及び高悪性度の双方の漿液性腫瘍の転写プロファイルが決定されているEtemadmoghadametal.2009のデータ集合からEtemadmoghadamD,deFazioA,BeroukhimR,MermelC,GeorgeJ,GetzG,TothillR,OkamotoA,RaederMB,AOCSStudyGroup,HarnettP,LadeS,AkslenLA,TinkerAV,LocandroB,AlsopK,ChiewYE,TraficanteN,FeredayS,JohnsonD,FoxS,SellersW,UrashimaM,SalvesenHB,MeyersonM,BowtellD.IntegratedGenome-WideDNACopyNumberandExpressionAnalysisIdentifiesDistinctMechanismsofPrimaryChemoresistanceinOvarianCarcinomas.ClinicalCancerResearch2009Feb.;154:1417-1427、FOXM1及びその標的のいくつかの対数発現を得た。この独立したデータから、FOXM1及びその標的のいくつかが、低悪性度卵巣癌と比較して漿液性卵巣癌において有意に上方制御されることが確認された図32。FOXM1転写因子ネットワークの25個の遺伝子が、高悪性度疾患において発現がより高い遺伝子を有意な割合で含んだかどうかを決定するため、本発明者らはEtemadmoghadamのデータを使用してスチューデントのt検定を実施した。ゲノム中723個の遺伝子5.4が、有意水準0.05で低悪性度癌と比べて高悪性度癌において有意に上方制御されることが分かったベンジャミニホッホバーグBenjamini-Hochbergの方法を用いて複数の試験に対して補正。FOXM1ネットワークは、その遺伝子のうち13個52が示差的に調節されることが分かった。これは、超幾何検定に基づけば有意な割合であるP<3.81012。従ってFOXM1ネットワーク遺伝子の高発現は、ゲノムにおける典型的な遺伝子の発現と比較したとき、高悪性度疾患と確かに特異的に関連するように思われる。乳癌及び肺癌を含め、多くの異なる癌におけるFOXM1の役割が十分に裏付けられているが、卵巣癌におけるその役割はいまだ研究がなされていない。FOXM1は、3つの既知のスプライシング型を有する多機能性転写因子であり、スプライシング型の各々は、細胞増殖及びDNA修復において様々な役割を有する個別的な遺伝子サブセットを調節する。この分析に関連するFOXM1の相互作用ネットワークの抜粋を図27に示す。FOXM1aアイソフォームは、AUKB、PLK1、CDC25、及びBIRC5を含めた、細胞増殖において既知の役割を有するいくつかの標的を直接調節する。対照的にFOXM1bアイソフォームは、DNA修復遺伝子BRCA2及びXRCC1を含む完全に異なる遺伝子サブセットを調節する。ATMの間接的な制御下にあるCHEK2は、FOXM1の発現レベルを直接調節する。卵巣患者のほとんどにおけるFOXM1の発現増加に加え、小さいサブセットはまた、CBSにより検出されるコピー数増幅の増加も有する19が、計測したゲノムにおける全遺伝子の上位5分位内のコピー数増加を伴う。従って、FOXM1の選択的スプライシング調節はDNA修復と細胞増殖との間の制御スイッチに関与し得る。しかしながら、個別のアイソフォーム活性を区別することは、アイソフォームを区別するエクソン構造及びエクソンアレイプローブの位置によって困難なため、現時点ではこの主張を裏付けるデータは不十分である。これらの試料のmRNAの将来のハイスループットシーケンシングが、示差的なFOXM1アイソフォームレベルの決定を促進し得る。PARADIGMがこの転写因子に集中している最高レベルの活性変化を検出したという観測は、FOXM1が細胞において重要な調節ポイントにあることを示唆している。診断本明細書に記載される方法を用いて、変化した遺伝子発現、mRNAの発現過剰に対する非存在/存在を検出及び定量化し、又は治療介入中のmRNAレベルをモニタし得る。発現の変化に関連する病態、疾患又は障害としては、特発性肺動脈高血圧症、二次性肺高血圧症、細胞増殖性障害、特に退形成乏突起膠腫、星状細胞腫、乏突起星細胞腫、膠芽腫、髄膜腫、神経節神経腫、神経細胞腫瘍、多発性硬化症、ハンチントン病、乳腺腺癌、前立腺腺癌、胃腺癌、転移性神経内分泌癌、非増殖性線維嚢胞性及び増殖性線維嚢胞性乳腺疾患、胆嚢炎及び胆石症、変形性関節症、及び関節リウマチ;後天性免疫不全症候群AIDS、アジソン病、成人呼吸窮迫症候群、アレルギー、強直性脊椎炎、アミロイドーシス、貧血症、喘息、アテローム性動脈硬化症、自己免疫性溶血性貧血、自己免疫性甲状腺炎、良性前立腺肥大症、気管支炎、チェディアック東症候群、胆嚢炎、クローン病、アトピー性皮膚炎、皮膚筋炎、真性糖尿病、気腫、胎児赤芽球症、結節性紅斑、萎縮性胃炎、糸球体腎炎、グッドパスチャー症候群、痛風、慢性肉芽腫性疾患、グレーブス病、橋本甲状腺炎、過好酸球増加症、過敏性腸症候群、多発性硬化症、重症筋無力症、心筋又は心膜炎症、変形性関節症、骨粗鬆症、膵炎、多嚢胞性卵巣症候群、多発筋炎、乾癬、ライター症候群、関節リウマチ、強皮症、重症複合型免疫不全症SCID、シェーグレン症候群、全身性アナフィラキシー、全身性エリテマトーデス、全身性硬化症、血小板減少性紫斑病、潰瘍性大腸炎、ブドウ膜炎、ウェルナー症候群、血液透析、体外循環、ウイルス、細菌、真菌、寄生虫、原虫、及び蠕虫感染症;プロラクチン産生障害、不妊症、例えば、卵管疾患、排卵異常、及び子宮内膜症、発情周期の乱れ、月経周期の乱れ、多嚢胞性卵巣症候群、卵巣過剰刺激症候群、子宮内膜腫瘍又は卵巣腫瘍、子宮筋腫、自己免疫障害、子宮外妊娠、及び催奇形性;乳癌、線維嚢胞性乳腺疾患、及び乳汁漏出症;精子形成の乱れ、精子の生理異常、良性前立腺肥大症、前立腺炎、ペイロニー病、インポテンス、女性化乳房;日光角化症、動脈硬化、滑液包炎、硬変、肝炎、混合結合組織病MCTD、骨髄線維症、発作性夜間ヘモグロビン尿症、真性赤血球増加症、原発性血小板血症、癌の合併症、癌、例えば、腺癌、白血病、リンパ腫、黒色腫、骨髄腫、肉腫、テラトカルシノーマ、特に、副腎癌、膀胱癌、骨癌、骨髄癌、脳癌、乳癌、子宮頸癌、胆嚢癌、ガングリオンの癌、消化管癌、心臓癌、腎癌、肝癌、肺癌、筋肉癌、卵巣癌、膵癌、副甲状腺癌、陰茎癌、前立腺癌、唾液腺癌、皮膚癌、脾臓癌、精巣癌、胸腺癌、甲状腺癌、及び子宮癌が挙げられる。別の態様において、本発明の核酸。本明細書に記載される方法を用いて、変化した遺伝子発現;mRNAの非存在、存在、又は過剰発現を検出及び定量化し;又は治療介入中のmRNAレベルをモニタし得る。発現の変化に関連する障害としては、静坐不能akathesia、アルツハイマー病、健忘、筋萎縮性側索硬化症、失調、双極性障害、緊張病、脳性麻痺、脳血管疾患クロイツフェルトヤコブ病、認知症、うつ病、ダウン症候群、遅発性ジスキネジー、ジストニー、癲癇、ハンチントン病、多発性硬化症、筋ジストロフィー、神経痛、神経線維腫症、ニューロパチー、パーキンソン病、ピック病、網膜色素変性症、統合失調症、季節性情動障害、老年痴呆、脳卒中、トゥレット症候群並びに腺癌、黒色腫、及びテラトカルシノーマを含む癌、特に脳の癌が挙げられる。遺伝子発現に関連する病態、疾患又は障害の診断の根拠を提供するため、正常な又は標準の発現プロファイルを確立した。これは、正常な対象動物又はヒトのいずれかから採取した生体試料を、ハイブリダイゼーション条件又は増幅条件下でプローブと組み合わせることにより達成され得る。標準的なハイブリダイゼーションは、正常な対象を用いて得られた値を、既知量の実質的に精製された標的配列が用いられる実験の値と比較することにより定量化され得る。このようにして得られた標準の値は、特定の病態、疾患、又は障害に症候を示す患者の試料から得られた値と比較され得る。標準の値から特定の病態に関連する値の方への偏差を使用して、当該の病態が診断される。また、かかるアッセイを用いて、動物試験及び臨床実験における特定の治療的処置レジメンの効力を評価し、又は個別の患者の治療をモニタし得る。病態の存在が確定し、治療プロトコルが開始されると、診断アッセイが定期的に繰り返され、患者の発現レベルが正常な対象で観測されるレベルに近付き始めるかどうかが決定され得る。継続的なアッセイから得られる結果を使用して、数日から数ヶ月の範囲の期間にわたる治療効力が示され得る。モデル系動物モデルは、ヒトのものと同様の毒性応答を示し、且つ曝露条件がヒト曝露に適切である場合に、バイオアッセイとして用いることができる。哺乳動物は最も一般的なモデルであり、低コスト、利用し易さ、及び豊富な毒性学文献から、ほとんどの毒性試験がラット又はマウスなどのげっ歯類で行われている。近交系げっ歯類系統は、対象遺伝子の過小発現又は過剰発現の生理学的影響の研究並びに疾患の診断及び治療方法の開発に好都合なモデルを提供する。特定の遺伝子例えば、乳中に分泌されるものを過剰発現するように同系交配した哺乳動物もまた、当該の遺伝子により発現する好都合なタンパク質供給源として機能し得る。毒性学毒性学は、薬剤の生物系に対する影響の研究である。毒性試験の大半はラット又はマウスで行われ、ヒトの健康に対するそれらの薬剤の影響予測を促進する。生理機能、行動、恒常性過程、及び致死性の定性的及び定量的変化の観測を用いて毒性プロファイルが生成され、薬剤への曝露後のヒトの健康に対する影響が評価される。遺伝毒性学では、遺伝子突然変異を生じさせる薬剤の能力が同定及び分析される。遺伝毒性薬剤は、通常、核酸との相互作用を促進する共通の化学的又は物理的特性を有し、且つ染色体異常が子孫に伝わるとき最も有害である。遺伝毒性学的研究では、受胎前の親、妊娠中の母親、又は発育中の生物体のいずれかに投与される場合に子孫における構造的又は機能的異常の頻度を増加させる薬剤が同定され得る。マウス及びラットは、統計学的要件を満足するのに必要な生物個体数を産生する生殖周期が短いため、これらの試験において最も高頻度で使用される。急性毒性試験は、薬剤を対象に単回投与して、その薬剤の症候学又は致死性を決定することに基づく。3つの実験が行われる:a初回用量範囲を判断する実験、b有効用量の範囲を絞り込む実験、及びc用量反応曲線を確定する最終実験。長期毒性試験は、薬剤の反復投与に基づく。これらの試験にはラット及びイヌが一般的に用いられ、異なるファミリーの種からのデータが提供される。発癌を除いて、薬剤を高用量濃度で3~4ヶ月間にわたり連日投与すると、成体動物におけるほとんどの毒性形態が明らかになるというエビデンスが多数ある。1年以上の継続期間を伴う慢性毒性試験は、薬剤の毒性の欠如又は発癌可能性のいずれかを実証するために用いられる。試験がラットで行われる場合、最低3つの試験群1つの対照群が用いられ、動物は、開始時及び実験中所定間隔で検査及び観察される。トランスジェニック動物モデル対象遺伝子を過剰発現又は過小発現するトランスジェニックげっ歯類を同系交配し、ヒト疾患のモデリング又は治療剤若しくは毒物の試験に用い得る参照により本明細書に援用される米国特許第4,736,866号明細書;同第5,175,383号明細書;及び同第5,767,337号明細書を参照のこと。ある場合には、導入された遺伝子が、胎児発生中又は出生後に特定の時点で特定の組織型において活性化され得る。実験的薬物治療の投与前、投与中、及び投与後のトランスジェニック動物における表現型又は組織特異的mRNA発現を分析することにより、導入遺伝子の発現が観察される。胚性幹細胞げっ歯類の胚から単離された胚性幹細胞ESは、胚を形成する可能性を保持している。ES細胞を担体の胚内に置くと、それらの細胞は正常な発生を再開し、生産する動物の全ての組織に寄与する。ES細胞は、実験的ノックアウト及びノックインげっ歯類系統の作製に使用される好ましい細胞である。マウス129/SvJ細胞株などのマウスES細胞は初期マウス胚から誘導され、当該技術分野において公知の培養条件下で成長される。ノックアウト系統用のベクターは、生体内で転写及び/又は翻訳を乱すマーカー遺伝子を含むように修飾された疾患遺伝子候補を含む。ベクターは、当該技術分野において公知の電気穿孔、リポソーム送達、マイクロインジェクションなどの形質転換法によってES細胞に導入される。内因性げっ歯類遺伝子は、細胞分裂中の相同組換え及び組込みによって破壊疾患遺伝子に置き換えられる。形質転換ES細胞が同定され、及び好ましくは、C57BL/6マウス系統由来のものなどのマウス細胞胚盤胞に顕微注入される。胚盤胞は偽妊娠の雌親に外科的に移植され、得られたキメラの子孫の遺伝子型同定及び繁殖により、ヘテロ接合又はホモ接合系統が産生される。ES細胞はまた、神経細胞、造血系、及び心筋細胞などのインビトロでの様々な細胞型及び組織の分化の研究にも用いられるBainetal.1995Dev.Biol.168:342-357;WilesandKeller1991Development111:259-267;及びKlugetal.1996J.Clin.Invest.98:216-224。最近の進展では、ヒト胚盤胞に由来するES細胞もまた、内胚葉、中胚葉、及び外胚葉ectodermnal細胞型を含む8つの別個の細胞系統に分化するようインビトロで操作し得ることが実証されているThomson1998Science282:1145-1147。ノックアウト解析遺伝子ノックアウト解析では、ヒト疾患遺伝子候補の領域が、ネオマイシンホスホトランスフェラーゼ遺伝子neo;例えば、Capecchi1989Science244:1288-1292を参照などの非哺乳動物遺伝子を含むように酵素的に修飾される。挿入されたコード配列が標的化された遺伝子の転写及び翻訳を阻止し、疾患候補タンパク質の生化学合成を妨げる。修飾された遺伝子は、培養胚性幹細胞上記に記載されるに形質転換され、形質転換細胞がげっ歯類胞胚に注入され、及び胞胚が偽妊娠の雌親に移植される。トランスジェニック子孫を異種交配させると、ホモ接合近交系が得られる。ノックイン解析胚発生初期に存在する全能性ES細胞を使用して、ノックインヒト化動物ブタ又はヒト疾患のトランスジェニック動物モデルマウス又はラットを作製することができる。ノックイン技術によってヒト遺伝子のある領域が動物ES細胞に注入され、ヒト配列が組換えによって動物細胞ゲノムに組み込まれる。組み込まれたヒト遺伝子を含む全能性ES細胞は、上記に記載したとおり操作される。近交系動物が試験及び処置され、類似のヒト病態に関する情報が得られる。このような方法は、いくつかのヒト疾患のモデル化に用いられている例えば、Leeetal.1998Proc.Natl.Acad.Sci.95:11371-11376;Baudoinetal.1998GenesDev.12:1202-1216;及びZhuangetal.1998Mol.CellBiol.18:3340-3349を参照のこと。非ヒト霊長類モデル動物試験の分野は、生理学、遺伝学、化学、薬理学及び統計学などの基礎科学のデータ及び方法論を取り扱う。こうしたデータは、非ヒト霊長類に対する治療剤の効果の評価において、それがヒトの健康に関係し得るため、最も重要である。サルは、ワクチン及び薬物評価においてヒトの代理として使用され、その応答は同様の条件下でのヒト曝露に関係する。カニクイザルマカカファシキュラリスMacacafascicularis、マカカムラタMacacamulata及びコモンマーモセットカリトリクスヤックスCallithrixjacchusが、これらの研究で使用される最も一般的な非ヒト霊長類NHPである。NHPコロニーの育成及び維持は巨費を伴うため、通常、初期研究及び毒性学的試験はげっ歯類モデルで行われる。薬物中毒などの行動測定を用いる試験では、NHPは第一選択の試験動物である。加えて、NHP及びヒト個体は、多くの薬物及び毒素に対して示差的感受性を示し、そうした作用物質のエクステンシブメタボライザー及びプアーメタボライザーとして分類される。本発明の例示的使用テーラーメイド医療は、主として利益を受ける可能性が高い患者に1つ又は複数の特異的な治療を送達するのに有望である。本発明者らは、治療化合物の約半数が、臨床的に関連する転写上又はゲノム上の乳癌サブタイプの1つ又は複数において選択的に有効であることを示した。これらの知見は、乳癌治療において応答に関連する分子サブタイプを定義することの重要性を支持している。本発明者らはまた、細胞株に関する転写及びゲノムデータのパスウェイ統合が、観測されるサブタイプ特異的応答についての機構的説明を提供するサブネットワークを明らかにすることも示す。細胞株と腫瘍との間のサブネット活性の比較分析は、サブタイプ特異的サブネットワークの大部分が細胞株と腫瘍との間で保存されていることを示す。これらの分析から、十分に特徴付けられた細胞株パネルにおける実験的化合物の前臨床スクリーニングにより、初期相の臨床試験における感受性濃縮に使用し得る応答に関連する分子シグネチャー候補を同定し得るという着想が裏付けられる。本発明者らは、このインビトロ評価手法が、化合物の臨床開発開始前に応答性の腫瘍サブタイプが同定される可能性を高め得るもので、それにより費用が低減され、最終的にFDA承認が下りる確率が高まり、及び場合により、応答性が低いと思われる患者を治療することに伴う毒性が回避されることを示唆する。この試験において、本発明者らは、転写サブタイプを定義する分子シグネチャーのみを評価し、反復されるゲノムCNAを選択した。本発明者らは、遺伝子突然変異、メチル化及び選択的スプライシングなどのさらなる分子的特徴を解析に含めると、本手法の検出力及び精度が増加し得ると予想する。同様に、細胞株パネルのサイズが大きくなるほど、パネル内でそれほど一般的でない分子パターンを評価する検出力が増加し、ヒト乳癌に存在する多様性のより完全な範囲を表現する確率が高まり得る。乳癌の発生は、生来の細胞と適応免疫細胞との双方の存在の有意な増加により特徴付けられ、B細胞、T細胞、及びマクロファージが、腫瘍性間質に存在する最も豊富な白血球に相当するDeNardoDG,CoussensLM.Inflammationandbreastcancer.Balancingimmuneresponse:crosstalkbetweenadaptiveandinnateimmunecellsduringbreastcancerprogression.BreastCancerRes.2007;94:212。腫瘍間質及び血清における高い免疫グロブリンIg値、及び過剰な濾胞性B細胞、調節性T細胞の存在量の増加、及び原発性腫瘍又はリンパ節におけるCD4/CD8又はTH2/TH1Tリンパ球の比の高さが、腫瘍悪性度、病期、及び全患者生存と相関することが示されているBates,G.J.etal.,2006,QuantificationofregulatoryTcellsenablestheidentificationofhigh-riskbreastcancerpatientsandthoseatriskoflaterelapse,24:5373-5380;細胞傷害性Tリンパ球CTL及びナチュラルキラーNK細胞を含め、白血球によっては抗腫瘍活性を呈し34Dunn,G.P.,Koebel,C.M.,andSchreiber,R.D.,2006,Interferons,immunityandcancerimmunoediting,6:836-848、マスト細胞、B細胞、樹状細胞、顆粒球、及びマクロファージなど、他の白血球は、腫瘍進行を妨害するか、或いは増強するその能力を通して、より二極化した役割を呈する35deVisser,K.E.andCoussens,L.M.,2006,Theinflammatorytumormicroenvironmentanditsimpactoncancerdevelopment,13:118-137。これらの研究における最も傑出した知見は、免疫応答TCR及びインターロイキンシグナル伝達、IL4、IL6、IL12及びIL23シグナル伝達における摂動の同定であり、これは予後判定値によるサブクラスの分類をもたらすものであった。ここで本発明者らは、これらのイベントがハイスループット分子データに反映され、乳房腫瘍の分子的な細分類に強力に干渉するというエビデンスを提供する。本開示はまた、HGS-OvCaにおける異常についての初めての大規模な統合的見解も提供する。概して、突然変異スペクトルは意外なほど単純であった。TP53における突然変異が優勢であり、HGS-OvCaの少なくとも96に起こった一方、BRCA1/2は、生殖細胞系突然変異と体細胞突然変異とを組み合わせたことで、腫瘍の22で変異した。7つの他の有意に変異した遺伝子が同定されたが、HGS-OvCaの2~6に過ぎなかった。対照的に、HGS-OvCaは著しい度合いのゲノム配列の乱れdisarrayを実証した。高頻度のSCNAは、膠芽腫による先のTCGAの知見46と顕著な対照をなし、TCGAの知見では、反復的に変異する遺伝子がより多く、染色体腕レベルの又は限局的なSCNAがはるかに少なかった図37A。HR成分を含む推定DNA修復遺伝子における突然変異及びプロモーターメチル化の高率発生は、SCNAの高率発生を説明し得る。突然変異スペクトルは、HGS-OvCaを他のOvCa組織学的サブタイプとは完全に異なるものとして特徴付ける。例えば、明細胞OvCaはTP53突然変異が少ないが、反復されるARID1A及びPIK3CA突然変異を有し47~49;類内膜OvCaは高頻度のCTTNB1、ARID1A、及びPIK3CA突然変異及びより低率のTP53を有する一方48、49、粘液性OvCaはKRAS変異の発生を有する50。卵巣癌サブタイプ間のこれらの差は、病因学的効果と系統効果との組み合わせを反映していると思われ、サブタイプ別のケアによって卵巣癌転帰を改善する機会を表す。新規治療手法の特定は、TCGAの主目的である。HR欠損を有するHGS-OvCaの約50はPARP阻害薬によって利益を受け得る。さらに、共通に調節解除されたパスウェイ、RB、RAS/PI3K、FOXM1、及びNOTCHは、治療的攻撃の機会を提供する。最後に、反復性の増幅領域における22個の遺伝子について阻害薬は既に存在し実施例XIII以下を参照のこと、これらの標的遺伝子が増幅される場合のHGS-OvCaにおける評価を正当化する。総じてこれらの発見は、異常な遺伝子又はネットワークを検出し、それらの特定の異常に有効となるよう選択される治療によって標的化するHGS-OvCa治療手法の布石となる。さらなる実施形態では、まだ開発されていない任意の分子生物学的技法において、そうした新しい技法が、限定はされないがトリプレット遺伝子コード及び特異的塩基対相互作用などの特性を含め、現在公知の核酸分子の特性に頼るものであるならば、ポリヌクレオチド核酸が用いられ得る。本発明は、以下の例を参照することにより、さらに容易に理解されるであろう。これらの例は、本発明の特定の態様及び実施形態を例示するために含まれるに過ぎず、限定として含まれるものではない。実施例I:データソースChin2007上記の乳癌コピー数データを、NCBIGeneExpressionOmnibusGEOから受託番号GPL5737に基づき、GSE8757からの関連するアレイプラットフォームアノテーションと共に入手した。プローブアノテーションを、UCSC癌ゲノミクスブラウザUCSCCancerGenomicsBrowserZhu:2009、上記で表示し、続いて分析するため、BED15フォーマットに変換した。アレイデータをプローブIDによってプローブアノテーションにマッピングした。Naderi2007、上記からの対応する発現データを、EBIのMIAMIExpressから受託番号E-UCon-1を使用して入手した。Human1AV2のプラットフォームアノテーション情報を、Agilent社ウェブサイトから入手した。発現データをプローブレベル中央値で正規化し、プローブIDによってHUGO遺伝子名にマッピングした。全てのデータを、全ての試料-プローブ値を含む順位付け手順を用いてノンパラメトリックに正規化し、順位に基づき各遺伝子-試料ペアについて符号付きp値を得た。0.05の最大p値を使用して、有意に変化した遺伝子-試料ペアを決定した。TCGAの膠芽腫データを、AffymetrixU133Aプラットフォームにおいて230例の患者試料及び10例の隣接する正常組織についての遺伝子発現を提供するTCGADataPortalから入手した。患者試料のプローブを、各プローブの正常値中央値を減じて正常組織に対して正規化した。加えて、同じ患者集合についてのCBSでセグメント化Olshen:2004上記1618頁されたコピー数データを入手した。双方のデータ集合を、乳癌データと同じ手順を用いてノンパラメトリックに正規化した。実施例II:パスウェイ概論本発明者らは、NationalCancerInstituteパスウェイ相互作用データベースNCIPIDから利用可能なキュレーションされたパスウェイの集合を集めたSchaefer:2009上記。各パスウェイは、内因性及び外因性の細胞下レベル、細胞レベル、組織レベル、又は生物体レベルのイベント及び表現型を表す高次の生体分子過程の周りに論理的にまとめ合わされる相互作用の集合を表す。BioPAXレベル2フォーマットのパスウェイをダウンロードした。全てのエンティティ及び相互作用を、RasqalRDFエンジンを使用してSPARQLクエリーで抽出した。本発明者らは、3つの物理的実体タンパク質コード遺伝子、小分子、及び複合体、遺伝子ファミリー、及び抽象的プロセスを含む5つの異なるタイプの生物学的実体エンティティを抽出した。BioPAXタンパク質のクロスリファレンスが異なる遺伝子からのタンパク質を挙げるときは常に、遺伝子ファミリーを作成した。遺伝子ファミリーは、いずれの単一の遺伝子も特定の機能を果たすのに十分である遺伝子のまとまりを表す。例えば、冗長な役割を有する相同体、及び互いを機能的に補うことが認められる遺伝子は、ファミリーに組み入れられる。抽出により、種々のタイプを説明するアノテーション付きのパスウェイで使用されるあらゆるエンティティ及び相互作用のリストが作られた。本発明者らは、アポトーシスなどの、NCIコレクションに存在し得る一般的なプロセスを指す抽象的プロセスもまた抽出した。例えば、p53腫瘍抑制遺伝子が関与する相互作用を詳述するパスウェイは、機械学習による分類用の特徴として利用することができるアポトーシス及び老化へのリンクを含む。予想どおり、C2E相関は中程度であったが、活性化相互作用間で、偶然で予想されるものと比べて正の相関に顕著な濃縮を有した図3。E2E相関はさらに強く、同様に濃縮した。従って、特徴付けしにくい癌のこの例であっても、パスウェイ相互作用の有意なサブセットがゲノム変化を遺伝子発現における調節と連結し、パスウェイレベルの手法が追及するに値するという着想が裏付けられる。実施例III:生物学的パスウェイのモデル化及び予測本発明者らは、初めに各NCIパスウェイを個別の確率モデルに変換した。p53アポトーシスパスウェイのごく一部の小例を図2に示す。NCIのパスウェイ図は、隠れ状態と観測状態との双方を含む因子グラフに変換した。この因子グラフは、エンティティ間の既知の相互作用を表す構造により、遺伝子関連及び生物学的プロセス関連の状態情報に関する観測を統合する。生物学的パスウェイを因子グラフで表現するため、本発明者らは、特定のmRNA又は複合体など、細胞中のエンティティの状態を変数を用いて記述し、及びそれらのエンティティ間の相互作用及び情報フローを、因子を用いて表現した。これらの変数は各エンティティの示差的状態\textit{differential}を、分子エンティティの直接的な濃度ではなく、対照又は正常レベルとの比較で表現する。この表現により、多くの場合に遺伝子の示差的状態を直接計測するか、或いは直接的な計測値をマッチングした対照との相対的な計測値に変換するかのいずれかであるDNAマイクロアレイによって検出される遺伝子発現などの、多くのハイスループットデータ集合をモデル化することが可能になる。また、遺伝子間の多様な調節関係についても可能になる。例えば、p53のMDM2媒介性ユビキチン依存性分解を表す相互作用は、p53タンパク質レベルを阻害する活性化MDM2としてモデル化することができる。因子グラフは、各エンティティのランダム変数X={x1,x1,....,xn}と、エンティティを制約して互いの関数として生物学的に意味のある値をとらせるm個の非負関数、すなわち因子の集合とを使用して細胞の状態を符号化する。j番目の因子φjは、エンティティの部分集合Xj⊂Xの確率分布を定義する。エンティティ及び因子のグラフ全体は、以下のとおりに全エンティティに関する同時確率分布を符号化する:式中、Z=ΠjΣsxjφjSは正規化定数であり、SXは、SがXにおける変数のセッティングであることを表す。各エンティティは、対照レベル例えば、正常組織で計測されるとおりのと比べて活性化、基準どおり、又は不活性化に対応する3つの状態のうちの1つをとることができ、それぞれ1、0、又は-1として符号化され得る。状態は、エンティティの種類例えば、遺伝子、タンパク質等に応じて様々に解釈され得る。例えば、活性化されたmRNAエンティティは過剰発現を表す一方、活性化されたゲノムコピーエンティティはゲノムに3コピー以上が存在することを表す。図2は、単一のタンパク質コード遺伝子の因子グラフの概念モデルを示す。パスウェイ中の各タンパク質コード遺伝子Gについて、エンティティを導入することによりゲノムのコピー数GDNA、mRNA発現GmRNA、タンパク質レベルGprotein、及びタンパク質活性Gprotein図2でDNA、mRNA、タンパク質、及び活性と表示された楕円を表現する。パスウェイの化合物、タンパク質複合体、遺伝子ファミリー、及び抽象的プロセス毎に、分子タイプが活性の単一変数を含める。図2の例は一つのプロセスアポトーシスのみを示すが、実際には多くのパスウェイが、遺伝子活性の出力例えば、アポトーシス及び老化から入力例えば、DNA損傷に至る全てを表現する複数のかかるプロセスを有する。因子の構成を単純にするため、本発明者らは初めにこのパスウェイを有向グラフに変換し、そのグラフの各エッジに正又は負のいずれかの影響ラベルを付与する。第一に、タンパク質コード遺伝子G毎に、GDNAからGmRNAに、GmRNAからGproteinに、及びGproteinからGproteinにラベル正を有するエッジを加え、そのコピー数からその活性化形態のタンパク質産物の存在に至る遺伝子の発現を反映する。パスウェイのいずれの相互作用も、有向グラフの単一のエッジに変換される。この有向グラフを使用して、次に本発明者らは、因子グラフを特定する因子のリストを作成する。変数xi毎に単一の因子φXiを加え、式中、Xi={xi}∪{Parents}xi}、及びParentsxiは、有向グラフにおけるxiの全ての親を指す。全ての値のあるセッティングに対する因子の値は、xiがParentsxiのセッティングによるその期待値に一致するかどうかに依存する。この研究について、期待値は親変数の過半数の票に設定される。親が正のエッジにより連結される場合、それは、因子の値に対してそれ自体の状態の1倍の票に寄与する。逆に、親が負のエッジにより連結される場合、そのとき変数はそれ自体の状態の-1倍を投票する。最小のラベルが付されたエッジによってxiに接続される変数は、単一の票を獲得し、その票の値はそれらの変数の最小値であり、AND状の連結をもたらす。同様に、最大のラベルが付されたエッジによってxiに接続される変数は、単一の票を獲得し、その票の値はそれらの変数の最大値であり、OR状の連結をもたらす。ゼロの票は棄権票として取り扱う。投票がない場合、期待される状態はゼロである。そうでない場合、過半数の票が期待される状態であり、及び1と-1とが同点であると、抑制因子及び欠失により重点が置かれて-1の期待される状態となる。期待される状態のこの定義を所与として、φixi,Parentsxiは以下のとおり特定される:ここに示す結果について、εは0.001に設定したが、イプシロンの選択における大きさの程度の違いは結果に有意な影響を及ぼさなかった。最後に、因子グラフに観測変数及び因子を加えて、パスウェイ及び多次元機能的ゲノミクスデータの統合を完了する図2。離散化した機能的ゲノミクスデータ集合の各々が、タンパク質コード遺伝子の分子タイプ一つと関連付けられる。コピー数変化のアレイCGH/SNP推定はゲノム型に関連付けられる。遺伝子発現データはmRNA型に関連付けられる。ここでの結果に提示しないが、将来の拡張には、mRNA型のDNAメチル化データ、並びにタンパク質型及び活性型のプロテオミクス及び遺伝子リシーケンシングデータが含まれ得る。各観測変数はまた3つ一組の値である。次に記載するとおり、各観測データ型に関連付けられる因子は全エンティティ間で共有され、データから学習される。実施例IV:推論及びパラメータ推定割り当ての集合D={x1=s1,x2=s2,x2,....,xk=sk}が、添え字1~kが付与された観測変数に関する患者のデータの完全集合を表すとする。{SDX}が、Dにおける割り当てと一致する変数Xの集合のあらゆる可能な割り当ての集合を表すとする;すなわち隠れ変数は変化し得るが、任意の観測変数x1はDにおけるその割り当てに固定される。患者データを所与として、特定の隠れエンティティx1が状態aであり得るかどうか、例えば、どの程度TP53のタンパク質活性が-1不活性化されるであり得るか、又はアポトーシスが1活性化されるであり得るかを推定しようとする。これを行うには、患者のデータを観測する前に事象の事前確率を計算しなければならない。Aiaが単集合の割り当て集合{x1=a}を表し、及びφが完全に規定された因子グラフである場合、この事前確率は:であり、式中、Zは方程式1に導入した正規化定数である。同様に、患者の全ての観測と共にx1が状態aである確率は:である。本発明者らは、パスウェイの大部分にHUGIN最新版によるジャンクションツリー推論アルゴリズムを使用した。患者当たりの推論が3秒間より長くかかるパスウェイについては、本発明者らは、逐次更新、収束許容値109、及び最大反復10,000回で信念伝搬法を使用する。全ての推論は、対数定義域に対するものとしての実数定義域で実施し、libDAIMooij:2009上記で実施した。観測因子のパラメータの学習に、本発明者らは期待値最大化EMアルゴリズムDempster1977上記を用いる。簡潔に言えば、EMは、隠れ変数の確率を推論することと、隠れ変数の確率を所与とする尤度を最大化するようにパラメータを変更することとを反復して行うことにより、隠れ変数を含むモデルのパラメータを学習する。本発明者らはコードを書いてlibDAIに投入し、EMを実行した。各パスウェイについて、患者毎の因子グラフを作成し、患者のデータを適用し、尤度の変化が0.1未満になるまでEMを実行した。各パスウェイから学習したパラメータを平均化し、次にそれらのパラメータを使用して、各変数についての最終的な事後信念を計算した。推論後、活性分子タイプを有する各変数について統合パスウェイ活性を出力した。方程式2及び3からの数量を使用して、エンティティiの活性が上昇又は低下するという信念が患者のデータによって増加する程度を反映する対数尤度比を計算した:次に、として、対数尤度比に基づき遺伝子iについての単一の統合パスウェイ活性IPAを計算した。直感的には、IPAスコアは、対数尤度比Lの符号付きの類似形を反映する。遺伝子が活性化される可能性が高い場合、IPAはLに設定される。或いは、遺伝子が不活性化される可能性が高い場合、IPAは対数尤度比の負値に設定される。遺伝子が変わらない可能性が最も高い場合、IPAはゼロに設定される。各パスウェイは、他のパスウェイとは独立して分析される。従って遺伝子は、一つの遺伝子が、それが現れる各パスウェイへと、複数の推論に関連付けられることができる。同じ遺伝子についての推論が異なることは、遺伝子のパスウェイ文脈に応じた代替的なデータ解釈と見ることができる。実施例V:有意性評価本発明者らは、データの2つの異なる順列によりIPAスコアの有意性を評価する。within順列について、初めにランダムな実サンプルを選択し、次に同じパスウェイの範囲内からランダムな遺伝子を選択することによって、新しいデータタプルすなわち対応した遺伝子発現と遺伝子コピー数の選択を、パスウェイの各遺伝子についてタプルが選択されるまで行うことにより、順列化したデータサンプルを作成する。any順列も手順は同じであるが、但しランダムな遺伝子を選択するステップは、ゲノムのどこからでも遺伝子を選択し得る。双方の順列タイプについて、1,000個の順列サンプルが作成され、順列化サンプル毎に摂動スコアが計算される。順列化サンプルからの摂動スコアの分布をヌル分布として用いて、真の試料の有意性を推定する。実施例VI:シグナル伝達パスウェイ影響解析SPIATarca2009,上記によるシグナル伝達パスウェイ影響解析SPIAをCで実行して実行時間を低減し、本発明者らの解析環境に適合するものとした。本発明者らはまた、より詳細な出力を提供する能力も追加して、SPIA出力とPARADIGM出力とを直接比較することができるようにした。本発明者らのSPIAバージョンは、パスウェイのエンティティ毎に累積摂動及び摂動係数を出力することができる。このコードは、要求に応じて利用可能である。実施例VII:デコイパスウェイ各癌データ集合についてデコイパスウェイ集合を作成した。各NCIパスウェイを使用してデコイパスウェイを作成し、これは同じ構造からなるが、パスウェイのあらゆる遺伝子をRefGeneのランダム遺伝子に代えて用いた。全ての複合体及び抽象的プロセスは同じままとし、実パスウェイ及びデコイパスウェイの双方を含むパスウェイの集合に対してPARADIGM及びSPIAの双方についての有意性分析を実行した。パスウェイは各方法内で順位付けし、全パスウェイに対する実パスウェイの割合を計算して視覚化した。実施例VIII:クラスタリング及びカプランマイヤー分析重心連結による非中心化相関階層的クラスタリングを、Eisen1998上記1621頁の方法を用いて膠芽腫データに対して実行した。75例の患者試料間で少なくとも0.25のシグナルを有するIPAのみをクラスタリングに使用した。目視検査により4つの明らかなクラスターが現れ、カプランマイヤー分析に使用した。Rを使用してカプランマイヤー曲線を計算し、ログランク統計によってp値を求めた。実施例IX:PARADIGMの検証EM訓練手順の品質を評価するため、遺伝子発現とコピー数とのタプルE,Cが遺伝子及び患者にわたり順列化されているヌルデータ集合と比べた実際の患者データを使用して、EMの収束を比較した。予想どおり、PARADIGMは、真のデータ集合ではヌルデータ集合と比べてはるかに急速に収束した。例として、遺伝子AKT1についてのIPAをEM反復の関数としてプロットした図4。活性が最初の数回の反復で急速に収束することが分かる。EMは実際の患者データで訓練したとき活性レベルに急速に収束したが、一方、ランダムデータを入れたときは活性不変に収束した。この収束は、パスウェイ構造及び推論が、統合された患者データにおける活性パターンを成功裏に特定可能であることを示唆する。本発明者らは次に、乳癌及びGBMの双方のコホートに対してPARADIGMを実行した。本発明者らは統計的シミュレーション手順を開発し、負の分布から予想され得るものと比べてどのIPAが有意に異なるかを決定した。負の分布は、パスウェイにおける全ての患者間及び遺伝子間で順列化することにより作成した。経験的に、各遺伝子がネットワークにより決定される異なるトポロジー的文脈を有するという事実の補正を促進するには、パスウェイの遺伝子間のみでの順列化が必要であることが分かった。乳癌データ集合では、56,172個のIPA全体の7が、対応する陰性対照と比べて有意に高い又は低いことが分かった。平均して、NCIパスウェイは患者当たり497個の有意なエンティティを有し、127個中103個のパスウェイが、20以上の患者で変化した少なくとも1つのエンティティを有した。GBMデータ集合では、141,682個のIPA全体の9が、対応する陰性対照と比べて有意に高い又は低いことが分かった。平均して、NCIパスウェイは患者当たり616個の有意なエンティティを有し、127個中110個のパスウェイが、20以上の患者で変化した少なくとも1つのエンティティを有した。別のコントロールとして、本発明者らは、NCIパスウェイの遺伝子と同じように連結される任意の遺伝子から統合活性を得ることができるかという問いを立てた。これを行うため、偽発見率を推定し、それをSPIATarca:2009上記と比較した。多くの遺伝的ネットワークが癌に関係していることが分かっているため、本発明者らは陰性対照集合として模擬デコイパスウェイを使用することを選択した。各NCIパスウェイについて、NCIパスウェイと同じネットワーク構造を使用してゲノム中のランダム遺伝子を互いに連結することにより、デコイパスウェイを作成した。次に本発明者らはPARADIGM及びSPIAを実行して、NCIパスウェイ及びデコイパスウェイの双方のIPAを導き出した。PARADIGMについては、パスウェイサイズで正規化した後に患者間で有意であることが分かったIPAの数により、各パスウェイを順位付けした。SPIAについては、パスウェイは、それらの計算された影響係数に従い順位付けした。本発明者らは、PARADIGMがSPIAと比べて最上位の活性化パスウェイからより多くのデコイパスウェイを除外することを見出した図5。例えば乳癌では、PARADIGMは、上位10に1個、上位30に2個、及び上位50に4個のデコイを順位付けする。それに対してSPIAは、上位10に3個、上位30に12個、及び上位50に22個のデコイを順位付けする。NCIIPAの順位の分布全体は、順位の累積分布をプロットすることにより観察すると、PARADIGMにおいてSPIAより高いP<0.009、K-S検定。実施例X:乳癌及びGBMにおける上位PARADIGMパスウェイ本発明者らは、我々の順列解析により検出された有意なIPAのエンティティ当たりの平均数に従いNCIパスウェイを並べ替え、乳癌表1及びGBM表2における上位15を計算した。上位15のうちのいくつかのパスウェイは、既にそのそれぞれの癌に関係付けられている。乳癌では、SPIA及びPARADIGMの双方が、エストロゲン関連パスウェイ及びErbB2関連パスウェイを検出することができた。最近の主なメタ分析研究においてWirapatiP,SotiriouC,KunkelS,FarmerP,PradervandS,Haibe-KainsB,DesmedtC,IgnatiadisM,SengstagT,SchuetzF,GoldsteinDR,PiccartM,DelorenziM.Meta-analysisofgeneexpressionprofilesinbreastcancer:towardaunifiedunderstandingofbreastcancersubtypingandprognosissignatures.BreastCancerRes.2008;104:R65、Wirapetietal.は、エストロゲン受容体及びErbB2の状態が、乳癌における僅か3つに過ぎない主要な予後判定シグネチャーのうちの2つであることを見出した。PARADIGMはまた、AKT1関連PI3Kシグナル伝達パスウェイを、いくつかの試料において有意なIPAを有する最上位のパスウェイとして特定することができた図6を参照のこと。抗アポトーシスAKT1セリン-スレオニンキナーゼは乳癌に関与することが知られ、ERBB2パスウェイと相互作用するJuX,KatiyarS,WangC,LiuM,JiaoX,LiS,ZhouJ,TurnerJ,LisantiMP,RussellRG,MuellerSC,OjeifoJ,ChenWS,HayN,PestellRG.Akt1governsbreastcancerprogressioninvivo.Proc.Natl.Acad.Sci.U.S.A.2007May;10418:7438-7443。GBMでは、FOXM1及びHIF-1-αの双方の転写因子ネットワークが広範に研究されており、高悪性度の膠芽腫においては悪性度が低い神経膠腫と比べて過剰発現することが示されているLiuM,DaiB,KangS,BanK,HuangF,LangFF,AldapeKD,XieT,PelloskiCE,XieK,SawayaR,HuangS.FoxM1Bisoverexpressedinhumanglioblastomasandcriticallyregulatesthetumorigenicityofgliomacells.CancerRes.2006Apr.;667:3593-3602;SemenzaGL.HIF-1andhumandisease:onehighlyinvolvedfactor.GenesDev.2000Aug.;1416:1983-1991。実施例XI:データ集合の視覚化PARADIGM推論の結果を視覚化するため、本発明者らは、パスウェイにおける各遺伝子を中心として複数のデータ集合を表示するCircleMap視覚化を開発した図7。この表示では、遺伝子の周りに同心円状の輪をプロットすることにより、各遺伝子がコホート間のその全てのデータと関連付けられ、ここで各輪は、一つの種類の計測値又は計算推論に対応する。輪の中の各度数刻みが一つの患者試料に対応し、一方、色が、活性化レベル赤色、不活性化レベル青色、又は未変化レベル白色の活性に対応する。本発明者らはErbB2パスウェイのサブセットについてCircleMapをプロットし、乳癌コホートからのER状態、IPA、発現、及びコピー数データを含めた。遺伝子発現データは、様々な癌の分子サブタイプを定義することに成功している。癌サブタイプは、薬物感受性及び全生存などの種々の臨床転帰と相関することが分かっている。本発明者らは、GBMについて、生の発現データではなくPARADIGMIPAを使用して情報提供性を有するサブタイプを同定することができるかという問いを立てた。IPAを使用することの利点は、IPAが、コピー数、発現、及び遺伝子間の既知の相互作用の要約を提供し、従って意味のある患者サブグループを明らかにするためのよりロバストなシグネチャーを提供し得ることである。本発明者らは、初めに、GBM試料間で少なくとも中程度に繰り返し活性化した全てのIPAを決定し、1,755個のエンティティが229例中少なくとも75例の試料で0.25のIPAを有したことを見出した。それらのエンティティのIPAを全て、活性行列に集めた。次に試料及びエンティティを、非中心化ピアソン相関及び重心連結による階層的クラスタリングを用いてクラスター化した図8。目視検査から、IPAに基づき4つの明らかなサブタイプが明らかとなり、第4のサブタイプは最初の3つと明らかに異なった。第4のクラスターは、明らかなHIF-1-α転写因子ネットワークの下方制御並びにE2F転写因子ネットワークの過剰発現を示す。HIF-1-αは、低酸素条件に対する応答の調節に関与する主転写因子である。対照的に、最初の3つのクラスターのうちの2つは、EGFRシグネチャーが高く、且つGATAインターロイキン転写カスケードを含むMAPキナーゼカスケードが不活性である。興味深いことに、EGFRの突然変異及び増幅が高悪性度神経膠腫並びに膠芽腫と関連付けられているKuanCT,WikstrandCJ,BignerDD.EGFmutantreceptorvIIIasamoleculartargetincancertherapy.Endocr.Relat.Cancer2001Jun.;82:83-96。増幅及び特定の突然変異により、二量体の自己刺激又はリガンド非依存性活性化のいずれかを介して構成的に活性なEGFRが生じ得る。EGFRが構成的に活性化すると、腫瘍形成及び固形腫瘍の進行が促進され得る。EGFRを標的化することが公知の分子ゲフィチニブが、現在、他のEGFR駆動型の癌におけるその効力について研究されている。従って、定性的にはこれらのクラスターは、患者を層別化し得る生物学的に有意義な論題に焦点を合わせているように思われる。これらの観測を定量化するため、本発明者らは、PARADIGMにより同定される異なるGBMサブタイプが異なる生存プロファイルと一致するかという問いを立てた。本発明者らは、4つのクラスターの各々について、初回診断後の月数に対して生存患者の割合をプロットすることにより、カプランマイヤー曲線を計算した。4つのクラスターの各々についてカプランマイヤー生存曲線をプロットし、個別的なIPAシグネチャーと関連付けられる任意のクラスターが生存転帰を予測し得るかどうかを確かめた図9。第4のクラスターは、他のクラスターと大きく異なるP<2.11×105;Cox比例ハザード検定。最初の3つのクラスターの患者の半数は18ヶ月を超えて生存する;クラスター4の患者について生存は有意に増加し、半数が30ヶ月超生存する。加えて、20~40ヶ月の範囲にわたり、クラスター4の患者は他のクラスターの患者の2倍生存する可能性が高い。実施例XII:クラスターについてのカプランマイヤー生存プロット生存分析から、クラスター4の患者が有意により良好な生存プロファイルを有することが明らかとなった。クラスター4は、網膜芽細胞腫抑制因子と共に作用するE2Fの上方制御を有することが分かった。従ってE2Fの上方制御は、クラスター4の患者の腫瘍試料における細胞周期進行の能動的抑制と整合する。加えて、クラスター4はHIF-1-α転写因子の不活性と関連付けられた。第4のクラスターにおけるこの不活性は、腫瘍がより酸素化されていて、より小さい、又はより新しい腫瘍であり得ることを示唆するマーカーとなり得る。従って、PARADIGMIPAは、顕著に異なる生存転帰を有するサブタイプを描出するのに意味のあるプロファイル集合を提供する。比較のため、本発明者らはまた、患者を発現データ又はCNAデータのみを用いたクラスタリングにより患者サブタイプを導き出すことを試みた。これらのデータソースのいずれを用いたクラスタリングからも明らかな群は見出されず、このデータ集合の元のTCGA解析における知見と一致したTCGA:2008図14を参照のこと。これは、遺伝子間の相互作用及びそこから得られる個々の遺伝子発現の組み合わせ出力が、患者転帰としてのかかる複雑な表現型のより優れた予測因子を提供し得ることを示唆している。実施例XIII:卵巣癌の統合ゲノム解析:試料及び臨床データこの報告は、489個の臨床的にアノテートされたステージII~IVのHGS-OvCa及び対応する正常なDNAの解析を対象にしている。患者は、HGS-OvCaと診断された個体の診断時の年齢、病期、腫瘍悪性度、及び手術転帰を反映した。臨床データは2010年8月25日現在で最新のものとした。HGS-OvCa標本は全身的治療の前に外科的に切除され、但し全ての患者が白金剤の投与を受け、及び94がタキサンの投与を受けた。このコホートの無進行生存及び全生存中央値は、既に公表されている治験と同様である11、12。患者の25が無病を維持し、及び45が最終経過観察時に生存していた一方、31は白金ベース療法の完了後6ヶ月以内に進行した。中央値の追跡は30ヶ月であった範囲0~179。TCGA解析用の試料は、70超の腫瘍細胞核を有し、且つ壊死が20未満のものを選択した。独立した部位に複数の分子アッセイを用いる協調的分子解析Coordinatedmolecularanalysesを、表4に掲載するとおりデータはhttp://tcga.cancer.gov/dataportalで入手可能二階層で実施した。階層1のデータ集合は公開されているが、階層2のデータ集合は個人を特定し得る臨床情報又はゲノム情報を含むため、http://tcga.cancer.gov/dataportal/data/access/closed/に記載されるとおり条件を満たす必要がある。実施例XIV:突然変異分析。316個のHGS-OvCa試料から単離したDNA及び各個体に対応する正常試料に対し、エキソームキャプチャー及びシーケンシングを実施した。キャプチャー試薬は、合計約33メガベースの非冗長性配列の約18,500個の遺伝子から約180,000個のエクソンを標的化した。IlluminaGAIIxプラットフォーム236試料ペア又はABISOLiD3プラットフォーム80試料ペアでの超並列シーケンシングにより、試料当たり約14ギガベース合計約9×109塩基が得られた。平均して、腫瘍及び対応正常試料の双方においてコード塩基の76が十分な深さで網羅され、確信的な突然変異検出が可能となった。19,356個の体細胞突然変異腫瘍当たり約61個をアノテートし、表4に分類した。HGS-OvCa病態生理学において重要であり得る突然変異を、aバックグラウンドと比べて有意に増加した頻度で存在する非同義的変異又はスプライス部位変異を調べ、bこの研究の突然変異をCOSMIC及びOMIMの突然変異と比較し、及びcタンパク質機能に対する影響を予測することにより、同定した。2つの異なるアルゴリズムにより、非同義的変異又はスプライス部位変異の数が突然変異分布モデルに基づき予想される数を有意に上回った9個の遺伝子が同定された表5。公表されている結果と一致して13、TP53は316個中303個の試料において変異し283個は自動化された方法による、及び20個は手動での精査後、BRCA1及びBRCA2は、それぞれ9及び8のケースで生殖細胞系列突然変異を有し、双方ともさらに3のケースで体細胞突然変異を示した。6つの他の統計的に繰り返し変異した遺伝子が同定された;RB1、NF1、FAT3、CSMD3、GABRA6、及びCDK12。CDK12はRNAスプライシング調節に関与するもので14、これまでに肺及び大腸の腫瘍に関係付けられた15、16。9個のCDK12突然変異のうちの5個はナンセンス又はインデルのいずれかであったことから、潜在的な機能喪失が示唆され、一方4個のミスセンス突然変異R882L、Y901C、K975E、及びL996Fがそのプロテインキナーゼドメインにクラスター化された。GABRA6及びFAT3は、双方とも有意に変異したと見られたが、HGS-OvCa又は卵管組織では発現しないものと思われたため、従ってこれらの遺伝子の突然変異がHGS-OvCaにおいて有意な役割を果たす可能性は低い。本研究による突然変異をCOSMIC17及びOMIM18データベースの突然変異と比較して、一般的にそれほど変異しないさらなるHGS-OvCa遺伝子を同定した。これにより、BRAFN581S、PIK3CAE545K及びH1047R、KRASG12D、及びNRASQ61Rにおける突然変異を含め、それぞれ477個及び211個のマッチが得られた。これらの突然変異は形質転換活性を呈することが示されているため、本発明者らは、これらの突然変異がまれであるものの、HGS-OvCaにおける重要なドライバーであると考える。本発明者らは、タンパク質ファミリー及び全脊椎動物ゲノムの配列アラインメントからの進化情報と、予測された局所的なタンパク質構造と、選択されたヒトSwissProtタンパク質特徴とを組み合わせることにより、既知の癌遺伝子及び腫瘍抑制因子における突然変異に関する訓練後、CHASMを使用して推定ドライバー突然変異を同定した19、20。CHASMにより、腫瘍形成性を有すると予測される122個のミスセンス突然変異が同定された。突然変異により駆動されるタンパク質機能の変化を、タンパク質ファミリー配列アラインメント及び既知の又は相同性に基づく三次元タンパク質構造における残基配置をMutationAssessorを使用して比較することによって、全ての確認済みの体細胞ミスセンス突然変異についての進化情報から推測した。ミスセンス突然変異の27パーセントがタンパク質機能に影響を与えると予測された。実施例XV:コピー数解析489個のHGS-OvCaゲノムに存在する体細胞性コピー数変化SCNAを同定し、図37Aにおいて多形膠芽腫データと比較した。SCNAは、広範な染色体領域に影響を及ぼす領域性の異常と、より小さい限局的な異常とに分けた。領域性異常の統計的解析により8個の繰り返し起こる獲得と22個の喪失が同定され、その全てが既に報告されている22図37B。獲得のうち5個及び喪失のうち18個は、50超の腫瘍で起こった。GISTICを使用して、繰り返し起こる限局性SCNAを同定した。これにより、8個以下の遺伝子をコードする26領域を含む63領域の限局性増幅図37Cが得られた。最も一般的な限局性増幅はCCNE1、MYC、及びMECOMをコードし図37C、各々20超の腫瘍で高度に増幅された。HGS-OvCaにおける新しい極めて限局化された増幅ピークは、活性化Cキナーゼ受容体ZMYND8;p53標的遺伝子IRF2BP2;DNA結合タンパク質阻害薬ID4;胚発生遺伝子PAX8;及びテロメラーゼ触媒サブユニットTERTをコードした。3つのデータソース:http://www.ingenuity.com/、http://clinicaltrials.gov及びhttp://www.drugbank.caを使用して、増幅された過剰発現遺伝子の潜在的治療阻害因子を同定した。この調査により、少なくとも10のケースで増幅されたMECOM、MAPK1、CCNE1及びKRASを含め、治療標的である22個の遺伝子が同定された。GISTICにより50個の限局的な欠失も同定された。腫瘍の少なくとも2において、既知の腫瘍抑制遺伝子PTEN、RB1、及びNF1がホモ接合性欠失領域にあった。重要なことには、RB1及びNF1はまた有意に変異した遺伝子のなかにもあった。一つの欠失は、5個の非同義的突然変異と2個のフレームシフト突然変異とを有する必須細胞周期制御遺伝子のCREBBPを含め、3つの遺伝子のみを含んだ。実施例XVI:mRNA及びmiRNA発現並びにDNAメチル化解析。3つの異なるプラットフォームAgilent、AffymetrixHuEx、AffymetrixU133Aからの11,864個の遺伝子についての発現計測値を、サブタイプの同定及び転帰の予測用に組み合わせた。個々のプラットフォーム計測値は、限定的な、しかし統計的に有意なバッチ効果を受けたが、一方、組み合わせたデータ集合はそれを受けなかった。組み合わせたデータ集合の解析から約1,500個の本質的に変化し易い遺伝子が同定され、それらをNMFコンセンサスクラスタリングに使用した。この解析により4つのクラスターが得られた図38a。同じ解析手法をTothilletal.の公的に利用可能なデータ集合に適用し、同様に4つのクラスターを得た。TothillのクラスターとTCGAクラスターとの比較は、明らかな相関を示した。従って本発明者らは、HGS-OvCaには少なくとも4つのロバストな発現サブタイプが存在すると結論付ける。本発明者らは、その4つのHGS-OvCaサブタイプを、クラスターに含まれる遺伝子及び先行する観測に基づき25、免疫反応型、分化型、増殖型及び間葉型と名付けた。T細胞ケモカインリガンドのCXCL11及びCXCL10、並びに受容体のCXCR3が、免疫反応型サブタイプを特徴付けた。HMGA2及びSOX11などの転写因子の高発現、卵巣腫瘍マーカーMUC1、MUC16の低発現並びにMCM2及びPCNAなどの増殖マーカーの高発現が、増殖型サブタイプを定義した。分化型サブタイプはMUC16及びMUC1の高発現並びに分泌性卵管マーカーSLPIの発現と関連付けられ、より成熟した発達段階が示唆された。HOX遺伝子の高発現並びに筋線維芽細胞FAP及び微小血管周皮細胞ANGPTL2、ANGPTL1などの間質成分の増加を示唆するマーカーが、間葉型サブタイプを特徴付けた。DNAメチル化の上昇及び腫瘍発現の低下により、168個の遺伝子が、卵管対照と比較してHGS-OvCaにおいて後成的にサイレンシングされたものとして関係付けられた26。DNAメチル化は全試料にわたり遺伝子発現の低下と相関した。AMT、CCL21及びSPARCL1は、大多数の腫瘍でプロモーター過剰メチル化を示したため、注目に値した。奇妙にも、卵巣癌で増幅及び過剰発現することが以前に報告されているRAB25もまた、一部の腫瘍で後成的にサイレンシングされるように見えた。BRCA1プロモーターは、既に報告されているとおり、489個中56個11.5の腫瘍で過剰メチル化及びサイレンシングされた。腫瘍間での様々なDNAメチル化のコンセンサスクラスタリングにより、年齢、BRCA不活性化イベント、及び生存の差と有意に関連付けられた4つのサブタイプが同定された。しかしながら、これらのクラスターは中程度の安定性しか示さなかった。TCGAデータ集合の転写サブタイプについて、生存期間に有意な差はなかった。増殖型群がMYC増幅及びRB1欠失の速度低下を示したのに対し、免疫反応型サブタイプは3q26.2MECOM増幅頻度の増加を示した。DNAメチル化クラスターと遺伝子発現サブタイプとの間の中程度の、しかし有意な重複が指摘されたp<2.21016、カイ二乗検定、調整ランド指数=0.07。215個の試料からの統合した発現データ集合を使用して、全生存を予測し得る193個の遺伝子の転写シグネチャーを定義した。一変量コックス回帰分析の後、108個の遺伝子が不良な生存率と相関し、85個が良好な生存率と相関した0.01のp値カットオフ。独立した255個のTCGA試料の集合並びに3個の独立した発現データ集合において予測力を検証した25、29、30。検証試料の各々に、その発現プロファイルと予後判定遺伝子シグネチャーとの間の類似性を反映して、予後判定遺伝子スコアを割り当てた31図38c。このシグネチャーのカプランマイヤー生存分析から、全ての検証データ集合で生存率との統計的に有意な関連性が示された図38d。miRNA発現データのNMFコンセンサスクラスタリングにより、3つのサブタイプが同定された。興味深いことに、miRNAサブタイプ1はmRNA増殖型サブタイプと重複し、及びmiRNAサブタイプ2はmRNA間葉性サブタイプ図38dと重複した。生存期間はiRNAサブタイプ間で有意に異なり、miRNAサブタイプ1腫瘍の患者の生存が有意に長かった図38e。実施例XVII:疾患に影響するパスウェイいくつかの分析により、316例の完全に解析されたケースからデータが統合され、HGS-OvCaに寄与する生物学が特定された。既知の癌関連パスウェイが1つ又は複数の突然変異、コピー数変化、又は遺伝子発現変化を含む頻度の解析から、RB1パスウェイ及びPI3K/RASパスウェイがそれぞれ67及び45のケースで調節解除されたことが示された図39A。HotNet33を使用して大規模タンパク質間相互作用ネットワーク32において変化したサブネットワークを調べることで、HGS-OvCa試料の23で変化したNotchシグナル伝達パスウェイを含め、いくつかの既知のパスウェイが同定された図39B。公表された研究では、変異した又はメチル化したBRCA1又は変異したBRCA2を含む細胞は相同組換えHR欠陥を有し、PARP阻害薬に高度な応答性を有することが示されている35~37。図39Cは、HGS-OvCaの20がBRCA1/2に生殖細胞系列突然変異又は体細胞突然変異を有し、11がDNA過剰メチル化によってBRCA1発現を失っており、BRCA1の後成的サイレンシングがBRCA1/2突然変異と相互排他的であるP=4.4×104、フィッシャーの正確検定ことを示している。BRCA状態の一変量生存分析から図39C、BRCA変異ケースはBRCA野生型ケースと比べて全生存OSがより良好であることが示された。興味深いことに、後成的にサイレンシングしたBRCA1ケースは、BRCA1/2WTHGS-OvCaと同様の生存を示したOS中央値41.5ヶ月対41.9ヶ月、P=0.69、ログランク検定。これは、BRCA1が相互排他的なゲノム機構及びエピゲノム機構により不活性化されること、及び患者生存が不活性化機構に依存することを示唆している。この研究で見出された、細胞をPARP阻害薬に対して感受性にし得る他のHR遺伝子におけるゲノム変化としては、EMSYの増幅又は突然変異8、PTENの限局的な欠失又は突然変異7;RAD51Cの過剰メチル化3、ATM/ATRの突然変異2、及びファンコニー貧血遺伝子の突然変異5が含まれる。概して、HR欠陥はHGS-OvCaの約半数に存在し得るもので、腫瘍これらのHR関連異常を標的化するPARP阻害薬の臨床試験の理論的根拠を提供する。BRCA不活性化イベントの完全集合を全ての繰り返し変化したコピー数ピークと比較することにより、BRCA不活性化を有するケースにおいて予想外に低いCCNE1増幅頻度が明らかとなったBRCA野生型ケースの26と比べて、BRCAが変化したケースの8がCCNE1増幅を有した、FDR調整P=0.0048。既に報告されているとおり39、全生存はCCNE1増幅を有する患者について、他の全てのケースと比較して短くなる傾向があったP=0.072、ログランク検定。しかしながら、BRCA野生型ケースのみを調べたとき、CCNE1が増幅したケースについての生存上の不利は見られなかったためP=0.24、ログランク検定、既に報告されているCCNE1の生存の差は、BRCAが変異したケースのより良好な生存により説明され得ることが示唆される。最後に、確率的グラフィカルモデルPARADIGM40により、NCIパスウェイ相互作用データベースNCIPathwayInteractionDatabaseにおける変化したパスウェイを調べ、FOXM1転写因子ネットワーク図39Dを87のケースで有意に変化したものとして同定した。FOXM1及びその増殖関連標的遺伝子;AURB、CCNB1、BIRC5、CDC25、及びPLK1は、一貫して過剰発現したが、転写調節の指標となるDNAコピー数変化によっては変化しなかった。TP53はDNA損傷後にFOXM1を抑制することから42、HGS-OvCaにおける高率のTP53突然変異がFOXM1過剰発現に寄与することが示唆される。他のデータ集合において、FOXM1パスウェイは隣接上皮組織と比べて腫瘍で有意に活性化され、HGS-OvCaと関連付けられる。実施例XVIII:卵巣漿液性癌において高頻度で変化したパスウェイコピー数及び遺伝子発現の双方の統合解析によって有意に変化したパスウェイを同定するため、本発明者らはPARADIGMを適用した。この計算モデルはコピー数変化、遺伝子発現データ、及びパスウェイ構造を組み込み、パスウェイデータベースに存在する遺伝子、複合体、及び遺伝的プロセス毎に統合パスウェイ活性IPAを生成する。本発明者らは用語エンティティを、それが遺伝子か、複合体か、又は小分子かを問わず、パスウェイの任意の分子を指して使用する。エンティティのIPAは、最終活性のみを指す。遺伝子について、IPAはタンパク質の活性状態の推論された活性のみを指し、これは、パスウェイにおける他の遺伝子のコピー数、遺伝子発現、及びシグナル伝達から推論される。本発明者らはPARADIGMを卵巣試料に適用し、NationalCancerInstitutesのパスウェイ相互作用データベースNCI-PIDに含まれるパスウェイに存在する多くの異なる遺伝子及び過程の変化を見出した。本発明者らは、1000回のランダムシミュレーションを用いて推論された変化の有意性を評価し、それらのシミュレーションでは、同じ構造を有するパスウェイが使用されるが、パスウェイの種々の点に任意の遺伝子を割り当てた。換言すれば、所与のパスウェイのあるランダムシミュレーションにおいて相互作用の集合は固定されたままであり、従って任意の遺伝子集合がそれらのパスウェイ相互作用によって互いに連結された。全ての試料のIPAの有意性は同じヌル分布に対して評価して、各試料におけるエンティティ毎の有意水準を得た。IPA及びそれらが有意である試料の割合及び標準偏差が少なくとも0.1のIPAを、図28にヒートマップとして表示する。表3は、PARADIGMにより見出された順列化試料に関して少なくとも3標準偏差だけ変化するパスウェイを示す。FOXM1転写因子ネットワークは、検証した全てのパスウェイのなかで最多数の試料において変化した-試料間で平均化したとき67のエンティティが活性の変化を有した。それに対して、卵巣コホートにおいて次に最も高いレベルの活性の変化を有したパスウェイとしては、PLK1シグナル伝達イベント27、オーロラBシグナル伝達24、及びトロンボキサンA2受容体シグナル伝達20が含まれた。従って、NCI-PIDのパスウェイのなかでは、卵巣試料に関して、FOXM1ネットワークが他のパスウェイと比べて有意により大きく変化する活性を含む。FOXM1転写因子ネットワークが、最も高い割合の患者試料において、正常対照と比較して腫瘍試料で示差的に変化することが分かった図29。FOXM1は、3つの既知の支配的なスプライシング型を有する多機能性転写因子であり、スプライシング型の各々は、細胞増殖及びDNA修復において様々な役割を有する個別的な遺伝子サブセットを調節する。FOXM1cアイソフォームは、AUKB、PLK1、CDC25、及びBIRC5を含めた、細胞増殖において既知の役割を有するいくつかの標的を直接調節する。他方で、FOXM1bアイソフォームは、DNA修復遺伝子BRCA2及びXRCC1を含む完全に異なる遺伝子サブセットを調節する。ATMの間接的な制御下にあるCHEK2は、FOXM1発現レベルを直接調節する。本発明者らは、FOXM1転写因子それ自体のIPAが、他の転写因子のIPAと比べてより高度に変化するかという問いを立てた。FOXM1の活性レベルをNCI-PIDにおける他の203個の転写因子の全てと比較した。NCI集合の他の転写因子と比較しても、FOXM1転写因子は有意に高いレベルの活性を有したp<0.0001;K-S検定ことから、それが重要なシグネチャーであり得ることがさらに示唆される図30。FOXM1はまた、上皮起源の多くの異なる正常組織においても発現するため、本発明者らは、PARADIGMにより同定されるシグネチャーが、他の組織では正常と見なされ得る上皮シグネチャーに起因するかという問いを立てた。これに答えるため、卵管上皮及び卵巣腫瘍組織がマイクロダイセクションされたGEOからの独立したデータ集合をダウンロードしGSE10971、遺伝子発現をアッセイした。FOXM1のレベルが正常試料と比較して腫瘍試料において有意に高いことが見出されたことから、実際に癌性組織においてはFOXM1調節が、正常上皮組織に見られるものを超えて亢進することが示唆される図31。TCGA卵巣コホート全体が高悪性度漿液性腫瘍に由来する試料を含んだため、本発明者らは、FOXM1シグネチャーが高悪性度漿液性に特異的であるのかという問いを立てた。低悪性度及び高悪性度の双方の漿液性腫瘍の転写プロファイルが決定されているEtemadmoghadametal.2009のデータ集合から、FOXM1及びその標的のいくつかの対数発現を得た。この独立したデータから、FOXM1及びその標的のいくつかが、低悪性度卵巣癌と比べて漿液性卵巣癌において有意に上方制御されることが確認された図32。FOXM1転写因子ネットワークの25個の遺伝子が、高悪性度疾患において発現がより高い遺伝子を有意な割合で含んだかどうかを決定するため、本発明者らはEtemadmoghadamのデータを使用してスチューデントのt検定を実施した。ゲノム中723個の遺伝子5.4が、有意水準0.05で低悪性度癌と比べて高悪性度癌において有意に上方制御されることが分かったベンジャミニホッホバーグBenjamini-Hochbergの方法を用いて複数の試験に対して補正。FOXM1ネットワークは、その遺伝子のうち13個52が示差的に調節されることが分かった。これは、超幾何検定に基づけば有意な割合であるP<3.81012。従ってFOXM1ネットワーク遺伝子の高発現は、ゲノムにおける典型的な遺伝子の発現と比較したとき、高悪性度疾患と確かに特異的に関連するように思われる。乳癌及び肺癌を含め、多くの異なる癌におけるFOXM1の役割が十分に裏付けられているが、卵巣癌におけるその役割はいまだ研究がなされていない。FOXM1は、3つの既知のスプライシング型を有する多機能性転写因子であり、スプライシング型の各々は、細胞増殖及びDNA修復において様々な役割を有する個別的な遺伝子サブセットを調節する。この分析に関連するFOXM1の相互作用ネットワークの抜粋を図27に示す。FOXM1aアイソフォームは、AUKB、PLK1、CDC25、及びBIRC5を含めた、細胞増殖において既知の役割を有するいくつかの標的を直接調節する。対照的にFOXM1bアイソフォームは、DNA修復遺伝子BRCA2及びXRCC1を含む完全に異なる遺伝子サブセットを調節する。ATMの間接的な制御下にあるCHEK2は、FOXM1の発現レベルを直接調節する。卵巣患者のほとんどにおけるFOXM1の発現増加に加え、小さいサブセットはまた、CBSにより検出されるコピー数増幅の増加も有する19が、計測したゲノムにおける全遺伝子の上位5分位内のコピー数増加を伴う。従って、FOXM1の選択的スプライシング調節はDNA修復と細胞増殖との間の制御スイッチに関与し得る。しかしながら、個別のアイソフォーム活性を区別することは、アイソフォームを区別するエクソン構造及びエクソンアレイプローブの位置によって困難なため、現時点ではこの主張を裏付けるデータは不十分である。これらの試料のmRNAの将来のハイスループットシーケンシングが、示差的なFOXM1アイソフォームレベルの決定を促進し得る。PARADIGMがこの転写因子に集中している最高レベルの活性変化を検出したという観測は、FOXM1が細胞において重要な調節ポイントにあることを示唆している。実施例XIX:データ集合及びパスウェイ相互作用コピー数及び発現データの双方をPARADIGM推論に組み入れた。8個の正常組織対照の集合が発現データにおける解析に利用可能であったため、患者の遺伝子値の各々を、正常な卵管対照で観察される遺伝子の中央値レベルを減じることにより正規化した。コピー数データは、腫瘍で検出された遺伝子レベルと、それに対する血液正常レベルとの間のコピー数の差を反映するように正規化した。PARADIGMへの入力のため、発現データは、サブタイプ解析に使用したものと同じ統合データ集合から取り、コピー数は、MSKCCAgilent1Mコピー数データのセグメント化コールから取った。131個のパスウェイ、11,563個の相互作用、及び7,204個のエンティティを含むNCI-PIDから一群のパスウェイを入手した。エンティティは、PARADIGMのグラフィカルモデルにおいてノードとして表現される分子、複合体、小分子、又は抽象概念である。抽象概念は、一般的な細胞プロセスアポトーシス又は光の吸収など、及びシグナルトランスデューサーのRASファミリーなどの機能的活性を共有する遺伝子のファミリーに対応する。本発明者らは、タンパク質間相互作用、転写調節相互作用、リン酸化及びユビキチン化などのタンパク質修飾相互作用を含む相互作用を収集した。実施例XX:パスウェイ文脈における統合分子活性の推論本発明者らは、コピー数、遺伝子発現、及び各エンティティのパスウェイ文脈を反映する統合パスウェイ活性IPAを割り当てるPARADIGMを使用した。データの遺伝子特異的及び患者特異的断面の順列を使用して、IPAの有意性を評価した。ゲノム中の各遺伝子についての値の遺伝子発現とコピー数とのペアを無作為に選択することにより、1000個のヌル患者のデータを作成した。PARADIGMIPAの有意性を評価するため、パスウェイ構造を維持しながらパスウェイにランダム遺伝子を割り当てることにより、ヌル分布を作成した。実施例XXI:FOXM1パスウェイの同定ランダムシミュレーションの間は、FOXM1ネットワーク内の全ての遺伝子を使用して統計的有意性を評価したが、FOXM1パスウェイの視覚化を可能にするため、図29により有意に変化したIPAを有するFOXM1に直接連結されたエンティティを、図27に含めるため選択した。これらのうち、DNA修復及び細胞周期制御において役割を有する遺伝子であって、FOXM1との相互作用について文献の裏付けがあると認められたものを表示した。元のNCI-PIDパスウェイに存在しなかったBRCC複合体メンバーを、NCI-PIDによればFOXM1の標的であるBRCA2と共にプロットに含めた。上流DNA修復標的を、他のNCIパスウェイにおけるCHEK2の上流直接因子を見つけることにより同定した例えば、PLK3シグナル伝達パスウェイにATMからの間接的リンクが見出された。実施例XXII:クラスタリング活性及び非活性の確率の変化を直接的に表す推論活性の使用は、それにより様々な種類のエンティティをまとめて一つのヒートマップにクラスター化することを可能にする。PARADIGM推論の結果を包括的に視覚化するため、EisenCluster3.0を使用して特徴フィルタリング及びクラスタリングを実施した。0.1の標準偏差フィルタリングにより7204個中1598個のパスウェイエンティティが残り、エンティティ及びサンプルの双方に対して平均連結法非中心化相関階層的クラスターを実施した。実施例XXIII:細胞株は多くの重要な腫瘍サブタイプ及び特徴をモデル化する臨床的に関連する分子的応答予測因子の同定に対する細胞株の有用性は、腫瘍における応答を決定する多様な分子機構が細胞株内でどの程度機能するかに依存する。本発明者らは、転写物及びゲノムコピー数の双方のレベルでの細胞株モデルと原発腫瘍との間の類似性に関して以前に報告しており9、ここでこの比較を、より高分解能のプラットフォーム及び解析技法を用いて改良する。具体的には、本発明者らは遺伝子発現プロファイルの階層的コンセンサスクラスタリングHCCを使用して、50個の乳癌細胞株及び5個の非悪性乳房細胞株を3つの転写サブタイプ:ルミナル、基底細胞及び新しく記載されるクローディン低に分類した図14A。これらのサブタイプは先述したものの改良版であり、基底細胞型及びクローディン低caludin-low型は、それぞれ先に指定した基底細胞A及び基底細胞Bサブタイプに対応する、表7。改良版高分解能SNPコピー数解析図14Bにより、細胞株パネルが、8q24MYC、11q13CCND1、17q12ERBB2、20q13STK15/AURKAの繰り返し生じる増幅領域、及び原発腫瘍に見られる9p21CDKN2Aのホモ接合性欠失をモデル化することが確認される。トラスツズマブ及びラパチニブ療法により決定されるとおりのERBB2腫瘍サブタイプの臨床的関連性を所与として、本発明者らは、ERBB2AMPと命名される特別なサブタイプとして、ERBB2のDNA増幅を伴う細胞株を調べた。概して、ルミナル型、基底細胞型、クローディン低型及びERBB2AMP型細胞株の本発明者らの同定は、臨床生物学と一致している。実施例XIX:細胞株はほとんどの治療化合物に示差的感受性を示す本発明者らは、77種の治療化合物に対する我々の細胞株パネルの感受性を調べた。細胞成長アッセイを用い、定量的エンドポイントは9つの濃度の各薬剤に3日間連続的に曝露した後に計測した。試験した抗癌化合物には、従来の細胞傷害剤例えば、タキサン、プラチノール、アントラサイクリンanthracylineと標的化薬剤例えば、SERM及びキナーゼ阻害薬との混合物が含まれた。多くの場合、いくつかの薬剤が同じタンパク質又は分子的作用機構を標的化した。本発明者らは、各化合物について応答の定量的尺度を、50だけ成長を阻害するのに必要な濃度GI50と称するとして決定した。基礎となる成長データは高品質であるものの、50阻害が達成されなかった場合、GI50は試験した最高濃度に設定した。GI50値は全ての化合物について表8に提供する。3つの化合物についてはPS1145、セツキシマブ及びバイカレイン、細胞株応答の変動が最小限であったためさらなる解析から除外した。SigmaAKT1-2阻害薬に対する応答の変動を、関連する転写サブタイプと共に示す代表的なウォーターフォールプロットを図10Aに示す。この化合物に対する感受性はルミナル型及びERBB2AMP型で最も高く、基底細胞型及びクローディン低型乳癌細胞株でより低かった。全ての化合物についての細胞株間のGI50値の分布を示すウォーターフォールプロットは、付録に載せる。本発明者らは、3回又は4回の反復による229個の化合物/細胞株の組み合わせについてのGI50値の絶対偏差中央値を計算することにより、データ集合全体の再現性を確立した。これらの反復にわたる平均偏差中央値は0.15であった図15。本発明者らは、GI50値の集合間のペアワイズなピアソン相関を計算することにより、8種の化合物に対する応答の一致を評価した図15B.同様の作用機構を有する薬物のペアに対する感受性は高度に相関したことから、同様の作用機構が示唆される。実施例XX:多くの化合物は細胞株のサブセットで選択的に有効である本研究の主な前提は、細胞株における予測性のある分子的特徴がヒト腫瘍に反映される場合に、前臨床細胞株分析で認められる応答と分子サブタイプとの間の関連性が診療室で繰り返し起こり得ることである。本発明者らは、ノンパラメトリックANOVAを用いて転写サブタイプ及びゲノミクスサブタイプ間のGI50値を比較することにより、応答-サブタイプの関連性を確立した。概して、試験した74個中33個の化合物が、転写サブタイプ特異的な応答を示したFDRp<0.2、表7及び表9。図10Cは、ルミナル、基底細胞、クローディン低及びERBB2AMPサブタイプの1つ又は複数と有意な関連性を有する34個の薬剤の階層的クラスタリングを示す。サブタイプと最も強い関連性を有する11個の薬剤は、受容体チロシンキナーゼシグナル伝達及びヒストンデアセチラーゼの阻害薬であり、ルミナル型及び/又はERBB2AMP型細胞株で最も高い効力を有した。次に最も高い3つのサブタイプ特異的薬剤-エトポシド、シスプラチン、及びドセタキセル-は、臨床的に観察されるとおり、基底細胞型及び/又はクローディン低型細胞株で選択的活性を示す。イクサベピロン、GSK461364ポロキナーゼ阻害薬及びGSK1070916オーロラキナーゼ阻害薬を含む分裂装置を標的化する薬剤もまた、基底細胞型及びクローディン低型細胞株に対する活性がより高かった。全てEGFR及び/又はERBB2を標的化するAG1478、BIBW2992及びゲフィチニブは、ERBB2増幅と正の関連性を有した。HSP90の阻害薬であるゲルダナマイシンもまた、ERBB2増幅と正の関連性を示した。興味深いことに、VX-680オーロラキナーゼ阻害薬及びCGC-11144ポリアミン類似体は、双方ともERBB2増幅と負の関連性を示したことから、これらはERBB2AMP型腫瘍には比較的不適切な療法であることが示される。本発明者らは、応答と、繰り返し起こる限局的な高レベルコピー数異常CNA;標本のt検定、FDRp<0.2、表10との間に、7つの関連性6つの固有の化合物を同定した。図10Dは以下のことを示す:a9p21CDKN2A及びCDKN2Bにおけるホモ接合性欠失が、ビノレルビン、イクサベピロン及びファスカプリシンfascalypsinに対する応答と関連性を有した。ファスカプリシンfascalypsinはCDK4を阻害した。この特異性は、CDK4の阻害におけるCDKN2Aのp16INK4A産物の役割と一致する20。b20q13AURKAをコードするにおける増幅は、AURKB及びAURKCを標的化するGSK1070916及びVX-680に対する感受性よりむしろ、耐性と関連性を有した23。これは、AURKAの増幅がAURKB及びAURKC阻害薬の迂回機構を提供することを示唆する。c11q13CCND1における増幅は、カルボプラチン及びAURKB/C阻害薬GSK1070916に対する感受性と関連性を有した。実施例XXI:サブタイプ特異性は成長速度効果を抑制する概して、本発明者らは、ルミナルサブタイプ細胞株が基底細胞型又はクローディン低型細胞より成長が遅くクラスカルワリス検定p=0.006、図16A及び表7、倍加時間の範囲が広かった18~300時間ことを見出した。これにより、最も感受性が高い細胞株は最も速く成長する細胞株であった可能性が生じた。そうであるならば、次に観測されたサブタイプとの関連性は共変量との関連性を表すことになり得る。この仮説を、共分散分析ANCOVAを用いてサブタイプ及び倍加時間の効果を同時に評価することにより検証したところ、33個中22個のサブタイプ特異的化合物が、倍加時間よりサブタイプとより良好な関連性を有したことが分かったp値の平均対数比=0.92、標準偏差1.11。これは、サブタイプの構成員が成長速度より良好な応答予測因子であるという着想を支持する。さらに、33個中15個のサブタイプ特異的化合物が、より緩徐に成長するルミナル型細胞株において有効性がより高かった表7。一つの薬剤、5-フルオロウラシル5-florouracilは、サブタイプの検証単独では有意でなかったが、ANCOVAモデルではクラス及び倍加時間の双方について強い有意性を示した。5-フルオロウラシル5-florouracilに対する応答は、ルミナル型細胞株及び基底細胞型細胞株の双方で倍加時間が増加するに従い低下した図16B。本発明者らは、ほとんどの場合に、3日間の成長阻害アッセイが、成長速度により強い影響を受けない分子シグネチャー特異的応答を検出していると結論付ける。実施例XXII:コピー数及び転写計測値の統合は、サブタイプ特異的応答のパスウェイを同定する本発明者らは、ネットワーク解析ツールPARADIGM24を使用して、細胞株パネルのサブタイプ間におけるパスウェイ活性の違いを同定した。この解析は、キュレーションされたパスウェイが部分的に重複しているという事実によって複雑化する。例えばEGFR、PI3キナーゼ及びMEKは、実際にはそれらが単一のより大きいパスウェイの構成要素であるときにも、別個のパスウェイとしてキュレーションされることが多い。この問題に対処するため、PARADIGMは約1400個のキュレーションされたシグナル伝達パスウェイ、転写パスウェイ及び代謝パスウェイを単一の重畳パスウェイSuperPathwayにまとめ合わせ、かかる冗長性を排除する。特定の細胞株についてのコピー数データ及び遺伝子発現データの双方を用いて、PARADIGMはパスウェイ相互作用を使用することにより、遺伝子毎、複合体毎、及び細胞プロセス毎に統合パスウェイレベルIPLを推論する。本発明者らは、PARADIGMIPLを用いて、細胞株を原発性乳房腫瘍とそれらのパスウェイ活性化により比較した。細胞株-腫瘍比較データは、癌ゲノムアトラスTheCancerGenomeAtlas:TCGAプロジェクトhttp://cancergenome.nih.govにより生成されたデータを使用して行われた。図11は、階層的クラスタリング後の腫瘍及び細胞株の各々についてのパスウェイ活性を示す。各サブタイプについて上位5つのパスウェイ特徴を表11に掲載する。概して、腫瘍及び細胞株サブタイプは同様のパスウェイ活性を示し、及び調節解除されたパスウェイは、元のサブタイプと比べて転写サブタイプとより良好な関連性を有した図13。しかしながら、クローディン低型細胞株サブタイプと関連付けられるパスウェイは腫瘍では十分に表現されない-おそらくは、クローディン低サブタイプが細胞株コレクションに過剰に見られ、且つルミナルAサブタイプが欠けているためである図12。実施例XXIII:サブタイプ特異的パスウェイマーカーの同定本発明者らは、サブタイプ間の違いの根底に固有のパスウェイ活性があるかという問いを立てた。そのため、本発明者らは、あるサブタイプの細胞株で他と比較して示差的に上方制御又は下方制御される遺伝子活性を含むSuperPathwayのサブネットワークを同定した。基底細胞型細胞株とコレクション内の他の全てとの間のパスウェイ活性の比較により、941本のエッジにより連結された965個のノードから構成されるネットワークを同定した。ここで、ノードはタンパク質、タンパク質複合体、又は細胞プロセスを表し、エッジはタンパク質リン酸化などの、それらの要素間の相互作用を表す図18~図22を参照のこと。図35Aは、増殖、血管新生、及び腫瘍形成に関連するMYC/MAXサブネットワークの上方制御;及び細胞周期、接着、浸潤、及びマクロファージ活性化を制御するERK1/2サブネットワークの上方制御を示す。FOXM1及びDNA損傷サブネットワークもまた、基底細胞型細胞株で顕著に上方制御された。クローディン低サブタイプを他の全てと比較することにより、基底細胞型細胞株におけるものと同じサブネットワークの多くの上方制御が示され、但し、基底細胞と比較したときのクローディン低型細胞株におけるβ-カテニンCTNNB1ネットワークの上方制御を含めて例外はあった図35B。β-カテニンは腫瘍発生に関与しているとされ、不良な予後と関連付けられる。ルミナル型細胞株を他の全てと比較することにより、黒色腫における腫瘍形成能を阻害するATF2ネットワークの下方制御、及びER調節性遺伝子の転写を制御し、且つ良好な予後のルミナル型乳癌に関与するとされるFOXA1/FOXA2ネットワークの上方制御が示された図35C。ERBB2AMP型細胞株を他の全てと比較することにより、ルミナル型細胞と共通する多くのネットワーク特徴が示された-ほとんどのERBB2AMP型細胞はまたルミナル型細胞として分類されるため、これは意外ではない。しかしながら図35Dは、ERBB2AMP型細胞株におけるRPS6KBP1を中心とする下方制御を示す。IPLを使用した細胞株間における示差的な薬物応答の比較分析から、応答機構に関する情報を提供するパスウェイ活性が明らかとなった。例えば、基底細胞型細胞株は、DNA損傷剤のシスプラチンに対して選択的な感受性を有し、また、シスプラチンに対する応答に関連する中心的存在であるATM、CHEK1及びBRCA1を含むDNA損傷応答サブネットワークの上方制御も示した34図36A。同様に、ERBB2AMP型細胞株はHSP90の阻害薬であるゲルダナマイシンに対して感受性を有し、またERBB2-HSP90サブネットワークにおける上方制御を示した図36B。この観測はゲルダナマイシンの作用機構と一致する:ゲルダナマイシンはERBB2に結合して、その分解degredationをもたらす。本発明者らは、ERBB2AMP型細胞株がオーロラキナーゼ阻害薬VX-680に耐性を示したこと図36C、上、さらに、この化合物に対する感受性が20q13AURKAにおける増幅と関連性を有しなかったことを見出した。これにより、この耐性が、FOXM1によりAURKBと同時制御されるCCNB1によって媒介され得る可能性が生じる。4つのサブタイプのうちERBB2AMPは、CCNB1の実質的な下方制御を示す唯一のものである図36C及び図22。この提案される機構は、原発腫瘍においてCCNB1遺伝子発現がAURKB遺伝子発現と有意に相関するという観察により裏付けられる。実施例XXIV:細胞成長阻害アッセイ及び成長速度本発明者らは、我々の55個の乳癌細胞株パネルにおいて77種の化合物の効力を評価した。このアッセイは既に記載されるとおり実施したKuo,W.L.etal.AsystemsanalysisofthechemosensitivityofbreastcancercellstothepolyamineanaloguePG-11047.BMCMed7,77,doi:1741-7015-7-77[pii]10.1186/1741-7015-7-772009。簡潔に言えば、細胞を1:5段階希釈の一組の9用量の各化合物で72時間処理した。CellTiterGloアッセイを用いて細胞生存度を決定した。未処理のウェルについての0hに対する72hの比から、倍加時間DTを推定した。本発明者らは、非線形最小二乗法を使用して、以下のパラメータでデータにゴンペルツ曲線をフィットした:上側及び下側漸近線、傾き及び変曲点。フィットした曲線は、NCI/NIHDTPHumanTumorCellLineScreenProcessにより記載され、且つ以前に記載された方法を用いて、GI曲線に変換したScreeningServices-NCI-60DTPHumanTumorCellLineScreen.http://dtp.nci.nih.gov/branches/btb/ivclsp.html.;Monks,A.etal.Feasibilityofahigh-fluxanticancerdrugscreenusingadiversepanelofculturedhumantumorcelllines.JNatlCancerInst83,757-7661991。本発明者らは、50だけ成長を阻害するのに必要な化合物濃度GI50、成長を完全に阻害するのに必要な濃度全成長阻害、TGI及び50だけ集団を減少させるのに必要な濃度50致死濃度、LC50を含め、様々な応答尺度を評価した。基礎をなす成長データが高品質でありながらエンドポイント応答GI50、TGI、LC50に達しなかった場合、値を試験した最高濃度に設定した。GI50は最初に達する閾値に相当し、従って最も正確な計測値集合を含む。薬物応答データを、以下の基準を満たすようフィルタリングした:19つのトリプリケートのデータ点間における標準偏差中央値<0.20;2特定の細胞株について、DT±2SDのDT中央値;3フィットした曲線の傾き>0.25;4明らかな応答がないデータ集合について、最高濃度での成長阻害<50。薬物プレートの約80が全てのフィルタリング要件を通った。本発明者らは、ロバストな形の標準偏差である絶対偏差中央値MADを使用して、我々のGI50の反復計測の信頼性を評価した。曲線の当てはめ及びフィルタリングは、特別に書いたRパッケージで実施した。実施例XXV:薬物スクリーニング統計的解析に含まれた各薬物は、データ品質についての以下のスクリーニング基準を満たした:1欠測値:細胞株の集合全体でGI50値の40以下が欠測していてもよい;2ばらつき:少なくとも3細胞株について、GI50>1.5.mGI50又はGI50<0.5.mGI50のいずれか、式中mGI50は所与の薬物についてのGI50中央値である。これらの基準を満たさない化合物は解析から除外した。実施例XXVI:SNPArray及びDNAコピー数解析AffymetrixゲノムワイドヒトSNPアレイ6.0を使用して、DNAコピー数データを測定した。アレイ品質及びデータ処理は、R統計フレームワークhttp://www.r-project.orgを使用してaroma.affymetrixに基づき実行した。乳癌細胞株SNPアレイを、記載のとおり20個の正常試料アレイを使用して正規化したBengtsson,H.,Irizarry,R.,Carvalho,B.Speed,T.P.Estimationandassessmentofrawcopynumbersatthesinglelocuslevel.BioinformaticsOxford,England24,759-7672008。データを、bioconductorパッケージのDNAcopyからサーキュラーバイナリセグメンテーションCBSを用いてセグメント化したOlshen,A.B.,Venkatraman,E.S.,Lucito,R.Wigler,M.Circularbinarysegmentationfortheanalysisofarray-basedDNAcopynumberdata.BiostatisticsOxford,England5,557-5722004。MATLABベースの、癌において有意な標的のゲノム同定GenomicIdentificationofSignificantTargetsinCancer:GISTICを使用して、有意なDNAコピー数変化を分析したBeroukhim,R.etal.Assessingthesignificanceofchromosomalaberrationsincancer:methodologyandapplicationtoglioma.ProcNatlAcadSciUSA104,20007-200122007。生データはEuropeanGenotypeArchiveEGAにおいて受託番号EGAS00000000059により利用可能である。コピー数の有意な変化を検出する確率が確実に最大となるようにするため、非悪性細胞株はGISTIC解析から除いた。各同質遺伝子細胞株ペアの一方のメンバーのGISTICスコアを使用して他方のゲノム変化を推論した:AU565はSKBR3から推論した;HCC1500はHCC1806から推論した;LY2はMCF7から推論した;ZR75BはZR751から推論した。実施例XXVII:エクソンアレイ解析細胞株の遺伝子発現データは、AffymetrixGeneChipHumanGene1.0STエクソンアレイから得た。発現の遺伝子レベルの要約を、aroma.affymetrixRパッケージを使用して、HuEx-_0-st-v2,coreチップタイプに基づく分位点正規化及び対数相加プローブレベルモデルPLMで計算した。転写物識別子を、BioMartRパッケージを使用してEnsemblデータベースに問い合わせることによりHGNC遺伝子記号に変換した。続いて得られた発現プロファイルをフィルタリングし、全細胞株においてlog2スケールで1.0より大きい標準偏差を表す遺伝子のみを捕捉した。生データはArrayExpressE-MTAB-181において利用可能である。実施例XXVIII:コンセンサスクラスタリング階層的コンセンサスクラスタリングを用いて細胞株サブタイプを同定したMonti,S.,Tamayo,P.,Mesirov,J.P.Golub,T.A.ConsensusClustering:AResampling-BasedMethodforClassDiscoveryandVisualizationofGeneExpressionMicroarrayData.MachineLearning52,91-1182003。コンセンサスは、細胞株の500サンプリング、試料当たり細胞株の80、凝集型階層的クラスタリング、ユークリッド距離計量及び平均連結法を用いて計算した。実施例XXIX:臨床的に関連するサブタイプと治療剤に対する応答との関連性本発明者らは、3つのスキームを使用してGI50を比較した:1ルミナル型対基底細胞型対クローディン低型;2ルミナル型対基底細胞型クローディン低型;及び3ERBB2-AMP型対非ERBB2-AMP型。群のGI50の違いを、順位尺度で、必要に応じてノンパラメトリックANOVA又はt検定により比較した。本発明者らは、3組の検定のp値を組み合わせ、偽発見率FDRを使用して複数の試験を補正した。3標本検定について、各群を他の全てと比較して、どの群が最も感度が高いかを決定することにより、有意なクラス効果を有する化合物に対して事後解析を実施した。事後検定のp値を合わせてFDR補正した。いずれの場合にも、FDRp<0.20を有意と見なした。スキーム2で基底細胞型クローディン低型の群が有意であると分かったが、スキーム1ではそれらの群のうち一つのみが有意であったケースであった場合、クラス特異度を割り当てる際にその3標本ケースに優先度を与えた。解析はRで実施した。実施例XXX:ゲノム変化と治療剤に対する応答との関連性本発明者らはt検定を用いて、繰り返し起こるコピー数変化8q24MYC、11q13CCND1、20q13STK15/AURKAにおけるものと薬物感受性との間の関連性を評価した。本発明者らは、増幅が低い、又はない細胞株を単一の群に組み合わせ、それを増幅が高い細胞株と比較した。欠失領域について同等の解析を行った。GI50が試験した最大濃度に等しかった細胞株は、解析から除外した。任意の群のサンプル数が5より小さい場合、化合物を除外した。実施例XXXI:成長速度と治療剤に対する応答との関連性細胞株クラス及び成長速度の薬物感受性に対する効果を評価するため、本発明者らは、上記に記載される3つの細胞株分類スキームの各々につき一つの、一組の二元配置共分散分析ANCOVA検定を実施した。これにより6組のp値が得られた2つの主効果×3分類スキーム;一回のFDR補正を用いて有意性を評価し、FDRp値<0.20が対象となることを宣言した。これらの分析はRにおいて関数lm及びANOVAで実施した。これらの関数はcarパッケージの一部として利用可能である。実施例XXXII:統合パスウェイ解析PARADIGMソフトウェアを使用して、コピー数、遺伝子発現、及びパスウェイ相互作用データの統合を実施した。簡潔に言えば、この手順は、単一の細胞株又は患者試料から、パスウェイ相互作用並びにゲノムデータ及び機能的ゲノムデータを使用して、遺伝子、複合体、及びプロセスについての統合パスウェイレベルIPLを推論する。詳細は実施例XLを参照のこと。実施例XXXIII:TCGA及び細胞株クラスタリング本発明者らは、細胞株について推論された活性が、TCGA腫瘍試料におけるそのそれぞれのサブタイプでクラスター化されるかという問いを立てた。高度に連結されたハブ遺伝子及び高度に相関した活性によって生じるバイアスを回避するため、相関分析により決定された2351個の非冗長な活性の集合を使用して、細胞株及び腫瘍試料をクラスター化した補足的方法を参照のこと。細胞株が同じサブタイプの腫瘍試料とクラスター化される程度を、コルモゴルフスミルノフ検定を用いて計算し、同じサブタイプの細胞株と腫瘍試料とのペア間の相関から計算されたt統計量の分布を、異なるサブタイプの細胞株ペアから計算された分布と比較した補足的方法を参照。詳細は実施例XLIを参照のこと。実施例XXXIV:サブタイプパスウェイマーカーの同定本発明者らは、まとめて特定のサブタイプに関する示差的活性を示す相互に連結された遺伝子を探索した。各サブタイプは、細胞株を2つの群に二分化するものとして扱った:一方の群は、そのサブタイプに属する細胞株を含み、第2の群は残りの細胞株を含んだ。本発明者らは、2クラスマイクロアレイ有意性解析SAMアルゴリズムのRインプリメンテーションを使用してTusher,V.G.,Tibshirani,R.Chu,G.Significanceanalysisofmicroarraysappliedtotheionizingradiationresponse.ProcNatlAcadSciUSA98,5116-5121,doi:10.1073/pnas.091062498[pii]2001、SuperPathwayの概念毎に示差的活性DAスコアを計算した。サブタイプについては、正のDAが、そのサブタイプで他の細胞株と比較してより高い活性に対応する。SuperPathwayにおける密接に関係する遺伝子の協調的な上方制御及び下方制御が、PARADIGMにより推論される活性を補強する。隣接遺伝子の活性もまた特定の表現型と相関する場合、高いDAスコアのサブネットワーク全体が見出されることが予想される。本発明者らは、2つの概念を連結するリンクであって、双方の概念が平均DA絶対値より高いDAスコアを有したリンクのみを保持することにより、高いDA絶対値の概念が相互連結されたSuperPathwayにおける領域を同定した。実施例XXXV:統合パスウェイ解析PARADIGMソフトウェアを使用して、コピー数、遺伝子発現、及びパスウェイ相互作用データの統合を実施した24。簡潔に言えば、この手順は、単一の細胞株又は患者試料から、パスウェイ相互作用並びにゲノムデータ及び機能的ゲノムデータを使用して、遺伝子、複合体、及びプロセスについての統合パスウェイレベルIPLを推論する。TCGABRCAデータは、2010年11月7日付けのTCGADCCから入手した。TCGA及び細胞株遺伝子発現データは、個別に、各データ集合内で中心となるプローブ中央値とした。データ集合全体細胞株又はTCGA腫瘍試料全体の値を全て順位変換し、-log10順位比に変換してから、PARADIGMに投入した。パスウェイは、http://pid.nci.nih.gov/からBioPaxレベル2フォーマットで得られ、NCI-PID、Reactome、及びBioCartaデータベースを含んだ。相互作用を組み合わせて、結合重畳パスウェイSuperPathwayとした。遺伝子、複合体、及び抽象的プロセス例えば細胞周期をパスウェイ概念として保持した。遺伝子概念を結合する前に、全ての遺伝子識別子をHUGO命名法に変換した。全ての相互作用が含まれ、矛盾する影響を解決することは試みなかった。P53最も多く連結された構成要素を始端とする幅優先無向探査を実行し、一つの単一構成要素を構築した。得られた結合パスウェイ構造は、3491個のタンパク質、4757個の複合体、及び520個のプロセスに相当する合計8768個の概念を含んだ。PARADIGM用の期待値最大化パラメータを細胞株データで訓練し、次にTCGA試料に適用した。次に細胞株及び腫瘍試料からのデータを組み合わせて単一のデータ行列とした。細胞株又は腫瘍試料のいずれかのデータにおいて0.5IPLを上回る少なくとも1つの値を有しないエントリーは全て、続く解析から除外した。実施例XXXVI:TCGA及び細胞株クラスタリングPARADIGMIPLを使用して細胞株をTCGA腫瘍試料と共にクラスター化し、細胞株が同じサブタイプの腫瘍試料と類似しているかどうかを決定した。SuperPathwayの十分に調べられた範囲は、多数の相互作用を有する遺伝子ハブと、直接的なデータは利用できない多数の中間的な複合体及び抽象的プロセスの大きいシグナル伝達鎖とを含む。ハブへのバイアスを回避するため、クラスタリング前に、細胞株及び腫瘍試料の双方に高度に相関するベクトルピアソン相関係数>0.9を含むパスウェイ概念を、単一のベクトルに統一した。この統一の結果、元の8939個のパスウェイ概念から2351個の非冗長性のベクトルとなった。得られた非冗長性概念の集合を使用して、試料をクラスター化した。47個の細胞株及び183個のTCGA腫瘍試料の双方についての推論されたパスウェイ活性の行列を、EisenClusterソフトウェアパッケージバージョン3.0に実装される完全連結法階層的凝集型クラスタリングを用いてクラスター化した。非中心化ピアソン相関をパスウェイ概念の計量として使用し、及び試料計量にユークリッド距離を使用した。細胞株が同じサブタイプの腫瘍試料とクラスター化される程度を定量化するため、本発明者らは、ピアソン相関から導き出されるt統計量の2つの分布を比較した。Csをサブタイプsの細胞株の集合とする。同様に、TsをサブタイプsのTCGA腫瘍試料の集合とする。例えば、Cbasal及びTbasalは、それぞれ、全ての基底細胞型細胞株及び基底細胞型腫瘍試料の集合である。第1の分布は、同じサブタイプの細胞株及び腫瘍試料を含むあらゆる可能なペア間のピアソン相関から導き出されたt統計量から構成された;すなわち全てのサブタイプsについて、a∈Cs及びb∈Tsであるようなペアa,b間のあらゆるペアワイズ相関t統計量を計算した。第2の分布は異なるサブタイプの細胞株間の相関t統計量から構成された;すなわち、a∈Cs及びb∈Cs及びs≠s’であるようなペアa,bに関して計算した。コルモゴルフスミルノフ検定を実施して分布を比較した。実施例XXXVII:統合パスウェイ解析PARADIGMソフトウェアを使用して、コピー数、遺伝子発現、及びパスウェイ相互作用データの統合を実施した24。簡潔に言えば、この手順は、単一の細胞株又は患者試料から、パスウェイ相互作用並びにゲノムデータ及び機能的ゲノムデータを使用して、遺伝子、複合体、及びプロセスについての統合パスウェイレベルIPLを推論する。TCGABRCAデータは、2010年11月7日付けのTCGADCCから入手した。TCGA及び細胞株遺伝子発現データは、個別に、各データ集合内で中心となるプローブ中央値とした。データ集合全体細胞株又はTCGA腫瘍試料全体の値を全て順位変換し、-log10順位比に変換してから、PARADIGMに投入した。パスウェイは、http://pid.nci.nih.gov/からBioPaxレベル2フォーマットで得られ、NCI-PID、Reactome、及びBioCartaデータベースを含んだ。相互作用を組み合わせて、結合重畳パスウェイSuperPathwayとした。遺伝子、複合体、及び抽象的プロセス例えば細胞周期をパスウェイ概念として保持した。遺伝子概念を結合する前に、全ての遺伝子識別子をHUGO命名法に変換した。全ての相互作用が含まれ、矛盾する影響を解決することは試みなかった。P53最も多く連結された構成要素を始端とする幅優先無向探査を実行し、一つの単一構成要素を構築した。得られた結合パスウェイ構造は、3491個のタンパク質、4757個の複合体、及び520個のプロセスに相当する合計8768個の概念を含んだ。PARADIGM用の期待値最大化パラメータを細胞株データで訓練し、次にTCGA試料に適用した。次に細胞株及び腫瘍試料からのデータを組み合わせて単一のデータ行列とした。細胞株又は腫瘍試料のいずれかのデータにおいて0.5IPLを上回る少なくとも1つの値を有しないエントリーは全て、続く解析から除外した。実施例XXXVIII:TCGA及び細胞株クラスタリングPARADIGMIPLを使用して細胞株をTCGA腫瘍試料と共にクラスター化し、細胞株が同じサブタイプの腫瘍試料と類似しているかどうかを決定した。SuperPathwayの十分に調べられた範囲は、多数の相互作用を有する遺伝子ハブと、直接的なデータは利用できない多数の中間的な複合体及び抽象的プロセスの大きいシグナル伝達鎖とを含む。ハブへのバイアスを回避するため、クラスタリング前に、細胞株及び腫瘍試料の双方に高度に相関するベクトルピアソン相関係数>0.9を含むパスウェイ概念を、単一のベクトルに統一した。この統一の結果、元の8939個のパスウェイ概念から2351個の非冗長性のベクトルとなった。得られた非冗長性概念の集合を使用して、試料をクラスター化した。47個の細胞株及び183個のTCGA腫瘍試料の双方についての推論されたパスウェイ活性の行列を、EisenClusterソフトウェアパッケージバージョン3.0に実装される完全連結法階層的凝集型クラスタリングを用いてクラスター化した45。非中心化ピアソン相関をパスウェイ概念の計量として使用し、及び試料計量にユークリッド距離を使用した。細胞株が同じサブタイプの腫瘍試料とクラスター化される程度を定量化するため、本発明者らは、ピアソン相関から導き出されるt統計量の2つの分布を比較した。Csをサブタイプsの細胞株の集合とする。同様に、TsをサブタイプsのTCGA腫瘍試料の集合とする。例えば、Cbasal及びTbasalは、それぞれ、全ての基底細胞型細胞株及び基底細胞型腫瘍試料の集合である。第1の分布は、同じサブタイプの細胞株及び腫瘍試料を含むあらゆる可能なペア間のピアソン相関から導き出されたt統計量から構成された;すなわち全てのサブタイプsについて、a∈Cs及びb∈Tsであるようなペアa,b間のあらゆるペアワイズ相関t統計量を計算した。第2の分布は異なるサブタイプの細胞株間の相関t統計量から構成された;すなわち、a∈Cs及びb∈Cs及びs≠s’であるようなペアa,bに関して計算した。コルモゴルフスミルノフ検定を実施して分布を比較した。実施例XXXIX:様々な遺伝分子レベルにおける腫瘍の分子サブタイプ乳房腫瘍に関して行われる全ゲノム遺伝子発現解析の先駆的研究から、最も顕著にはエストロゲン受容体ER陰性基底細胞様サブグループ及びER陽性ルミナルサブグループに属する種々のサブクラスが同定されておりPerou,C.M.etal.,2000,Molecularportraitsofhumanbreasttumours,406:747-752、その臨床転帰には違いがある14Sorlie,T.etal.,2001,Geneexpressionpatternsofbreastcarcinomasdistinguishtumorsubclasseswithclinicalimplications,98:10869-10874。いくつかの分子サブタイプの存在はまた、DNAコピー数解析2Russnesetal.2007上記、DNAメチル化Ronnebergetal.2011上記及びmiRNA発現解析Enerlyetal.2011上記によっても観察されている。しかしながら、問題は、分子解析によって様々な新しい分子レベルで得られたこれらの新しいプロファイルが、mRNA発現による当初見出されたサブクラスをどの範囲まで利用するか、及び臨床上重要な新規患者サブグループを同定するこれらの新しい分類の潜在力はどのようなものかである。これらの問いに取り組むため、本発明者らは初めにMicMaデータ集合の乳癌患者を、偏りのない教師なし方法を用いて、調べる各分子レベルに従いクラスター化した図23。分子レベル毎の患者のクラスタリングのヒストグラムを個別に、及び患者サブグループ毎の生存KMプロットを、図23に示す。興味深いことに、このクラスタリング手順により、Pam50分類から導き出されるクラスターと高度に相関した7つのmRNA発現クラスターの同定がもたらされる。これはPam50と一致し、但しルミナルAクラスターはexp1~4mRNAクラスターの間に、基底細胞及びERBB2は最後の3つexp5~7のクラスターの間に分かれた。miRNAレベルでは、Enerlyetal.2011上記に既に記載されるとおり、3つの異なるクラスターが得られた;メチル化レベルでは、記載されるとおりの3つの主クラスターと、はるかに小さい1つが認められ、この第4のクラスターはRonnebergetal.2011,上記でも観察されたものであるが、それ以上は考察されていない。CNAレベルでは6つの異なるクラスターが出現した。明らかに、あらゆるレベルで個別的な患者クラスターが特定の生存パターンと関連付けられた図23。次に、同じ患者が異なる分子レベルで対応するクラスターを形成したかどうかを評価した。実際、異なるレベルのクラスタリング間に、最も顕著にはDNAメチル化及びmRNA発現及びDNAコピー数の間に、良好な一致が大いにあった表12。しかしながら、試料によってはいかなるレベルでも常に共にクラスター化され、別の試料は研究における特定の各分子的エンドポイントに応じて異なる群にクラスター化される。一つの分子レベルに由来する一つのサブクラスを、もう一つによるクラスタリングにより一貫して分けることで、重要な生物学的関係が明らかとなり得る。例えば、3で考察されるとおり、メチル化とmRNA発現との間の良好な相関に基づく分類が認められたがp=2.29106、ルミナル-AクラスmRNA発現によるはさらになお、2つの異なるメチル化クラスターの間に分けられた。同じことが基底細胞様腫瘍にも該当したことから、mRNA発現クラスターとの強い一致にもかかわらず、DNAメチル化に従うクラスタリングによってさらなる情報が提供されたことが示唆される。異なるDNAメチル化プロファイルを有するルミナルA試料は、生存が異なる3Ronneberg,J.A.etal.,2011,Methylationprofilingwithapanelofcancerrelatedgenes:associationwithestrogenreceptor,TP53mutationstatusandexpressionsubtypesinsporadicbreastcancer,5:61-76。本発明者及び他の研究者の双方からの新しいデータ集合の数が増えることで、将来、これらのクラスターがいくつかの最も高頻度及び多数の低頻度の組み合わせに収束するかどうかが明らかとなるであろう。異なる分子レベルでの再分類は、異なるレベルで影響を受ける新しい有益な生物学的パスウェイが指示され得るため、さらに調べるに値するものの、このようにクラス同士でサンプルを水平に入れ替えるなかでの情報量は限定的であり得る。パスウェイ毎にこれらのクラスターの範囲内で示差的に発現/変化する遺伝子を調べることは、先験的知識及び既知の相互作用の選択に依存し、新規のパスウェイを同定することはできない。さらに、これらの手法は異なるデータ集合における遺伝子及び計測値を独立変数として扱い、パスウェイにおける遺伝子の位置、又はその相互に作用するパートナーの数すなわちパスウェイのトポロジーは考慮せず、遺伝子集合中の1つ又は少数の遺伝子の発現における大きい変動に対して脆弱であり得る。一般に、癌における多くの腫瘍で特定のパスウェイが調節解除され得るが、詳細な遺伝子及び調節解除の方法は腫瘍によって異なることが観察されているCancerGenomeAtlasResearchNetwork.Comprehensivegenomiccharacterizationdefineshumanglioblastomagenesandcorepathways.Nature2008Oct.;4557216:1061-1068。従って、次に本発明者らは、パスウェイ及び関連する臨床データの文脈における腫瘍内での各遺伝子の活性レベルを特徴付けるため、単一の遺伝子に関する異なるデータ型の計測値間の相互作用並びに既知の遺伝子間相互作用をモデル化するパスウェイベースのモデリング法を適用した。本発明者らは各遺伝子の統合パスウェイレベルIPLを用いて患者を直接同定して、これらの調節解除されたパスウェイに従い分子データ型間で分類し、次に様々な分子レベルにおけるこの新しいクラスターの先述のクラスとの関係を調べた。実施例XL:予後判定有意性により浸潤癌を分類するためのPARADIGM腫瘍表現型を説明し、且つ腫瘍を標的化治療され易いものにし得る個別の生物学的機能を、ゲノム変化がどのように妨げるかを理解するためには、パスウェイレベルの摂動を理解する必要がある。PARADIGMは、遺伝子を単一のレベルで調べた場合には区別できない患者サブセットにおいて、一貫した活性パスウェイを同定する。この方法は、既知のパスウェイ構造に対して、確率的グラフィカルモデルPGMから統合ゲノミクスデータに至るまでの技法を用いる。これは既に、TCGA膠芽腫及び卵巣データ集合のコピー数及びmRNA発現データの解析に適用されている。PARADIGM解析はまた、DNAメチル化又はコピー数、mRNA及びmiRNA発現などの、複数のレベルでゲノム変化を関係付けるためにも用いることができ、従って個々の各試料における任意の数のオミクスデータ層を統合することができる。DNAメチル化及びmiRNA発現は、ここで観察された調節解除パスウェイに寄与し、且つMicMaコホートにおいて各々それ自体で乳癌の予後判定及び分子プロファイルに明確な寄与を有するように思われるが図23、本発明者らは、これらの2つの分子プロファイルタイプを追加することによるPARADIGMクラスターの予後判定値の向上は認めなかった。これについての一つの説明は、miRNA及びDNAメチル化解析の予後判定値が、その高い相関から、mRNA発現によっても利用されるというものである。しかしながら、そのように結論付けるには、例えば、解析プラットフォームの選択メチル化についての限られたIllumina1505CpG癌パネル及び真のmiRNA標的についての我々の限られた知識が、総合的に計測し、且つmiRNA及びDNAメチル化情報を有効にモデル化する我々の能力を制限する要因となり得るかどうかに関して、さらなる分析を要する。MicMaコホートのmRNA発現及びコピー数変化に基づくPARADIGM解析により、5つの異なるクラスターの存在が同定され図24A、mRNA発現とDNAコピー数とを組み合わせると、個別に研究される分子レベルのいずれと比べても、予後に関して患者のより優れた判別につながることが示された図24B及び図23。その摂動がこの分類に最も強力に寄与したパスウェイは、アンジオポエチンAngiopoientin受容体Tie2媒介性シグナル伝達のもの、及び最も顕著には、免疫応答TCR及びインターロイキンシグナル伝達のものであり、ここではパスウェイにおけるほぼ全ての遺伝子又は複合体が正常値から外れた図25A。最も著明に認められたのは、IL4、IL6、IL12及びIL23シグナル伝達であった。他の著明なパスウェイは、エンドセリン、同様に卵巣で調節が解除されるFOXM1転写、並びに同様に乳癌及び卵巣癌で調節が解除されることが既に分かっている膠芽腫TCGAデータ集合及びERBB4であった。この解析に基づき、本発明者らは、有意に異なる予後を有する以下の患者群を同定した。これらの群は、おおまかに以下のとおり特徴付けることができる:pdgm.1=高FOXM1、高免疫シグナル伝達、pdgm.2=高FOXM1、低免疫シグナル伝達、マクロファージ優勢、pdgm.3=低FOXM1、低免疫シグナル伝達、pdgm.4=高ERBB4、低アンギオポエチンシグナル伝達、pdgm.5=高FOXM1、低マクロファージシグネチャー。2つの既に発表されているデータ集合でParadigmクラスターの同定を検証した。一つはChinetal2007Chin,S.F.etal.,2007,Usingarray-comparativegenomichybridizationtodefinemolecularportraitsofprimarybreastcancers,26:1959-1970によるもので、これはMicMaデータ集合と比較して、ER及び高悪性度腫瘍の頻度がより高かったとともに、さらに一層興味深いことに、別の集合では非悪性DCIS非浸潤性乳管癌が濃縮された12Muggerud,A.A.etal.,2010,MoleculardiversityinductalcarcinomainsituDCISandearlyinvasivebreastcancer,4:357-368図25B、図25C。純粋なDCIS腫瘍のヒートマップを図25Dに示す27。MicMaにおいて最も不良な予後のクラスターpdgm.2では、IL4シグナル伝達がSTAT6と共に強く下方制御され、STAT6は、ヒト乳癌細胞において成長阻害を妨げることが示されている16Gooch,J.L.,Christy,B.,andYee,D.,2002,STAT6mediatesinterleukin-4growthinhibitioninhumanbreastcancercells,4:324-331。IL4シグナル伝達の下方制御はまた、より高い腫瘍成長を支持し得るマスト細胞の活性化も促進している17deVisser,K.E.,Eichten,A.,andCoussens,L.M.,2006,Paradoxicalrolesoftheimmunesystemduringcancerdevelopment,6:24-37。逆にpdgm.5では、IL23シグナル伝達により、マクロファージ活性化が低下し、ナチュラルキラー細胞活性が増加する。免疫応答が、一方の側面でTh-2細胞及びB細胞の動員に、他方でTh-1の増殖に向かい、癌依存的に二極化することが考察されている1Ursini-Siegel,J.etal.,2010,Receptortyrosinekinasesignalingfavorsaprotumorigenicstateinbreastcancercellsbyinhibitingtheadaptiveimmuneresponse,70:7776-7787。一定の条件下でTh1/CTL免疫応答が、マウスにおける過形成の腺腫への移行を妨げ得る一方、Th2応答が、癌腫への移行を促進する慢性炎症状態を付与することにより得るという仮説が立てられている。IL4はTh-2由来のサイトカインであり、B細胞分化及び癌細胞での慢性炎症を刺激する。さらにTh-2細胞は、これらの癌における免疫抑制を媒介するIL10を分泌する。この免疫抑制は、主に基底細胞癌及びERBB2癌で起こることが示されている。この裏付けとして、最近、腫瘍促進性微小環境で抗腫瘍性獲得免疫プログラムが奪われ、代わりに、上皮細胞挙動の調節に機能的に関与する先天性免疫系の細胞成分が関わることにより悪性腫瘍が促進され得ることが示されているDeNardo,D.G.etal.,2009,CD4Tcellsregulatepulmonarymetastasisofmammarycarcinomasbyenhancingprotumorpropertiesofmacrophages,16:91-102。ここに提案されるこの免疫分類と、十分に確立されたmRNA発現による分類ルミナルA、B、基底細胞、ERBB2、正常様との間には、少なからぬ一致があった図24.基底細胞及びERBB2クラスターに属するサンプルは、主に、prgm1最も不良な予後、ルミナルA-prgm3最も良好な予後のものであった。しかしながら、Paradigmクラスタリングは、ルミナルAprgm3とルミナルBprgm4とのクラスター間についてのかなり有意な区別、並びに予後が極めて不良な基底細胞型腫瘍サブセットprgm2の同定を提供する。実施例XLI:摂動がPARADIGMクラスタリングに特異的に影響を与える同定されたパスウェイFOXM1転写FOXM1は、細胞周期進行の主要調節因子であり、その内因的FOXM1発現は細胞周期の相に従い振動する。ヒト癌原遺伝子として確認されたFOXM1は、肝癌、乳癌、肺癌、前立腺癌、子宮頸癌、結腸癌、膵癌、脳癌並びに最も一般的なヒト癌の基底細胞癌を含め、ヒト固形癌の大多数で上方制御されることが分かっている。FOXM1は、細胞周期及び染色体/ゲノム維持におけるその複数の役割によって腫瘍形成を促進すると考えられるWonsey,D.R.andFollettie,M.T.,2005,LossoftheforkheadtranscriptionfactorFoxM1causescentrosomeamplificationandmitoticcatastrophe,65:5181-5189。初代ヒト皮膚ケラチノサイトにおけるFOXM1の異常な上方制御は、ヘテロ接合性の喪失LOH及びコピー数異常の形でのゲノム不安定性を直接引き起こし得るTehM,GemenetzidisE,ChaplinT,YoungBD,PhilpottMP.UpregulationofFOXM1inducesgenomicinstabilityinhumanepidermalkeratinocytes.Mol.Cancer2010;9:45。最近の報告では、成人ヒト上皮幹細胞におけるFOXM1の異常な上方制御が3D器官型組織再生系において前癌表現型-ヒト過形成と同様の状態-を引き起こすことが示されたGemenetzidis,E.etal.,2010,Inductionofhumanepithelialstem/progenitorexpansionbyFOXM1,70:9515-952。この著者らは、FOXM1の過剰発現が、分化経路に干渉し、それにより前駆細胞コンパートメントを拡大することによって幹細胞の固有の自己再生増殖能を利用することを示した。従って、FOXM1が幹細胞/前駆細胞の増殖によって癌の発生を引き起こすという仮説が立てられた。本発明者らは、主にインターロイキンシグナル伝達活性に従い分けられた、このパスウェイの活性が高い群及び低い群の2つの乳癌患者群を明らかに認めている。図26は、クラスターpdgm3最良の予後についての、このパスウェイの逆の活性化のし方を示し活性化されたとき赤色、対して不活性化されたとき青色、これはより不良な生存及びそれに寄与する分子レベル図の形に従いmRNA、CNA、miRNA又はDNAメチル化を伴うその他のクラスターとは対照的である。pdgm3におけるMMP2の下方制御がDNAメチル化による一方、その他の腫瘍ではDNA欠失によることが分かる。miRNAのなかで、has-let7-bはpgm3において上方制御され、その他では下方制御され、その標的AURKBと相補的であった。DNA増幅及びmRNA発現の双方が発現の調節解除の原因と見られた。アンギオポエチン受容体tie2媒介性シグナル伝達Angファミリーは、ヒト癌の発症及び成長の間の血管新生において重要な役割を果たす。Ang2の血管新生における役割は、概してAng1の拮抗体であって、血管成熟及び安定化に重要な、Ang1が促進するTie2シグナル伝達を阻害すると考えられる23。Ang2は、別の重要な血管新生因子である血管内皮増殖因子AVEGFAと共に、協働的に血管新生を調整するHashizume,H.etal.,2010,Complementaryactionsofinhibitorsofangiopoietin-2andVEGFontumorangiogenesisandgrowth,70:2213-2223。新しいデータは、ヒト癌が進行する間の浸潤性の表現型の癌細胞の血管新生におけるAng2についてのより複雑な役割を示唆する。特定のアンギオポエチンAngファミリーメンバーはTie1を活性化することができ、例えばAng1は内皮細胞においてTie1リン酸化を誘発する2Yuan,H.T.etal.,2007,ActivationoftheorphanendothelialreceptorTie1modifiesTie2-mediatedintracellularsignalingandcellsurvival,21:3171-3183。しかしながら、Tie2が内皮細胞において下方制御されるとき、Ang1はTie1リン酸化を誘発できず、及びAng1が存在しない場合に、Tie1リン酸化は構成的に活性な形態のTie2又はTie2作動性抗体のいずれかにより誘発されるため、Tie1リン酸化はTie2依存性である25Yuanetal.2007上記。Ang1媒介性AKT及び42/44MAPKリン酸化は主にTie2媒介性であり、Tie1はこのパスウェイを下方制御する。従ってTie1の主な役割は、Tie2駆動性のシグナル伝達を下方制御するその能力及び内皮生存により、血管形態形成を調整することである。Tie2媒介性シグナル伝達並びにVEGFR1及び2媒介性シグナル伝達の双方及び特異的シグナルが、このデータ集合で観察された。ERBB4ERBB4は、乳房形態形成における増殖及び細胞運動並びにErbb4を発現する乳房原始上皮の方向性のある細胞運動に寄与する一方で、乳房細胞運命を促進する。Nrg3/Erbb4シグナル伝達の候補エフェクターは同定されており、ここでは初期の乳腺発生及び癌に関連する他のシグナル伝達パスウェイと相互作用することが示されている。ErbB4の生体内での主要な機能の一つは、妊娠及び乳汁分泌誘導中の乳腺の成熟にある。妊娠及び長期の乳汁分泌期間は乳癌リスクの低下と関連付けられており、従って腫瘍抑制におけるErbB4の役割は、乳汁分泌におけるその役割と結び付けられ得る。ほとんどの報告は、思春期に他のErbBファミリーメンバーにより惹起される成長刺激を反転させることにおけるErbB4の役割と一致するが、しかしながらERBB4発現との生存の有意な関連性は確認されていない2Sundvall,M.etal.,2008,RoleofErbB4inbreastcancer,13:259-268。実施例XLII:非浸潤性乳管癌DCISにおける分類のためのPARADIGMマウスモデルでの前癌性過形成腺における免疫応答の関与を所与として18Ursini-Siegel,J.etal.,2010,Receptortyrosinekinasesignalingfavorsaprotumorigenicstateinbreastcancercellsbyinhibitingtheadaptiveimmuneresponse,70:7776-7787、本発明者らは、DCISケースから構成される既発表のデータ集合を分析し、浸潤性腫瘍で観察される強い免疫応答及びインターロイキンシグナル伝達が前癌期においても同様に存在するかどうかを見出した。非浸潤性乳管癌DCISは非浸潤型の乳癌であり、ここでは一部の病変が急激に浸潤性乳管癌IDCに移行する一方、他の病変は変化しないままであると考えられる。本発明者らは、過去に、31例の純粋なDCIS、36例の純粋な浸潤癌及び42症例の混合診断浸潤癌で非浸潤性要素を伴うものの遺伝子発現パターンを研究し1Muggerudetal.2010上記、及び組織学的悪性度が高いDCIS間でのトランスクリプトームの多様性を観測して、進行性腫瘍により類似した遺伝子発現特徴を有する個別的なDCISサブグループを同定している。このコホート全体IDC及びILCを含むについてのPARADIGM結果のヒートマップは図25C、及び純粋なDCIS試料についてのヒートマップは図25D。純粋なDCIS腫瘍には、高マクロファージ活性に典型的なシグナル伝達によって特徴付けられるprgm2型のものはなかった図25。それと一致して、実験研究から、原発性乳腺腺癌におけるマクロファージが、その血管新生促進特性の結果として後期発癌を調節するとともにLin,E.Y.andPollard,J.W.,2007,Tumor-associatedmacrophagespresstheangiogenicswitchinbreastcancer,67:5064-5066;Lin,E.Y.etal.,2007,Vascularendothelialgrowthfactorrestoresdelayedtumorprogressionintumorsdepletedofmacrophages,1:288-302、悪性乳房上皮細胞に上皮成長因子EGFを提供することにより肺転移を助長することが実証されている。ここでも、DCISにおいてPARDIGM解析により同定された上位の調節解除パスウェイには、IL2、4、6、12、23、及び23シグナル伝達に関与するものがあった。双方のデータ集合DCIS、MicMaにおいて、ナイーブCD8T細胞におけるTCRシグナル伝達は、CD8T細胞を動員することが知られる多数のケモカインと共にリストの上位にあった。一つはIL-12であり、これは、NK細胞及びT細胞からのIFN-γ産生を刺激することが示された抗原提示細胞により産生される。IFN-γパスウェイは、調節解除されたパスウェイの一つで、DCISにおけるリストの上位であった。IFNγはTh1細胞及びNK細胞から産生され、抗腫瘍免疫応答を惹起することが示された。第I相臨床試験では、トラスツズマブハーセプチンherceptinの臨床効果が、HER2過剰発現腫瘍を有する患者にIL-12を同時投与すると増強されることが示されており、この効果はNK細胞におけるIFNγ産生の刺激により媒介される29。DCISでは、他の最も強力な寄与因子表8は84_NOX4であった。酸素感受性NAPHDオキシダーゼのNOX4、及び食細胞型Aオキシダーゼは、好中性顆粒球における多量の活性酸素種ROSの産生、一次免疫応答に関与するものと同様である。また、FN1フィブロネクチン、及び血小板由来成長因子受容体であるPDGFRBも、COL1A2、IL12/IL12R/TYK2/JAK2/SPHK2、ESR1及びKRT14と共に、DCISにおいて合わせて特異的に繰り返し現れた。これらの遺伝子/パスウェイは全て、細胞外マトリックスにおける機能、細胞間相互作用、並びに線維症及び角質化に寄与しているように見える。例えば、フィブロネクチン-1、FN1は、細胞表面上、細胞外液、結合組織、及び基底細胞膜に存在する高分子量糖タンパク質ファミリーに属する。フィブロネクチンは他の細胞外マトリックスタンパク質及び細胞性リガンド、例えば、コラーゲン、フィブリン、及びインテグリンと相互作用する。フィブロネクチンは細胞の接着及び遊走プロセスに関与する。血小板由来成長因子受容体PDGFRは、上皮成長因子EGFと共に、重要な受容体チロシンキナーゼRTKであるEGF及びPDGF受容体を介してシグナルを伝達する。重要なことに、ここで特定のDCISにおいて過剰発現することが分かったPDGFRは、スニチニブの標的であり30Fratto,M.E.etal.,2010,Newperspectives:roleofsunitinibinbreastcancer,161:475-482、及びメシル酸イマチニブグリベックGleevecの二次標的であるWeigel,M.T.etal.,2010,Invitroeffectsofimatinibmesylateonradiosensitivityandchemosensitivityofbreastcancercells,10:412。INFγ産生を増加させることにより媒介される上記のトラスツズマブハーセプチンherceptinの免疫賦活の役割とは逆に、イマチニブはTCR活性化CD4T細胞によるインターフェロン-γ産生を阻害することが示された。これらの観察は、DCIS及び悪性細胞の表面に提示される成長因子受容体間の相互作用並びに免疫構成を明らかにする点で、本発明者らの議論にとって興味深い。PDGFRに対する刺激性自己抗体が、I型コラーゲンの発現増加をもたらすRas、ERK1/ERK2、及び活性酸素種ROSを含む細胞内ループを惹起すると見られることが示された。これは、本発明者らの研究で同様にDCISにおいて調節解除されるものとして観察されたCOL1A2発現と一致する。実施例XLIII:材料及び方法本解析は、約110例の乳癌から収集したデータに適用し、mRNA発現はAgilent全ヒトゲノム4×44K1色法オリゴアレイにより解析した。コピー数変化CNAは、IlluminaHuman-1109KBeadChipを使用して解析した。このSNPアレイは遺伝子中心的であり、平均物理距離が30kbのゲノム全体を網羅するマーカーを含み、15,969個の固有の遺伝子を表す2004年5月作成、hg17、NCBIBuild35。各試料を全ゲノム増幅に供した。BeadStudiov.2.0、Illuminaを使用して、dbSNPbuild125の順方向の対立遺伝子の向きを参照して遺伝子型レポート及びlogR値を抽出し、CNAについてlogR値を調整した。AgilentTechnologiesHumanmiRNAMicroarrayKitV2を使用して、製造者のプロトコルに従い全RNAからのmiRNAプロファイリングを実施した。AgilentScannerG2565Aでのスキャニング及びFeatureExtractionFEv9.5を使用してシグナルを抽出した。デュプリケートのハイブリダイゼーションを用いて99試料、種々のアレイ及び時間点で実験を行った。2つの試料は1回のみプロファイル決定した。複製プローブのmiRNAシグナル強度をプラットフォームで平均化し、log2変換して、75パーセンタイルに対して正規化した。FEv9.5の初期設定により、miRNA発現状態を各試料中における各遺伝子の存在又は非存在としてスコア化した。DNAメチル化。EpiTect96重亜硫酸キットBisulfiteKitQiagenGmbH、独国を使用して、1マイクログラムのDNAを重亜硫酸処理した。500ngの重亜硫酸処理したDNAを、807個の癌関連遺伝子における1505ヶ所のCpG部位を同時に解析するGoldenGateメチル化癌パネルIMethylationCancerPanelIIlluminaInc、CA、米国を使用して解析した。遺伝子当たり少なくとも2つのCpG部位が解析され、ここで1つのCpG部位はプロモーター領域にあり、1つのCpG部位は最初のエクソンにあり、メチル化データの初期処理にはBeadstudioソフトウェアを、製造者のプロトコルに従い使用した。各CpG部位の検出p値を使用して標本パフォーマンスを検証し、及び検出p値に基づきデータ集合をフィルタリングし、ここで検出p値>0.05のCpG部位は以降の解析から除外した。データ前処理及びParadigmパラメータ。コピー数をCBSを用いてセグメント化し、次に、hg18のRefSeq遺伝子座標にわたる全てのセグメントの中央値を取ることにより、遺伝子レベルの計測値に対応させた。mRNA発現について、計測値は、各プローブについて発現値中央値を減じることにより、第1のプローブで正規化した。各プローブの製造者のゲノム位置を、UCSCのliftOverツールを使用してhg17からhg18に変換した。次に、RefSeq遺伝子に重なる全てのプローブの中央値を取ることにより、遺伝子毎の計測値を得た。製造者の説明を用いてメチル化プローブを遺伝子と対応させた。各データ集合を個別に分位点変換することにより、先のとおりParadigmを実行し10、しかしデータは、5及び95分位点においてではなく、等しいサイズのビンに離散化した。パスウェイファイルは、先に解析したとおりPID36からであった。図26は、各データ型における上のビン又は下のビンのいずれかの観測割合をカウントし、次に任意のデータ型において観測割合が最大のビンで各ノードをラベル化することによる、離散化した入力データの及びIPL値ではない要約を示す。HOPACH教師なしクラスタリング。Rバージョン2.12上で動くHOPACHRインプリメンテーションバージョン2.1037を使用して、クラスターを導き出した。全てのデータ型で相関距離計量を使用し、但しParadigmIPLは例外で、正規分布でなく、且つゼロ値が見られたため、cosangleを使用した。含まれるサンプル数が5より小さい任意のサンプルクラスターについては、各サンプルを、より大きいクラスター中の最も類似するサンプルと同じクラスターにマッピングした。MicMaデータ集合のParadigmクラスターは、MicMaデータ集合における各クラスターのメドイドmediod中央値関数を使用を決定し、次に別のデータ集合の各サンプルを、cosangle距離により最も近かったいずれかのクラスターメドイドmediodに割り当てることにより、他のデータ型にマッピングした。カプランマイヤーKaplain-Meier、クラスター濃縮。Rバージョン2.12を使用して、カプランマイヤー統計、プロット、及びクラスター濃縮を決定した。coxph比例ハザードモデルからワルド検定を使用してコックスp値を決定し、及びsurvdiff関数からカイ二乗検定によりログランクp値を決定した。クラスタリングについての遺伝子の値又はパスウェイメンバーの値の全体的な濃縮をANOVAにより決定し、及び特定のクラスターラベルの遺伝子の濃縮を、特定のクラスターにおける遺伝子の値と、他の全てのクラスターにおける遺伝子の値とのT検定により決定した。p値調整のベンジャミニとホッホバーグBenjaminiHochbergの方法を用いてFDRを決定した。実施例XLIV:データ集合及びパスウェイ相互作用コピー数及び発現データの双方をPARADIGM推論に組み入れた。8個の正常組織対照の集合が発現データにおける解析に利用可能であったため、患者の遺伝子値の各々を、正常な卵管対照で観察される遺伝子の中央値レベルを減じることにより正規化した。コピー数データは、腫瘍で検出された遺伝子レベルと、それに対する血液正常レベルとの間のコピー数の差を反映するように正規化した。PARADIGMへの入力のため、発現データは、サブタイプ解析に使用したものと同じ統合データ集合から取り、コピー数は、MSKCCAgilent1Mコピー数データのセグメント化コールから取った。131個のパスウェイ、11,563個の相互作用、及び7,204個のエンティティを含むNCI-PIDから一群のパスウェイを入手した。エンティティは、PARADIGMのグラフィカルモデルにおいてノードとして表現される分子、複合体、小分子、又は抽象概念である。抽象概念は、一般的な細胞プロセスアポトーシス又は光の吸収など、及びシグナルトランスデューサーのRASファミリーなどの機能的活性を共有する遺伝子のファミリーに対応する。本発明者らは、タンパク質間相互作用、転写調節相互作用、リン酸化及びユビキチン化などのタンパク質修飾相互作用を含む相互作用を収集した。実施例XLV:パスウェイ文脈における統合分子活性の推論本発明者らは、コピー数、遺伝子発現、及び各エンティティのパスウェイ文脈を反映する統合パスウェイ活性IPAを割り当てるPARADIGMを使用した。データの遺伝子特異的及び患者特異的断面の順列を使用して、IPAの有意性を評価した。ゲノム中の各遺伝子についての値の遺伝子発現とコピー数とのペアを無作為に選択することにより、1000個のヌル患者のデータを作成した。PARADIGMIPAの有意性を評価するため、パスウェイ構造を維持しながらパスウェイにランダム遺伝子を割り当てることにより、ヌル分布を作成した。実施例XLVI:FOXM1パスウェイの同定ランダムシミュレーションの間は、FOXM1ネットワーク内の全ての遺伝子を使用して統計的有意性を評価したが、FOXM1パスウェイを視覚化するため、図29により有意に変化したIPAを有するFOXM1に直接連結されたエンティティを、図27に含めるため選択した。これらのうち、DNA修復及び細胞周期制御において役割を有する遺伝子であって、FOXM1との相互作用について文献の裏付けがあると認められたものを表示した。元のNCI-PIDパスウェイに存在しなかったBRCC複合体メンバーを、NCI-PIDによればFOXM1の標的であるBRCA2と共にプロットに含めた。上流DNA修復標的を、他のNCIパスウェイにおけるCHEK2の上流直接因子を見つけることにより同定した例えば、PLK3シグナル伝達パスウェイにATMからの間接的リンクが見出された。実施例XLVII:クラスタリング活性及び非活性の確率の変化を直接的に表す推論活性の使用は、それにより様々な種類のエンティティをまとめて一つのヒートマップにクラスター化することを可能にする。PARADIGM推論の結果を包括的に視覚化するため、EisenCluster3.0を使用して特徴フィルタリング及びクラスタリングを実施した。0.1の標準偏差フィルタリングにより7204個中1598個のパスウェイエンティティが残り、エンティティ及びサンプルの双方に対して平均連結法非中心化相関階層的クラスターを実施した。実施例XLVIIIゲノムDNAの単離血液試料2~3mlを患者から採取し、使用時まで、EDTAが入った試験管に-80°Cで保存する。この血液試料から、DNA単離キットPUREGENE,GentraSystems,MinneapolisMNを使用して製造者の指示に従いゲノムDNAを抽出する。DNA純度は、ベックマン分光光度計で計測した260nm及び280nmにおける吸光度の比1cm光路;A260/A280として計測される。実施例XLIX:SNPの同定患者のDNA試料からの遺伝子領域は、PCRにより、その領域用に特別に設計されたプライマーを使用して増幅する。PCR産物は、上記に開示されるとおりの、当業者に公知の方法を用いて配列決定する。配列トレースで同定されるSNPを、Phred/Phrap/Consedソフトウェアを使用して検証し、NCBISNPデータバンクに寄託されている既知のSNPと比較する。実施例L:統計的解析値は平均値±SDとして表現される。χ2分析WebChiSquareCalculator,GeorgetownLinguistics,GeorgetownUniversity,WashingtonDCを使用して、正常な対象と障害を有する患者とにおける遺伝子型頻度の差を評価する。事後解析を伴う一元配置ANOVAを示されるとおり実施して、異なる患者群間の血行動態を比較する。当業者は、上述の実施形態の様々な適応例及び改良例を本発明の範囲及び精神から逸脱することなく構成し得ることを理解するであろう。当該技術分野において公知の他の好適な技法及び方法が、当業者によって数多くの具体的な様式で、且つ本明細書に記載される本発明の説明をふまえて適用され得る。従って、本発明は本明細書に具体的に記載される以外に従い実施され得ることが理解されるべきである。上記の説明は例示的であることを意図し、限定的であることを意図するものではない。上記の説明を考察することで、当業者には多くの他の実施形態が明らかであろう。従って本発明の範囲は、かかる特許請求の範囲に適用される均等物の全範囲と共に、添付の特許請求の範囲を参照して決定されるべきである。"