(57)【要約】生物システムの状態のプロファイルを作製するための方法およびシステムが開示される。このプロファイルは、複数のデータセット間の類似性、相違、および／または相関の識別に基づく。このデータセットは、１つ以上の生体分子成分のタイプ、１つ以上の生物学的サンプルのタイプ、および／または１つ以上の測定のタイプから算出される。本発明の１つの局面において、生物システムの状態（例えば、疾患状態）のプロファイルを作製するために、複雑な生物学的サンプルに関する複数の測定が実施される。その後、相関分析およびネットワークモデリングと組み合わせた、遺伝子、遺伝子転写物、タンパク質、および／または代謝産物についての包括的なプロファイリングが、システムレベルでの生物システムへの洞察を提供する。
2007502992.tif 000002
【特許請求の範囲】
【請求項１】
哺乳動物における生物システムの状態をプロファイリングする方法であって、該方法は、
（ａ）生物システムの複数のデータセットを統計分析で評価し、そして、該複数のデータセットの少なくとも１部分の間の差異を１セット以上決定するために、該複数のデータセット間で特徴を比較する工程；ならびに
（ｂ）工程（ａ）の結果に基づいて該生物システムの状態についてのプロファイルを作製する工程
を包含し、該複数のデータセットは、１つより多い生物学的サンプルのタイプ、１つより多い測定技術のタイプ、１つより多い生体分子成分のタイプ、または、生物学的サンプルのタイプ、測定技術、および生体分子成分のうちの少なくとも２つとの組合せから算出される測定値を含む、方法。
【請求項２】
請求項１に記載の方法であって、前記生物システムは、ヒトにおけるシステムである、方法。
【請求項３】
請求項１に記載の方法であって、前記統計分析は、多変量解析を含む、方法。
【請求項４】
請求項１に記載の方法であって、前記生物学的サンプルのタイプは、血液、血漿、血清、脳脊髄液、胆汁、唾液、滑液、胸膜液、心膜液、腹水、汗、糞便、鼻液、眼液、細胞内液、細胞間液、リンパ液、尿、肝細胞、上皮細胞、内皮細胞、腎臓細胞、前立腺細胞、血液細胞、肺細胞、脳細胞、皮膚細胞、脂肪細胞、腫瘍細胞、および乳房細胞を含む群より選択される、方法。
【請求項５】
請求項１に記載の方法であって、複数のデータセットは、別々に処理される１つの生物学的サンプルのタイプから算出されるか、あるいは、異なる時間に収集または分析される１つの生物学的サンプルのタイプから算出される、方法。
【請求項６】
請求項１に記載の方法であって、前記測定技術は、液体クロマトグラフィー、ガスクロマトグラフィー、高速液体クロマトグラフィー、キャピラリー電気泳動、質量分析、液体クロマトグラフィー質量分析、ガスクロマトグラフィー質量分析、高速液体クロマトグラフィー質量分析、キャピラリー電気泳動質量分析、核磁気共鳴分析、並行ハイブリダイゼーションアッセイ、並行サンドイッチアッセイ、および、競合アッセイを含む群より選択される、方法。
【請求項７】
請求項１に記載の方法であって、複数のデータセットは、単一のタイプの測定技術の異なる機器構成からの測定値を含む、方法。
【請求項８】
請求項１に記載の方法であって、前記生体分子成分のタイプは、遺伝子、遺伝子転写物、タンパク質、または代謝産物である、方法。
【請求項９】
請求項１に記載の方法であって、該方法は、前記生物システムの状態についての前記プロファイルと、プロファイルデータベースとを比較する工程を包含する、方法。
【請求項１０】
請求項１に記載の方法であって、該方法は、前記生物システムの状態についての前記プロファイルと、生物システムの別の状態についてのプロファイルとを比較する工程を包含する、方法。
【請求項１１】
コンピューターが読み取り可能な媒体を有する製品であって、請求項１に記載の方法を実施するために、コンピューターが読み取り可能な命令が該媒体に具体化されている、製品。
【請求項１２】
哺乳動物における生物システムの状態をプロファイリングする方法であって、該方法は、
（ａ）生体分子成分のタイプについての複数のデータセットを統計分析で評価し、そして、該複数のデータセットの少なくとも１部分の間の差異を１セット以上決定するために、該複数のデータセット間で特徴を比較する工程；
（ｂ）別の生体分子成分のタイプについての複数のデータセットを統計分析で評価し、そして、該複数のデータセットの少なくとも１部分の間の差異を１セット以上決定するために、該複数のデータセット間で特徴を比較する工程；ならびに
（ｃ）該生物システムの状態のついてのプロファイルを作製するために、工程（ａ）および工程（ｂ）の結果を関連付ける工程
を包含する、方法。
【請求項１３】
請求項１２に記載の方法であって、生体分子成分のタイプまたは別の生体成分のタイプについての前記複数のデータセットが、１つより多い生物学的サンプルのタイプ、１つより多い測定技術のタイプ、または、生物学的サンプルのタイプと測定技術との組合せから算出される測定値を含む、方法。
【請求項１４】
請求項１２に記載の方法であって、前記生体分子成分のタイプは、タンパク質であり、そして、他の生体分子成分のタイプは、代謝産物である、方法。
【請求項１５】
請求項１２に記載の方法であって、前記生体分子成分のタイプは、遺伝子転写物であり、そして、他の生体分子成分のタイプは、代謝産物である、方法。
【請求項１６】
哺乳動物における生物システムの状態をプロファイリングする方法であって、該方法は、
（ａ）少なくとも２つの生体分子成分のタイプからの測定値を含む複数のデータセットを統計分析で評価し、そして、該複数のデータセットの少なくとも１部分の間の差異を１セット以上決定するために、該複数のデータセット間で特徴を比較する工程；および
（ｂ）工程（ａ）の結果に基づいて該生物システムの状態についてのプロファイルを作製する工程
を包含する、方法。
【請求項１７】
請求項１６に記載の方法であって、前記複数のデータセットは、１つより多い生物学的サンプルのタイプ、１つより多い測定技術、または生物学的サンプルのタイプと測定技術との組合せから算出される測定値を含む、方法。
【請求項１８】
請求項１６に記載の方法であって、前記評価する工程は、
生体分子成分のタイプについての複数のデータセットを評価し、そして、該複数のデータセットの少なくとも１部分の間の差異を１セット以上決定するために、該複数のデータセット間で特徴を比較する工程；ならびに
別の生体分子成分のタイプについての複数のデータセットを評価し、そして、該複数のデータセットの少なくとも１部分の間の差異を１セット以上決定するために、該複数のデータセット間で特徴を比較する工程
を包含する、方法。
【請求項１９】
請求項１６に記載の方法であって、前記少なくとも２つの生体分子成分のタイプは、タンパク質および代謝産物を含む、方法。
【請求項２０】
請求項１６に記載の方法であって、前記少なくとも２つの生体分子成分のタイプは、遺伝子転写産物および代謝産物を含む、方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本出願は、米国仮特許出願第６０／４９６，６５７号（２００３年８月２０日出願）の利益およびそれに対する優先権を主張する。この開示の全体は、本明細書中において参考として援用される。
【０００２】
（発明の分野）
本発明は、データの処理および評価の分野に関する。より具代的には、本発明は、生物システム（例えば、哺乳動物（例えば、ヒト））の状態をプロファイリングするための方法およびシステムに関する。
【背景技術】
【０００３】
（背景）
生物学を理解するための現在のアプローチ（例えば、ゲノミクスおよびプロテオミクス）は、代表的に、常に生物システムの単一の局面に集中する。「オミクス」の技術革命、特にゲノミクスの技術革命は、単細胞生物（例えば、酵母）および単純な多細胞のシステム（例えば、ウニ胚）の両方において、単一のタイプの生体分子の研究のための土台を提供した。両方のタイプの研究において、上記システムは、環境の変化によってか、および／または、数多くの異なるシナリオにおいて遺伝子発現の変化と相関することを可能にする遺伝子操作によって、摂動を受ける。インシリコ相互作用ネットワークの構築は、いくつかの異なる観点から、遺伝子間の相互依存性を考察することにより、容易にされる。しかしながら、現代の定量的なゲノミック技術は容易に利用可能ではあるが、結果として得られた情報は、精度および有用性が低くあり得る。例えば、１つのウニ胚の研究において、摂動は、それが遺伝子発現において３倍以上の大きな変化を生じる場合のみ有意であると判断された。多くの実験的要因が、システムの最終的な変動性に寄与し得、そして精度を損ない得るが、有意な生物学的効果は、３倍での切捨てのもとでよく生じる変化によって明示され得る。
【０００４】
複雑な多細胞生物体（例えば、哺乳動物）を分析し、かつ理解することは、はるかにずっと複雑である。複雑な生物システムの状態を研究する場合、独特の遺伝子発現ならびにタンパク質レベルおよび代謝産物レベルを有する細胞および組織のタイプの多様性はもちろんのこと、複数のコンパートメントに分かれたシステムの特徴にも気を配られなければならない。生物システムの単一の局面（例えば、分子または単一のタイプの標的）の分析に依存する現在の研究は、通常、特定の分子経路または特定の疾患に関与し得る生物システム全体または生物サブシステムを理解する程、十分に強いものではない。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
哺乳動物の生物システムを理解し、かつ複雑な多因性疾患のための新薬を開発するという重要な挑戦とは、生物マーカー／代理マーカーの同定および評価である。さらに、生物システムの状態を示している単一の生物マーカーの代わりに、生物マーカーのパターンまたは生物マーカーのセットが、生物システムについての恒常性または疾患状態を特徴付け、かつ診断することが必要であり得るようである。ここで、生物システムの複数のレベルは、上記分析において同時に考慮される。従って、生物システムを全体として考慮し、そして、ヒトの疾患の研究、ならびに治療薬の発見および開発を前進させることを可能とする方法およびシステムが必要である。
【課題を解決するための手段】
【０００６】
（発明の要旨）
本特許出願の出願人は、「システム生物学」として公知の分野における開拓者である。生物システムの個々の局面の分析とは対照的に、システム生物学とは、統合した生物システム（遺伝子成分、タンパク質成分、および代謝成分、ならびに、流動的かつ相互依存的なそれらの経路を含む）としての生物学の研究である。複雑な生物の生物学の土台となる生物学的プロセス（例えば、ヒトの疾患に関与する生物学的プロセスか、または薬物応答を支配する生物学的プロセス）の固有の複雑さを人工的に単純化するのではなく、本明細書中に記載される方法およびシステムは、生物システム内に含まれる複雑さおよび相互依存性を包括的にとらえる。生物システムの複雑さを適切に可視化し、かつ考慮することによって、当業者は、全体としての生物システムへの洞察を提供する、生物システムの状態についてのプロファイルを作製して、システムレベルで、生物学的研究を行い得る。
【０００７】
本出願は、哺乳動物（ヒトを含む）の複雑な臨床サンプルを、生物システムレベルで分析する方法およびシステムを記載する。この方法およびシステムは、伝統的な化学またはゲノミクス単独ではこれまでに入手不可能であった生物システムの状態に関する新たな情報を提供する。本明細書中に記載される方法およびシステムを使用して、疾患および薬物応答に関する生物学的経路および機構への洞察を得ることが可能である。より具体的には、上記方法およびシステムは、生体分子成分のタイプのレベル（すなわち、遺伝子／遺伝子転写物レベル、タンパク質レベル、および代謝産物レベル）で、データを分析し得、かつそのデータを統合し得、健康および疾患に関する分子機構への新たな洞察を提供することで、薬学的な研究開発を進歩させる知識を作製し、そして、さらに、ヒトの疾患を処置するための新規治療を開発および発見する。
【０００８】
生物システムの状態（例えば、疾患状態）のプロファイルを作製するために、複雑な生物学的サンプルに関する複数の測定が実施される。その後、相関分析およびネットワークモデリングと組み合わせた、遺伝子、遺伝子転写物、タンパク質、および／または代謝産物についての包括的なプロファイリングが、システムレベルでの生物システムへの洞察を提供する。その結果、何千もの多様で測定可能な分子成分との間の関連、相関、および関係が、得られ得る。その後、このような知識が、治療剤または生物マーカーの開発のために直接的に使用され得、臨床情報と合わせて使用され得、そして／あるいは、病態生理学的機構をさらに解明するために設計された直接仮説主導タイプ実験のための土台として役立ち得る。さらに、生物システムのプロファイルの変化を追跡することは、薬学的な発見および開発の多くの局面（薬物安定性および薬物効力、薬物応答、ならびに疾患の病因を含む）を改善し得る。
【０００９】
本出願は、複数のデータセット（２つ以上の生体分子成分のタイプを含み得るデータセット）を統合する能力を有する、方法およびシステム、または「技術プラットフォーム」を提供することにより、現在のプロファイリング技術における限界に取り組み、成分間の関連または成分間の相互作用ネットワークについての情報を解明する。上記方法およびシステムは、複数のデータセット（例えば、分光測定データ）の統計学的分析に利用し、生物システム（例えば、哺乳動物（例えば、ヒト））の状態に関するプロファイルを作製する。上記データセットは、生物システムの複数の測定値を含み、かつ３つの主な供給源（生物学的サンプルのタイプ、測定技術、および生体分子成分のタイプ）から算出される。本出願はさらに、サンプルまたは生物システム内の単一の生体分子成分のタイプ内だけでなく、２つ以上の生体分子成分のタイプ間の類似性、差異、および／または相関の識別を容易にする技術プラットフォームを記載する。
【００１０】
広範な局面において、生物システムの状態のプロファイリング方法は、生物システムの複数のデータセットを統計分析で評価し、そして、この複数のデータセットの少なくとも１部分の間の差異を１セット以上決定するために、この複数のデータセット間で特徴を比較する工程を包含する。上記複数のデータセット間で特徴を比較する行為は、別のデータセット中の対応する特徴と、第一のデータセット中の特徴の１つを直接比較する工程を包含し得る。上記特徴を比較する行為はまた、データセット間で、特徴を相関させるかまたは特徴を関連付ける工程を包含し得、これは、例えば、多変量解析のような統計分析と関連するか、および／または、その統計分析から得られる相関である。上記評価および比較の結果に基づき、生物システムの状態に関するプロファイルが作製され得る。
【００１１】
哺乳動物のおける生物システムの状態をプロファイルする別の方法は、生体分子成分のタイプについての複数のデータセットを統計分析で評価し、そして、この複数のデータセットの少なくとも１部分の間の差異を１セット以上決定するために、この複数のデータセット間で特徴を比較する工程；別の生体分子成分のタイプについての複数のデータセットを統計分析で評価し、そして、この複数のデータセットの少なくとも１部分の間の差異を１セット以上決定するために、この複数のデータセット間で特徴を比較する工程；ならびに上記生物システムの状態のついてのプロファイルを作製するために、上記分析の結果を関連付ける工程を包含する。
【００１２】
哺乳動物における生物システムの状態をプロファイリングするさらなる方法は、少なくとも２つの生体分子成分のタイプからの測定値を含む複数のデータセットを統計分析で評価し、そして、この複数のデータセットの少なくとも１部分の間の差異を１セット以上決定するために、この複数のデータセット間で特徴を比較する工程；ならびに上記分析の結果に基づいて生物システムの状態についてのプロファイルを作製する工程を包含する。
【００１３】
本明細書中に記載される方法およびシステムの中核を成すものは、複数のデータセットの分析である。この複数のデータセットは、１つより多い生物学的サンプルのタイプ、１つより多い測定技術のタイプ、１つより多い生体分子成分のタイプ、または、生物学的サンプルのタイプと、測定技術と、生体分子成分のタイプとの組合せのうちの少なくとも２つから算出される測定値を含む。生物システムは、好ましくは、哺乳動物（例えば、ヒト）におけるシステムである。生体分子成分のタイプとしては、タンパク質、糖タンパク質、遺伝子、遺伝子転写物、および代謝産物が挙げられる。
【００１４】
生物学的サンプルのタイプとしては、特に、血液、血漿、血清、脳脊髄液、胆汁、唾液、滑液、胸膜液、心膜液（ｐｅｒｉｃａｒｄｉａｌ  ｆｌｕｉｄ）、腹水、汗、糞便、鼻液、眼液、細胞内液、細胞間液、リンパ液、尿、肝細胞、上皮細胞、内皮細胞、腎臓細胞、前立腺細胞、血液細胞、肺細胞、脳細胞、皮膚細胞、脂肪細胞、腫瘍細胞、および乳房細胞（ｍａｍｍａｒｙ  ｃｅｌｌ）が挙げられる。データセットは、別々に処理される１つの生物学的サンプルのタイプからか、あるいは、別々の時間で収集されるかまたは分析される１つの生物学的サンプルのタイプからの測定値を含み得る。
【００１５】
測定技術としては、特に、液体クロマトグラフィー、ガスクロマトグラフィー、高速液体クロマトグラフィー、キャピラリー電気泳動、質量分析、液体クロマトグラフィー質量分析、ガスクロマトグラフィー質量分析、高速液体クロマトグラフィー質量分析、キャピラリー電気泳動質量分析、核磁気共鳴分析、並行ハイブリダイゼーションアッセイ（ｐａｒａｌｌｅｌ  ｈｙｂｒｉｄｉｚａｔｉｏｎ  ａｓｓａｙ）、並行サンドイッチアッセイ（ｐａｒａｌｌｅｌ  ｓａｎｄｗｉｃｈ  ａｓｓａｙ）、および、競合アッセイが挙げられる。データセットは、単一タイプの測定技術の異なる機器構成からの測定値を含み得る。
【００１６】
生物システムの状態についてのプロファイルを作製した後に、このプロファイルは、生物システムの別の状態のプロファイルと比較され得る。この生物システムは、同じシステムかまたは違うシステムである。プロファイルはまた、この生物システムの状態が、公知の状態と一致するかまたは類似しているかどうかを評価するために、プロファイルのデータベースと比較され得る。本明細書中に記載される方法は、上記方法を実施するために、コンピューターが読み取り可能な命令を具体化されたコンピューターが読み取り可能な媒体を有する製品によって、実行され得る。
【００１７】
本発明の他の局面および利点は、以下の図面、以下の詳細な説明、および上記特許請求の範囲から明確となり、これらの全ては、ほんの一例として本発明の原理を例示する。
【発明を実施するための最良の形態】
【００１８】
本明細書中で開示される方法およびシステムは、生物学的サンプルの複数の測定値（代謝産物、タンパク質、遺伝子および転写遺伝子物の分析を含む）に依存し、これらの因子のうちの１つだけを検査するアプローチよりも、当業者が、生物システムをより深く理解することを可能にする。生物システムを全体として理解することは、薬学的な発見および開発の複数の局面（薬物安全性および薬物効力、薬物応答、ならびに疾患の病因）を改善し得る。本明細書に記載されるように、システム生物学のプラットフォームは、ゲノミクス、プロテオミクスおよびメタボロミクス、ならびにバイオインフォマティクスを統合し得、そして、測定可能な何千もの分子成分間の関連、相関、および関係を生成して生物システムの状態のプロファイルを作製する、データ統合プラットフォームおよび知識管理プラットフォームをもたらす。結果として作製されるプロファイルは、生物システムの状態に関する知識を増大するために、臨床情報と合せられ得る。
【００１９】
生物システムの「プロファイル」とは、生物システム（例えば、ヒトのような哺乳動物）の特有な特徴または特性を表すデータの概要または分析である。このデータは、生物学的サンプルのタイプ、測定技術のタイプ、および生体分子成分のタイプから算出される測定値、またはそれらのタイプの特徴を含み得る。多くの場合、上記データは、グラフ、表、またはこれらに多少類似するデータ編集形態における、スペクトルの特徴かまたはクロマトグラムの特徴である。代表的に、プロファイルとは、生物システムの状態の特徴付けを可能にする特徴のデータセットである。
【００２０】
プロファイルは、生物システムの１つ以上の「生物マーカー」を含むことを考慮され得る。一般に、生物マーカーとは、生物学的成分のタイプ（例えば、遺伝子、遺伝子転写物、タンパク質、または代謝産物）を呼ぶ。生物システムにおける、生物マーカーの定性的および／または定量的な存在または欠如は、哺乳動物の生物学的状態の指標である。従って、プロファイルは、生物システムの状態の特徴付けを可能にする、特有の生物マーカー（スペクトルの特徴またはクロマトグラムの特徴）のセットとして考慮され得る。プロファイルはまた、相関および上記データセット分析の他の結果（例えば、因果関係）を含むことが考慮され得る。従って、プロファイルは、上記の複数の異なる要素を含むか、またはこれらの要素のうち１つ（例えば、生物マーカー）のみを含み得る。
【００２１】
「生物システムの状態」とは、生物システムが自然に存在する状態か、または摂動を与えた後の状態に依存するうちのいずれかを呼ぶ。生物システムの状態の例としては、標準状態または健康状態、疾患状態、薬理学的因子応答、毒性学的状態、生化学的調節（例えば、アポトーシス）、年齢応答、環境応答、およびストレス応答が挙げられるが、限定はされない。この生物システムは、好ましくは、ヒトおよび非ヒト哺乳動物（例えば、マウス、ラット、モルモット、イヌ、ネコ、サルなどを含む）を含む哺乳動物におけるシステムである。
【００２２】
生物システムの状態のプロファイルは、１つのプロファイルと別のプロファイルとを比較し、そのプロファイルが同じ状態（例えば、健康状態または罹患状態）にあるかどうかを決定することを可能にする。生物システムは、同じ変数の複数の測定値を使用するよりも、多変量分析を使用するほうが、より良く特徴付けられる。なぜなら、多変量解析は、生物システムを全体として想定するからである。複数の異なる供給源からの不同性データは、多次元ではなく、一次元であるかのように、処理される。結果として、データの分析は、より情報的であり、そして代表的には、複数の成分を別個に、体系的に評価することにより作製されたか、または特定の生体分子成分の１つのタイプに依存する、プロファイルよりも、より頑強かつ予測的なプロファイルを提供する。
【００２３】
「生体分子成分のタイプ」とは、生物システムのレベルと一般的に関連する生体分子の分類を指す。例えば、遺伝子および遺伝子転写物（本明細書中において交換可能に言及され得る）は、生物システム中の遺伝子発現と一般的に関連する生体分子成分のタイプの例であり、ここでの生物システムのレベルとは、ゲノミクスまたは機能的ゲノミクスとして呼ばれる。タンパク質およびそれらの構成ペプチド（本明細書中において交換可能に言及され得る）は、タンパク質の発現および修飾と一般的に関連する生体分子成分のタイプの別の例であり、ここでの生物システムのレベルは、プロテオミクスとして呼ばれる。糖タンパク質もまた、生体分子成分のタイプであると考慮される。生体分子成分のタイプの別の例は、代謝産物（低分子としてもまた言及され得る）であり、これは、メタボロミクスとして呼ばれる生物システムのレベルと一般的に関連する。代謝産物としては、脂質、ステロイド、アミノ酸、有機酸、胆汁酸、エイコサノイド、神経ペプチド、ビタミン、神経伝達物質、炭水化物、イオン性有機物、ヌクレオチド、無機物、生体異物、ペプチド、微量元素、ならびにファーマコフォアおよび薬物分解産物が挙げられるが、限定はされない。
【００２４】
本明細書中に記載される方法は、任意の単一の生体分子成分のタイプに基づいて、ならびに２つ以上の生体分子成分のタイプに基づいて、生物システムの状態のプロファイルを作製するために使用され得る。生体分子成分のタイプのプロファイルは、生物システムの異なるレベルでの包括的なプロファイル（例えば、ゲノムプロファイル、トランススクリプトームプロファイル、プロテオームプロファイル、およびメタボロームプロファイル）の作製を容易にし、それらの統合および分析を可能にする。すなわち、上記方法は、１つ以上の生物学的サンプルのタイプ、１つ以上の測定技術のタイプ、または各生物学的サンプルのタイプのうちの少なくとも１つと測定技術との組合せから算出される測定値を分析するために使用され得る。その結果、上記方法は、単一の生体分子成分のタイプかまたは２つ以上の生体分子成分のタイプ間の類似性、差異、および／または相関の評価を可能にする。これらの測定値から、根底にある生物学的機構のより良い洞察が得られ得、新規バイオマーカー／代理マーカーが検出され得、そして、介入経路（ｉｎｔｅｒｖｅｎｔｉｏｎ  ｒｏｕｔｅ）が開発され得る。
【００２５】
「生物学的サンプルのタイプ」としては、血液、血漿、血清、脳脊髄液、胆汁酸、唾液、滑液、胸膜液、心膜液（ｐｅｒｉｃａｒｄｉａｌ  ｆｌｕｉｄ）、腹水、汗、糞便、鼻液、眼液、細胞内液、細胞間液、リンパ液、尿、組織、肝細胞、上皮細胞、内皮細胞、腎臓細胞、前立腺細胞、血液細胞、肺細胞、脳細胞、脂肪細胞、腫瘍細胞、および乳房細胞（ｍａｍｍａｒｙ  ｃｅｌｌ）が挙げられるが、限定はされない。生物学的サンプルのタイプの供給源は、異なる被験体；異なる時間における同じ被験体；異なる状態（例えば、薬物処置前、および薬物処置後）における同じ被験体；異なる性別；異なる種（例えば、ヒトおよび非ヒト哺乳動物）；ならびに種々の他の順列である。さらに、生物学的サンプルは、例えば、異なる精密検査プロトコルを使用する評価の前に、別々に処理され得る。
【００２６】
「測定技術」とは、生物システムの状態の分析において有用なデータを生成するかまたは提供する、任意の分析技術を指す。例えば、測定技術としては、質量分析（「ＭＳ」）、核磁気共鳴分析（「ＮＭＲ」）、液体クロマトグラフィー（「ＬＣ」）、ガスクロマトグラフィー（「ＧＣ」）、高速液体クロマトグラフィー（「ＨＰＬＣ」）、キャピラリー電気泳動（「ＣＥ」）、ゲル電気泳動（「ＧＥ」）、および、低分解能様式かまたは高分解様式のハイフンでつないだ質量分析の任意の公知な形態（例えば、ＬＣ／ＭＳ、ＧＣ／ＭＳ、ＣＥ／ＭＳ、ＭＳ／ＭＳ、ＭＳｎ、および他の改変形）が挙げられるが、限定はされない。測定技術としては、生物学的画像化（例えば、磁気共鳴画像法（「ＭＲＩ」）、映像信号、および蛍光アレイ（例えば、空間内の点からの光強度、および／または色）、ならびに他のハイスループットデータ収集技術かまたは高度な並行データ収集技術が挙げられる。
【００２７】
測定技術としてはまた、光学分光学、デジタル画像、オリゴヌクレオチドアレイハイブリダイゼーション、タンパク質アレイハイブリダイゼーション、ＤＮＡハイブリダイゼーションアレイ（「遺伝子」チップ）、免疫組織化学的分析、ポリメラーゼ連鎖反応、核酸ハイブリダイゼーション、心電図検査、コンピューター横断断層撮影、ポジトロン放出断層撮影、およびテキストベースの臨床的なデータ報告に見られるような主観的な分析が挙げられる。特定の分析に関して、異なる測定技術とは、同じ測定技術に関する、異なる機器構成かまたは異なる設定を含み得る。
【００２８】
「測定値」とは、測定技術によって得られたデータセットの要素を指す。「データセット」とは、１つ以上の供給源から算出された測定値を含む。例えば、測定技術から算出されたデータセットは、同じ技術（すなわち、関連する測定値の収集またはデータセット）によって収集された一連の測定値を含む。さらに、データセットは、多様なデータ（例えば、タンパク質発現データ、遺伝子発現データ、代謝物質濃度データ、核磁気共鳴画像データ、心電図データ、遺伝子型データ、一塩基多型データ、および他の生物学的データ）の収集を、より広範に表し得る。すなわち、研究されている生物システムに関する測定可能または数値化可能な任意の局面は、所定のデータセットを生成するための土台として役立ち得る。
【００２９】
データセットの「特徴」とは、別のデータセットと比較され得るデータセットと関連する特定の測定値を指す。例えば、プロファイルは、代表的には、生物システムの状態の特徴付けを可能にする特徴を有するデータセットである。
【００３０】
データセットとは、１つ以上の測定技術と関連するデータの実質的に全てのセットかまたは部分セットを指し得る。例えば、異なるサンプル供給源の分光測定の測定値と関連するデータは、異なるデータセットへグループ化され得る。結果として、第１のデータセットは、実験グループのサンプルの測定値を指し得、第２のデータセットは、コントロールグループのサンプルの測定値を指し得る。加えて、データセットは、他の関連すると考えられる分類に基づいてグループ化したデータを指し得る。例えば、単一のサンプル供給源の分光測定の測定値と関連するデータは、測定を実施するために使用される機器に基づいて、異なるデータセット（サンプルを採取した時間、サンプルの外観、または他の同定可能な変数および特性）へグループ化され得る。
【００３１】
従って、１つのデータセットは、別のデータセットの部分セットを含み得る。例えば、サンプルの外観に基づくグループ化は、実験グループのデータセットを１つ以上含み得る。測定技術がＮＭＲである場合、データセットは、１つ以上のＮＭＲスペクトルを含み得る。測定技術が紫外線（ＵＶ）分光法である場合、データセットは、ＵＶ放出スペクトルまたはＵＶ吸収スペクトルを１つ以上含み得る。同様に、測定技術がＭＳである場合、データセットは、１つ以上の質量スペクトルを含み得る。測定技術がＬＣ／ＭＳまたはＧＣ／ＭＳのようなクロマトグラフ−ＭＳ技術である場合、データセットは、１つ以上のクロマトグラムを含み得る。あるいは、クロマトグラフ−ＭＳ技術のデータセットは、全イオン電流（「ＴＩＣ」）クロマトグラムまたは再配列したＴＩＣクロマトグラムを１つ以上含み得る。加えて、用語「データセット」が、生の分光測定データと前処理されたデータの両方を含むことが、理解されるべきである。前処理とは、例えば、ノイズを除去すること、基線を補正すること、データを平滑化すること、ピークを検出すること、および／またはデータを標準化することである。
【００３２】
「分光測定データ」とは、グラフ、表、ベクトル、配列、または多少類似するデータ編集の形態で表され得る、任意のデータを指す。そして、分光測定データは、任意の分光測定技術またはクロマトグラフ技術からのデータを含み得る。用語「分光測定の測定値」は、任意の分光測定技術かまたはクロマトグラフ技術によって生成された測定値を含む。
【００３３】
本明細書中に開示される方法の中核を成すものは、複数のデータセットの統計分析である。「統計分析」としては、パラメトリック分析、ノンパラメトリック分析、単変量分析、多変量分析、線形分析、非線形分析、および当業者に公知の他の統計学的方法が挙げられる。多変量分析は、見かけ上は無秩序なデータにおいてパターンを決定する分析であり、この分析としては、主成分分析（「ＰＣＡ」）、判別分析（「ＤＡ」）、ＰＣＡ−ＤＡ、正準相関（「ＣＣ」）、クラスター分析、部分最小二乗法（「ＰＬＳ」）、予測的線形判別分析（「ＰＬＤＡ」）、ニューラルネットワーク、およびパターン認識技術が挙げられるが、限定はされない。
【００３４】
多変量分析を行う前には当然、生データが、別のデータセットとの比較に役に立つように前処理され得る。特に、異なる生体分析成分のタイプを横断してデータを比較するためには、適切な前処理が行われるべきである。データの前処理は、
（ｉ）異なるサンプルのスペクトルのピークを揃えるために、データセット間のデータポイントを（例えば、部分的な線形フィット技術を使用して）揃えること；
（ｉｉ）ピークの高さを調整するために、（例えば、各測定において標準物質を使用して）データセットのデータを標準化すること；
（ｉｉｉ）ノイズを減らし、そして／またはピークを検出すること（例えば、潜在的な基線ノイズから実際の種の存在を識別するために、ピークの閾値を設定すること）；ならびに／あるいは
（ｉｖ）当該分野において公知の他のデータ処理技術
を包含し得る。データ処理は、米国特許第６，７４３，３６４号において開示されるようなエントロピーベースのピーク検出、および部分的な線形フィット技術（例えば、Ｊ．Ｔ．Ｗ．Ｅ．Ｖｏｇｅｌｓら、「Ｐａｒｔｉａｌ  Ｌｉｎｅａｒ  Ｆｉｔ：Ａ  Ｎｅｗ  ＮＭＲ  Ｓｐｅｃｔｒｏｓｃｏｐｙ  Ｐｒｏｃｅｓｓｉｎｇ  Ｔｏｏｌ  ｆｏｒ  Ｐａｔｔｅｒｎ  Ｒｅｃｏｇｎｉｔｉｏｎ  Ａｐｐｌｉｃａｔｉｏｎｓ」、Ｊｏｕｒｎａｌ  ｏｆ  Ｃｈｅｍｏｍｅｔｒｉｃｓ、第１０巻、１９９６年、ｐｐ．４２５〜３８において見られる）を含み得る。
【００３５】
組成物が特定の成分を有するか、含むか、または包含するものとして記載されるか、あるいは、プロセスが特定のプロセス工程を有するか、含むか、または包含すると記載される、本明細書全体を通して、本発明の成分がまた、列挙した成分から、本質的になるかまたは単に構成されること、ならびに、本発明のプロセスがまた、列挙したプロセス工程から、本質的になるかまたは単に構成されることが、企図される。
【００３６】
工程の順番または特定の行為を行う順番は、本発明が実施可能である（すなわち、生物システムのプロファイルが作製される）限りは重要でないことが理解されるべきである。さらに、２つ以上の工程または行為は、同時に行われ得る。
【００３７】
本明細書中に記載される方法は、一般に、生物システムの複数のデータセットを統計分析で評価し、そして、そのデータセット間の特徴を比較し、比較に基づいて生物システムの状態についてのプロファイルを作製するために、データセットの少なくとも１部分間の差異を１つ以上決定する工程を包含する。いくつかの実施形態において、上記データセットは、１以上の生物学的サンプルのタイプから算出され、そして、このデータセットは、１つ以上の測定技術から算出される測定値を含む。他の実施形態において、上記データセットは、２つ以上の生物学的サンプルのタイプから算出され、そして、このデータセットは、生物システムのサンプルの分光測定における１つ以上の異なるタイプを含む。
【００３８】
特定の実施形態において、データセットは、多変量分析を使用して前処理および評価される。他の実施形態において、１つより多くの統計分析は、複数のデータセットに関して、複数のデータセットの種々の順列に関して、および／または特定の統計分析の結果に関して実施される。例えば、プロファイルは、生物システム中のタンパク質から算出された測定値を含む複数のデータセット、および生物システム中の代謝産物から算出された測定値を含む複数のデータセットを、別々に評価することによって作製され、その後、タンパク質と代謝産物の両方を含む生物システムについてのプロファイルを作製するために、個々の分析の結果を統計分析で評価し得る。あるいは、生物システムのタンパク質および代謝産物に関する上記複数のデータセットは、統計分析で同時に評価され得る。
【００３９】
同じように、プロファイルは、タンパク質および遺伝子；タンパク質および遺伝子転写物；遺伝子および遺伝子転写物；遺伝子および代謝産物；ならびに遺伝子転写物および代謝産物から算出される測定値を含むデータセットから作製され得る。プロファイルはまた、タンパク質、遺伝子、および遺伝子転写物；タンパク質、遺伝子、および代謝産物；タンパク質、遺伝子転写物、および代謝産物；ならびに遺伝子、遺伝子転写物、および代謝産物；ならびにタンパク質、遺伝子、遺伝子転写物、および代謝産物から算出される測定値を含むデータセットから作製され得る。加えて、上記順列の各々は、追加または代用として糖タンパク質を含み得る。
【００４０】
特定の生体分子成分のタイプについての測定値は、通常、特定の生体分子成分のタイプについて当該分野において頻繁に使用されかつ公知な測定技術によって生成される。例えば、代謝産物の分析は、ＮＭＲ（例えば、１Ｈ−ＮＭＲ；ＬＣ／ＭＳ；ＧＣ／ＭＳ；およびＭＳ／ＭＳ）を使用し得る。他の生体分子成分のタイプの分析は、ＬＣ／ＭＳ；ＧＣ／ＭＳ；およびＭＳ／ＭＳを使用し得る。
【００４１】
１つの実施形態において、一般に、上記方法は、
生物学的サンプルを選択する工程；
調査されるべき生化学的成分、および利用されるべき分光測定技術に基づいて生物学的サンプルを調製する工程；
生物学的サンプル中の成分を、分光測定技術およびクロマトグラフ技術を使用して測定する工程；
選択された分子のサブクラスを、化合物を研究するためにＮＭＲおよびＭＳアプローチを使用して測定する工程；
生データを処理する工程；
以下により詳細が記載される統計分析を使用して、前処理されたデータを分析し、分子の単一のサブクラスの測定においてかまたは成分の測定において、ＮＭＲまたはＭＳを使用してパターンを同定する工程；ならびに
統計分析を使用して、別の実験からのデータセットを組み合わせ、そしてこのデータ中の対象となるパターンを同定する工程
を包含する。
【００４２】
この技術のプラットフォームはまた、生体分子成分のタイプ間のデータ比較を容易にするために、複数のデータセットを標準化することを含み得る。本発明はまた、適切なデータセットの生体分子成分のタイプ間の関連性／相関を、線形、非線形、または他の数学的ツールを使用して決定する技術を提供する。さらに、生体分子成分の相互作用のネットワークを予測するためにこれらの関連性および／または相関を使用し、これらの関係の間の因果関係を決定し、かつ、データセットを生じる観察の根底にある生物学的プロセスについての仮説を証明することは、本明細書中に記載される方法およびシステムのなお別の局面である。
【００４３】
本出願はまた、本明細書中に開示される方法の機能性が、コンピューターが読み取り可能な媒体（例えば、フロッピー（登録商標）ディスク、ハードディスク、光ディスク、磁気テープ、ＰＲＯＭ、ＥＰＲＯＭ、ＣＤ−ＲＯＭ、またはＤＶＤ−ＲＯＭが挙げられるが、限定はされない）に組み込まれた製品を提供する。上記方法の機能性は、かなり多数のコンピューターが読み取り可能な命令または言語（例えば、ＦＯＲＴＲＡＮ、ＰＡＳＣＡＬ、Ｃ、Ｃ＋＋、ＢＡＳＩＣ、およびアセンブリ言語）で、コンピューターに読み取り可能な媒体に組み込まれ得る。さらに、上記コンピューターが読み取り可能な命令は、スクリプトでか、マクロで書かれ得るか、または市販のソフトウェア（例えば、ＥＸＣＥＬまたはＶＩＳＵＡＬ  ＢＡＳＩＣ）中に機能的に組み込まれ得る。他の局面において、本出願は、本明細書中に記載される方法を実施するために適合されたシステムを提供する。
【００４４】
データ処理装置は、本明細書中において開示される１つ以上の方法の機能性を、分光測定機器によって提供された情報の少なくとも１部を使用して、実行するのに適合したアナログ回路よび／またはデジタル回路を含み得る。いくつかの実施形態において、上記データ処理装置は、本明細書中に記載される方法の機能性を、多目的コンピューター上のソフトウェアとして実行し得る。加えて、このようなプログラムは、分光測定の測定値の取得、データセットの統計分析、および／または、生物システムについてのプロファイルの作製に影響する制御論理を提供するために、コンピューターのランダムアクセスメモリ部分を取っておくことができる。このような実施形態において、上記プログラムは、多数の高級言語（例えば、ＦＯＲＴＲＡＮ、ＰＡＳＣＡＬ、Ｃ、Ｃ＋＋、またはＢＡＳＩＣ）のうちのいずれか１つで書かれ得る。さらに、上記プログラムは、スクリプトでか、マクロで書かれ得、あるちは、登録商標権を有するソフトウェアまたは市販のソフトウェア（例えば、ＥＸＣＥＬまたはＶＩＳＵＡＬ  ＢＡＳＩＣ）中に機能的に組み込まれ得る。その上、上記ソフトウェアは、コンピューター上で、常駐マイクロプロセッサを対象とするアセンブリ言語で実行され得る。例えば、上記ソフトウェアは、このソフトウェアがＩＢＭ  ＰＣかまたはＩＢＭ  ＰＣクローン上で実行されるように設定される場合、Ｉｎｔｅｌ  ８０×８６のアセンブリ言語で実行され得る。上記ソフトウェアは、コンピューターが読み取り可能な媒体（例えば、フロッピー（登録商標）ディスク、ハードディスク、光ディスク、磁気テープ、ＰＲＯＭ、ＥＰＲＯＭ、またはＣＤ−ＲＯＭ）を含むが、限定はされない、製品に組み込まれ得る。
【００４５】
図１に示されるように、いくつかの実施形態において、上記方法は、罹患集団および健康な集団の両方から抽出された複合サンプルから算出される遺伝子転写産物（ｍＲＮＡ）、タンパク質、および代謝産物の定量的プロファイルを並行分析することから始まる。測定された全ての化合物についての、平均量は、分散および範囲と同様に、遺伝子応答、タンパク質活性、および代謝産物のダイナミクスに関連する分子を同定するためにパターンを認識するような方法を使用して、まとめて分析される。本明細書中において開示される方法は、ＢｉｏＳｙｓｔｅｍａｔｉｃｓＴＭと名付けられたものであり、次いで、これは、共変する遺伝子のセット（遺伝子転写物、タンパク質、および代謝産物を含む）を翻訳して、必要に応じて臨床情報とともに、生物システムのプロファイル、および標的の情報を解明するためのそれらの生化学的相互作用を理解するために利用され得る。この情報は、共変する分子の特定のグループから、既存の経路の知識にまでわたり、分子ネットワークを構築し、かつそれらの生物学的背景を踏まえて化合物を配置して、生物システムの状態のプロファイルを作製するために使用される。
【００４６】
図２は、分析方法（２００）の１つの実施形態についてのフローチャートを示す。以下の記載される１つ以上の工程が、省略され得、そして／または、工程の順番は、この実施形態が実施可能である（すなわち、生物システムの状態のプロファイルを作製し得る）限りは変更され得ることが、理解されるべきである。２つ以上の生体分子成分のタイプから得られた１つ以上のデータセット（２０５）は、さらなるデータ分析の前の最初の前処理工程（２１０）に供される。好ましい実施形態において、上記最初の処理工程は、代表的に、複数のデータセットの１つ以上を連結することを含む。この最初の前処理工程はまた、適切なスキーマまたは適切なデータ階層に基づいて、データセットを一緒に統合することを含み得る。いくつかの実施形態において、上記最初の処理工程は、連結工程と統合工程の両方を含む。上記最初の処理工程は、必要に応じて、種々の形態の前処理工程（データの平滑化、ノイズの減少、基線の補正、およびピークの検出を含むが、限定はされない）を、含んでも、その工程の後に続けても、その工程の前に先行させてもよい。
【００４７】
上記最初の前処理工程の対象であるデータセットは、研究される生物システムの測定可能かまたは定量可能な局面のうちの任意のものを含み得る。例えば、上記データセットは、例えば、タンパク質発現データ、遺伝子発現データ、代謝産物濃度データ、核磁気共鳴画像データ、心電図データ、遺伝子型データ、および／または一塩基多型データの収集物を表し得る。主成分分析のような統計的手法は、データセットをスペクトル因子（これらは、単に生データの処理された形態である）に変換するために利用され得る。
【００４８】
特に、広範囲にわたるデータセットが利用され得ると仮定すると、完全に無関係な現象についてのデータセットと異なる測定単位とを比較する手段が、必要となる。図２を参照すると、以下でより詳細に記載される標準化工程（２１５）が、このような異なるデータセットについて実行され得る。一般に、個々のデータセットは、最大尤度推定器を使用して計算された最適なスケール化パラメーターを用いてデータセットをスケール化することによって、標準化される。標準化は、１つ以上の生体分子成分のタイプから得られるデータセットの比較を容易にする。
【００４９】
抽出工程（２２０）は、代表的に、処理されたデータに関して行われる。抽出工程において、統計学的に有意な変化を示す成分の１つ以上のリストが、抽出される。上記成分は、代表的に、生物学的成分のタイプか、または、より具体的には、生体分子成分のタイプである。さらに、これらの変化はまた、抽出工程の１部分として数値化される。上記抽出工程は、代表的に、データセット間の差異、および／または類似性を識別するための統計分析を必要とする。上記抽出工程、および関連した差異の数値化は、調査中の生物学的サンプルについて、２つ以上の生体分子成分のタイプ間の類似性、差異、および／または相関を識別することを容易にする。
【００５０】
成分タイプ間の変化を数値化するのに適切な統計分析に適した形態としては、例えば、主成分分析（「ＰＣＡ」）、判別分析（「ＤＡ」）、ＰＣＡ−ＤＡ、正準相関（「ＣＣ」）、部分最小二乗（「ＰＬＳ」）、予測線形判別分析（「ＰＬＤＡ」）、ニューラルネットワーク、およびパターン認識技術が挙げられる。１つの実施形態において、ＰＣＡ−ＤＡは、スコアプロット（すなわち、２つの主成分に関するデータのプロット）を生成する第一の相関レベルで、行われる。その後、同じかまたは異なる統計分析が、このデータセットに関して、前の分析から識別された差異、および／または類似性に基づいて、行われる。
【００５１】
例えば、処理されたデータセットが、ＰＣＡ−ＤＡスコアプロットを含む、１つの実施形態において、統計処理の次のレベルは、ＰＣＡ−ＤＡ分析によって生成された負荷プロットであり得る。この第二の相関レベルは、上記第一のレベルへの階層的関係を有する。なぜなら、負荷プロットが、ＰＣＡ−ＤＡに対する個々の入力ベクトルの寄与の情報を提供し、次いでこの情報が、スコアプロットを生成するために使用されるからである。例えば、各データセットが、複数の質量クロマトグラムを含む場合、スコアプロットの点は、１つのサンプル源から生成される質量クロマトグラムを表す。対照的に、負荷プロット上の点は、データセット間の相関に対する特定の質量または質量範囲の寄与を表す。同様に、各データセットが、複数のＮＭＲスペクトルを含む場合、スコアプロットの点は、１つのＮＭＲスペクトルを表す。対照的に、対応する負荷プロット上の点は、データセット間の相関に対する特定のＮＭＲの化学シフト値または値の範囲の寄与を表す。
【００５２】
図２はまた、抽出工程（２２０）に続く、相関ネットワーク生成工程（２２５）を図示する。相関ネットワークの公式化により、前記工程で前もって作製された、成分に関する抽出されたリスト間での、潜在的な関連性が示される。相関ネットワークは、１つ以上のサンプルグループ間で存在量が異なる、システムの生体分子成分のタイプについての（図式的、数学的、または他の）表示である。２つの成分が、いくらか同期した様式で変動する場合、この２つの成分は、「相関性がある」。例えば、グループ２と比較してグループ１において遺伝子およびタンパク質の両方が、アップレギュレートしており、このアップレギュレーションが、グループ１を含む全ての生物学的サンプルの間で一貫してアップレギュレートしている場合、この遺伝子とタンパク質とは、「相関性がある」と考慮される。同じように、生体分子成分のタイプはまた、反相関でもあり得る。さらに、異なる「相関強度」が存在し、これは、２つ以上の生体分子のタイプ間の関係がいかに密接に同期的であるかに依存する。
【００５３】
比較工程（２３０）は、相関ネットワークが構築された後に、実施される。相関ネットワークの関連性（相関および反相関の両方を含む）は、調査中の成分または生物システムに関する既存の知識に基づいて比較され、そして評価される。この知識は、確立された情報源（例えば、研究論文、および／または実験研究）から確認され得る関連性に関する。
【００５４】
続いて、摂動工程（２３５）は、代表的には、より大きな分析の１部分として実施される。調査の対象となる生物システムは、代表的には、実験パラメーターを変化させて、規定の時間にわたって、システムをモニターすることによって、摂動が与えられる。摂動の例としては、薬物を誘導するか、遺伝子を変更するか、環境条件を変化させるか、または別の適切な変化を与えることが挙げられるが、限定はされない。摂動はまた、種間で比較するというアイデア（すなわち、動物システムについてのワークフローを実施し、そして、ヒトシステムについて同じワークフローを実施して、種間の類似性および／または差異を調査すること）も含む。
【００５５】
摂動工程（２３５）の後、新たなデータセットおよび新たな相関ネットワークが、生成される（２４０）。このようにして、所与の生物システムまたは生物学的サンプルへ摂動を与えた結果、測定可能な新たなデータセットが生じる。同様に、工程（２４０）の１部分として、新たな相関ネットワークが、それらの摂動後の新規データセットに基づいて作製され得る。上記新たなデータセットにおける統計学的に有意な変化は、摂動前のデータセットと比較して決定されるので、これらは、上記新たなデータセットにおける統計学的に有意な生物学的成分のタイプと、前の実験結果の成分のタイプとを比較することにより、識別される（２４５）。システムに摂動を与える前と後で、生体分子成分のタイプ間の統計学的変化を考察すること（２４５）に加えて、相関ネットワークが、同様の方法により分析され得る。従って、上記相関ネットワークの関連性が、摂動前と摂動後で比較され得る（２５０）。これら２つのレベルの比較（２４５、２５０）が実施された後、成分と関連性との間の変更または変化が、同定され得る（２５５）。
【００５６】
その後、調査されるシステムへの摂動は、繰り返され得る（２６０）。フィードバックループは、システムへの最初の摂動、システム自体、新たなデータセットの生成、有意な成分と前の実験との比較、新たな相関ネットワークの関連性と前の関連性との比較、および変化の同定の間で生じる。フィードバックループは、複数の生体分子成分のタイプと、生物システムに対するそれらの影響を特徴付ける相関およびネットワークとの間で因果関係が同定され得る（２６５）まで、繰り返され得る。
【００５７】
再度、図２中の標準化工程（２１５）および上記の工程を参照して、遺伝子発現データ、タンパク質データ、および代謝産物レベルのデータについての標準化の方法が、次に記載される。サンプルの多様な効果、アレイの効果、および染色効果は、対数線形モデルに導入され、そして、最大尤度を極大化する技術が、上記モデルの全てのパラメーターを計算するために適用されて、各アレイおよび色素について適切な倍率を決定する。上記標準化方法は、一般的であり、種々のデータ、実験装置、および実験設計に適用され得る。以下に記載されるモデルでは、遺伝子発現分析の専門用語を使用する。例えば、プロテオミクス実験における「アレイ」とは、１つの質量スペクトルの試行であり得、そして「色素」とは、単一の試行の間に使用される全てのサンプルを表し得る。そうはいうものの、他の生体分子成分のタイプは、以下に記載されるモデルを使用して分析され得る。
【００５８】
（標準化モデル）
データマトリックス（ｘ）は、遺伝子インデックス（ｇ（ｇ＝１．．．Ｎｇ））、アレイインデックス（ｉ（ｉ＝１．．．Ｎｉ）、色素インデックス（ｋ（ｋ＝１．．．Ｎｋ））、および変種インデックス（ｖ（ｖ＝１．．．Ｎｖ））によって特徴付けられる。各々の変種（ｖ）について、それに対応するＣｖサンプルが存在するので、従って、
【００５９】
【数１】
2007502992.tif 000003
となる。変種の割り当ては、アレイインデックスおよび色素インデックスの関数なので、各データ点は、インデックスｇ、ｉ、およびｋにより一意的に記載される。便宜上、上記マトリックスは、対数に変換される：
ｙｇｉｋ＝ｌｏｇ（ｘｇｉｋ）      （１）
データは、以下のモデルにより記載される：
ｙｇｉｋ＝μｇｖ＋Ａｉ＋Ｄｋ＋εｇｉｋ      （２）
ここで、遺伝子の効果および変種の効果は、μｇｖ、アレイの効果は、Ａｉ、染色の効果はＤｋ、および誤差関数はεｇｉｋで記載される。この誤差関数は、通常、ゼロ平均および分散σ２ｇｖで分布されると仮定される。すなわち、この分散は、各遺伝子および変種によって異なることを許容している。上記変種インデックス（ｖ）は、ｉおよびｋの独特な関数であり、｛ｉ，ｋ｝∈Ｖのように書かれ得る。上記の遺伝子の効果、変種の効果、および染色の効果は、固定されると仮定されるので、発現レベルの分散は、以下：
【００６０】
【数２】
2007502992.tif 000004
のように記載され得る。
最大尤度の推定は、データを適切に標準化するために使用される最適なスケール化パラメーターを算出するために使用される。パラメーターμｇｖ、Ａｉ、Ｄｋ、およびσｇｖを解くことで、以下の式：
【００６１】
【数３】
2007502992.tif 000005
が導かれる。従って、各アレイおよび色素についての最適な倍率は、以下：
Ｓｉｋ＝−Ａｉ−Ｄｋ      （５）
となり、上記標準化された発現レベルは、以下：
【００６２】
【数４】
2007502992.tif 000006
となる。
【００６３】
（有意検定およびブートストラップ法）
上記標準化されたデータは、帰無モデルと比較され得、そして、この帰無モデルからのデータの偏差が、ランダム誤差に寄与し得る確率を測定するｐ値が、算出され得る。比較に使用されるパラメーターは、２つの選択された変種の間の倍率比である。上記方法を評価するため、ｔ検定が、上記２つの選択された変種を比較するために実施される（Ｓｈｅｓｋｉｎ，Ｈａｎｄｂｏｏｋ  ｏｆ  Ｐａｒａｍｅｔｒｉｃ  ａｎｄ  Ｎｏｎｐａｒａｍｅｔｒｉｃ  Ｐｒｏｃｅｄｕｒｅｓ，Ｃｈａｐｍａｎ  ＆  Ｈａｌｌ／ＣＰＣ，Ｂｏｃａ  Ｒａｔｏｎ，ＦＬ（２０００））。上記対応するｐ値は、各遺伝子について算出された。各遺伝子の倍率変化の統計的有意性を評価する場合、
【００６４】
【数５】
2007502992.tif 000007
でのいくつかのｐ値として算出された全Ｎｇ個のｐ値が期待されることを、考慮に入れる必要がある。これを説明するため、Ｎｇ個の遺伝子のうちのいずれかについてｐ値≦ｐを観察する、全体の尤度、Ｐ（ｐ）が使用される。全遺伝子の独立性を仮定して、全体の尤度は、
【００６５】
【数６】
2007502992.tif 000008
で推定される。
【００６６】
遺伝子の独立性を仮定することは、明らかに過度の単純化であり、ｐ値およびＰ（ｐ）値を算出するための正しい方法は、一般的な無作為データセットに対して使用される帰無モデルのパラメーター（μｇｖ、Ａｉ、Ｄｋ、σｇｖ）とともに、ブートストラップ法を使用する方法である。
【実施例】
【００６７】
（実施例１．ＡＰＯＥ＊３−Ｌｅｉｄｅｎトランスジェニックマウスの肝臓からの遺伝子発現データの標準化）
上記標準化方法を例示するために、ＡｐｏＥ３−Ｌｅｉｄｅｎトランスジェニックマウスの研究を、実施した。合計９，５９６個の遺伝子を、１０枚のｃＤＮＡマイクロアレイを使用して分析した。サンプルは、４匹のＡｐｏＥ３−Ｌｅｉｄｅｎトランスジェニック（ＴＧ）マウス、および４匹の野生型（ＷＴ）マウスから集めた。最適化した実験設計を、図３に示す。従って、変種のベクトルは、
Ｖａｒｓ＝「１１１２２１１２２１１２２１１２２２２１」      （８）
であった。
【００６８】
トランスジェニックマウスおよび野生型マウスの標準化した値を比較するｔ検定を、適用した。図４は、ｔ検定からのｐ値および倍率比に基づいたデータの有意プロットを示す。上部の水平線は、Ｐ（ｐ）＝０．０５カットオフの全体の尤度を示すが、一方で、下方の線は、ｐ＝０．０５カットオフを示す。１６個の遺伝子だけが、最も厳しい上記の基準を満たすが、一方で、ｐ＜０．０５の範囲中に、７１３個の遺伝子が存在する。
【００６９】
（肝臓からのタンパク質データ）
８匹の異なる哺乳動物（４匹のトランスジェニック、および４匹の野生型）からの８サンプルを、８つの実験において分析した。従って、変種のベクトルは、
Ｖａｒｓ＝「１１１１２２２２」      （９）
となる。
【００７０】
質量分析（ＭＳ）のスペクトルを、各々が１６００個のピークを含む、合計４つのフラクションから選択した。上記ＭＳスペクトルを、ＩＭＰＲＥＳＳアルゴリズム（Ｌｅｉｄｅｎ大学で開発され、かつ米国特許６，７４３，３６４号明細書中に記載される）を使用して処理した。ＩＭＰＲＥＳＳピークの特徴付けソフトウェアを、ピークの有意性（０と１の間）を決定するために理論尺度（ＩＱ）情報を使用する。ＩＱ＞０．５の上記データセット中のピークは、上記サンプルの大多数（すなわち、８サンプルのうち５サンプル以上）について保持された。合計１０５９個のピーク（フラクション１から５個、フラクション３から２７１個、フラクション４から４５４個、およびフラクション５から３２９個）を、選択した。上記有意プロットは、図５に示される。Ｐ（ｐ）＝０．０５カットオフの基準を満たすピークはないが、ｐ＜０．０５では８４個のピークが、存在する。この場合、より多くのデータが、異なるフラクションに関して別々に標準化が実施されるべきであるかどうかを決定する必要がある。
【００７１】
（合成「ＧＩＳＴ」データ）
多数の色素を有するデータについて標準化する方法の検定を実施するために、２０００個のピーク、５つの色素、３つの変種、および６つの実験を有する合成データについての実験を、実施した。これは、潜在的に、Ｇｌｏｂａｌ  Ｉｎｔｅｒｎａｌ  Ｓｔａｎｄａｒｄｓ  Ｔｅｃｈｎｏｌｏｇｙ（「ＧＩＳＴ」）（Ｃｈａｒｋｒａｂｏｒｔｙ，Ａ．およびＲｅｇｎｉｅｒ，Ｆ．，Ｊ．Ｃｈｒｏｍａｔｏｇ．Ａ  ９４９，１７３〜８４（２００２））を使用して実施されたプロテオミクス実験に対応し得る。上記実験設計は、図６に示され、そしてまた、以下の変種のベクトル：
Ｖａｒｓ＝「１１２２３２２１１３３１１２２３２２１１１１３２２２２３１１」      （１０）
によって記載され得る。
【００７２】
各ピークのバックグラウンドは、Ｇａｕｓｓｉａｎ乱数生成器（同じ平均値および偏差に設定）を使用して選択されている。その後、３つの大きなピークが、各々の変種１および変種２について、それぞれ加えられているが、一方で、変種３は、コントロールとして残されている。図７〜図９は、各変種についての散布図および正規分布プロットを示す。３つのはずれ値が、明らかに変種１および変種２についてみられる。以下の倍率比：
【００７３】
【数７】
2007502992.tif 000009
を、各ピークについて算出し、そして、ｔ検定を、上記２つの変種を比較するために使用した。上記有意プロットは、図１０に示される。予想どおり、６つのはずれ値だけが、Ｐ（ｐ）＝０．０５カットオフの基準を満たすが、一方で、各ピーク（６つのはずれ値を除く）が、各サンプルについて、独立して無作為に生成されたという事実にもかかわらず、ｐ＜０．０５を満たす合計９４個のピークが存在する。
【００７４】
（図２におけるワークフローの例示的な実施例）
３つのさらなる実施例が、図２において例示されたフローチャート中に概説された上記実験方法、技術、および分析的アプローチをさらに例示するために、本明細書中に開示される。より詳細なフローチャートは、図１１、図１２、および図１３中に提示される。このフローチャートは、生物学的サンプルからデータセットを作製し、その後、上記閾値を超える存在量の変化を示す、遺伝子、タンパク質、または代謝産物のうちのいずれかのリストを抽出することを記載する。図１１、図１２、および図１３は、図２のより詳細な図として理解され得、特に、図２中の工程（２０５）から工程（２２０）までの工程に焦点を合わせている。図１４は、ネットワークの関連性と文献中の既知の関連性との比較に使用され得る、相関ネットワークを生成するために、上記抽出された成分リストを統合することを例示する（図２中の工程（２２０）、工程（２２５）、および工程（２３０））。上記例示した実施形態のより詳細な状況を提供するために、個々の図１５〜図２９が、提示される。これらの図は、図２、図１１、図１２、図１３、および図１４中に示される個々の工程に関して、直接的に描写したものである。
【００７５】
（実施例２．ＡＰＯＥ＊３−Ｌｅｉｄｅｎトランスジェニックマウスのシステム生物学的分析）
哺乳動物システムに対するシステム生物学的分析の適用についてのテストケースとして、アポリポタンパク質Ｅ３−Ｌｅｉｄｅｎ（ＡＰＯＥ＊３−Ｌｅｉｄｅｎ、ＡＰＯＥ＊３）トランスジェニックマウスを、選択した。ＡｐｏＥは、超低密度リポタンパク質（ＶＬＤＬ）およびＶＬＤＬレムナントの成分であり、そしてこのＡｐｏＥは、肝臓によるレセプター媒介性のリポタンパク質の再取り込みに必要とされる（ＧｌａｓｓおよびＷｉｔｚｔｕｍ，Ｃｅｌｌ  １０４，５０２（１９８９））。上記ＡＰＯＥ＊３−Ｌｅｉｄｅｎの変異は、コドン１２０〜１２６の直列重複によって特徴付けられ、そしてこの変異は、ヒトにおける家族性異常βリポタンパク質血症と関連する（ｖａｎ  ｄｅｎ  Ｍａａｇｄｅｎｂｅｒｇら、Ｂｉｏｃｈｅｍ．Ｂｉｏｐｈｙｓ．Ｒｅｓ．Ｃｏｍｍｕｎ．１６５，８５１（１９８６）；およびＨａｖｅｋｅｓら、Ｈｕｍ．Ｇｅｎｅｔ．７３，１５７（１９８６））。ヒトＡＰＯＥ＊３−Ｌｅｉｄｅｎを過剰発現するトランスジェニックマウスは、肝臓のＬＤＬレセプター認識の低下が原因で、食餌性高リポタンパク質血症およびアテローム性動脈硬化症に対する感受性が高いが、標準的な固形飼料の食餌を与えた場合、そのマウスは、９ヶ月目で、中程度Ｉ型（マクロファージ泡沫細胞）およびＩＩ型（細胞内脂質蓄積を有する脂肪線条）の病変を示すのみである（Ｊｏｎｇら、Ａｒｔｅｒｉｏｓｃｌｅｒ．Ｔｈｒｏｍｂ．Ｖａｓｃ．Ｂｉｏｌ．１６，９３４（１９９６））。
【００７６】
ＡＰＯＥ＊３−Ｌｅｉｄｅｎトランスジェニックマウス系統を、ヒトＡＰＯＥ＊３−Ｌｅｉｄｅｎ遺伝子、ＡＰＯＣ１遺伝子、およびＡＰＯＣ１およびＡＰＯＥ＊３の間に肝臓制御領域と名付けられた調節性エレメントを含む、２７ｋｂのゲノムＤＮＡ構築物をマウスの受精卵の雄性前核中にマイクロインジェクションすることによって、作製した。卵の供給源は、過剰排卵させた（Ｃ５７Ｂｌ／６Ｊ  ｘ  ＣＢＡ／Ｊ）Ｆ１雌であった。初代トランスジェニックマウスをさらに、トランスジェニック系統を確立するために、Ｃ５７Ｂｌ／６Ｊマウスとさらに交配させた。Ｆ２１世代〜Ｆ２２世代のトランスジェニック同腹仔および非トランスジェニック同腹仔を、これらの実験において使用した。全てのマウスに、標準的な固形飼料の食餌（ＳＲＭ−Ａ，Ｈｏｐｅ  Ｆａｒｍｓ，Ｗｏｅｒｄｅｎ，Ｔｈｅ  Ｎｅｔｈｅｒｌａｎｄｓ）を与え、そして、９週間目で屠殺した。この時に、血漿サンプル、尿サンプル、および肝臓組織サンプルを採取し、そして液体窒素中で凍結した。その後、各個体からのサンプルを、個々の遺伝子発現分析、タンパク質分析、および代謝産物分析用に細分した。標準的な固形飼料の食餌が与えられ、そして９週齢で屠殺された、野生型マウスおよびＡＰＯＥ＊３−Ｌｅｉｄｅｎマウスから採取した肝臓組織、血漿、および尿に対して適用された、ｍＲＮＡ発現プロファイル分析、可溶性タンパク質プロファイル分析、および脂質の差異的なプロファイル分析を組み合わせた結果が、以下に提示される。野生型マウスは、トランスジェニックマウスの特徴を比較するためのツールとして、すなわち、言い換えると、コントロールマウスとして使用される。
【００７７】
図１１〜図１３を参照して、調査されるべき生物学的状態（１１０５、１２０５、１３０５）は、トランスジェニック哺乳動物システムにおける脂質代謝であり、具体的には、ＡＰＯＥ＊３−Ｌｅｉｄｅｎトランスジェニックマウスにおけるアテローム性動脈硬化症および高脂血症である。（１１１０、１２１０、１３１０）で収集したサンプルは、トランスジェニックマウスから採取した肝臓組織、血漿、および尿からのサンプルである。
【００７８】
（肝臓遺伝子発現）
図１１を参照して、全ｍＲＮＡを、購入したＲＮＡｅａｓｙ  ｋｉｔ（Ｑｉａｇｅｎ，Ｇｅｒｍａｎｔｏｗｎ，Ｍａｒｙｌａｎｄ）を使用して、ホモジェナイズした肝臓組織から抽出した。その後、ｍＲＮＡを、購入したＯｌｉｇｏｔｅｘキット（Ｑｉａｇｅｎ，Ｇｅｒｍａｎｔｏｗｎ，Ｍａｒｙｌａｎｄ）を使用して、上記全ｍＲＮＡ調製物から、抽出した（１１１５）。遺伝子発現マイクロアレイデータを、マウスのＵｎｉＧｅｎｅ  １をスポットしたｃＤＮＡアレイ（Ｉｎｃｙｔｅ  Ｇｅｎｏｍｉｃｓ，Ｓｔ．Ｌｏｕｉｓ，Ｍｉｓｓｏｕｒｉ）を使用して取得した。１つ実施形態において、分散（ＡＮＯＶＡ）モデルの分析を、その技術固有の多様性を最適に減少する、サンプルペアリングの設計のために選択した。
【００７９】
ｍＲＮＡ量の実験（１１２０）を、肝臓組織について実施した。１つの実施形態において、この実験は、ｍＲＮＡのハイブリダイゼーションを含む。遺伝子発現、および／またはパターン認識の連続的な分析が、実施され得る。１つの実施形態において、ＰＡＲＣパターン認識プログラムを使用する。図１５は、ｍＲＮＡ量の実験を例示する。特に、遺伝子発現分析が、ＡＰＯＥ＊３トランスジェニックマウス対野生型マウスについての、マウスの肝臓ｍＲＮＡ発現比のプロットによって例示される。遺伝子発現データセット（１１２５）の例は、図１５に例示される肝臓遺伝子発現分析のみを含むのではなく、図１６に例示される遺伝子発現データ、および図１７に例示される遺伝子発現量の結果も含む。
【００８０】
（肝臓および血漿から抽出されたタンパク質のプロファイリング）
タンパク質を、凍結した肝臓組織サンプルおよび血漿サンプル（１２１０）から抽出した（１２１５）。クロマトグラフィー工程（１２２０）は、上記サンプルをさらに特徴付けるために利用され得る。１つの実施形態において、上記タンパク質は、上記クロマトグラフィー工程（１２２０）の後に、化学修飾される（１２２５）。別の実施形態において、上記タンパク質を、上記クロマトグラフィー工程（１２２０）または化学修飾工程（１２２５）のいずれかの後に、ペプチドへフラグメント化する（１２３０）。１つの実施形態において、フラグメント化（１２３０）は、上記タンパク質の部分的な加水分解によって実施される。第二のクロマトグラフィー工程（１２３５）は、フラグメント化工程（１２３０）の後であり得、そして、質量分光測定工程（１２４０）は、上記クロマトグラフィー工程（１２３５）の後であり得る。１つの実施形態において、ＰＡＲＣパターン認識プログラムを、上記タンパク質を数値化するために、使用する。ＧＩＳＴ同位体標識法もまた、利用され得る。タンパク質の同定は、質量分析またはＢｉｏＳｙｓｔｅｍａｔｉｃｓのいずれかで実施し得る。
【００８１】
タンパク質由来のデータセット（１２４５）の例が、図１８〜図２０中に示される。図１８は、ＡＰＯＥ＊３トランスジェニックマウス対野生型マウスからの、血漿のＬＣ／ＭＳ全イオンクロマトグラム（ＴＩＣ）の強度プロットを例示する。図１９において、かすかに検出可能な差異を明らかにし得る、ＬＣ／ＭＳプロファイリングからのＴＩＣが、示される。図１８および図１９の両方は、それらが１０００より多くのペプチドピークを含むものとして、データセット（１２４５）の複雑さを例示する。図２０は、５匹のトランスジェニックマウスおよび５匹の野生型マウスの上記消化した肝臓タンパク質から取得したＬＣ／ＭＳクロマトグラムを例示する。１つの実施形態において、ＬＣ／ＭＳは、エレクトロスプレーイオン化（ＥＳＩ）プローブを備える、ＬＣＱ  ＤｅｃａＸＰ（ＴｈｅｒｍｏＦｉｎｎｉｇａｎ，Ｓａｎ  Ｊｏｓｅ，ＣＡ）４重極イオントラップ型質量分析システムを使用して実施される。
【００８２】
（尿および血漿から抽出した代謝産物のプロファイル）
代謝産物を、尿サンプルおよび血漿サンプル（１３１０）から抽出した。上記尿サンプルを、１次元の１Ｈ  ＮＭＲを使用してプロファイルした（１３１５）。ＮＭＲスペクトルは、データセット（１３４０）のうちの１例である。データセット（１３４０）はまた、クロマトグラフィー工程（１３２０）によって血漿データから生成され得る。そして、その後に、代謝産物の化学修飾（１３２５）が続く。上記修飾された代謝産物（１３２５）は、一連のクロマトグラフィー（１３３０）工程、および質量分析（１３３５）工程によって、特徴付けられ、データセット（１３４０）を生成し得る。１つの実施形態において、上記血漿サンプルは、ＥＳＩによってイオン化され、かつＬＣ／ＭＳを使用して特徴付けられる。
【００８３】
代謝産物データセット（１３４０）の例が、図２１および図２２中に示される。図２１は、ＡＰＯＥ＊３マウスおよび野生型マウスについての、血漿から抽出した代謝産物の１Ｈ  ＮＭＲスペクトルを例示する。ＭｅＯＤ（δ＝３．３０）の−ＣＨ３シグナルを参照した後で、一覧表を、標準的なＶａｒｉａｎ  ＮＭＲソフトウェアを使用して準備した。これらの一覧表を得るために、閾値（ＳＮ比の約３倍に対応する）を上回るスペクトル中の全ての共鳴を収集し、そして統計学的分析の適用にとって適切なデータファイルフォーマットに変換した。図２２は、ＡＰＯＥ＊３マウスおよび野生型マウスについての、ＬＣ／ＭＳを使用して記録した血漿脂質の質量クロマトグラムを例示する。
【００８４】
（データセットの結合）
再度、図１１〜図１３を参照して、１つの実施形態において、遺伝子のデータセット（１１２５）、タンパク質のデータセット（１２４５）、および代謝産物のデータセット（１３４０）は、分子機能を決定し、そして細胞機構を明らかにするために、並行して分析される。多くのバイオインフォマティクスツールが、遺伝子の応答、タンパク質の活性、および代謝産物の運動性と結びつけるために利用され得る。データセット（１１２５、１２４５、１３４０）は、データ処理工程（１１３０、１２５０、１３４５（または、図２では２１０））に供される。ＩＭＰＲＥＳＳアルゴリズムは、ＬＣ／ＭＳクロマトグラムおよびＮＭＲスペクトルの両方において、バックグラウンドのノイズを減らすために使用され得る。別の実施形態において、上記ＩＭＰＲＥＳＳアルゴリズムを、ＰＡＲＣアルゴリズムに入力するためのＩＱファイルを生成するために、使用する。
【００８５】
１つの実施形態において、上記データを前処理する工程（１１３０、１２５０、１３４５）から算出されるデータは、統計分析工程（１１３５、１２５５、１３５０）で処理される。統計分析の適切な形態は、上記により詳細に記載される。上記前処理されたデータは、ＡＮＯＶＡアルゴリズムを使用して標準化され得る。別の実施形態において、標準化は、ＰＡＲＣアルゴリズムを使用して、データセットに関して実施され得る統計分析工程の後に、起こる。１つの実施形態において、分化するスペクトル成分は、上記統計分析によって生成された上記因子スペクトル中において同定される。
【００８６】
図２３は、標準化工程（２１５）によって処理されたスペクトルを示す。個々の遺伝子のスペクトル、タンパク質のスペクトル、および代謝産物のスペクトルは、上記のモデルを使用して標準化され、そしてその後に、上記個々の標準化されたスペクトルは、単一の因子スペクトルへ連結される。図２３において、データは、マウスの肝臓から抽出された生物学的サンプルに関して測定したデータである。上記連結したスペクトルを使用して、生体分子成分のタイプ間での直接比較が、実施され得る。
【００８７】
図２４〜図２５は、上記統計分析工程（１１３５、１２５５、１３５０）、およびその後の上記検査工程（１１４０、１２６０、１３５５）の例示的な実施形態を提供する。単純にするために、タンパク質血漿分析のみが提示されるが、この方法は、遺伝子および代謝産物の両方まで拡大適用され得る。図２４は、野生型マウスデータおよびＡＰＯＥ＊３トランスジェニックマウスのデータを、ペプチドイオン質量データに関するＰＣ−ＤＡ（１２５５）を使用して、クラスタリングすることを例示する。図２４中に示される２つの別のクラスターの検査（１２６０）は、イオンの質量が、２つのクラスターを区別することを明らかにする。図２５は、差因子スペクトルにおいてプロットされた有意差を示すペプチドイオンの質量を示す。１つの実施形態において、ｔ検定を、各々の上記分化するイオンに適用して、それらの有意性を検定する。別の実施形態において、負荷プロット（ｌｏａｄｉｎｇ  ｐｌｏｔ）が、因子スペクトルの代わりに使用される。
【００８８】
さらなる質量分析の分析工程（１２６５、１３６０）は、存在量レベルの閾値を上回る変化を示すさらなるタンパク質、ペプチド、または代謝産物を分析するために、実施され得る。１つの実施形態において、ＭＳ／ＭＳは、上記タンパク質、ペプチド、または代謝産物を分析し、かつ同定するために使用される。別の実施形態において、統計的に有意な変化を示す遺伝子、タンパク質、ペプチド、または代謝産物は、手動の検査工程（１１４０、１２６０、１３３５）の間に同定される。全ての遺伝子、タンパク質、ペプチド、および代謝産物の同定（１１４５、１２７０、１３６５）に続いて、これらの遺伝子、タンパク質、ペプチド、および代謝産物の一覧表が、今後の比較のために、抽出され、かつ保存される（１１５０、１２７５、１３７０）。
【００８９】
図２６は、マウスの血漿から抽出した上記タンパク質の加水分解によって生成された、上記ペプチドのＭＳ／ＭＳスペクトルを示す（図１２中の工程（１２６５）に対応する）。これらのペプチドフラグメント（ｂ７〜ｂ１７およびｙ５〜ｙ１６と標識される）は、データベースと比較され、その結果、フラグメント化されたタンパク質が、同定され得、かつ配列決定され得る。これは、図１２中の同定工程（１２７０）に対応する。この特定の場合において、同定されたタンパク質は、トランスジェニック操作によって誘導されたタンパク質である、ヒトＡｐｏＥ３である。
【００９０】
表Ｉは、遺伝子、タンパク質、および代謝産物の一覧から抽出した、重要な差次的に発現した成分を列挙している。この一覧表は、図１１〜図１３中に例示される工程（１１５０、１２７５、１３７０）に従って生成した。この抽出された成分の一覧はまた、図２中の成分の一覧を抽出する工程（２２０）に対応する。
【００９１】
【表１】
2007502992.tif 000010
１つの実施形態において、表Ｉに列挙された上記個々の生体分子成分は、標準化され、その結果、生体分子成分のタイプ間でより意味のある比較が、実施され得る。別の実施形態において、表Ｉに列挙された生体分子成分の一覧は、図２中の工程（２２５）と図１４中の工程（１４２０）とを踏まえて相関ネットワークを生成するために使用される。図２７は、生体分子成分のタイプ間の相関ネットワークを例示する。このネットワークは、非線形ＰＣＡ特徴の相関で生成されものであり、個々の生体分子成分間の潜在的な関連性を例示する。その後、上記相関ネットワークの関連性は、文献かまたは別の公的な情報源からの既存の知識と比較され得る。これは、図２中の工程（２３０）、または図１４中の工程（１４２５）に対応する。図２８は、相関ネットワークの関連性と刊行された情報との間の既知の関係のマップを例示する。
【００９２】
再度、図１４を参照すると、行為（１４３０）の生物マーカーまたは機構を決定するために分析される、相関ネットワークの関連性の例示的な実施形態を示す。上記既知の関係は、行為（１４３０）の生物マーカーまたは機構を決定するために分析され得る。１つの実施形態において、この相関ネットワークの関連性は、生体分子成分のタイプ間で関連し、かつ原因となる関係を決定するために、使用される（１４３５）。上記既知の関係はまた、生体分子成分のタイプ間で関連し、かつ原因となる関係を決定するために、使用され得る（１４３５）。
【００９３】
再度、図２を参照して、１つの実施形態において、上記システムは、摂動が与えられる（２３５）。上述のとおり、その後、この摂動が与えられたシステムは、上記摂動の原因となる機構を推論する前に、新たなデータセット、新たな相関ネットワーク、および新たな相関ネットワークの関連性を生成するために使用され得る。上記システムに対する摂動は、原因となる関係が、複数の生体分子成分のタイプ間で決定されるまで、繰り返され得る。
【００９４】
システム生物学的分析により決定された上記生物マーカーによって、上述のものと同じように、病気の集団と健康な集団とを区別するマーカーが、得られ得る。その後、この情報は、例えば、マーカーが、脱制御された経路の原因因子としてか、または下流の生成物のいずれかとして同定され得るときを決定するために、適切な生物学的背景に位置付けられ得る。上述のとおり、相関分析とネットワークモデリングとを組み合わせた、包括的な遺伝子プロファイリング、タンパク質プロファイリング、および代謝産物プロファイリングは、生物学的背景への洞察を提供する。そして、この知識のレベルは、治療剤を開発するために使用され得るか、または、病態生理学的機構をさらに明らかするために設計される有向の仮説主導型実験の土台として役立ち得る。
【００９５】
図２９は、生物マーカーまたは治療剤の観点から、代表的な、システム生物学的分析より生じ得る「提供物」または「送達可能物」を例示する。以下に記載される２つの例は、代表的なシステム生物学的分析を例示するだけでなく、これらのシステム生物学的分析から算出された情報が使用されるべき治療剤だけでなく、さらなる研究が必要な病態生理学的機構を決定するために、どのように使用されるのかという、より詳細な説明を例示する。
【００９６】
（実施例３．ＡＰＯＥ＊３−Ｌｅｉｄｅｎトランスジェニックマウスのシステム生物学的分析）
標準的な固形飼料の食餌を与え、そして９週齢にて屠殺した野生型マウスおよびＡＰＯＥ＊３−Ｌｅｉｄｅｎマウスから採取した肝臓組織、血漿、および尿に対して適用される、ｍＲＮＡ発現、可溶性タンパク質、および脂質の差次的なプロファイリング分析を組み合わせた結果が、以下に示される。各生体分子成分のタイプの分類分析の結果から、疾患素因を早期に発見するマーカーの存在を明らかにする。加えて、相関分析の結果は、協奏的な変化を起こす分子（遺伝子、タンパク質、および脂質まで及ぶ）のネットワークを示唆させる。
【００９７】
（動物）
ＡＰＯＥ＊３−Ｌｅｉｄｅｎトランスジェニックマウス系統を、ヒトＡＰＯＥ＊３−Ｌｅｉｄｅｎ遺伝子、ＡＰＯＣ１遺伝子、およびＡＰＯＣ１およびＡＰＯＥ＊３の間に肝臓制御領域と名付けられた調節性エレメントを含む、２７ｋｂのゲノムＤＮＡ構築物をマウスの受精卵の雄性前核中にマイクロインジェクションすることによって、作製した。卵の供給源は、過剰排卵させた（Ｃ５７Ｂｌ／６Ｊ  ｘ  ＣＢＡ／Ｊ）Ｆ１雌であった。初代トランスジェニックマウスをさらに、トランスジェニック系統を確立するために、Ｃ５７Ｂｌ／６Ｊマウスとさらに交配させた。Ｆ２１世代〜Ｆ２２世代のトランスジェニック同腹仔および非トランスジェニック同腹仔を、これらの実験において使用した。全てのマウスに、標準的な固形飼料の食餌（ＳＲＭ−Ａ，Ｈｏｐｅ  Ｆａｒｍｓ，Ｗｏｅｒｄｅｎ，Ｔｈｅ  Ｎｅｔｈｅｒｌａｎｄｓ）を与え、そして、９週間目で屠殺した。この時に、血漿サンプル、尿サンプル、および肝臓組織サンプルを採取し、そして液体窒素中で凍結した。その後、各個体からのサンプルを、個々の遺伝子発現分析、タンパク質分析、および代謝産物分析用に細分した。
【００９８】
（肝臓遺伝子発現）
全ｍＲＮＡを、購入したＲＮＡｅａｓｙ  ｋｉｔ（Ｑｉａｇｅｎ，Ｇｅｒｍａｎｔｏｗｎ，Ｍａｒｙｌａｎｄ）を使用して、ホモジェナイズした肝臓組織から抽出した。その後、ｍＲＮＡを、購入したＯｌｉｇｏｔｅｘキット（Ｑｉａｇｅｎ，Ｇｅｒｍａｎｔｏｗｎ，Ｍａｒｙｌａｎｄ）を使用して、上記全ｍＲＮＡ調製物から、抽出した。遺伝子発現マイクロアレイデータを、マウスのＵｎｉＧｅｎｅ  １をスポットしたｃＤＮＡアレイ（Ｉｎｃｙｔｅ  Ｇｅｎｏｍｉｃｓ，Ｓｔ．Ｌｏｕｉｓ，Ｍｉｓｓｏｕｒｉ）を使用して取得した。１つ実施形態において、分散（ＡＮＯＶＡ）モデルの分析を、その技術固有の多様性を最適に減少する、サンプルペアリングの設計のために選択した。
【００９９】
（肝臓タンパク質プロファイリング）
凍結肝臓組織を、液体窒素を加えることで冷たく保った予冷乳鉢中で、粉末化した。その後、Ｔ−ＰＥＲタンパク質抽出試薬（Ｐｉｅｒｃｅ  Ｃｈｅｍｉｃａｌ  Ｃｏ．，Ｒｏｃｋｆｏｒｄ，Ｉｌｌｉｎｏｉｓ）を、８μＬ／ｍｇの組織に加え、そのサンプルをさらに、超音波処理によってホモジェナイズした。その後、サンプルを１０，０００×ｇで５分間遠心分離し、上清を収集した。全タンパク質相対濃度を、Ｓｕｐｅｒ  ＳＷ３０００  ＴＳＫゲルカラム（Ｔｏｓｏｈ  Ｂｉｏｓｅｐ，Ｔｏｋｙｏ）およびＬＣ  Ｐａｃｋｉｎｇｓ  Ｕｌｔｉｍｔｅポンプ（Ｄｉｏｎｅｘ，Ｍａｒｌｔｏｎ，ＮＪ）から構成される、サイズ排除クロマトグラフィーシステムへ注入されたアリコートの全体統合クロマトグラム（ｉｎｔｅｇｒａｔｅｄ  ｗｈｏｌｅ−ｃｈｒｏｍａｔｏｇｒａｍ）から決定した。サンプルの複雑さを減少するために、タンパク質の上清を、０．１％のトリフルオロ酢酸（ＴＦＡ）の存在下で、水／アセトニトリル（ＭｅＣＮ）勾配で溶出したＰＯＲＯＳ  Ｒ２／Ｈカラム（４．６×１００ｍｍ）（Ａｐｐｌｉｅｄ  Ｂｉｏｓｙｓｔｅｍｓ，Ｆｏｓｔｅｒ  Ｃｉｔｙ，Ｃａｌｉｆｏｒｎｉａ）から構成される、ＶＩＳＩＯＮワークステーション（Ａｐｐｌｉｅｄ  Ｂｉｏｓｙｓｔｅｍｓ，Ｆｏｓｔｅｒ  Ｃｉｔｙ，Ｃａｌｉｆｏｒｎｉａ）上の逆相クロマトグラフィーを介して、分画した。タンパク質を、（１００ｍＭの重炭酸アンモニウム、５ｍＭの塩化カルシウム、および１０ｍＭのジチオスレイトール）中で、７５℃にて３０分間、切断し、熱変性し、そして還元させ、２５ｍＭのヨードアセトアミドで、７５℃にて３０分間、アルキル化し、その後、０．３％（ｗ／ｗ  トリプシン／タンパク質）で、３７℃にて２４時間、切断した。
【０１００】
（タンパク質のＬＣ／ＭＳ分析）
液体クロマトグラフィータンデム質量分析（ＬＣ／ＭＳ）を、エレクトロスプレーイオン化プローブから構成される、ＬＣＱ  ＤｅｃａＸＰ（ＴｈｅｒｍｏＦｉｎｎｉｇａｎ，Ｓａｎ  Ｊｏｓｅ，ＣＡ）４重極イオントラップ型質量分析システムを使用して、実施した。上記ＬＣの構成要素は、Ｓｕｒｖｅｙｏｒ自動サンプラーと４つ一組の勾配ポンプ（ＴｈｅｒｍｏＦｉｎｎｉｇａｎ，Ｓａｎ  Ｊｏｓｅ，ＣＡ）からなった。サンプルを移動相に懸濁し、Ｖｙｄａｃ  ｌｏｗ−ＴＦＡ  Ｃ１８カラム（１５０×１ｍｍ、５μｍ）（ＧｒａｃｅＶｙｄａｃ，Ｈｅｓｐｅｒｉａ，ＣＡ）を通して溶出した。このカラムを、溶媒Ａ（水／ＭｅＣＮ／酢酸／ＴＦＡ、９５／４．９５／０．０４／０．０１、ｖｏｌ／ｖｏｌ／ｖｏｌ／ｖｏｌ）で、無勾配的に２分間、５０μＬ／分で、溶出した。その後、４３分間にわたり、７５％の溶媒Ｂ（水／ＭｅＣＮ／酢酸／ＴＦＡ、２０／７９．９５／０．０４／０．０１、ｖｏｌ／ｖｏｌ／ｖｏｌ／ｖｏｌ）まで線形勾配が続く。エレクトロスプレーイオン化の電圧を、４．２５ｋＶに設定し、加熱した移動用キャピラリーを、２００℃に設定した。窒素シースおよび補助ガスの設定は、それぞれ、２５単位および３単位である。トリプシンペプチドの数値化のために、単一全スキャン質量スペクトルから構成されるスキャン周期は、陽イオンモードにおいて、４００〜２０００ｍ／ｚを越えて得られた。データ依存性生成物イオン質量スペクトル（ＭＳ／ＭＳ）はまた、ＴｕｒｂｏＳＥＱＵＥＳＴアルゴリズム（ＴｈｅｒｍｏＦｉｎｎｉｇａｎ，Ｓａｎ  Ｊｏｓｅ，ＣＡ）を使用するペプチド同定のために、得た。
【０１０１】
（肝臓脂質プロファイリング）
肝臓組織を、凍結肝臓し、微粉砕し、その後、超音波槽中で２時間、組織１ｍｇあたり２０μＬのイソプロパノールで、抽出した。その後、このサンプルを、遠心分離し、上清を収集した。その後、サンプルを４容量の水で希釈し、ＬＣ／ＭＳ分析のために取った。ＬＣ／ＭＳデータを、エレクトロスプレーイオン化プローブから構成される、ＬＣＱ（ＴｈｅｒｍｏＦｉｎｎｉｇａｎ，Ｓａｎ  Ｊｏｓｅ，ＣＡ）４重極イオントラップ型質量分析計を使用して、取得した。このＬＣの構成要素は、Ｗａｔｅｒｓ  ７１７シリーズ自動サンプラーおよび６００シリーズ単一勾配形成ポンプ（Ｗａｔｅｒｓ，Ｍｉｌｆｏｒｄ，Ｍａｓｓａｃｈｕｓｅｔｔｓ）からなった。サンプルを、２連で、順序不同で、Ｒ２ガードカラム（Ｃｈｒｏｍｐａｃｋ）によって保護されるＩｎｅｒｔｓｉｌカラム（ＯＤＳ  ３．５ｍｍ、１００×３ｍｍ）に、注入した。３つの移動層を、以下：
（１）（水／ＭｅＣＮ／酢酸アンモニウム／酢酸、９３．９／５／１／０．１、ｖｏｌ／ｖｏｌ／ｖｏｌ／ｖｏｌ）
（２）（アセトニトリル／イソプロパノール／酢酸アンモニウム／酢酸、６８．９／３０／１／０．１、ｖｏｌ／ｖｏｌ／ｖｏｌ／ｖｏｌ）および、
（３）（イソプロパノール／ジクロロメタン／酢酸アンモニウム／酢酸、４８．９／５０／１／０．１、ｖｏｌ／ｖｏｌ／ｖｏｌ／ｖｏｌ）
の溶出において使用した。上記カラムを、０．７ｍＬ／分で、以下：
工程（１）０分〜１５分間、７０％のＡ、３０％のＢ、０％のＣで開始し、５％のＡ、９５％のＢ、および０％のＣで終了し、そして
工程（２）Ａにおいては変更なし、Ｂは９５％〜３５％、Ｃは０％〜６０％で２０分間の勾配の
２工程勾配工程を使用して、溶出した。エレクトロスプレーイオン化電圧を、４．０ｋＶに設定し、加熱した移動用キャピラリーを２５０℃に設定した。窒素シースおよび補助ガスの設定は、それぞれ、７０単位および１５単位である。代謝産物の数値化のために、単一全スキャン（１秒／スキャン）質量スペクトルから構成されるスキャン周期は、陽イオンモードにおいて２５０〜１２００ｍ／ｚにわたり得られた。
【０１０２】
（ＬＣ／ＭＳデータの前処理）
ＬＣ／ＭＳデータセットを、Ｘｃａｌｉｂｅｒ機器制御ソフトウェア（ＴｈｅｒｍｏＦｉｎｎｉｇａｎ，Ｓａｎ  Ｊｏｓｅ，Ｃａｌｉｆｏｒｎｉａ）に機能的に組み込まれているファイル変換器を使用して、ＡＮＤＩ（．ｃｄｆ）フォーマットに変換した。その後、ＩＭＰＲＥＳＳアルゴリズム（ＴＮＯ  Ｐｈａｒｍａ，Ｚｅｉｓｔ，Ｔｈｅ  Ｎｅｔｈｅｒｌａｎｄｓ）を、自動ピーク検出およびピークデータの質の評価のために、上記変換したファイルに適用した。このプログラムは、各々の質量トレースのクロマトグラフィーの質について、それらの情報量を見積もることによって、各質量トレースを評価する。各質量電荷比でのＬＣ／ＭＳクロマトグラムを、ノイズスパイクを除去するために平滑化し、その後、上記トレースのエントロピーを、式１２を使用して算出した。Ｈの逆数値を取り、そして、全ての結果を最大値までスケール化することで、各々の質量トレースに、スケール化したクロマトグラムの質番（Ｉｍｐｒｅｓｓ  Ｑｕａｌｉｔｙ（ＩＱ）と呼ばれる）を与えた：
【０１０３】
【数８】
2007502992.tif 000011
その後、ＩＱの閾値を選択した。そして、ピークのＩＱがこの閾値を下回った場合は、このピークは、質が良くないと判断され、かつ以下に記載されるクラスター分析には使用しなかった。
【０１０４】
（マイクロアレイデータの標準化）
上に記載されるように、上記データは、以下のモデルによって表され得る：
ｙｇｉｋ＝μｇｖ＋Ａｉ＋Ｄｋ＋εｇｉｋ      （１３）
ここで、遺伝子および変種の効果は、μｇｖによって、アレイの効果は、Ａｉによって、色素の効果は、Ｄｋによって、およびエラーは、εｇｉｋによって記載される。このエラーは、通常、ゼロ平均で分布され、そして、分散（σ２ｇｖ）は、各遺伝子および変種について異なることを許容しない。上記モデルの最適なパラメーターは、最大尤度推定器を使用して計算される。従って、上記各々の特定のアレイおよび色素について、上記サンプルは、以下：
【０１０５】
【数９】
2007502992.tif 000012
のようにスケール化される。
【０１０６】
（有意性の統計的検定）
トランスジェニックサンプルおよび野生型サンプルからの、異なる平均値が標準化された強度の統計的有意性を評価するために、ｔ検定を、Ｎ個の遺伝子の各々に適用し、そして、対応するｐ値を算出した。各遺伝子についての倍率変化の統計的有意性を評価する場合、全部Ｎ個のｐ値を集め、その結果、いくつかのｐ≦０．０５でのｐ値を、予測した。これを考慮するために、Ｎ個の遺伝子のうちのいずれかについてｐ値≦ｐを観察する全体尤度Ｐ（ｐ）を、使用した。全ての遺伝子の独立性を仮定して、上記全体尤度を、以下：
【０１０７】
【数１０】
2007502992.tif 000013
で推定した。
【０１０８】
（ＰＣＤＡ分析および収集プロット）
主成分および判別分析（ＰＣＤＡ）を、ＩＭＰＲＥＳＳアルゴリズムにより上述のように前処理された、トリプシンペプチドおよび液体ＬＣ／ＭＳのプロファイルに適用した。これは、ＷＩＮＬＩＮ統計ソフトェア（ＴＮＯ  Ｐｈａｒｍａ，Ｚｅｉｓｔ，Ｔｈｅ  Ｎｅｔｈｅｒｌａｎｄｓ）を使用して行った。
【０１０９】
（肝臓遺伝子発現のマイクロアレイ）
マウスの肝臓ｍＲＮＡサンプルを、ＵｎｉＧｅｎｅ  １  ｃＤＮＡをスポットしたマイクロアレイでのハイブリダイゼーションのために、図３０Ａに示される「ループ設計」に従って、組み合わせた。このペアリング方法は、ＡＮＯＶＡモデルに基づいており、これは、遺伝子発現データの最適な標準化の土台を提供し、かつ因子から生じ得る可変性（例えば、核酸の効果または色素の効果の間の、不揃いのハイブリダイゼーション率）の寄与を最小化するために、設計された。ｍＲＮＡサンプルを、示されるように、二重のハイブリダイゼーションのために、Ｃｙ３およびＣｙ５で標識した。
【０１１０】
図３０Ｂに示されるｃＤＮＡマイクロアレイデータの分散プロットから明らかなように、比較的わずかな遺伝子が、９５％の信頼水準で差次的に発現した。値を、野生型マウスおよびＡＰＯＥ＊３−Ｌｅｉｄｅｎトランスジェニックマウスにおける発現の平均値として、プロットした。データ点を、統計学的意義に基づいて色分けした。より厳格な全体尤度Ｐ（ｐ）を満たすものはかなり少ないので、データが無作為に、しかし誤って、ｐ値＜０．０５を有し得るという起こりうる事象を除外するために、評価を試みる。
【０１１１】
表ＩＩは、トランスジェニックと野生型コントロールとの間の倍率比が、０．８よりも小さいかまたは１．２よりも大きいかのいずれかであった、遺伝子のサンプルセットを列挙する。発現における差異がより狭い限度であるにもかかわらず、観察された比較的低いｐ値が、ＡＮＯＶＡモデルの統計的な利点を反映する。トランスジェニック動物におけるアポリポタンパク質ＡＩおよび、アポリポタンパク質Ｂのアナログの発現がより低いレベルであることは注目すべきであるが、アポリポタンパク質Ｆのアナログが、より高かった。興味深いことに、ＡＰＯＥ＊３−Ｌｅｉｄｅｎマウスから得られた血漿の前分析は、タンパク質レベルで約２倍のダウンレギュレートを明らかにした。加えて、ペルオキシソーム増殖因子活性化レセプターα（ＰＰＡＲα）の発現は、２つの集団間で差異はなかったが、肝臓脂肪酸結合タンパク質（Ｌ−ＦＡＢＰ）は、トランスジェニックにおてい４３％高かった。ＰＰＡＲαの重要な役割は、脂質の代謝に関与するタンパク質の遺伝子発現を開始することであるが、実験的な事実が、Ｌ−ＦＡＢＰが、リガンドの活性化を示す速度を制御することによって、転写因子の活性を制御し得ることを示唆する。脂質プロファイリング分析は、導入遺伝子の存在、およびＰＰＡＲαレベルにおいて変化がないことにより、脂質代謝に実際に影響が与えられることを示す。これらのデータは、Ｌ−ＦＡＢＰについての制御的な役割を支持する。
【０１１２】
【表２】
2007502992.tif 000014
（肝臓タンパク質の定量的プロファイリング）
サンプルの複雑さを約２０の因子まで減らすために、可溶性肝臓タンパク質のオフラインの逆相分離を、最初に行った。ＥＳＩ−ＬＣの構成を、何百という継続的な注入を操作することが可能である質量分析計と組み合わせた。次に、データを、連続的なＭＳ／ＭＳスキャンを取得することなく、ＭＳのみのスキャン周期を使用して取得し、スキャン間でカラムが溶出する間に、周期時間を減少し、かつ情報の損失の最小化した。図３１Ａに示されるように、ＬＣ／ＭＳクロマトグラムを、５匹のＡＰＯＥ＊３−Ｌｅｉｄｅｎマウスおよび５匹の野生型マウス由来の消化した肝臓タンパク質フラクションについて取得した。その後、ＩＭＰＲＥＳＳアルゴリズムを各データセットに対して適用し、ピーク強度およびシグナルのクオリティーに関する情報を抽出した。０．５のＩＭＰＲＥＳＳの質値は、下限の閾値として選択された。この閾値は、低いクオリティーのシグナルデータをさらなる分析から排除する。その後、クラスタリングを、ＷＩＮＬＩＮソフトウェアに組み込まれた主成分判別分析（ＰＣＤＡ）ツールを使用して、実施した。図３１Ｂに示されるように、２つの明確なクラスターが、一方はトランスジェニックマウスに関して、他方は野生型マウスに関して観察された。図３１Ｃに図解される、因子スペクトルの検査は、２つのクラスターを区別したイオンの質量を提供した。ｔ検定を、区別するイオンの各々に対して適用し、有意性を検定し、そして各ペプチドについてのＬＣ／ＭＳ／ＭＳスペクトルを取得した。６個のトリプシンペプチドは各々、Ｌ−ＦＡＢＰの消化から生じ、質量電荷比４４６、５９９、７０６、８９２、８９５、および１０５８を有し、それらは、図３１Ｃにおいて標識される。因子スペクトルは、自然状態では半定量的であるので、ＩＭＰＲＥＳＳによって収集したピーク強度情報を、相対的差異を計算するために使用した。このプロファイリング分析の結果は、Ｌ−ＦＡＢＰが、野生型コントロールに対して、トランスジェニックマウスにおいて４４％までアップレギュレートされることを示した。これは、基本的に、上述のｍＲＮＡ発現の観察に１対１で相関した。表ＩＩＩは、タンパク質分析の結果をまとめる。
【０１１３】
【表３】
2007502992.tif 000015
（肝臓脂質の定量的プロファイリング）
脂質を、タンパク質分析について使用したものと同様のストラテジーを使用してプロファイルした。２組のデータセットを、各動物について取得した。抽出プロトコルおよびＬＣシステムは、より大きな非極性脂質（例えば、ジアシルグリセロール（ＤＧ）およびトリアシルグリセロール（ＴＧ））を分画するために、設計された。この取得において得られたものはまた、ホスファチジルコリン（ＰＣ）およびリゾホスファチジルコリン（ＬｙｓｏＰＣ）脂質の定量的プロファイルであった。ピーク情報を取得するためのＩＭＰＲＥＳＳでのデータ前処理の後に、ＰＣＤＡクラスタリング分析を、ＷＩＮＬＩＮを使用して実施した。図３２Ａに示されるように２つの集団のマウスが、２つの明確なクラスターを形成した。図３２Ｂに図解されるＰＣＤＡ因子スペクトルは、多くの脂質が、２つの集団間における差異に寄与することを示す。リゾホスファチジルコリン（ＬｙｓｏＰＣ）、ジアシルグリセロール（ＤＧ）、ホスファチジルコリン（ＰＣ）、およびトリアシルグリセロール（ＴＧ）の大部分を含む質量電荷比の範囲が示される。
【０１１４】
表ＩＶに要約されるように、多数のトリアシルグリセロールが、トランスジェニックマウス中ではより高い存在量であった一方で、より低い存在量のものは見出されなかった。同様に、２つの、リゾホスファチジルコリン、１−パルミトイル−２−ヒドロキシ−ｓｎ−グリセロ−３−ホスホコリン（ＬｙｓｏＰＣ  Ｃ１６：０）、および１−ステアロイル−２−ヒドロキシ−ｓｎ−グリセロ−３−ホスホコリン（ＬｙｓｏＰＣ  Ｃ１８：０）は、ＡＰＯＥ＊３−Ｌｅｉｄｅｎマウス中で、より高いレベルにて発見された一方で、他のＬｙｓｏＰＣについては、有意な差異が観察されなかった。興味深いことに、ジアシルグリセロールのサブクラスとホスファチジルコリンのサブクラスとの間では、トランスジェニック動物におけるより高い存在量に関する全体的な傾向は、観察されなかった。このことは、導入遺伝子の挿入によって負わされた脂質代謝の崩壊が、脂質レベルの制御において複雑な多因子性変化をもたらすことを示唆した。
【０１１５】
【表４】
2007502992.tif 000016
（議論）
図３３Ａ〜図３３Ｃに強調されるように、差次的なゲノムプロファイリング、プロテオミックプロファイリング、メタボロミックプロファイリングに基づく包括的なシステムの分析は、マウスが基本的に疾患の臨床的徴候を何も示さない状況のもとで、ＡＰＯＥ＊３−Ｌｅｉｄｅｎトランスジェニックマウスと野生型コントロールとを区別する多くの新規知見をもたらした。ＰＣＤＡクラスター分析、および区別する因子の同定の後で、生体分子成分のタイプ、ｍＲＮＡ、タンパク質、および脂質の各々の相対存在量が計算され、それぞれが、図３３Ａ、図３３Ｂ、および図３３Ｃに示される。値は、ｎ＝４〜５の別々の動物についての±ＳＥＭの平均（＊ｐ＜０．０５）を表す。個々に関して、これら実体の各々は、被験体を高脂血症およびアテローム性動脈硬化症に罹りやすくする異常代謝状態の生物マーカーとして、役立ち得る。
【０１１６】
ＡＰＯＥ＊３−Ｌｅｉｄｅｎマウスにおける疾患の早期マーカーとして同定された、アテローム性動脈硬化症における主要な種類を、図３４に例示する。ヒトにおいて、ＡＰＯＥ＊３−Ｌｅｉｄｅｎの変異は、低密度リポタンパク質レセプター（ＬＤＬＲ）に関する親和性が減少した機能障害性アポタンパク質Ｅ改変体を生じる。同様に、ＡＰＯＥ＊３−Ｌｅｉｄｅｎトランスジェニックマウスはまた、高脂血症を発症し、食餌誘発性アテローム性動脈硬化症に患りやすい。標準的な固形飼料の食餌で飼育された若年マウスにおけるシステム生物学を介して見出された、病理学における早期マーカーを、矢印で示す（上向きの矢印は、トラスジェニックにおけるアップレギュレートを意味し、一方で下向きの矢印は、トラスジェニックにおけるダウンレギュレートを意味する）。これらのマーカーとしては、Ａｐｏ  ＡＩおよびＬ−ＦＡＢＰのｍＲＮＡおよびタンパク質、ならびに種々の脂質分子が挙げられる。例えば、リポタンパク質関連ホスホリパーゼＡ２（これは、血小板活性化因子アセチル加水分解酵素としても記載される）は、循環においてＰＣからのＬｙｓｏＰＣの生成を触媒する酵素であり、心臓疾患に関する危険因子として同定されている（Ｐａｃｋａｒｄら、Ｎ．Ｅｎｇｌ．Ｊ．Ｍｅｄ．３４３  １１４８（２０００））。ＬｙｓｏＰＣは、病因に寄与する早期の炎症促進性現象に寄与し、この現象において、ＬｙｓｏＰＣは，脂肪線の発生の間に単球の接着および化学走性を増加させる。本研究において、ＡＰＯＥ＊３−Ｌｅｉｄｅｎトランスジェニックマウスの肝臓中において増加した、２つのＬｙｓｏＰＣ化合物を、同定し、肝臓中の早期の炎症後現象が、アテローム性動脈硬化症の病因となり得ることを示唆した。
【０１１７】
アポリポタンパク質およびＬ−ＦＡＢＰは、生物マーカーの第２の高分子グループを構成する。アポリポタンパク質ＡＩ（ＡｐｏＡＩ）は、ＡＰＯＥ＊３−Ｌｅｉｄｅｎマウスの血漿中において、野生型コントロールと比べて、有意により低い。ここで、このアポリポタンパク質のｍＲＮＡ転写物は、肝臓中でより低くなることが見出され、これは、以前の観察を支持し、従って、疾患に対する素因に寄与する因子としての低下したＡｐｏＡＩレベルおよびＨＤＬレベルに関する役割を支持する。
【０１１８】
上昇したＬ−ＦＡＢＰの証拠はまた、ゲノム分析およびプロテオミック分析の両方によって提供された。ＡｐｏＥ欠損マウスはまた、脂肪細胞の脂肪酸結合タンパク質が欠損しており、障害性マクロファージ機能に関与する機構を介するアテローム性動脈硬化症に対して保護された（Ｍａｋｏｗｓｋｉら、Ｎａｔ．Ｍｅｄ．７，６９９（２００１））。Ｌ−ＦＡＢＰは、同じ細胞内脂肪酸結合タンパク質のファミリーの一員である。Ｌ−ＦＡＢＰは、ＰＰＡＲαのリガンドに関するシャトルとして作用することによって転写制御において役割を果たすと考えられている（Ｗｏｌｆｒｕｍら、Ｐｒｏｃ．Ｎａｔｌ．Ａｃａｄ．Ｓｃｉ．ＵＳＡ  ９８，２３２３（２００１））。ヒトにおいて、ＡｐｏＡＩの発現は、ＰＰＡＲαによって転写制御される。特に興味深いのは、本研究の結果は、Ｌ−ＦＡＢＰとＰＰＡＲα媒介性ＡｐｏＡＩ発現との間の関係において脱共役を示す。なぜなら、Ｌ−ＦＡＢＰレベルが上昇し、ＰＰＡＲαレベルが変化せず、そしてＡｐｏＡＩ発現が低下したからである。従って、これらの結果は、追加であるが必須の、存在しないかまたはダウンレギュレートされたことを示唆する。この因子が、ＰＰＡＲαに関する特定のリガンドであり得ることが推測されることが興味深い。
【０１１９】
まとめると、本出願人らは、ｍＲＮＡレベル、タンパク質レベル、および脂質レベルでのプロファイリングのシステム生物学的アプローチの結果が、アテローム性動脈硬化症の発症に関するＡＰＯＥ＊３−Ｌｅｉｄｅｎトランスジェニックマウスの初期の素因に関する多くの新規生物マーカーを明らかにしていることを示してきた。まとめると、このような実体の収集は、多因子性疾患を区別する際により高い精度を有する独特の複合的な生物マーカーを構成し得る。このシステム生物学的アプローチは、いくつかのこれらの生物マーカー間の相関関係の解明を可能にし、疾患の機構、および治療処置のための手段の両方に関しての洞察を提供した。
【０１２０】
（実施例４．システム生物学的アプローチ：ＡｐｏＥ３−Ｌｅｉｄｅｎトランスジェニックマウスモデルの多並行分析）
複雑な哺乳動物の高脂血症モデルおよびアテローム性動脈硬化症モデルの病因学的プロセスにおけるシステム生物学的分析の結果が、以下に示される。プロテオミック分析とメタボロミック分析とを統合し、そしてトランスジェニックシステムに内在する疾患因子を、定量的に区別するプラットフォームが記載される。多因子疾患（例えば、高脂血症およびアテローム性動脈硬化症）への洞察を得るために、ＡｐｏＥ＊３−Ｌｅｉｄｅｎトランスジェニックマウスの全血漿中のタンパク質および代謝産物の成分をプロファイルするためのシステム生物学的アプローチを、使用した。この結果を、公知の脂質代謝プロセスであると確認し、トランスジェニック疾患モデルにおけるリポタンパク質レベルおよび脂質レベルでの新規差異を解明する。
【０１２１】
本研究に適用された、システム分析に対するアプローチの全体（全血漿並行タンパク質−代謝プロファイリングスキーム）は、図３５において模式的に概説される。ＡｐｏＥ＊３−Ｌｅｉｄｅｎマウスおよびコントロールマウス由来の全血漿、脂質、およびタンパク質のフラクションを、ＮＭＲおよびＭＳによって分析した。代謝産物データセットおよびタンパク質データセットの両方を、ＩＭＰＲＥＳＳアルゴリズムを通してフィルタ処理し、テキストに記載されるように、ＷＩＮＬＩＮ統計ソフトェアを使用して同時にクラスタリングした。分離分析方法および分光学的分析方法（例えば、ＨＰＬＣ、ＮＭＲ、およびＬＣ／ＭＳ）を、強力な統計学的パターン認識アルゴリズム（例えば、判別分析）と組み合わせ、コントロール  対  遺伝的に摂動を与えた動物の血漿中において、生化学的成分を迅速にクラスタリングし、かつ同定した。結果はこの主要な差異（＞２倍）および目立たないが、統計学的に有意な差異（ｔ検定でｐ＜０．０５）を、タンパク質レベルおよび代謝産物レベルで示す。
【０１２２】
（動物）
ＡＰＯＥ＊３−Ｌｅｉｄｅｎトランスジェニックマウス系統を、ヒトＡＰＯＥ＊３−Ｌｅｉｄｅｎ遺伝子、ＡＰＯＣ１遺伝子、およびＡＰＯＣ１およびＡＰＯＥ＊３の間に肝臓制御領域と名付けられた調節性エレメントを含む、２７ｋｂのゲノムＤＮＡ構築物をマウスの受精卵の雄性前核中にマイクロインジェクションすることによって、作製した。卵の供給源は、過剰排卵させた（Ｃ５７Ｂｌ／６Ｊ  ｘ  ＣＢＡ／Ｊ）Ｆ１雌であった。初代トランスジェニックマウスをさらに、トランスジェニック系統を確立するために、Ｃ５７Ｂｌ／６Ｊマウスとさらに交配させた。Ｆ２１世代〜Ｆ２２世代のトランスジェニック同腹仔および非トランスジェニック同腹仔を、これらの実験において使用した。全てのマウスに、標準的な固形飼料の食餌（ＳＲＭ−Ａ，Ｈｏｐｅ  Ｆａｒｍｓ，Ｗｏｅｒｄｅｎ，Ｔｈｅ  Ｎｅｔｈｅｒｌａｎｄｓ）を与え、そして、９週間目で屠殺した。この時に、血漿サンプル、尿サンプル、および肝臓組織サンプルを採取し、そして液体窒素中で凍結した。その後、各個体からのサンプルを、個々の遺伝子発現分析、タンパク質分析、および代謝産物分析用に細分した。
【０１２３】
（血漿リポタンパク質プロファイリング）
通常の固形飼料の食餌（ＳＲＭ−Ａ，Ｈｏｒｅ  Ｆａｒｍｓ，Ｗｏｅｒｄｅｎ，Ｔｈｅ  Ｎｅｔｈｅｒｌａｎｄｓ）を与え続けられた９週齢のマウス由来の血漿を、ＬＣパッキングクロマトグラフィーシステム（Ｄｉｏｎｅｘ，Ｍａｒｌｔｏｎ，ＮＪ）上のＳｕｐｅｒ  ＳＷ３０００  ＴＳＫゲルカラム（Ｔｏｓｏｈ  Ｂｉｏｓｅｐ，Ｔｏｋｙｏ）を介する、サイズ排除クロマトグラフィーによって、分画した。各サンプルについてのタンパク質の濃度を、Ｂｒａｄｆｏｒｄアッセイによって決定し、最も低い濃度に標準化した１０μＬの全血漿を、２０ｍＭのビス−トリスプロパン（ｐＨ６．９）；１００ｍＭのＮａＣｌ中に、５０μＬ／分で注入し、かつ無勾配敵に溶出した。３００ｋＤより大きな分子量の範囲に対応する基底分離（ｂａｓｅ−ｒｅｓｏｌｖｅｄ）ピークを、分散したフラクションとして収集した。タンパク質を、１００ｍＭの重炭酸アンモニウム、５ｍＭの塩化カルシウム、および１０ｍＭのジチオスレイトール中で、７５℃にて３０分間、消化し、熱変性し、そして還元させ、２５ｍＭのヨードアセトアミドで、７５℃にて３０分間、アルキル化し、その後、０．３％（ｗ／ｗ  トリプシン／タンパク質）で、３７℃にて２４時間、消化した。
【０１２４】
（タンパク質のＬＣ／ＭＳ分析）
液体クロマトグラフィー質量分析（ＬＣ／ＭＳ）を、エレクトロスプレーイオン化プローブを備える、ＬＣＱ  ＤｅｃａＸＰ（ＴｈｅｒｍｏＦｉｎｎｉｇａｎ，Ｓａｎ  Ｊｏｓｅ，ＣＡ）４重極イオントラップ型質量分析システムを使用して、実施した。上記ＬＣの構成要素は、Ｓｕｒｖｅｙｏｒ自動サンプラーと４つ一組の勾配ポンプ（ＴｈｅｒｍｏＦｉｎｎｉｇａｎ，Ｓａｎ  Ｊｏｓｅ，ＣＡ）からなった。サンプルを移動相に懸濁し、Ｖｙｄａｃ  ｌｏｗ−ＴＦＡ  Ｃ１８カラム（１５０×１ｍｍ、５μｍ）（ＧｒａｃｅＶｙｄａｃ，Ｈｅｓｐｅｒｉａ，ＣＡ）を通して溶出した。このカラムを、溶媒Ａ（水／アセトニトリル／酢酸／トリフルオロ酢酸、９５／４．９５／０．０４／０．０１、ｖｏｌ／ｖｏｌ／ｖｏｌ／ｖｏｌ）で、無勾配的に２分間、５０μＬ／分で溶出し、その後、４３分間にわたり、７５％の溶媒Ｂ（水／アセトニトリル／酢酸／トリフルオロ酢酸、２０／７９．９５／０．０４／０．０１、ｖｏｌ／ｖｏｌ／ｖｏｌ／ｖｏｌ）まで線形勾配が続く。エレクトロスプレーイオン化の電圧を、４．２５ｋＶに設定し、加熱した移動用キャピラリーを、２００℃に設定した。窒素シースおよび補助ガスの設定は、それぞれ、２５単位および３単位であった。トリプシンペプチドの数値化のために、単一全スキャン質量スペクトルから構成されるスキャン周期を、陽イオンモードにおいて、４００〜２０００ｍ／ｚにわたって得た。データ依存性生成物イオン質量スペクトル（ＭＳ／ＭＳ）はまた、ＴｕｒｂｏＳＥＱＵＥＳＴアルゴリズム（ＴｈｅｒｍｏＦｉｎｎｉｇａｎ，Ｓａｎ  Ｊｏｓｅ，ＣＡ）をＭＡＳＣＯＴ検索アルゴリズム（Ｍａｔｒｉｘ  Ｓｃｉｅｎｃｅ）をＮＣＢＩｎｒ、ＳｗｉｓｓｐｒｏｔおよびＭＳＤＢのデータベース検索と組み合わせて使用するペプチド同定のために、得た。
【０１２５】
（代謝産物の分析）
マウスの血漿サンプルを、脂質全体および代謝産物全体の分析のために、０．６ｍＬのイソプロパノールを１５０μＬの全血漿に加えることで調製し、その後、タンパク質を沈殿させ、かつ除去するために遠心分離した。上清の５００μＬのアリコートを、濃縮して乾燥させ、ＮＭＲ分析の前に、７５０μＬのＭｅＯＤ中に再溶解させた。ＬＣ／ＭＳ用のサンプルを調製するために、４００μＬの水を、１００μＬの上清に加え、２００μＬのこの混合液を、ＬＣ／ＭＳ用の自動サンプラーに移動した。
【０１２６】
（ＮＭＲ分析）
ＮＭＲスペクトルを、Ｖａｒｉａｎ  ＵＮＩＴＹ  ４００ＭＨｚ分光計上にて、２９３Ｋの温度で作動する陽子ＮＭＲ設定を使用し、完全自動様式にて、３連で記録した。自由誘導減衰（ＦＩＤ）を、８．０００ＨＺのスペクトル幅で、６４Ｋのデータポイントとして収集し、４５度のパルスを、４．１０秒の収集時間、および２秒のリラグゼーションディレイ（ｒｅｌａｘａｔｉｏｎ  ｄｅｌａｙ）で、使用した。このスペクトルは、５１２ＦＩＤの累積によって得られた。上記スペクトルを、標準的なＶａｒｉａｎソフトウェアを使用して処理した。０．５Ｈｚのラインを広幅化し、かつ手動で基線を補正する指数ウィンドウ関数を、全てのスペクトルに適用した。ＣＤ３ＯＤ（δ＝３．３０）の−ＣＤ３シグナルを参照した後で、線の一覧を、標準的なＶａｒｉａｎ  ＮＭＲソフトウェアを使用して作製した。これらの一覧を得るために、約３倍の信号雑音比に対応する閾値を上回る、スペクトル中の全ての線を、収集し、そして、統計分析のアプリケーションに適するデータファイルに変換した。
【０１２７】
（ＬＣ／ＭＳ分析）
ＬＳＱ  Ｃｌａｓｓｉｃ（ＴｈｅｒｍｏＦｉｎｎｉｇａｎ，Ｓａｎ  Ｊｏｓｅ）を、血漿脂質および代謝産物の成分のＭＳスペクトルを取得するために使用した。このＬＣの構成要素は、Ｗａｔｅｒｓ  ７１７シリーズ自動サンプラーおよび６００シリーズ単一勾配形成ポンプ（ＷａｔｅｒｓＣｏｒｐｏｒａｔｉｏｎ，Ｍｉｌｆｏｒｄ，ＭＡ）からなった。サンプルを、Ｒ２ガードカラム（Ｃｈｒｏｍｐａｃｋ）によって保護されるＩｎｅｒｔｓｉｌカラム（ＯＤＳ  ３．５μｍ、３ｍｍ×１００ｍｍ）上に、注入した。マウス血漿抽出物の７５μＬのアリコートを、無作為な順番で、２度注入した。ランダムシーケンスを、統計学的統計値から得た結果についての分析の間、可能性のあるドリフトの有害な影響を阻止するために、適用した。溶出勾配を、３つの移動層：
（１）（水／アセトニトリル／酢酸アンモニウム（１Ｍ／Ｌ）／ギ酸、９３．９：５：１：０．１、ｖｏｌ／ｖｏｌ／ｖｏｌ／ｖｏｌ）
（２）（アセトニトリル／イソプロパノール／酢酸アンモニウム（１Ｍ／Ｌ）／ギ酸、６８．９：３０：１：０１、ｖｏｌ／ｖｏｌ／ｖｏｌ／ｖｏｌ）および、
（３）（イソプロパノール／ジクロロメタン／酢酸アンモニウム（１Ｍ／Ｌ）／ギ酸、４８．９：５０：１：０．１、ｖｏｌ／ｖｏｌ／ｖｏｌ／ｖｏｌ）
を使用して形成した。上記サンプルを、４段階の勾配：
（１）３０％〜９５％の緩衝液Ｂで１５分間以上の勾配；
（２）９５％〜３５％のＢおよび６０％のＣで２０分間の勾配、この工程で５分間保持される；
（３）３５％のＢおよび６０％のＣの、それぞれ、９５％および０％までの迅速な１分間勾配；ならびに
（４）９５％の緩衝液Ｂから３０％まで、５分間以上の時間で、戻す  
によって、０．７ｍＬ／分で、分画した。
【０１２８】
エレクトロスプレーイオン化電圧を、４．０ｋＶに設定し、加熱した移動キャピラリーを２５０℃に設定した。窒素シースおよび補助ガスの設定は、それぞれ、７０単位および１５単位であった。代謝産物の数値化のために、単一全スキャン（１秒／スキャン）質量スペクトルから構成されるスキャン周期を、陽イオンモードにおいて２００〜１７００ｍ／ｚにわたって得た。
【０１２９】
（ＮＭＲデータの前処理）
ＮＭＲスペクトルを、ＷＩＮＬＩＮ統計ソフトウェアパッケージ（ＴＮＯ  Ｐｈａｒｍａ，Ｚｅｉｓｔ，Ｔｈｅ  Ｎｅｔｈｅｒｌａｎｄｓ）を用いて、手動で整列化した。
【０１３０】
（ＬＣ／ＭＳデータの前処理）
ＬＣ／ＭＳデータファイルを、Ｘｃａｌｉｂｕｒソフトウェア（ＴｈｅｒｍｏＦｉｎｎｉｇａｎ）を使用して、ＮｅｔＣＤＦフォーマットに変換した。この変換されたファイルを、ＩＭＰＲＥＳＳ取得後ノイズ減少および標準化ソフトウェア（ＴＮＯ  Ｐｈａｒｍａ，Ｚｅｉｓｔ，Ｔｈｅ  Ｎｅｔｈｅｒｌａｎｄｓ）で、評価して、各ＬＣ／ＭＳファイルについてのフィンガープリントスペクトルを取得した。このプログラムは、その情報量を評価することで、そのクロマトグラフィーの質についての各質量トレースを評価する。これは、スパイクを除去するための平滑化の後に、そして式１２に従ってトレースのエントロピーを各質量について計算することによって、実施される。Ｈの逆数値を取り、そして、最大値に対して全ての結果をスケール化することで、各質量トレースに、スケール化されたクロマトグラフの質、またはＩＱを与える。
【０１３１】
（ＰＣＡ分析およびＰＣ−ＤＡ分析）
主成分分析（ＰＣＡ）および主成分判別分析（ＰＣ−ＤＡ）を、整列化された血漿ＮＭＲスペクトルのフィンガープリントスペクトル、およびＩＭＰＲＥＳＳ前処理されたＬＣ／ＭＳスペクトルに適用した。これは、ＷＩＮＬＩＮ統計ソフトウェア（ＴＮＯ  Ｐｈａｒｍａ，Ｚｅｉｓｔ，Ｔｈｅ  Ｎｅｔｈｅｒｌａｎｄｓ）を使用して行った。
【０１３２】
（差次的な代謝産物ＮＭＲ分析）
代謝産物の分析についてのパターン認識およびクラスタリング法を評価するために、２連のアプローチを使用した。ここで、ＮＭＲを、ＬＣ／ＭＳの後に続く最初のスクリーニング法（種々の生物システムにおけるメタボロームプロファイリングのための、ベンチマーク分析法として、確立されている）として、利用した（Ｒａａｍｓｄｏｎｋら、Ｎａｔｕｒｅ  Ｂｉｏｔｅｃｈ．１９，４５（２００１）；Ｎｉｃｈｏｌｓｏｎら、Ｘｅｎｏｂｉｏｔｉｃａ  ２９，１１８１（１９９９）；Ｆｉｅｎら、Ａｎａｌ．Ｃｈｅｍ．７２，３５７３（２０００））。ＮＭＲデータの処理を容易にするため、上記ＷＩＮＬＩＮソフトウェアパッケージを、上記野生型のデータセットとトランスジェニックのデータセットとの間での分散の程度をクラスター化し、かつ見積もる為に適用した。予備ＮＭＲスクリーンに基づく、十分な差異は、ＭＳおよびＭＳ／ＭＳを使用する、より詳細な分析に値することを明らかにした。
【０１３３】
２０匹のマウス（各グループについて、ｎ＝１０）からの全血漿サンプルを、全域代謝産物ＮＭＲ分析のために使用した。代表的な、４００ＭＨｚ  １Ｈ  ＮＭＲについて、ＭｅＯＤ中の７５０μＬのタンパク質を取り除いたサンプルを、野生型マウス血漿サンプル（ＷＴ）およびＬｅｉｄｅｎマウス血漿サンプル（ＴＧ）の両方について、３連のスペクトルを生成するために使用した（これらは、図３６に例示される）。ＭｅＯＤ（δ＝３．３０）の−ＣＤ３シグナルを参照した後で、線の一覧を、標準的なＶａｒｉａｎ  ＮＭＲソフトウェアを使用して作製した。これらの一覧を得るために、約３倍の信号雑音比に対応する閾値を上回る、スペクトル中の全ての共鳴を、収集し、そして、統計分析のアプリケーションに適するデータファイルに変換した。血漿代謝産物成分の最初の分析のためのＮＭＲフィンガープリンティングを使用する意図は、特定の化合物に対するシグナルを与えることではなく、十分なクラスター化を示し、従って、より詳細な分析に値するサンプルであるかを確立することである。上記ＮＭＲデータの厳密な検査は、匹敵する線の共鳴位置における小さな違いを明らかにした。線の上記位置の差は、上記サンプル中の化合物の相対濃度、および装置の不安定さ（例えば、温度および磁場の均一性）に起因し、これらは、手動で補正した。このように処理されたスペクトルを、判別成分分析（ＰＣ−ＤＡ）クラスタリングのために、ＷＩＮＬＩＮ統計分析ツールにインポートした。
【０１３４】
図３７は、Ｌｅｉｄｅｎマウス（三角で表される）およびコントロールマウス（丸で表される）についてのＮＭＲデータのクラスタリングを示す、ＰＣ−ＤＡスコアプロットを例示する。ＷＩＮＬＩＮは、データが標準化され、そして主成分分析（ＰＣＡ）を受けた後の、図式的なクラスタリングの結果を与える。上記クラスター内の各点は、上記前処理されたスペクトルの３連のセットのうちの１つを表すように、空間的に位置される。各々の上記３連のスペクトルからの濃度強度を、ＰＣ−ＤＡクラスターセットを構築するために使用した。主成分分析の第一の工程は、多くの直交する新たな変数のセット（主成分と呼ぶ）を得るために、分散／共分散マトリックスから、固有ベクトルを抽出することである。これらのベクトルは、オリジナルデータ中の分散の最大量を説明するそれらの能力において、最適化される。強く相関したデータにおいて、少数の上位に入る主成分は、上記オリジナルセット中で、有意な分散を再生成するのに十分である。ＰＣＡを、適用して、コントロールマウスおよびＡＰＯＥ＊３  Ｌｅｉｄｅｎマウスの、ＮＭＲスペクトルを整列化した部分的な線形一致（ＰＬＦ）を調査するのに必要な特徴の数を減少した。その後、上記サンプルの１／１５の主成分軸への投影を、線形判別分析についての起始点として使用した。
【０１３５】
因子スペクトルを、負荷ベクトルの図式的に回転させることによって、上記スコアプロット中のクラスターの位置と上記スペクトル中のオリジナルの特徴とを相関させるために使用した（Ｗｉｎｄｉｇら、Ａｎａｌ．Ｃｈｅｍ．５６、２２９７（１９８４））。因子スペクトルプロットの差異（図３８に示される）は、種々の代謝産物成分を表す多くの線によって特徴付けられる。これは、寄与因子の範囲、より具体的には、トランスジェニックマウスおよびコントロールマウスの集団のクラスタリングを容易にするイオンｍ／ｚによって定義される。上記プロットの軸の上下の線の高さは、全体の分散への寄与の幅に、直接的に関連し、軸の下に伸びる因子は、トランスジェニック動物におけるより高いスペクトル強度に対応する。ＰＣ−ＤＡが、単一の特定の方向におけるクラスターに分かれるので、中央軸の下に投影する線は、トランスジェニックマウスの血漿におけるより高い強度のＮＭＲスペクトルパターン成分を表す。中央軸の上に伸びる線は、上記コントロールグループに対して、より高い絶対濃度で存在する因子を象徴する。
【０１３６】
２つのカテゴリーの最下分離の方向において作製される因子スペクトルを、上記観察されたカテゴリーの分離の原因である代謝産物のタイプへの洞察を与えるために、使用した。仮の結果は、Ｌｅｉｄｅｎサンプルとコントロールサンプルとの間の定量的な分散に対する主要な寄与物としての、δ３．８ｐｐｍ〜δ４．２ｐｐｍの領域および脂質の領域（δ１．２ｐｐｍ〜δ０．８ｐｐｍ）に対するＰＣ−ＤＡ負荷プロット点に基づいた。
【０１３７】
ＮＭＲ分光測定の制限は、上記技術の低い固有の感受性から、ならびにＮＭＲスペクトルの高い複雑さおよび情報の内容から、生じる。上記技術の感受性はまた、検出される化合物の最小閾値濃度により影響される。それらの制限にかかわらず、パターン認識技術と合わせたメタボロームプロファイリングに基づいたＮＭＲが、メタボリックデータを包括的なシステムレベルでの分析へ統合するための、強力な分析的アプローチであることは、明らかである。この研究において、しかしながら、ＮＭＲスクリーンの目的は、特定の分子を同定することではなく、むしろ、サンプル集団間での定性的な差異の程度が存在するかどうかを決定する方法を使用することである。
【０１３８】
（代謝産物およびタンパク質の成分の同期的な分析が、予測されたパターンおよび新規のパターンを生じる）
トランスジェニック（ｎ＝４）マウスおよびコントロール（ｎ＝４）マウスの血漿からの代謝産物抽出物を、イソプロパノール沈降法により調製した。１００μＬの抽出物に４００μＬの水を加えることにより、上記サンプルを、ＬＣ／ＭＳ分析に供した。図３９は、４００〜１７００ｍ／ｚの質量範囲にわたる単一のスキャンモードを使用して収集されたＴＩＣを示す。ＬＣ／ＭＳスペクトルに統計分析を適用するために、上記生データファイルを、最初にＮｅｔＣＤＦフォーマットに変換し、そしてＩＭＰＲＥＳＳノイズ減少および標準化ソフトウェアを使用して処理した。このプログラムは、各々の質量トレースのクロマトグラフィーの質について、それらの情報量を見積もることによって、各質量トレースを評価する。各質量電荷比でのＬＣ／ＭＳクロマトグラムを、スパイクを除去するために平滑化し、その後、上記トレースのエントロピーを、式１２を使用して算出した後に行った。ＩＭＰＲＥＳＳによって標準化された質量強度は、スケール化されたクロマトグラフィーの質の番号かまたはＩＱを割り当てる。主成分分析を実施するために、図３９中のクロマトグラムに基づいた上記ＩＱを、ＷＩＮＬＩＮにインポートし、そして、判別分析の分離を、２つの最初の主成分ベクトルに基づいて得た。
【０１３９】
プロテオミクス全血漿分析を、リポタンパク質複合体を含むフラクションの方へ偏らせた。これは、上記Ｌｅｉｄｅｎの変異と関連する統計的に最も明らかな変化がこのクラスのタンパク質で生じるという、選択された上記トランスジェニックモデルに基づいた予測と一致した。上記トランスジェニック（ｎ＝４）およびコントロール（ｎ＝４）の動物からの全血漿サンプルを、分析的サイズ排除クロマトグラフィーによって分画し、そして、高分子量の血漿タンパク質成分と対応するフラクションを、実験プロトコルに記載されるように単離した。２３分および２７分において溶出された２つの主要な初期ピーク（それぞれ、全血漿の成分のＶＬＤＬ（フラクション１）およびＨＤＬ（フラクション２）と対応する）を、次に続く全ての操作に使用した。フラクション１とフラクション２とを含むタンパク質をトリプシンで処理して、タンパク質分解性ペプチドを生成した。
【０１４０】
ＭＳ分析からのＶＬＤＬフラクションからのＴＩＣが、野生型マウス（ＷＴ）およびＬｅｉｄｅｎマウス（ＴＧ）に関して図４０中に示される。８つの代表的なサンプル全てについて収集したＭＳ／ＭＳスペクトルを、ＮＣＢＩの非重複の、ヒトおよびマウスのデータベースに対するヒットを得るために、ＴｕｒｂｏＳＥＱＵＥＳＴによって分析した。これらの最初のヒットの同一性を、ＭＡＳＣＯＴデノボ配列決定およびデータベース検索ツールを使用して、さらに検証した。タンパク質の同一性を指定するための閾値は、全残基数のうちの２０％で設定した最小限の配列補償範囲に基づいた。上記タンパク質ＭＳデータを、ＩＱ値スペクトルを生成することによって、上記代謝産物成分と同じような方法で、クラスター化し、その後、判別分析を行った。
【０１４１】
血漿の代謝産物およびタンパク質の成分間での定量的な関係性を観察するために、連結した異種のデータセットの集合体を使用した。オリジナルの個々のデータセットを、別々に統合し、そして、これらのセットのＩＭＰＲＥＳＳの質ｍ／ｚの値を合計し、そして統計クラスタリング分析に供した。結果として生じるスコアプロット（図４１中に例示される）は、野生型（ＷＴ）およびトランスジェニック（ＴＧ）動物についてのＰＣ−ＤＡクラスターを示す。これは、Ｄ１における最大分離を達成するために、回転した２つの主成分に基づき生成した。各点は、上記個々の動物についての、代謝産物およびタンパク質分散因子（オリジナルデータセットの６０％）の一次結合を示す。
【０１４２】
代謝産物およびペプチドのスペクトルからのフィルタ処理されたｍ／ｚ強度を、因子プロット中に線形様式でまとめた（図４２中に示される）。中央軸に沿う線形分布は、コントロールグループとトランスジェニックグループとの間における、２方向の、算出された分散の寄与をとともに、タンパク質および代謝産物の成分を表す。積極的に寄与する主要な因子は、５０の名目カットオフ量を上回る突出が見られる。全体の分散に対して消極的な寄与因子は、設定境界の−５０を下回って突き出ている。
【０１４３】
１６０１および３４０１の名目値を、第二のタンパク質および代謝産物の成分における各ｍ／ｚ値に、それぞれ加えることによって、異種の実験のデータを、図４２中に示されるように、並行して分析した。有意な寄与強度は、因子プロットの特異的な閾値パラメーター（この例では、５０に設定した）に基づいて、スコアを付けた。上記ＷＴとＴＧのデータセット間で、主要な差別化要因であると見出された質量を、抽出して、そして、ＬＣ／ＭＳ／ＭＳによって同定した。差別化因子の強度の組合せ（生データおよびＩＱスコア）を、統計的有意性（Ｐ＜０．０５）および倍率変化の計算のためにＬＣ／ＭＳクロマトグラムにおいて、直接、測定した。
【０１４４】
上記結果は、上記ＡＰＯＥ＊Ｌｅｉｄｅｎの表現型と関連するリポタンパク質および脂質の異常度に関するこれまでの知見を実証する複合プロファイルを指す（Ｍｅｎｓｅｎｋａｍｐら、Ｊ．Ｈｅｐａｔ．３３，１８９（２０００）；ｖａｎ  ｄｅｎ  Ｍａａｇｄｅｎｂｅｒｇら、Ｊ．Ｂｉｏｌ．Ｃｈｅｍ．２６８，１０５４０（１９９３）；Ｗｉｌｌｉａｍｓ  ｖａｎ  Ｄｉｊｋら，Ａｒｔｅｒｉｏｓｃｌｅｒ．Ｔｈｒｏｍｂ．Ｖａｓｃ．Ｂｉｏｌ．１９，２９４５（１９９９）；およびＭｅｎｓｅｎｋａｍｐら、Ｊ．Ｂｉｏｌ．Ｃｈｅｍ．２７４，３５７１１（１９９９））。具体的には、ヒトＡＰＯＥ*３Ｌｅｉｄｅｎの対立遺伝子改変体が、上記トランスジェニック動物において、発現しかつ機能的に活性であることを、タンパク質レベルで、本発明者らが示すことが可能であった。このことは、リポタンパク質由来の血漿の、ＶＬＤＬ（図４２中のタンパク質成分１）フラクション、およびＬＤＬ／ＨＤＬ（図４２中のタンパク質成分２）フラクションへの、取り込によって証明された。あるいは、マウスのＡｐｏＡ１は、トランスジェニックマウスの血漿中の量が２倍少ないことが見出されている。このことは、これらの動物において、上記ＬＤＬ／ＨＤＬ複合体への上記アポリポタンパク質の取り込みのより低い程度を示している。
【０１４５】
ＨＤＬ代謝を支配する根底にあるプロセスは、完全に定義されてはないが、血漿中のＨＤＬレベルは、アテローム性動脈硬化症の罹患率に反比例することが示されている（Ｃａｌｌｏｗら、Ｇｅｎｏｍｅ  Ｒｅｓ．１０，２０２２（２０００）；ならびにＧｌａｓｓおよびＷｉｔｚｔｕｍ）。多くの異なる機構が、ＨＤＬ血漿を制御し得る。血漿ＨＤＬの低下に寄与する、マウスモデル中で同定された最も突出した因子としては、ａｐｏＡ１、ａｐｏＥ、リン脂質転移タンパク質（ＰＬＴＰ）における欠陥、および、コレステロールエステル転移タンパク質（ＣＥＴＰ）またはスキャベンジャー受容体ＳＲＢの過剰発現が挙げられる（Ｃａｌｌｏｗら；Ｗｉｌｌｉａｍｓｏｎら、Ｐｒｏｃ．Ｎａｔｌ．Ａｃａｄ．Ｓｃｉ．ＵＳＡ  ８９，７１３４（１９９２）；およびＷａｎｇら、Ｊ．Ｂｉｏｌ．Ｃｈｅｍ．２７３，３２９２０（１９９８））。Ｌｅｉｄｅｎの変異が、欠陥ＡＰＯＥ対立遺伝子と機能的に類似していると仮定すると、Ｌｅｉｄｅｎモデルとの関連において、より低いＨＤＬレベルは、少なくとも部分的には、上記ＡｐｏＥ＊３の導入遺伝子の機能の結果である可能性が非常に高い。全ての内在性ＡｐｏＡ１が減少する可能性の１つとしては、ｈＡｐｏＥ３の成分の過剰発現、およびＬＤＬ／ＨＤＬ集合体のための優先的な補充に起因する、化学量論的不均衡である。
【０１４６】
本研究は、高度に複雑なシステムを特徴付けるためのマルチレベルのアプローチの有用性を実証する。高度な内容の分析用出力を生成し、かつ複合データセットから算出される統合された主要成分因子を比較することによって、ＡｐｏＥ＊３−Ｌｅｉｄｅｎの表現型を定義する主要なリポタンパク質代謝媒介物における同一性および相対存在量の迅速な解明が可能であった。生物流体の分析のみを踏まえると、この試みは、疾患を説明するために、定量的なプロテオミクスデータとメタボロームデータとを統合する方法で、システム生物学の原理を適用する最初の試みであることを示す。将来的に、このアプローチは、複数の組織の差次的な転写分析の形態にゲノム成分を含めることによって、また、遺伝子の摂動の多面的（ｐｌｅｏｔｒｏｐｉｃ）効果の理解を、真に全域的にすることで向上させられる可能性がある。
【０１４７】
（実施例５．システム生物学アプローチ：代謝疾患研究）
（要旨）
本実施例の最終的な目標は、本発明に従った分子分析およびデータ統合能力を実証することである。医学的な関心の一般的な領域は、代謝疾患であり、そして分析されるべき物質は、２つの動物種（げっ歯類および非ヒト霊長類）およびヒト被験体由来の血清サンプルであった。げっ歯類の各グループのサブセット（罹患およびコントロール）を、薬物処置した。プロジェクトの初期相（第Ｉ相）の間、試験者は、３種のサンプル源（げっ歯類、非ヒト霊長類、およびヒト）が存在することを知っていたが、各々の種におけるサンプルの分類の詳細については、情報を伏せた（ｂｌｉｎｄ）。
【０１４８】
本研究の具体的な目的は、以下の通りであった。
（第Ｉ相）
・動物およびヒト被験体由来の情報が伏せられた血清サンプルの代謝産物分析およびタンパク質分析を行うこと；そして
・血清の代謝産物プロファイルおよびタンパク質プロファイルに基づいて、サンプルの分類を行うこと。
（第ＩＩ相）
・情報を明らかにした後、決定されたサンプルの分類と、実際のサンプルグループとを比較すること；
・サンプルタイプの各々に対して、サンプルの１つのグループを別のグループから差別化するために使用され得る分子成分（生物マーカー）を定義すること；
・罹患表現型または薬物処置された表現型の根底にある生化学プロセスの洞察を得るために、その生物マーカーに対する相関ネットワークを構築すること；および
・罹患げっ歯類をコントロールのげっ歯類から差別化する分子成分が、罹患ヒト患者をコントロールのヒト被験体から差別化する分子成分と類似するか否かを決定すること。
【０１４９】
ラット血清サンプルに対する、代謝産物プロファイルおよびタンパク質プロファイルの情報を伏せた分析は、情報を明らかにした際に、実際のサンプルグループ（罹患＋ビヒクル、罹患＋薬物、コントロール＋ビヒクル、コントロール＋薬物）に正確に対応する４つの明確に異なるグループを示した。非ヒト霊長類サンプルに対するプロファイルの情報を伏せた分析は、情報を明らかにした際に、罹患グループおよびコントロールグループに正確に対応する２つの異なるグループを示した。ヒトサンプルに対して、代謝産物プロファイルおよびタンパク質プロファイルの情報を伏せた分析は、使用された分析プラットフォームに依存して、異なる数（４つまたは２つ）のグループを示した。脂質プロファイルのみに基づいた分析は、情報を明らかにした際に、８６％の精度で罹患患者に対応し、８９％の精度でコントロールの被験体に対応する２つのグループを示した。
【０１５０】
動物血清サンプルグループとヒト血清サンプルグループとの間を差別化する多数の代謝産物およびタンパク質を同定した。サンプルにおけるこれらの生物マーカーの相対レベルは、疾患および薬物応答の根底にある生化学プロセスへの洞察を提供した。注目すべき知見の１つは、血清タンパク質レベルにおける、罹患げっ歯類での薬物処置の影響であった。２つ目の異なる知見は、罹患げっ歯類および罹患患者において、対応するコントロール被験体におけるレベルに対して、１５０個を超える血清脂質のレベルが、ほとんど同一に広範に変化したことであった。ヒト疾患のモデルとしてのげっ歯類モデルの検証として、試験者はまた、優れた精度で罹患患者をコントロールのヒト被験体から区別するために、罹患げっ歯類対コントロールのげっ歯類を正確に分類することが見出されている血清脂質生物マーカーのセットも使用可能であった。
【０１５１】
（導入）
本実施例の最終的な目標は、前臨床および臨床の血清サンプルの比較研究に適用する場合の、プロテオミクス、メタボロミクスおよびインフォマティクス技術の統合されたプラットフォームを評価するための基礎を提供することであった。血清サンプルは、代謝疾患のげっ歯類モデルにおける薬物処置研究、ヒト被験体における代謝疾患の比較研究、および非ヒト霊長類における関連状態の研究から提供された。このプロジェクトを、２つの相に分けた。第Ｉ相において、試験者をサンプル情報に対して情報を伏せ、ＮＭＲ技術とＭＳ技術との組み合わせを用いた代謝産物およびタンパク質の比較定量プロファイリングを実施した。教師なしクラスタリング分析のようなインフォマティクス法をデータに適用し、実験グループが正確に区別されているか否かを決定した。第Ｉ相の最後に、データの情報を明かし、そして使用した方法が、高い精度でグループを決定したことを示した。第二の相の重点は、げっ歯類薬物処置／疾患研究における４つの実験グループの差別化に寄与する代謝産物およびタンパク質の同定、および個々の分子種が互いに相関する程度の決定であった。さらに、罹患ヒト被験体およびコントロールのヒト被験体と、それらのげっ歯類モデルの対応物との間の相関を探索し、ヒト疾患と動物モデルとの間の類似点および相違点を明らかにした。本発明およびその技術を例示するために、本実施例は、特定の結果のみを強調する。
【０１５２】
（サンプル情報）
本研究の第Ｉ相において、試験者を、サンプルが影響されていない（正常）かまたは影響されている（疾患および／または薬物処理）かのいずれかに関して、情報を伏せた。第ＩＩ相の前に、サンプルの情報を明らかにした。実験グループおよびサンプル数を、以下に列挙する。
【０１５３】
（Ａ．代謝疾患のげっ歯類モデルにおける薬物処置研究）
罹患のげっ歯類および非罹患のげっ歯類（コントロール）に治療薬物を投与した薬物処置研究からの合計３２個の血清サンプル（各６００μＬ）を、以下のように細分した。
【０１５４】
ｎ＝８  ビヒクルで処置されたコントロール
ｎ＝８  薬物で処置されたコントロール
ｎ＝８  ビヒクルで処置された罹患グループ
ｎ＝８  薬物で処置された罹患グループ。
【０１５５】
（Ｂ．ヒト被験体における代謝疾患の比較研究）
代謝疾患と診断された個体およびコントロールからの合計４２個の血清サンプル（サンプル当たり３００〜４００μＬ）を、以下のように細分した。
【０１５６】
ｎ＝１４  代謝疾患と診断された被験体
ｎ＝２８  コントロール。
【０１５７】
（Ｃ．非ヒト霊長類の疾患研究）
非ヒト霊長類からの合計２４個の血清サンプル（サンプル当たり３００〜８５０μＬ）を、プロファイルした。
【０１５８】
ｎ＝１３  正常の非ヒト霊長類のサル
ｎ＝１２  罹患した非ヒト霊長類のサル。
【０１５９】
（利用方法−分析的プロファイリング）
差示的なプロテオミクスおよびメタボロミクスへの実施例におけるアプローチとして、広範な分子成分の定量的プロファイリングを可能にするいくつかの異なる分析方法を使用する。これらの方法は、分析のエンドポイントとして、ＮＭＲまたはＭＳのいずれかを利用する。プロファイリングのプラットフォームは、頑強性、再現性、選択性、およびダイナミックレンジを考慮に入れて最適化されており、量の規模のオーダーおよび生化学的分類の範囲にかかり得る分子を調査するように設計されている。各プラットフォームは、１回の分析で多数の成分（数百〜数千）をプロファイルするための容量を有し、そしてソフトウェアツールを、計算的分析およびインフォマティクス分析への統合のための定量的情報の抽出を容易にするために使用した。本研究において適用された方法を、以下に列挙する。
１．タンパク質ＬＣ／ＭＳ：ペプチドおよびタンパク質のプロファイリングおよび同定を可能にする。
２．ＣＰＭＧ  ＮＭＲ：低分子量代謝産物の増強されたＮＭＲ測定。
３．拡散を校正したＮＭＲ（ｄｉｆｆｕｓｉｏｎ−ｅｄｉｔｅｄ  ＮＭＲ）：リポタンパク質関連代謝産物の増強された測定。
４．脂質ＬＣ／ＭＳ：脂質および非極性代謝産物のプロファイリングのために最適化されている。
【０１６０】
（利用方法−データ処理）
プロファイリング実験から得られたＮＭＲスペクトルおよびＬＣ／ＭＳクロマトグラムの結果は、数百の分子の相対量を示す数百ものピークを含み得る。データ処理ソフトウェアツールは、各データファイルからのこの情報の抽出ならびにサンプルセットをまたがって測定されたピーク強度の比較を可能にするために使用される。上記のように、代表的には、データ処理工程は、ピーク検出および相対強度（ピーク統合）の測定、１つのサンプル分析と別のサンプル分析とで発生し得るピーク位置における小さな差異（すなわち、特定のピークに対するＮＭＲ化学シフトまたはＬＣ／ＭＳ保持時間における小さな差異）の補正のための「整列」工程、ならびに各ピークに対する識別番号（またはインデックス数）の割り当て（この結果、サンプルにまたがって比較され得る）を含む。
【０１６１】
（利用方法−データ分析）
データをいくつかの異なる統計学的アプローチを使用して分析した：（１）サンプルの教師なし（ｕｎｓｕｐｅｒｖｉｓｅｄ）クラスタリング（ＣＯＳＡ階層型クラスタリングを含む）、（２）サンプルグループの間で異なるピークを決定するための一変量統計、および（３）全てのサンプルに対する代謝産物の個々の成分とタンパク質のセットとの間の相関を同定するための相関ネットワーク分析。さらに、分類の目的のためのサポートベクトルマシン（ＳＶＭ）分類器を用いるいくつかの予備的データ分析を行った。図４３は、データ分析ワークフローの模式図である。このデータ分析プロセスの要素を、実施された順に以下に列挙する。
１．データの標準化：データセット内のプラットフォーム特異的変動に対する調整。
２．探索性教師なしクラスタリング法の適用：
−ＣＯＳＡ
−主成分分析
−Ｋ−ｍｅａｎｓクラスタリング（ヒトサンプルのみ）
−ニューラルネットワーク（ヒトサンプルのみ）。
３．同定のためのピーク選択：一変量統計法（ペアワイズ、両側ｔ検定）を用いて、有意な、識別力のあるピークを同定し、同定に対して優先順位を付ける。
４．相関ネットワーク：ピークのペアの間の統計的相関を決定する。
５．データの可視化：ソフトウェアツールを使用し、データベースの情報と実験から生成されたデータとを統合する。
【０１６２】
（血清サンプルの分析に関する代謝疾患のげっ歯類モデルに対する結果および考察−教師なしクラスタリング）
最初の分析は、情報が伏せられたげっ歯類血清サンプルから収集されたデータの教師なしクラスタリングに焦点を当てた。教師なしクラスタリングは、サンプルのクラス分類のまたはサンプル収集における異なるグループの数の予備知識なしで、サンプルを分類することを試みる統計学的方法である。そのワークフローの概略を、図４４に提供する。概して、複数の分析プラットフォーム由来の複数のデータセットを、標準化およびクラスタリングした。個々のデータセットの範囲では正確または明確にクラスタリングされないが、複数のデータセットは、さらなるクラスタリング分析に連結（すなわち、組み合わせおよび／または相関）され得る。本実施例において、特定の個々のデータセットは適切なクラスタリングを示したが、このデータセットを、さらにより頑強な分析を得るために、連結および／または統合および／または相関させた。連結されたデータを標準化およびクラスタリングし、その結果を生物システムのプロファイルとして記録した。
【０１６３】
全ての個々のプラットフォームから収集されたデータが、情報を伏せられた血清サンプルから異なるグループへクラスタリングされた。プラットフォーム間の差異は、形成されたクラスターの数のみであった。４グループへのクラスタリングを、タンパク質プラットフォームおよび脂質プラットフォームの両方で観察した。最終的に同定されたこれら４グループは、サンプル１〜８、９〜１６、１７〜２４、２５〜３２からなった。
【０１６４】
ＬＣ／ＭＳプロテオームのデータのクラスタリング（すなわち、単一分析的プラットフォーム）を、図４４Ａに示す。図４４Ａは、データの整列および標準化後の、げっ歯類の血清のプロテオームＬＣ／ＭＳ分析のＣＯＳＡクラスタリング分析の一例である。この分析において、少なくとも２８／３２のげっ歯類（サンプルの８７％を超える）で現れた２，９７７個のピークを、クラスタリングに使用した。他の代謝産物プラットフォームであるＣＰＭＧ  ＮＭＲおよび拡散を校正したＮＭＲから得られたデータは、より少ないグループにサンプルをクラスタリングしたが、区分は、脂質分析およびタンパク質分析の間に見られたグループと一致した。
【０１６５】
図４４Ｂは、より頑強に提示された４つのグループ（上記の通り）を示す。図４４Ｂは、全てのプラットフォーム由来の組み合わせデータに適用されたＣＯＳＡクラスタリングの結果である。ＣＰＭＧ  ＮＭＲデータを使用したクラスタリングは、３つのクラスターのみを示したが、ＤＥ  ＮＭＲデータのみを使用したクラスタリングは、２つのクラスターを示した（示さず）。プロテオミクス、脂質ＬＣ／ＭＳ、ＣＰＭＧ  ＮＭＲおよびＤＥ  ＮＭＲからのデータ（合計４８５１個の変量）を組み合わせて、４つの明らかなグループを得た。この分類は、プロテオミクスデータおよび脂質プロファイリングデータの個々の処理の結果と一致した。
【０１６６】
サンプルの情報を明らかにして、これらの方法を用いて区切られたグループが、以下の表Ｉにまとめられるように、異なるげっ歯類のコホートに正確に対応したことを示した。
【０１６７】
表Ｉ：クラスター分析後に提供されるサンプルの識別
【０１６８】
【表５】
2007502992.tif 000017
（血清サンプルの分析に関する代謝疾患のげっ歯類モデルに対する結果および考察−代謝産物およびペプチドピークの同定）
引き続いて同定するために、一変量統計法を、第Ｉ相においてプロファイルされたピークに適用し、げっ歯類の４つのグループの間で異なる量を示したピークを選択した。この第一次統計学的分析は、有意なレベルであるα＝０．０５を有するペアワイズｔ検定からなった。この分析のためのワークフローを、図４５にまとめる。概して、複数の分析プラットフォーム由来の複数のデータセットを、連結、統合および相関させ、次いで標準化した。疾患グループとコントロールグループとの間の統計学的に異なる成分を抽出し、そしてこの差異を定量化した。次いで、罹患グループに薬物を投与することによりシステムに摂動を与え、そして同様の分析を行い、処置グループとコントロールグループとの間の差異を決定した。最後に、同定された全ての成分を、２つの実験間で比較し、生物システムのプロファイルを得た。
【０１６９】
代謝産物およびペプチドの間で観察された差異を示す代表例の抜粋を、表４５Ａに示す。（これらの成分はまた、相関ネットワーク分析においても観察され得る（図４６）。この相関ネットワークにおいて、これらの成分は、それらの間で、ならびに同定された他のペプチドおよび代謝産物と相関を示す。）この提示されたデータを見ると、例えば、２種の血清タンパク質（タンパク質１およびタンパク質２）が、罹患げっ歯類と（ビヒクル処置された）コントロールげっ歯類との間で、差示的かつ対比的に制御されることが見出されたこと、ならびに薬物を用いた処置は、基本的に、罹患タンパク質１のレベルを、コントロール動物のタンパク質１のレベルまで低下させるが、タンパク質２をコントロールよりも約２倍高いレベルまで増加させることが分かる。別の興味深い知見は、脂質レベルの選択における薬物処置の差示的な効果である。
【０１７０】
各分子成分に対して、結果は以下の順で示されることに注意されたい。
１．  罹患＋ビヒクル／コントロール＋ビヒクル．．．．．疾患の影響
２．  罹患＋薬物／罹患＋ビヒクル．．．．．疾患状態に対する薬物処置の効果
３．  罹患＋薬物／コントロール＋薬物．．．．．薬物処置された疾患と処置されたコントロールとの比較
４．  罹患＋薬物／コントロール＋ビヒクル．．．．．薬物処置された疾患と未処置のコントロールとの比較
５．  コントロール＋薬物／コントロール＋ビヒクル．．．．．薬物の「副作用」
これは、５種全ての比較がなされた実例に対する実施例を通じての、げっ歯類血清サンプルの全ての分析に対する表示の順である。
【０１７１】
（血清サンプルの分析に関する代謝疾患のげっ歯類モデルに対する結果および考察−相関ネットワーク分析）
グループの間の成分量レベルの変化に加えて、個々の成分間および個々の成分にわたる相関の試験は、研究される種々の成分の間での重要な関係を明らかにするのに有用である。このような相関の分析は、量レベルの情報に相補的であり、しばしば、疾患または薬物応答の根底にある生化学的なプロセスについての情報を提供する。
【０１７２】
図４６は、８匹の罹患した薬物処置げっ歯類および８匹の罹患したビヒクル処置げっ歯類のペアワイズ比較における、プロテオーム、メタボロームおよび臨床化学データに由来する代表的な相関ネットワーク（疾患状態における薬物の効果）である。凡例において見ることができるように、そのネットワークの成分（または「節点（ｎｏｄｅ）」）は、種々のプラットフォームにより測定された種々のタンパク質、代謝産物または臨床的化学である。この図およびこの図に類似した図における全ての節点は、（ｉ）同定され、そして（ｉｉ）ｐ＜０．０５で±１５％より大きい倍率変化（ｆｏｌｄ−ｃｈａｎｇｅ）を示した、成分である。
【０１７３】
多数の独立したレベルの情報が、このタイプの相関ネットワークにおいて示される。第一に、節点の特定の形状は、成分の測定に使用されたプラットフォームを示す。例えば、図４６において、正方形の形状の節点は、質量分析により測定および同定された（すなわち、配列決定および確認された）ペプチドである。第二に、所定の節点の影は、比較された２グループの血清における量の差異を反映する；これは、標準化されたグループ平均の差異である。第三に、節点のペア間の線は、Ｐｅａｒｓｏｎ係数が０．８０と１．００との間または−０．８０と−１．００との間である相関を示す。負の相関値は、薄い線として示され、一方、正に相関した成分は、図示において、濃い線で視覚的に結合される。一般的な話をすると、正に相関する２つの成分は、統計学的に有意な相互の挙動を反映する。この挙動は、グループにおける全てのサンプルにわたって、第一の成分における変化が、第二の成分における類似の変化に同時に関連することにより特徴付けられる。些細な例は、類似して挙動する同じタンパク質由来のペプチド成分のペア、または同じ分子由来の２つのＮＭＲ共鳴成分であり得る。生物化学的に関連する相関はまた、例えば、同じ生合成経路の一部である代謝産物間で、または同じ高分子構造の成分である実体間でも観察され得る。相関のこのタイプの例を、図４６に示す。ここで、タンパク質２のペプチドは、血清中の多数の脂質成分と高度に正に相関する；この高度の相関は、これらの脂質が、血清中でタンパク質２と同じリポタンパク質起源を共有し得ることを示唆する。負の相関は、例えば同じ経路の一部である成分間で生じ得るが、それらの成分は、酵素の阻害および基質の限界の点により、区別され得る。さらに、過去に関係付けられた生合成の枝分かれ点に入る成分は、互いに負の相関を示し得る。
【０１７４】
構造の全体的なトポロジーは、自己構築（ｓｅｌｆ−ａｓｓｅｍｂｌｉｎｇ）といわれるものであり、そして高度に内部相関している成分のクラスターを反映する。互いに近接する節点は、特に高密度の相互相関を反映している。そのトポロジーは、教師なしの自動様式で生成される。
【０１７５】
このような構造を調査することにより、多数の興味深い知見が明白になる。例えば、脂質２は、処置に際して量がより多い（節点は、最も大きな円構造におけるほぼ４時の方向にある）こと、そしてさらに脂質２が、多くの他の脂質成分と負に相関されることが、見られる。この図は本発明の原理および技術を図示していることが、理解されるべきである；これは、可能な多くのこのような相関のうちの１つである。
【０１７６】
（血清サンプルの分析に関する代謝疾患のげっ歯類モデルに対する結果および考察−ヒートプロット（Ｈｅａｔ  Ｐｌｏｔ）分析）
罹患薬物処置グループおよび罹患ビヒクル処置グループの比較に対する相関情報の代替の見方を、図４７に示す。この「ヒートプロット」は、同定された代謝産物およびペプチドのピークの各ペアに対して計算された相関係数のアレイを示す。成分ピークのペアに対する対角線外のスポットは、そのピーク間の相関係数の符号に対応する（正か負のいずれか）が、その色の強度は、相関の規模に比例する。
【０１７７】
複雑ではあるが、この可視化は、相関の完全なアレイの高速な検査を可能にする。図４７に示したような分析方法に従って成分が分類される場合、異なる成分クラス間の相関は、明白である。例えば、インデックス数２２〜３２のペプチドおよびインデックス数１１０〜１４０の脂質が並ぶ対角線外の領域は、高度に正および高度に負の両方の領域を示す。この場合、正に相関されたペプチド（２２〜２６）は、タンパク質１由来であるが、脂質はトリグリセリドである。倍率変化情報は、図４７には示されないことに留意されたい；影のスケールは、Ｐｅａｒｓｏｎ相関係数を示す。
【０１７８】
（血清サンプルの分析に関する代謝疾患のげっ歯類モデルに対する結果および考察−げっ歯類タンパク質比）
特定のタンパク質は、脂質の代謝に欠くことのできない役割を果たす。このため、これらのタンパク質のいくつかに関連するペプチドのレベルの差異が、本研究の一部として試験された異なるサンプルのコホートにおいて見られることは、驚くことではない。図４８は、異なるグループの間の比として表される、４個のこのようなタンパク質である、タンパク質Ａ（タンパク質１）、タンパク質Ｂ、タンパク質Ｃおよびタンパク質Ｄ（タンパク質２）における差異を示す。６個のトリプシンペプチドを、タンパク質Ａから、１個をタンパク質Ｂから、１個をタンパク質Ｃから、そして２個をタンパク質Ｄから観察した。図４８におけるプロットは、（標準化およびスケール化後の）各グループの中のピーク強度値の平均に基づく、グループの間の比を示す。異なるグループの間で、有意な倍率の変化が存在することが、明白である。特に著しいのは、薬物で処置された罹患げっ歯類とビヒクルで処置された罹患げっ歯類との間、ならびにビヒクルで処置された罹患げっ歯類とビヒクルで処置されたコントロールのサブグループとの間の、タンパク質Ｄの比の変化である。
【０１７９】
（ヒト血清サンプルの分析に関する代謝症候群の研究に対する結果および考察−教師なしクラスタリング）
教師なしクラスタリングを、全ての別個のプラットフォーム（タンパク質、脂質およびＮＭＲ）を用いて誘導されたヒトデータに適用した。代謝疾患のげっ歯類モデルに対して上述したように、これは、サンプルの分類または異なるグループの数の予備知識なしでサンプルを分類させる。ペプチドデータのＣＯＳＡ分析は、サンプルを４つの弱いクラスターに分類した。ＮＭＲの大域的な（ｇｌｏｂａｌ）代謝産物データを用いるクラスタリングは、サンプルを２グループに分割した。一旦サンプル情報が伏せられると、これらの分類は、罹患対コントロールのコホートに対応しないことが明白であった。
【０１８０】
対照的に、脂質データのＣＯＳＡ分析は、２つのクラスターを示す（図４９）。このＣＯＳＡ距離クラスタリングは、７７９個のヒトＬＣ／ＭＳ脂質ピークを使用した。これらのクラスターは、罹患患者に８６％の精度（１２／１４）で対応し、コントロール被験体に８９％の精度（２５／２８）で対応する。多変量分析は、脂質が、罹患サンプルとコントロールサンプルとの間の最も強い識別子であることを示した。
【０１８１】
３つのプラットフォームのうち２つにおいて強いクラスタリングを欠いたことは、クラスタリングが他の因子（例えば、投薬法、性別、年齢または環境）により支配されることを示す。いくつかのプラットフォームに対してＣＯＳＡを用いて誘導されたこのような弱いクラスターを考慮して、他のクラスタリング技術（例えば、Ｋ−Ｍｅａｎｓおよびニューラルネットワーク）について、同じデータセットを用いて調査した。これらの技術は、グループの間の境界にある少数のサンプルを除いては、ＣＯＳＡに類似した結果を与えた。
【０１８２】
（ヒト血清サンプルの分析に関する代謝症候群の研究に対する結果および考察−代謝産物およびペプチドピークの同定）
げっ歯類の研究において見られたように、サンプルタイプ間のレベルにおいて有意に異なるピークを強調することにより、潜在的に興味深いピークが発見され得る。この研究の目的のために、ヒトサンプルを最初に２グループ（１４人の罹患患者および２８人のコントロール被験体）に分割した。２つのサンプルのｔ検定を、この２グループの間の平均差異を試験するために、各ピークに対して行った。この結果を、同定のために入力されるピークに対する一覧表にした。
【０１８３】
脂質プラットフォームに対して、罹患患者とコントロールの被験体との間で差異を示したピークのサブセットを、参照データベースおよび標的化ＭＳ／ＭＳ法を用いて同定した。一般的に、ピーク同定の際に、罹患患者における特定の脂質分子のレベルが、コントロールの被験体におけるこれらの脂質のレベルから有意に異なることを発見した。興味深いことに、以下のげっ歯類／ヒト比較研究において見られるように、多くのこれらの脂質レベルもまた、コントロールのげっ歯類と比較して罹患げっ歯類において有意に異なる。
【０１８４】
さらに、ヒトタンパク質の一覧表を、「ショットガン」タンデム型質量分析（ＭＳ／ＭＳ）法を用いて、本研究の一部として同定した。ショットガンＭＳ／ＭＳにより配列決定するための、ＭＳプロファイリング段階の間に選択されたピークのセットと、２グループのヒトサンプルの血清間で有意なレベルの差異を示したピークのセットとの間に、オーバーラップは存在しなかった。
【０１８５】
（げっ歯類サンプルとヒトサンプルとの比較に対する結果および考察）
本研究のこの部分において、目的は、罹患のビヒクル処置げっ歯類およびコントロールのビヒクル処置げっ歯類由来の血清における脂質成分と、罹患したヒトおよびコントロールのヒト由来の血清における対応する脂質とを比較することであった。薬物処置しないグループを、これらの分析に含めた。ＬＣ／ＭＳ血清脂質プラットフォーム由来のデータを使用し、具体的には、両方の種に共通の５７１個のＬＣ／ＭＳピークを使用した。図５０は、この分析のワークフローを示す。
【０１８６】
このフレームワークにおいて、２つの問題点を扱った。第一の問題点は、げっ歯類の測定に基づくヒトサンプルのクラスタリングおよび分類における精度に関する。第二の問題点は、２つの種にわたる脂質の量の変化および相関の比較に関する。
【０１８７】
（げっ歯類サンプルとヒトサンプルとの比較に対する結果および考察−クラスタリングおよび分類）
２つの種に共通であった５７１個のピークの間で、３６６個において、２つのげっ歯類グループの間で有意な平均の変化（両側のペアワイズｔ検定を用いて０．０５の有意性レベル）が存在した。探索工程として、罹患したビヒクル処置されたげっ歯類と一緒になった罹患したヒト、およびコントロールのビヒクル処置されたげっ歯類と一緒になったコントロールのヒトからなるデータにおいて、天然のクラスターが存在するか否かを決定するために、このセットの３６６個のピークを使用した。この分析の結果を、図５０Ａに示す。具体的には、ヒト血清サンプルのＣＯＳＡ分析の結果を示す。この結果において、分類のために使用されたインプットのデータセットは、罹患げっ歯類モデルから選択された３６６個の脂質ピークからなった。この図は、罹患サンプルおよびコントロールサンプルによく対応した２つの主要なグループを示す：２８個のコントロールのヒトのうちの２７個、および８個全てのコントロールのげっ歯類が、１つのグループに属し、そして１４個のうち１１個の罹患したヒトおよび全ての罹患したげっ歯類が、第二のグループに属する。ヒトの診断が未知である場合、２つのげっ歯類グループにおいて形成されるクラスターを検査することにより、高精度で推論され得ることが、この分析から結論付けられる。
【０１８８】
分類する目的のために、サポートベクトルマシン（ＳＶＭ）線形分析器を使用した。ここで、３３６個のげっ歯類脂質測定がモデル構築セットの機能を果たし、そして対応する３６６個のヒト脂質測定が独立した試験セットとしての機能を果たした。正確に分類されたヒトサンプルのパーセンテージは、図５１に見られるように、７６％（４２個のサンプルのうちの３２個）と９３％（４２個のサンプルのうちの３９個）の間を変動した。図５１は、脂質ピーク数の関数としての、ＳＶＭ線形分析器の成功率を示す。この分析において、げっ歯類のデータをモデル構築のために使用し、成功率は、ｌｅａｖｅ−ｏｎｅ−ｏｕｔ法において正確に分類されたげっ歯類のパーセンテージである。またこの分析において、ヒトデータは試験セットとして使用され、そして成功率は、げっ歯類モデルにより正確に分類されたヒトのパーセンテージである。分類およびピーク削減手順のさらなる調査が、罹患げっ歯類モデルがヒトにおける代謝疾患に対する優れたモデルであるという確証をもたらし得る。
【０１８９】
（げっ歯類サンプルとヒトサンプルとの比較に対する結果および考察−共通成分）
両方の種に共通であった５７１個のＬＣ／ＭＳ脂質ピークの比較は、この５７１個の脂質ＬＣ／ＭＳピークのうちの１９５個に対して、両方の種において、罹患グループとコントロールグループとの間の有意な平均の差異（両側のペアワイズｔ検定を用いて０．０５の有意性レベル）が存在したことを示した。これらの１９５個のピークのうち、１８５個が両方の種において同じ傾向（罹患対コントロールにおける、より多い血清量またはより少ない血清量）を示した。さらに、０．７より大きいＰｅａｒｓｏｎ相関係数の絶対値を用いると、脂質ピークのペア間の多数の相関がヒトサンプルおよびげっ歯類サンプルの両方に存在した。これは、量の差異が保存されているだけでなく、これらの脂質レベルの制御に関わる根底にある機構が、種間で存在されているようであり得ることを示す。この結果の抜粋を、図５２にまとめる。
【０１９０】
より具体的には、図５２は、ヒトおよびげっ歯類種にまたがる脂質量の変化および相関の比較を示す。この図において、大きな円は要素からなり、この要素の各々が異なるＬＣ／ＭＳ脂質ピークを示す。その要素の影は、罹患サンプル対コントロールサンプルにおける脂質の相対量に対応する。相対量は、標準化されたグループの平均の差異である。１９５個のこのような要素が存在し、全てがｐ＜０．０５の脂質を示す。外側の大きな円は、罹患げっ歯類グループ対コントロールげっ歯類グループの比較を示し、一方内側の同心円は、罹患ヒトグループ対コントロールヒトグループの比較を示す。図において要素のペアを結ぶ線は、両方の種において存在する、Ｐｅａｒｓｏｎ係数が｜Ｃｉｊ｜＞０．７０の相関である。
【０１９１】
（要旨および結論）
動物被験体またはヒト被験体由来の、情報が伏せられた血清サンプルの代謝産物分析およびタンパク質分析を行った。これは、それらの血清の代謝産物プロファイルおよびタンパク質プロファイルに基づいたサンプルの分類を可能にした。クラスタリング分析を用いて同定されたグループは、１００％の精度で動物被験体の表現型カテゴリーを反映し、そして高度の精度（＞８０％）でヒト被験体の表現型カテゴリーを反映した。その後の分析は、被験体を差別化する多くの分子成分を同定した。
【０１９２】
これらの独立した方法は、元来、参考になる。さらに、相関ネットワークを用いて結合された場合、疾患または薬物応答の根底にある生化学プロセスの詳細が理解され始める。より興味深い結果のうちの１つは、コントロールのげっ歯類から罹患げっ歯類を差別化する分子成分が、コントロールの被験体から罹患ヒトを差別化する分子成分と非常に類似していることである。本研究により生成されたデータの財産は、プロテオミクス、メタボロミクスおよびインフォマティクス技術の統合されたプラットフォームを利用したシステム生物学的アプローチの長所を示す。
【０１９３】
（本実施例において使用される専門用語／用語）
（略語および用語）
ＣＯＳＡ：属性のサブセットに関するオブジェクトのクラスタリング（Ｃｌｕｓｔｅｒｉｎｇ  Ｏｂｊｅｃｔｓ  ｏｎ  Ｓｕｂｓｅｔｓ  ｏｆ  Ａｔｔｒｉｂｕｔｅｓ）
ＣＰＭＧ  ＮＭＲ：Ｃａｒｒ−Ｐｕｒｃｅｌｌ−Ｍｅｉｂｏｏｍ−ＧｉｌｌスピンエコーＮＭＲ
ＤＥ  ＮＭＲ：拡散を校正したＮＭＲ（Ｄｉｆｆｕｓｉｏｎ−ｅｄｉｔｅｄ  ＮＭＲ）
ＬＣ：液体クロマトグラフィー
ＭＳ／ＭＳ：タンデム型質量分析
ＭＳ：質量分析
ＮＭＲ：核磁気共鳴。
【０１９４】
（タンパク質の専門用語）
ショットガン配列決定：「データ依存性」の機器様式において取得されたタンデム型質量スペクトル（ＭＳ／ＭＳ）を使用して、ペプチドの配列情報を得る方法であって、この方法により、機器は、可能な限り多くのペプチドピークについてのＭＳ／ＭＳスペクトルを測定するように構成される。この様式において、機器は、ペプチドピーク信号の最初の調査的スキャンから構成される繰返しスキャンサイクルを実行し、最も強い３〜４つを選択し、その後、選択したピークの各々についてＭＳ／ＭＳスキャンする。
標的化配列決定：特定のペプチドピークについて取得されたタンデム型質量スペクトル（ＭＳ／ＭＳ）を使用して、ペプチドの配列情報を得る方法。
【０１９５】
（実施例６  システム生物学的アプローチ：ヒトの心血管系疾患）
この実施例の目的は、ヒト心脈管系疾患の患者を、健康な被験体と区別する、血漿代謝産物を明らかにすることであった。研究に先立って、被験体のサンプルを、罹患カテゴリーまたはコントロールカテゴリーのいずれかに分類した（心脈管系疾患、および対応するコントロール被験体からの血漿サンプル）。ＮＭＲ、ＬＣ／ＭＳおよびＧＣ／ＭＳ技術、ならびにデータ前処理ソフトウェアを使用する、いくつかのメタボロミクスプラットフォームを、８０の血漿サンプルの比較研究に適用した。メタボロミクスプロファイリングプラットフォームは、最初は同定されていない、数百のスペクトルのピークを含むデータベースを生成する。その代わりに、統計的に有意なピークを決定した。これらの全体を、分析の第２相において、データベース、追加のＭＳ／ＭＳデータ、および専門家の解釈を使用して、同定のために目印を付けた。メタボロミクスデータセットの一変量および多変量の統計分析は、研究被験体の２つのグループ間で有意に異なった測定特徴を明らかにした。このプロジェクトの第２相の開始前に、疾患の重篤度の臨床的インデックスに基づく、罹患した被験体のさらなる分類を使用し、そして、任意の測定特徴が、罹患グループにおいて心脈管系疾患の重篤度と相関する場合に、さらなる統計分析を行なった。多数の特徴が、１つ以上の分析において有意性を示し、そして、多数の特徴が同定された。次いで、相関ネットワークを構築して、同定された、有意な代謝産物の間の統計的かつ生物学的な関係性を可視化した。
【０１９６】
（目的）
この研究の目的は、心脈管系疾患の患者、および対応するコントロール被験体から採取した血漿サンプル間での分子の差異として生物マーカー分子を同定することであった。
【０１９７】
（研究設計）
研究を、２つの相で実行した：
第Ｉ相：メタボロミクスプラットフォームを用いて、男性の心脈管系疾患患者（４０サンプル、平均年齢５３．４歳）、または年齢が対応するコントロール被験体（４０サンプル、平均年齢５１．６歳）のいずれかに由来するものと説明される８０の血漿サンプルを比較的にプロファイリングした。分析プラットフォームは、ＣＰＭＧ  ＮＭＲ、拡散を校正したＮＭＲ、ＧＣ／ＭＳ、脂質ＬＣ／ＭＳ、およびアミノ酸／大域  ＬＣ／ＭＳであった。ソフトウェアアルゴリズムを使用して、生データから、スペクトルおよびクロマトグラフィーのピーク情報を抽出した。さらなる前処理を行なって、比較統計分析のために、各プラットフォームからのデータベース間でピークを揃えた（すなわち、ＬＣ−ＭＳおよびＧＣ／ＭＳについてのクロマトグラフィー保持時間の整列）。ピークは、統計的有意性に基づく識別のために印を付けられるまで、同定されないままであった。識別行為は、２つの実験グループ間で異なるレベルの量を有したピークについて開始した。
【０１９８】
第ＩＩ相：このプロジェクトの第２相の開始前に、疾患の重篤度の臨床的インデックスに基づく、罹患した被験体のさらなる分類を行い、そして、任意の測定特徴が、罹患グループにおいて疾患の重篤度と相関した場合に、さらなる統計分析を行なった。可能な場合、有意と思われる特徴についてさらなる識別情報を得た。次いで、相関ネットワークを構築して、同定された有意な代謝産物間の統計的かつ生物学的な関係性を可視化した。
同定された、有意な代謝産物の間の統計的かつ生物学的な関係性を可視化した。
【０１９９】
（方法の概要）
広範囲の代謝産物の比較的プロファイリングを可能にする多数の分析的手法を使用した。このサンプルを、いくつかの分析的手法を使用して分析し、そして、未識別のピークについて統計を実施した。使用した方法が、以下に列挙され、簡単に説明される。
【０２００】
（ｉ）ＣＰＭＧ  ＮＭＲ：１００μＭより多い濃度における低分子量代謝産物の増強されたＮＭＲ測定（例えば、アミノ酸、アミノ酸代謝産物、有機酸、糖）
（ｉｉ）ＧＣ／ＭＳ：広範囲の代謝産物の分類のプロファイリングのために設計された大域的方法（ｇｌｏｂａｌ  ｍｅｔｈｏｄ）（例えば、アルコール、アルデヒドおよびシクロヘキサノール、アミノ酸、アシルアミノ酸、スクシニルアミノ酸、アミン、芳香族化合物、脂肪酸（Ｃ６より大きい）、有機酸、ホスホ有機酸（ｐｈｏｓｐｈｏ−ｏｒｇａｎｉｃ  ａｃｉｄ）、糖、糖酸、糖アミン、糖リン酸）
（ｉｉｉ）脂質ＬＣ／ＭＳ：脂質および非極性代謝産物のプロファイリングに最適化されている（例えば、リゾリン脂質、リン脂質、コレステロールエステル、ジアシルグリセロール、トリアシルグリセロール）
（ｉｖ）アミノ酸／大域  ＬＣ／ＭＳ：アミノ酸および極性代謝産物のプロファイリングに最適化されている。血液抗凝固因子として使用されるクエン酸の存在に起因して、このプラットフォームは、有用なデータを得られず、そして、第ＩＩ相において使用しなかった。
【０２０１】
（ｖ）拡散を校正したＮＭＲ：リポタンパク質に関連する代謝産物の増強された測定。プロファイルされたピークは、多くの脂質部分からの信号の複合物であり、それゆえ、非特異的である。固有に同定された分子実体が、生物マーカーとして好ましかったので、この方法は、第ＩＩ相において実行しなかった。
【０２０２】
上記の分析の各々は、サンプルあたり数百〜数千のピークを含む生のデータセットを得た。全サンプルセットにわたって代謝産物のピーク情報の比較的な分析を可能にするために、いくつかのアルゴリズムを、ピークの検出および信号の統合のために、各生データファイルに適用した。次に、ＬＣ／ＭＳ技術およびＧＣ／ＭＳ技術についての保持時間の観点、または、ＮＭＲ技術についての化学シフトの重要でない差異の点から生じ得る、ピーク位置の重要ではないシフトについて比較するために、アルゴリズムを、これらのピークを「整列する」ために使用した。この処理の結果として、プロファイル中の各代謝産物のピークに、ピーク識別番号（または、インデックス数）を割り当てた。この同じ識別番号を使用して、全ての他のサンプルからのプロファイルにおいて見られる類似のピークを説明し、そして、それゆえ、統合したピーク強度の比較分析が可能となった。
【０２０３】
各プラットフォームからのデータの一変量および多変量の統計分析の後に、罹患と健康な被験体とを区別した代謝産物を、適用された統計学によりランク付けしたものとしての第ＩＩ相における識別のために列挙した。
【０２０４】
（一変量の結果）
データの整列および標準化に続いて、誤った発見率のためのコントロールを用いた一変量の等分散ｔ検定を、この研究において使用した全ての生物分析学的プラットフォームからの同定された代謝産物の分析物について、行なった。結果は、Ｂｅｎｊａｍｉｎｉ−Ｈｏｃｈｂｅｒｇアプローチを使用して、１０％の誤った発見コントロールに基づき０．０５未満の調整されたｐ値を有する２４の分析物を示した。
【０２０５】
（多変量の結果）
罹患サンプルとコントロールサンプルとを分類し得る、スペクトルのピークのセットを見出すための複数の分析物のアプローチをまた実行した。文献において、サンプルのグループを分離し得る１つ以上の分子成分から構成される生物マーカーを見出すというこの問題は、「分類上の問題」と呼ばれる。この場合、確信的かつ固有に識別された分析物のみを使用した；９４のこのような分析物が、分析時に存在した。この数は、同位体、付加体、冗長的な１ＮＭＲ共鳴ピークなどを含まず、これらはまた、同定されていてもよい。簡単に述べると、分類の課題は、最も多くの情報を与える分析物の最少数から構成される、複数の分析物の生物マーカーを決定することである。
【０２０６】
１つ以上の成分から構成される生物マーカーを考慮する際、多数の点が考慮された。これらは、どの分析物のサブセットがマーカー中に含むための最適なものでるか；最終的な生物マーカーが、目前にあるサンプルセットをいかに良く正確に分類するのか；そして、最終的な生物マーカーが、独立したサンプルセットからいかに良く正確にサンプルを分類するのか、を決定することを含む。上記の項目に加えて、生物マーカーを構成する成分の生化学的関連性がまた重要であり、最終的な生物マーカーについての実用的な診断アッセイを開発する実行可能性もまた重要である。後者を心に留めて、最良の予測性能基準を達成する分析物の最小の最適な数を決定した。図５３は、この分析の工程の概略を図示する。一般に、複数の分析的プラットフォームからの複数のデータセットが、連結、統合、および相関され、次いで、標準化される。このデータは、さらに、生物システムのプロファイルを得るために、教師なしクラスタリング分析を通じて分析される。複数の分析物の生物マーカーを構成する方法論の概要は、以下に示される。
【０２０７】
疾患サンプルとコントロールサンプルとを最良に分離するスペクトルピークの最小の最適なサブセットを決定するために、回帰特徴排除（Ｒｅｃｕｒｓｉｖｅ  Ｆｅａｔｕｒｅ  Ｅｌｉｍｉｎａｔｉｏｎ）として知られるアプローチを使用する。このアプローチは、以下のとおりに進行する。
【０２０８】
１．Ｎ個の成分の入力（すなわち、Ｎ個のスペクトルピーク）として受け入れる、「分類アルゴリズム」を選択し、そして、（ｉ）Ｎ個の成分の線形組み合わせにより達成される、コントロールサンプルと疾患サンプルとの分離の成功（特異性および感受性により測定される）、および（ｉｉ）分類へのその寄与に基づいたＮ個の入力成分のランク付けを返す。
【０２０９】
２．分類アルゴリズムへの入力として、（整列され、標準化され、かつ前処理された）全ての分析物を許可する。
【０２１０】
３．入力としてこれらの成分を使用して、アルゴリズムを実行して、使用される入力分析物の線形組み合わせ上に集中させ、コントロールサンプルと疾患サンプルとを分類する。
【０２１１】
４．各分析物についてのランク付け基準（「重み」）を記録する。この重みは、アルゴリズムにより決定されるような、入力成分の線形組み合わせにおける係数である（最終的な重みは、実際の平均重みであり；複数回の相互検証の繰返しにわたって平均化される）。
【０２１２】
５．相互検証法（以下に議論される）、ならびに、この相互検証検定についての標準誤差を使用して、コントロールサンプルと疾患サンプルとを分類する際に、スペクトルピークのこの組み合わせの「相互検証」性能を算出する。
【０２１３】
６．最も低い重みを有する分析物を除外する。
【０２１４】
７．１つの分析物のみが残るまで、工程３〜工程６を繰り返す。
【０２１５】
８．コントロールサンプルと疾患サンプルとの分離において最高の成功を達成するために必要とされる分析物の最小数を決定する；この生物マーカーは、分析物値の線形組合せから構成され、この組み合わせの係数は、各分析物に対応する重みである。
用語「回帰特徴排除（Ｒｅｃｕｒｓｉｖｅ  Ｆｅａｔｕｒｅ  Ｅｌｉｍｉｎａｔｉｏｎ）」は、工程３〜６の各繰り返しについて１つのスペクトルピークによるスペクトルピークの一覧表の連続的な剪定を反映する。
【０２１６】
本研究において、１つの分類アルゴリズムを適用した。このアルゴリズムは、「Ｌｏｇｉｓｔｉｃ  Ｃｌａｓｓｉｆｉｅｒ」（Ａｎｄｅｒｓｏｎ，１９８２）と呼ばれる、当該分野の技術水準であるアプローチを含む。この方法は、手書きおよび生物測定のパターンの認識にその起源がある。この方法は、冗長性を避け、生物マーカーの大きさを最小にするための望ましい特性である。低い相互相関を有する成分を含む最終的な生物マーカーを選択するために設計される。この技術の一般的な原理は公知であるが、現行の分析は、以前に議論された特定の生物分析学的プロファイリングのプラットフォームから算出されたデータを用いて機能するように最適化している。
【０２１７】
この節において概説されるプロセスに適用された、２つの異なる性能試験がある。
【０２１８】
１．「相互検証性能」は、利用可能なサンプルのサブセットに基づいて構築され、演繹的かつ意図的に残された残りのサンプルについて試験された、生物マーカーの分類の成功例である（Ｈａｓｔｉｅ，２００１）。本研究についての代表的な状況は、無作為に選択された３４の罹患サンプルと、３４のコントロールサンプルについてのみ基づいて生物マーカーを構築し、そして、除外された残りの６の罹患サンプルおよび６のコントロールサンプルを分類するのに得られた生物マーカーの性能（分類の成功）を試験することである。この処理は、何度も首尾よく反復され、そして、無作為に選択された６＋６のサンプルの異なるセットが「残される」。報告された生物マーカーについての「相互検証性能」は、多くのこのような順列の平均的な性能であり；代表的に、１０回の相互検証ラウンドが使用される。
【０２１９】
相互検証の目的は、比較的制限された数の個々のサンプルの利用可能性により提示される制限の範囲内で、生物マーカーの一般化を評価することである、ということに注意することが重要である。異なる集団の患者からの独立したサンプルがなければ、相互評価性能は、サンプルの独立した試験セットについての生物マーカーの性能の推定である。このような外挿は、多くの順列および利用可能なサンプルのサブセットの組み合わせについて、生物マーカーの性能を測定することによって可能となる；この処理は、より多くのサンプルが利用可能である状況を効果的に模倣する。
【０２２０】
２．「順列性能」は、サンプル標識が無作為に順列された場合の、多変量の生物マーカーの選択アルゴリズムの性能である。これは、多くのこのような無作為な順列にわたって生じ、そして、その平均性能が報告される。訓練セットに対してオーバーフィットしない頑強な分類器は、約５０％の順列性能（すなわち、機会性能）を得るはずである。
【０２２１】
（結果および考察）
これらの分類法の結果は、図５４に図式的に示される。１５個の分子成分の生物マーカーセットを、ヒト心脈管系疾患のプロファイルの一部として同定した。生物マーカーセットのこれらの分子成分は、図５６に示されるような多変量統計分析法、ならびに１つ以上のタイプの測定技術のためのデータベース、および１つ以上の生体分子成分についてのデータベースを含む、複数のデータベースの統合を使用することよって発見された。この方法論的なアプローチは、８０個のサンプルを分類し得る生物マーカーセットを生成するために首尾よく使用された。図５５は、これらの生物マーカーを使用して、疾患グループおよびコントロールグループのメンバーとしての各被験体の分類を示す。９３％の選択性、および９４％の特異性が得られた。
【０２２２】
この実施例で使用される略語は、適切な場合、実施例５において使用したものと同じである。
【０２２３】
本明細書中で開示される、特許文献および科学技術文献の各々は、あらゆる目的のために、本明細書中で参考として援用される。
【０２２４】
本発明は、特定の実施形態を参照して具体的に示され、記載されてきたが、本発明の精神、本質的な特徴または範囲から逸脱することなく、その形態および細部に種々の変更がなされ得ることは、当業者により理解されるべきである。従って、上述の実施形態は、全て例示的な観点であり、本明細書中に記載される発明を限定するものではないと考慮されるべきである。本発明の範囲は、従って、上記の説明ではなく、添付の特許請求の範囲により示され、そして、特許請求の範囲の意味および等価物の範囲内に含まれる全ての変更は、従って、その中に包含されることが意図される。
【図面の簡単な説明】
【０２２５】
本発明における上記および他の目的、特徴、および利益は、上記の種々の例示的な実施形態の記載を添付の図面と共に読むと、より完全に理解され得る。図面おいて、一般に類似の参照文字は、異なる視野においても同じ部分を指す。図面は、必ずしも本発明の原理を拡大縮小および強調するものではなく、その代わりに、本発明の原理を例示する。
【図１】図１は、生物システムのプロファイルを作製するための、ゲノムデータセット、プロテオミクスデータセット、メタボロミクスデータセット、および臨床データセットの統合を例示する、概略的なフローチャートである。
【図２】図２は、本発明の例示的な実施形態に従って複数のデータセットに適用される種々の分析工程および処理工程のフローチャートである。
【図３】図３は、ＡｐｏＥ３−Ｌｅｉｄｅｎトランスジェニックマウスの遺伝子発現実験の実験設計を例示する。
【図４】図４は、遺伝子発現実験についての有意なプロットを例示する。
【図５】図５は、４つの肝臓フラクションから選択された１０５９個のペプチドのピークの有意なプロットを例示する。
【図６】図６は、合成ＧＩＳＴデータの実験のためのブロック設計を例示する。
【図７】図７は、合成ＧＩＳＴデータセットの変種１について、散布図および正規分布のプロットを例示する。
【図８】図８は、合成ＧＩＳＴデータセットの変種２について、散布図および正規分布のプロットを例示する。
【図９】図９は、合成ＧＩＳＴデータセットの変種３について、散布図および正規分布のプロットを例示する。
【図１０】図１０は、合成ＧＩＳＴデータセットについての有意なプロットを例示する。
【図１１】図１１は、生物学的サンプルから算出される遺伝子発現データの処理を記載するフローチャートを例示する。
【図１２】図１２は、生物学的サンプルから算出されるタンパク質データの処理を記載するフローチャートを例示する。
【図１３】図１３は、生物学的サンプルから算出される代謝産物データの処理を記載するフローチャートを例示する。
【図１４】図１４は、２つ以上の生体分子成分のタイプから算出される複数のデータセットの統合を記載するフローチャートを例示する。
【図１５】図１５は、ｍＲＮＡの存在量を明らかにする遺伝子発現分析を例示する。
【図１６】図１６は、遺伝子発現分析から選択されたグループについての結果を例示する。
【図１７】図１７は、遺伝子発現分析から選択されたグループについての結果を例示する。
【図１８】図１８は、血漿サンプルからのタンパク質のＬＣ／ＭＳ全イオンクロマトグラムの強度プロットを例示する。
【図１９】図１９は、血漿サンプルからのタンパク質のＬＣ／ＭＳプロファイリングからの全イオンクロマトグラムを例示する。
【図２０】図２０は、５匹のトランスジェニックマウスおよび５匹の野生型マウスの消化した肝臓タンパク質から得た、ＬＣ／ＭＳクロマトグラムを例示する。
【図２１】図２１は、トランスジェニックマウスおよび野生型マウスからの血漿から抽出した代謝産物の１Ｈ  ＮＭＲスペクトルを例示する。
【図２２】図２２は、トランスジェニックマウスおよび野生型マウスについてのＬＣ／ＭＳを使用して記録した、血漿脂質のマスクロマトグラムを例示する。
【図２３】図２３は、個々の生体分子成分のタイプ間の比較のための単一因子スペクトルを形成するために、標準化され、かつ連結された、個々の遺伝子スペクトル、タンパク質スペクトル、および代謝産物スペクトルを例示する。
【図２４】図２４は、主成分および判別（「ＰＣ−ＤＡ」）統計分析の結果として生じる、野生型マウスおよびトランスジェニックマウスのデータをクラスタリングすること例示する。
【図２５】図２５は、有意差を示すペプチドの差異因子スペクトルを例示する（ｍ／ｚ値  １３６６に注目）。
【図２６】図２６は、ＬＣ／ＭＳ／ＭＳを使用して記録したマウスの血漿からのペプチド（ｍ／ｚ値  １３６６）の質量スペクトルおよび配列を例示する。ＭＳ／ＭＳのスペクトルから推定したこのペプチドは、ヒトのアポリポタンパク質Ｅ３の配列中の残基５７〜７９として同定される。
【図２７】図２７は、生体分子成分のタイプ間の相関ネットワークを例示する。
【図２８】図２８は、相関ネットワークの関連性と公開された情報との間の既知の関係を示すマップを例示する。
【図２９】図２９は、システム生物学分析から生じ得る生物マーカー（「マーカー」）または治療剤の観点から、代表的な「提供物」または「送達可能物」を例示する。
【図３０Ａ】図３０Ａは、ＡｐｏＥ３−Ｌｅｉｄｅｎトランスジェニックマウス実験の実験設計を例示する。
【図３０Ｂ】図３０Ｂは、ｃＤＮＡのマイクロアレイデータの散布図を例示する。
【図３１Ａ】図３１Ａは、１０個のサンプルについて消化した肝臓タンパク質のフラクションのＬＣ／ＭＳクロマトグラムを例示する。
【図３１Ｂ】図３１Ｂは、トリプシンペプチドのプロファイルのクラスタリング分析を例示する。
【図３１Ｃ】図３１Ｃは、肝臓タンパク質データの因子スペクトルを例示する。
【図３２Ａ】図３２Ａは、肝臓脂質データセットの主成分分析から得られるクラスタリングを例示する。
【図３２Ｂ】図３２Ｂは、肝臓脂質データセットの因子スペクトルを例示する。
【図３３Ａ】図３３Ａは、３つの生体分子成分のタイプからのデータに基づく包括的なシステム分析を例示する。１．０の相対存在量は、１００％である（図３３Ａ−ｍＲＮＡ）。
【図３３Ｂ】図３３Ｂは、３つの生体分子成分のタイプからのデータに基づく包括的なシステム分析を例示する。１．０の相対存在量は、１００％である（図３３Ｂ−タンパク質）。
【図３３Ｃ】図３３Ｃは、３つの生体分子成分のタイプからのデータに基づく包括的なシステム分析を例示する。１．０の相対存在量は、１００％である（図３３Ｃ−脂質）。
【図３４】図３４は、血管内の高脂血症およびアテローム性動脈硬化症を例示する概要図である。
【図３５】図３５は、全血漿の並列タンパク質代謝産物プロファイリング（ｐａｒａｌｌｅｌｌ  ｐｒｏｔｅｏ−ｍｅｔａｂｏｌｉｃ  ｐｒｏｆｉｌｉｎｇ）のスキームを例示する。
【図３６】図３６は、野生型マウス血漿サンプル（ＷＴ）およびトランスジェニックマウス血漿サンプル（ＴＧ）についてのＮＭＲスペクトルを例示する。
【図３７】図３７は、トランスジェニックマウス（三角で表示）、および野生型（またはコントロール）マウス（丸で表示）に関するＮＭＲデータのクラスタリングを示すＰＣ−ＤＡスコアのプロットを例示する。
【図３８】図３８は、種々の代謝成分を示す多数の線によって特徴付けられた差異スペクトルを例示する。
【図３９】図３９は、２００〜１７００ｍ／ｚ以上の質量範囲にて得た質量スペクトルを有するＬＣ次元において４段階勾配によって分析された、トランスジェニック（ＴＧ）マウスおよび野生型（ＷＴ）マウスからのタンパク質を除いた脂質フラクションの全イオンクロマトグラム（ＴＩＣ）を例示する。
【図４０】図４０は、トリプシンペプチドから得られたトランスジェニック（ＴＧ）マウスおよび野生型（ＷＴ）マウスのタンパク質フラクションの全イオンクロマトグラムを例示する。
【図４１】図４１は、野生型（ＷＴ）マウスおよびトランスジェニック（ＴＧ）マウスについてのＰＣ−ＤＡクラスターを示す、スコアのプロットを例示する。
【図４２】図４２は、タンパク質成分および代謝産物成分の差異因子スペクトルを例示する。
【図４３】図４３は、データ分析のワークフローの概略図を例示する。
【図４４】図４４は、複数のプラットフォームの教師なしクラスタリング分析についてのワークフローを例示する。
【図４４Ａ】図４４Ａは、４つの異なるクラスターを明らかにするＬＣ／ＭＳのプロテオミクスデータのＣＯＳＡ教師なしクラスタリングを例示する。
【図４４Ｂ】図４４Ｂは、結び付けられた複数のデータセットのＣＯＳＡ教師なしクラスタリングを例示する。
【図４５】図４５は、別のサンプルとは異なる１つのサンプルの成分の選択および比較についてのワークフローを図解する。
【図４５Ａ】図４５Ａは、単変量統計的手法を使用して同定されたラットのグループ間での選択されたタンパク質、脂質、および代謝産物の差異に関する代表的なグラフを例示する。
【図４６】図４６は、薬物処置された罹患げっ歯類とビヒクル処置された罹患げっ歯類との間の比較のための相関ネットワーク（疾患に対する薬物効果）を例示する。
【図４７】図４７は、薬物処置された罹患げっ歯類、およびビヒクル処置された罹患げっ歯類において、成分のペア間の相関に関する強度プロットの可視化（疾患に対する薬物効果）を例示する。
【図４８】図４８は、特定のタンパク質由来のペプチドと関連する、各グループ内のピーク強度値の平均値に基づいて（標準化およびスケール化後）、グループ間の比率を示すプロットを例示する。
【図４９】図４９は、ヒトのＬＣ／ＭＳ脂質ピークを使用するＣＯＳＡ距離クラスタリングを例示する。
【図５０】図５０は、非ヒトサンプルデータとのヒトサンプルデータの比較および相関についてのワークフローを例示する。
【図５０Ａ】図５０Ａは、ヒト血清サンプルのＣＯＳＡ分析の結果を例示し、分類のために使用される入力データセットは、ヒト疾患のげっ歯類モデルから選択された３６６個の脂質のピークから構成された。
【図５１】図５１は、ＳＶＭ線形分類器の成功率を、脂質ピークの数の関数として例示する。
【図５２】図５２は、ヒト種およびげっ歯類種間の、脂質の存在量の変化および相関の比較を例示する。
【図５３】図５３は、いくつかのデータセットの分析に関するワークフローを例示する。
【図５４】図５４は、生物マーカーについての分析物の選択をグラフ表示で例示する。
【図５５】図５５は、グループ分けしたサンプル中の１５個の分析物の生物マーカーの性能を例示する。
【図５６】図５６は、図５５からの分析物のリストを例示する。
【図１】
2007502992.tif 000018
【図２】
2007502992.tif 000019
【図３】
2007502992.tif 000020
【図４】
2007502992.tif 000021
【図５】
2007502992.tif 000022
【図６】
2007502992.tif 000023
【図７】
2007502992.tif 000024
【図８】
2007502992.tif 000025
【図９】
2007502992.tif 000026
【図１０】
2007502992.tif 000027
【図１１】
2007502992.tif 000028
【図１２】
2007502992.tif 000029
【図１３】
2007502992.tif 000030
【図１４】
2007502992.tif 000031
【図１５】
2007502992.tif 000032
【図１６】
2007502992.tif 000033
【図１７】
2007502992.tif 000034
【図１８】
2007502992.tif 000035
【図１９】
2007502992.tif 000036
【図２０】
2007502992.tif 000037
【図２１】
2007502992.tif 000038
【図２２】
2007502992.tif 000039
【図２３】
2007502992.tif 000040
【図２４】
2007502992.tif 000041
【図２５】
2007502992.tif 000042
【図２６】
2007502992.tif 000043
【図２７】
2007502992.tif 000044
【図２８】
2007502992.tif 000045
【図２９】
2007502992.tif 000046
【図３０Ａ】
2007502992.tif 000047
【図３０Ｂ】
2007502992.tif 000048
【図３１Ａ】
2007502992.tif 000049
【図３１Ｂ】
2007502992.tif 000050
【図３１Ｃ】
2007502992.tif 000051
【図３２Ａ】
2007502992.tif 000052
【図３２Ｂ】
2007502992.tif 000053
【図３３Ａ】
2007502992.tif 000054
【図３３Ｂ】
2007502992.tif 000055
【図３３Ｃ】
2007502992.tif 000056
【図３４】
2007502992.tif 000057
【図３５】
2007502992.tif 000058
【図３６】
2007502992.tif 000059
【図３７】
2007502992.tif 000060
【図３８】
2007502992.tif 000061
【図３９】
2007502992.tif 000062
【図４０】
2007502992.tif 000063
【図４１】
2007502992.tif 000064
【図４２】
2007502992.tif 000065
【図４３】
2007502992.tif 000066
【図４４】
2007502992.tif 000067
【図４４Ａ】
2007502992.tif 000068
【図４４Ｂ】
2007502992.tif 000069
【図４５】
2007502992.tif 000070
【図４５Ａ】
2007502992.tif 000071
【図４６】
2007502992.tif 000072
【図４７】
2007502992.tif 000073
【図４８】
2007502992.tif 000074
【図４９】
2007502992.tif 000075
【図５０】
2007502992.tif 000076
【図５０Ａ】
2007502992.tif 000077
【図５１】
2007502992.tif 000078
【図５２】
2007502992.tif 000079
【図５３】
2007502992.tif 000080
【図５４】
2007502992.tif 000081
【図５５】
2007502992.tif 000082
【図５６】
2007502992.tif 000083
【国際調査報告】
2007502992.tif 000084
2007502992.tif 000085
2007502992.tif 000086