特許 5725274 乳癌の予後の検査方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】5725274

(24)【登録日】2015年4月10日

(45)【発行日】2015年5月27日

(54)【発明の名称】乳癌の予後の検査方法

(51)【国際特許分類】

   C12Q 1/68 20060101AFI20150507BHJP

   G01N 33/53 20060101ALI20150507BHJP

   G01N 37/00 20060101ALI20150507BHJP

   C12N 15/09 20060101ALN20150507BHJP

【ＦＩ】

   C12Q1/68 A

   G01N33/53 M

   G01N37/00 102

   !C12N15/00 A

【請求項の数】7

【全頁数】31

(21)【出願番号】特願2010-98935(P2010-98935)

(22)【出願日】2010年4月22日

(65)【公開番号】特開2011-223957(P2011-223957A)

(43)【公開日】2011年11月10日

【審査請求日】2013年2月26日

【国等の委託研究の成果に係る記載事項】（出願人による申告）平成２１年度、文部科学省、「関西広域バイオメディカルクラスター構想（大阪北部（彩都）地域）に伴う研究委託業務」に係る再委託契約、産業技術力強化法第１９条の適用を受ける特許出願

(73)【特許権者】

【識別番号】504176911

【氏名又は名称】国立大学法人大阪大学

(73)【特許権者】

【識別番号】390014960

【氏名又は名称】シスメックス株式会社

(74)【代理人】

【識別番号】110000280

【氏名又は名称】特許業務法人サンクレスト国際特許事務所

(72)【発明者】

【氏名】野口 眞三郎

(72)【発明者】

【氏名】直居 靖人

(72)【発明者】

【氏名】岸 和希

(72)【発明者】

【氏名】後藤 健吾

【審査官】 森井 文緒

(56)【参考文献】

【文献】 Nature (2002) vol.415, no.6871, p.530-536

【文献】 実験医学 (2007) vol.25, no.17, p.58-62

【文献】 がん研究に関わる特定領域研究研究報告書録 平成２０年度 (2009) p.248

【文献】 Bioinformatics (2002) vol.18, no.12, p.1600-1608

【文献】 Expert Opin. Biol. Ther. (2005) vol.5, no.8, p.1069-1083

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｃ１２Ｎ １５／００−１５／９０

Ｃ１２Ｑ １／６８

ＰｕｂＭｅｄ

ＣＡｐｌｕｓ／ＭＥＤＬＩＮＥ／ＥＭＢＡＳＥ／ＢＩＯＳＩＳ（ＳＴＮ）

ＪＳＴＰｌｕｓ／ＪＭＥＤＰｌｕｓ／ＪＳＴ７５８０（ＪＤｒｅａｍＩＩＩ）

Ｔｈｏｍｓｏｎ Ｉｎｎｏｖａｔｉｏｎ

(57)【特許請求の範囲】

【請求項1】

（Ａ）被験者から採取された検体からＲＮＡを抽出する工程、
（Ｂ）前記工程（Ａ）で抽出されたＲＮＡを用いて測定用試料を調製する工程、
（Ｃ）前記工程（Ｂ）で得られた測定用試料を用いて、表１−１−１および表１−２−１に記載の遺伝子群の各遺伝子の発現量を測定する工程、
（Ｄ）前記工程（Ｃ）で測定された前記各遺伝子の発現量を解析する工程、および
（Ｅ）前記工程（Ｄ）で得られた解析結果に基づいて、乳癌の予後を検査する工程
を含む、乳癌の予後の検査方法。

【表1】

【表2】

【請求項2】

前記工程（Ｄ）において、前記発現量を、クラス分け手法を用いて解析する、請求項１に記載の乳癌の予後の検査方法。

【請求項3】

前記クラス分け手法が、Ｂｅｔｗｅｅｎ−ｇｒｏｕｐ ａｎａｌｙｓｉｓである、請求項２に記載の乳癌の予後の検査方法。

【請求項4】

前記工程（Ｄ）において、前記発現量と、下記式（１）：

【数1】

｛式（１）中、ｉは表１−１および表１−２に記載の核酸に付与された遺伝子番号を示し、ｗ_iは表１−１および表１−２に記載された遺伝子番号ｉの核酸に対応する重み係数を示し、Ｘ_iは下記式（２）：

【数2】

〔式（２）中、ｊは各検体に付与された検体番号を示し、ｙ_ijは遺伝子番号ｉの核酸に対応する遺伝子の検体番号ｊの検体での標準化された発現量を示し、ｍｉｎは括弧内の値の最小値を示し、ｒｏｕｎｄは括弧内の値の小数点以下第一位を四捨五入した値を示し、ａｂｓは括弧内の値の絶対値を示し、ｙ_iは、下記式（３）：

【数3】

（式（３）中、ｘ_iは遺伝子番号ｉの核酸に対応する遺伝子の発現量を示し、ｕ_iは遺伝子番号ｉの核酸に対応する遺伝子の発現量の検体に渡る平均値を示す。）
に示される遺伝子番号iの核酸に対応する遺伝子の標準化された発現量を示す。〕
により標準化され正規化された発現量を示し、Σ_iは各核酸に渡る総和を示す。｝
で表される判別式とを用いて前記判別式の解Ｄを算出し、
前記工程（Ｅ）において、前記判別式の解Ｄが正の値であるとき、予後が不良であると予測され、解Ｄがゼロまたは負の値のとき、予後が良好であると予測される、請求項１に記載の乳癌の予後の検査方法。

【表3】

【表4】

【請求項5】

前記工程（Ｄ）において、前記発現量を、階層的クラスター分析により解析する、請求項１に記載の乳癌の予後の検査方法。

【請求項6】

前記工程（Ｄ）において、前記発現量を、スコア化手法により解析する、請求項１に記載の乳癌の予後の検査方法。

【請求項7】

前記工程（Ｃ）において、前記発現量を、前記表１−１−１および表１−２−１に記載の核酸を少なくとも有するマイクロアレイを用いて測定する請求項１〜６のいずれかに記載の乳癌の予後の検査方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、乳癌の予後の検査方法に関する。

【背景技術】

【0002】

原発性乳癌の症例のうち約２／３の症例では、乳癌細胞にエストロゲン受容体（ＥＲ）が存在している（「ＥＲ陽性」という）。ＥＲ陽性の乳癌細胞では、ＥＲへのエストロゲンの結合が、細胞増殖の一因となっている。
したがって、リンパ節転移陰性・ＥＲ陽性の乳癌患者に対する治療においては、ＥＲを標的とするホルモン療法は、重要な役割を果たしている。

【0003】

前記リンパ節転移陰性・ＥＲ陽性の乳癌患者に対するホルモン療法では、例えば、当該患者に対して、タモキシフェンなどの抗エストロゲン剤を投与し、乳癌細胞のＥＲへのエストロゲンの結合をブロックして乳癌細胞の増殖を抑制することによって、乳癌の転移や再発の抑制が図られている。かかるホルモン療法を前記患者に施した場合、当該患者のうちのほとんどが、比較的良好な予後を示す。
しかしながら、前記患者の約２０％において、乳癌が再発することがある。

【0004】

そのため、前記リンパ節転移陰性・ＥＲ陽性乳癌患者は、ほとんどの場合、化学療法が不要であると考えられるにもかかわらず、乳癌の再発率を減らすために、前記リンパ節転移陰性・ＥＲ陽性乳癌患者のほとんどに対して、ホルモン療法だけでなく、アジュバント化学療法が施されているのが現状である。
したがって、高い再発リスクがある患者に対してのみアジュバント化学療法を施すために、前記リンパ節転移陰性・ＥＲ陽性乳癌患者の予後を予測することが重要であると考えられている。

【0005】

近年、網羅的遺伝子発現プロファイルの解析に基づき、乳癌患者における乳癌の予後を予測することが試みられている（例えば、特許文献１、２、非特許文献１などを参照）。

【0006】

前記特許文献１には、１１７人の乳癌患者の腫瘍サンプルを用いて同定された遺伝子マーカーを用い、乳癌患者の細胞サンプルにおける前記遺伝子マーカーの発現と対照における当該遺伝子マーカーの発現とを比較した差により、乳癌患者を「初期診断から５年以内に遠位転移なし」の患者または「初期診断から５年以内に遠位転移あり」の患者に分類する方法が記載されている。また、前記特許文献１には、前記遺伝子マーカーとして、ＥＲの有無を識別することができる遺伝子マーカー、ＢＲＣＡ１遺伝子の変異を有する腫瘍と散発性腫瘍とを識別することができる遺伝子マーカーおよび「初期診断から５年以内に遠位転移なし」の患者と「初期診断から５年以内に遠位転移あり」の患者とを識別することができる遺伝子マーカーを用いることが記載されている。

【0007】

また、前記特許文献２には、２８６人のリンパ節転移陰性の乳癌患者の腫瘍サンプルを用いて同定された予後診断の指標となる７６個の遺伝子を用い、乳癌患者の生体サンプルにおける前記遺伝子の発現プロファイルを得、当該発現レベルを所定のカットオフレベルと対比することで予後診断をする方法が記載されている。

【0008】

さらに、前記非特許文献１には、１８９例の侵襲性乳癌患者と３つの公知の乳癌の遺伝子発現データセットとを用いて同定された９７個の遺伝子に基づくＧｅｎｏｍｉｃ Ｇｒａｄｅ Ｉｎｄｅｘ（ＧＧＩ）を用い、組織学的悪性度２として従来分類されていた乳癌の症例を、さらに、再発高リスク群と再発低リスク群とに分け、予後を予測する方法が記載されている。

【先行技術文献】

【特許文献】

【0009】

【特許文献1】特表２００９−１３１２６２号公報

【特許文献2】特開２００７−５２８２１８号公報

【非特許文献】

【0010】

【非特許文献1】クリストス・ソテリウ（Ｓｏｔｉｒｉｏｕ Ｃｈｒｉｓｔｏｓ）ら、「乳癌における遺伝子発現プロファイル：予後を改善するための組織学的悪性度の理解（Ｇｅｎｅ ｅｘｐｒｅｓｓｉｏｎ ｐｒｏｆｉｌｉｎｇ ｉｎ ｂｒｅａｓｔ ｃａｎｃｅｒ： ｕｎｄｅｒｓｔａｎｄｉｎｇ ｔｈｅ ｍｏｌｅｃｕｌａｒ ｂａｓｉｓ ｏｆ ｈｉｓｔｏｌｏｇｉｃ ｇｒａｄｅ ｔｏ ｉｍｐｒｏｖｅ ｐｒｏｇｎｏｓｉｓ．）」、ジャーナル・オブ・ザ・ナショナル・キャンサー・インスティチュート（Ｊｏｕｒｎａｌ ｏｆ ｔｈｅ Ｎａｔｉｏｎａｌ Ｃａｎｃｅｒ Ｉｎｓｔｉｔｕｔｅ）、２００６年２月１５日発行、第９８巻、第４号、ｐ．２６２−２７２

【発明の概要】

【発明が解決しようとする課題】

【0011】

しかしながら、これらの方法は、調べられた施設間の差や人種間の差などの影響を受け、必ずしも適切に予後を予測することができるとは限らないのが現状である。

【0012】

本発明は、前記従来技術に鑑みてなされたものであり、予後を適切に予測することができる乳癌における予後の検査方法を提供することを目的とする。

【課題を解決するための手段】

【0013】

すなわち、本発明は、
〔１〕（Ａ）乳癌に罹患した被験者から採取された検体からＲＮＡを抽出する工程、
（Ｂ）前記工程（Ａ）で抽出されたＲＮＡを用いて測定用試料を調製する工程、
（Ｃ）前記工程（Ｂ）で得られた測定用試料を用いて、表１−１および表１−２に記載の遺伝子群の各遺伝子の発現量を測定する工程、
（Ｄ）前記工程（Ｃ）で測定された前記発現量を解析する工程、および
（Ｅ）前記工程（Ｄ）で得られた解析結果に基づいて、乳癌の予後を検査する工程
を含む、乳癌の予後の検査方法、
〔２〕前記工程（Ｄ）において、前記発現量を、クラス分け手法を用いて解析する、前記〔１〕に記載の乳癌の予後の検査方法、
〔３〕前記クラス分け手法が、Ｂｅｔｗｅｅｎ−ｇｒｏｕｐ ａｎａｌｙｓｉｓである、前記〔２〕に記載の乳癌の予後の検査方法、
〔４〕前記工程（Ｄ）において、前記発現量と、下記式（１）：

【0014】

【数1】

【0015】

｛式（１）中、ｉは表１−１および表１−２に記載の核酸に付与された遺伝子番号を示し、ｗ_iは表１−１および表１−２に記載された遺伝子番号ｉの核酸に対応する重み係数を示し、Ｘ_iは下記式（２）：

【0016】

【数2】

【0017】

〔式（２）中、ｊは各検体に付与された検体番号を示し、ｙ_ijは遺伝子番号ｉの核酸に対応する遺伝子の検体番号ｊの検体での標準化された発現量を示し、ｍｉｎは括弧内の値の最小値を示し、ｒｏｕｎｄは括弧内の値の小数点以下第一位を四捨五入した値を示し、ａｂｓは括弧内の値の絶対値を示し、ｙ_iは、下記式（３）：

【0018】

【数3】

【0019】

（式（３）中、ｘ_iは遺伝子番号ｉの核酸に対応する遺伝子の発現量を示し、ｕ_iは遺伝子番号ｉの核酸に対応する遺伝子の発現量の検体に渡る平均値を示す。）
に示される遺伝子番号iの核酸に対応する遺伝子の標準化された発現量を示す。〕
により標準化され正規化された発現量を示し、Σ_iは各核酸に渡る総和を示す。｝
で表される判別式とを用いて前記判別式の解Ｄを算出し、
前記工程（Ｅ）において、前記判別式の解Ｄが正の値であるとき、予後が不良であると予測され、解Ｄがゼロまたは負の値のとき、予後が良好であると予測される、前記〔１〕に記載の乳癌の予後の検査方法、
〔５〕前記工程（Ｄ）において、前記発現量を、階層的クラスター分析により解析する、前記〔１〕に記載の乳癌の予後の検査方法、
〔６〕前記工程（Ｄ）において、前記発現量を、スコア化手法により解析する、前記〔１〕に記載の乳癌の予後の検査方法、ならびに
〔７〕前記工程（Ｃ）において、前記発現量を、表１−１および表１−２に記載の核酸を少なくとも有するマイクロアレイを用いて測定する前記〔１〕〜〔６〕のいずれかに記載の乳癌の予後の検査方法
に関する。

【発明の効果】

【0020】

本発明の乳癌の予後の検査方法によれば、予後を適切に予測することができるという優れた効果を奏する。

【図面の簡単な説明】

【0021】

【図1】実施例１において、プローブ数と、ＲＯＣ曲線の曲線下面積との関係を調べた結果を示すグラフである。

【図2】実施例２において、１０５例の乳癌患者について、判別式による予測結果と、観察結果とを比較した結果を示す図である。

【図3】実施例２において、術後からの期間と、非再発生存率との関係を調べた結果を示すグラフである。

【図4】実施例３において、１０５例の乳癌患者それぞれにおける核酸（プローブセット）に対応する遺伝子の発現量のデータの階層的クラスター分析を行なった結果を示す樹形図である。

【図5】実施例４において、１０５例の乳癌患者それぞれの発現量のデータに基づき算出した第一主成分スコアおよび第二主成分スコアの散布図である。

【発明を実施するための形態】

【0022】

本発明の乳癌の予後の検査方法は、
（Ａ）被験者から採取された検体からＲＮＡを抽出する工程、
（Ｂ）前記工程（Ａ）で抽出されたＲＮＡを用いて測定用試料を調製する工程、
（Ｃ）前記工程（Ｂ）で得られた測定用試料を用いて、表１−１および表１−２に記載の遺伝子群の各遺伝子の発現量を測定する工程、
（Ｄ）前記工程（Ｃ）で測定された前記発現量を解析する工程、および
（Ｅ）前記工程（Ｄ）で得られた解析結果に基づいて、乳癌の予後を検査する工程
を含むことを特徴としている。

【0023】

【表1】

【0024】

【表2】

【0025】

本明細書において、「プローブセット．ＩＤ」は、アフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社製のマイクロアレイ〔商品名：ＧｅｎｅＣｈｉｐ（登録商標） Ｓｙｓｔｅｍ〕において、基材上に固定されたプローブの１１〜２０個をまとめたプローブセットそれぞれにつけられているＩＤ番号を示す。前記プローブセット．ＩＤで示された核酸（プローブセット）の塩基配列は、ウェブページｈｔｔｐｓ：／／ｗｗｗ．ａｆｆｙｍｅｔｒｉｘ．ｃｏｍ／ａｎａｌｙｓｉｓ／ｎｅｔａｆｆｘ／ｉｎｄｅｘ．ａｆｆｘにより容易に入手することができる（２００９年６月３０日更新のデータベース）。「ＵｎｉＧｅｎｅ．ＩＤ」は、ＮＣＢＩが公開しているデータベースであるＵｎｉＧｅｎｅのＩＤ番号を示す。ＧｅｎＢａｎｋアクセッション番号は、前記アフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社製のマイクロアレイ〔商品名：ＧｅｎｅＣｈｉｐ（登録商標） Ｓｙｓｔｅｍ〕において、基材上に固定されたプローブそれぞれの配列の設計に用いられた公開データベースＧｅｎＢａｎｋのアクセッション番号を示す。

【0026】

なお、本明細書において、「表１−１および表１−２に記載の遺伝子群の各遺伝子の発現量」の用語は、表１−１および表１−２に記載のプローブセット．ＩＤに対応する表１−１および表１−２に記載のＧｅｎＢａｎｋアクセッション番号で示される核酸を保持する遺伝子の発現量を意味する。ＧｅｎＢａｎｋは、米国生物工学情報センター（Ｎａｔｉｏｎａｌ Ｃｅｎｔｅｒ ｆｏｒ Ｂｉｏｔｅｃｈｎｏｌｏｇｙ Ｉｎｆｏｒｍａｔｉｏｎ）により提供されているデータベースであり、ウェブページｈｔｔｐ：／／ｗｗｗ．ｎｃｂｉ．ｎｌｍ．ｎｉｈ．ｇｏｖ／ｓｉｔｅｓ／ｅｎｔｒｅｚ？ｄｂ＝ｎｕｃｌｅｏｔｉｄｅなどにより利用可能である。そして、前記表１−１および１−２に記載のＧｅｎＢａｎｋアクセッション番号が付された配列は、上記データベースより入手することができる。また、前記ＧｅｎＢａｎｋアクセッション番号は、２００９年６月３０日時点での最新リリースでの番号を示す。なお、本明細書において、「遺伝子」とは、遺伝子転写産物としてのＲＮＡが抽出される塩基配列の単位またはその一部であればよく、ＥＳＴ（ｅｘｐｒｅｓｓｅｄ ｓｅｑｕｅｎｃｅ ｔａｇ）も含む概念である。

【0027】

本発明の検査方法では、まず、被験者から採取された検体からＲＮＡを抽出する〔工程（Ａ）〕。

【0028】

前記被験者とは、乳癌患者および乳癌の疑いのある患者をいう。ここで、乳癌患者としては、特に限定されないが、例えば、リンパ節転移陰性・ＥＲ陽性の乳癌患者などが挙げられる。かかるリンパ節転移陰性・ＥＲ陽性の乳癌患者は、抗エストロゲン剤を投与するホルモン療法が施されたリンパ節転移陰性・ＥＲ陽性の乳癌患者であってもよい。

【0029】

本発明の検査方法では、被験者が、例えば、抗エストロゲン剤を投与するホルモン療法が施されたリンパ節転移陰性・ＥＲ陽性の乳癌患者である場合、患者が予後良好であるという予測を高い精度で行なうことができる。なお、本明細書において、「予後良好」とは、術後１０年間に再発がないことをいう。

【0030】

前記検体としては、例えば、手術時に切除された腫瘍組織、生検により被験者から採取された検体などが挙げられる。

【0031】

検体からのＲＮＡの抽出は、公知の方法によって行なうことができる。また、検体からのＲＮＡの抽出は、ＲＮＡを抽出するための市販のキットを用いることもできる。ここで、市販のキットとしては、キアゲン（Ｑｉａｇｅｎ）社製、商品名：Ｑｉａｇｅｎ ＲＮｅａｓｙ ｋｉｔ（登録商標）などが挙げられる。

【0032】

つぎに、前記工程（Ａ）で抽出されたＲＮＡを用いて測定用試料を調製する〔工程（Ｂ）〕。

【0033】

本工程（Ｂ）では、遺伝子の発現量、すなわち、当該遺伝子に対応する転写産物（ｃＲＮＡ、ｃＤＮＡ、ｍＲＮＡなど）の産生量などを測定するのに適した測定用試料を調製する。具体的には、測定用試料は、例えば、前記工程（Ａ）で抽出されたＲＮＡを用いて、対応するｃＲＮＡまたはｃＤＮＡを増幅すること、前記工程（Ａ）で抽出されたＲＮＡよりｍＲＮＡを精製することなどによって調製することができる。また、本発明においては、遺伝子の発現量を測定することが可能であるのであれば、前記工程（Ａ）で抽出されたＲＮＡをそのまま測定用試料として用いてもよい。

【0034】

前記ｃＲＮＡの増幅は、公知の方法を用いて行なうことができる。ｃＲＮＡの増幅には、ｃＲＮＡの増幅するための市販のキットを用いることもできる。ここで、前記市販のキットとしては、例えば、アフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社製、商品名：Ｏｎｅ−Ｃｙｃｌｅ Ｔａｒｇｅｔ Ｌａｂｅｌｉｎｇ ａｎｄ Ｃｏｎｔｒｏｌ Ｒｅａｇｅｎｔｓなどが挙げられる。また、前記ｃＤＮＡの増幅は、公知の方法を用いて行なうことができる。ｃＤＮＡの増幅には、ｃＤＮＡの増幅するための市販のキットを用いることもできる。前記ｍＲＮＡの精製は、公知の精製方法を用いて行なうことができる。また、ｍＲＮＡの精製には、市販の精製キットを用いてもよい。

【0035】

つぎに、前記工程（Ｂ）で得られた測定用試料を用いて、前記表１−１および表１−２に記載の遺伝子群の各遺伝子の発現量を測定する〔工程（Ｃ）〕。

【0036】

本工程（Ｃ）において、前記発現量の測定には、例えば、マイクロアレイ、定量的ＲＴ-ＰＣＲ、定量的ＰＣＲなどを用いることができる。これらのなかでは、前記発現量を迅速、かつ簡便に測定することができることから、マイクロアレイを用いて測定することが好ましい。この場合、前記発現量は、以下の工程において、マイクロアレイにおける蛍光強度のまま用いてもよい。マイクロアレイによる前記発現量の測定は、公知の方法を用いて行なうことができる。

【0037】

表１−１および表１−２に記載の遺伝子群の各遺伝子の発現量は、例えば、かかる表１−１および表１−２に記載の核酸（プローブセット．ＩＤで示されたプローブセット）を利用することで測定することができる。本発明の検査方法においては、前記表１−１および表１−２に記載の核酸（プローブセット）は、乳癌の予後の検査に際して、予後因子として用いられる。かかる核酸（プローブセット）は、乳癌の多くの症例における予後の検査において、効果量が大きいものとして本発明者らにより見出されたものである。また、本発明の検査方法における予後因子に用いる前記核酸（プローブセット）の数は、９５個であり、検査の精度が最も高くなる数とされている。したがって、本発明の検査方法によれば、種々の症例について、予後を適切に予測することができる。

【0038】

本発明の検査方法では、マイクロアレイを用いて前記表１−１および表１−２に記載の遺伝子群の各遺伝子の発現量を測定する場合、前記マイクロアレイとして、例えば、前記表１−１および表１−２に記載の核酸（プローブセット）を少なくとも有するマイクロアレイなどを用いることができる。かかるマイクロアレイとしては、例えば、ヒトゲノムの発現解析が可能なマイクロアレイであるアフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社製の商品名：Ｈｕｍａｎ Ｇｅｎｏｍｅ Ｕ１３３ Ｐｌｕｓ ２．０ Ａｒｒａｙなどが挙げられる。例えば、前記アフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社製の商品名：Ｈｕｍａｎ Ｇｅｎｏｍｅ Ｕ１３３ Ｐｌｕｓ ２．０ Ａｒｒａｙを用いた場合、表１−１および表１−２に記載のプローブセット．ＩＤで示される９５個の核酸（プローブセット）によって、前記表１−１および表１−２に記載の遺伝子群の各遺伝子の発現量を一度に測定することができる。

【0039】

つぎに、前記工程（Ｃ）で測定された前記発現量を解析する〔工程（Ｄ）〕。その後、前記工程（Ｄ）で得られた解析結果に基づいて、乳癌の予後を検査する〔工程（Ｅ）〕。

【0040】

前記工程（Ｄ）において、前記発現量は、例えば、クラス分け手法、階層的クラスター分析およびスコア化手法を用いて、解析することができる。ここで、前記発現量は、測定された発現量の生データを、例えば、ＲＭＡアルゴリズム、ＭＡＳ５アルゴリズム、ＰＬＩＥＲアルゴリズムなどによって正規化したものを用いることができる。前記ＲＭＡアルゴリズムは、例えば、解析ソフトウェア〔アフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社製、商品名：Ａｆｆｙｍｅｔｒｉｘ Ｅｘｐｒｅｓｓｉｏｎ Ｃｏｎｓｏｌｅソフトウェア〕上で利用可能である。

【0041】

前記クラス分け手法として、公知の方法を用いることができる。かかるクラス分け手法としては、例えば、Ｂｅｔｗｅｅｎ−ｇｒｏｕｐ ａｎａｌｙｓｉｓ（ＢＧＡ）〔カルヘイン（Ｃｕｌｈａｎｅ，Ａ．Ｃ．）ら、バイオインフォマティックス（Ｂｉｏｉｎｆｏｒｍａｔｉｃｓ）、２００２年、第１８巻、ｐｐ．１６００−１６０８参照），「マイクロアレイデータのＢｅｔｗｅｅｎ−ｇｒｏｕｐ ａｎａｌｙｓｉｓ（Ｂｅｔｗｅｅｎ−ｇｒｏｕｐ ａｎａｌｙｓｉｓ ｏｆ ｍｉｃｒｏａｒｒａｙ ｄａｔａ）」を参照〕、サポートベクターマシン（ＳＶＭ）、対角線形判別（ＤＬＤＡ）、ｋ最近傍分類（ｋＮＮ）、決定木、Ｒａｎｄｏｍ Ｆｏｒｅｓｔ、ニューラルネットなどが挙げられる。これらのなかでは、予後良好であることが予測されるものと、予後不良であることが予測されるものとを良好にクラス分けをすることができる観点から、ＢＧＡが好ましい。かかるクラス分け手法を用いて発現量を解析する場合、前記発現量に基づいて、予後良好であることが予測されるものと、予後不良であることが予測されるものとにクラス分けがされる。したがって、この場合、前記工程（Ｅ）においては、かかるクラス分けの結果によって、乳癌の予後を予測することができる。前記ＢＧＡを用いる場合には、判別式を構築して、この判別式の解Ｄに基づいて、乳癌の予後を予測してもよい。

【0042】

前記判別式としては、下記式（１）

【0043】

【数4】

【0044】

｛式（１）中、ｉは前記表１−１および表１−２に記載の核酸に付与された遺伝子番号を示し、ｗ_iは前記表１−１および表１−２に記載された遺伝子番号ｉの核酸に対応する重み係数を示し、Ｘ_iは下記式（２）：

【0045】

【数5】

【0046】

〔式（２）中、ｊは各検体に付与された検体番号を示し、ｙ_ijは遺伝子番号ｉの核酸に対応する遺伝子の検体番号ｊの検体での標準化された発現量を示し、ｍｉｎは括弧内の値の最小値を示し、ｒｏｕｎｄは括弧内の値の小数点以下第一位を四捨五入した値を示し、ａｂｓは括弧内の値の絶対値を示し、ｙ_iは、下記式（３）：

【0047】

【数6】

【0048】

（式（３）中、ｘ_iは遺伝子番号ｉの核酸に対応する遺伝子の発現量を示し、ｕ_iは遺伝子番号ｉの核酸に対応する遺伝子の発現量の検体に渡る平均値を示す。）
に示される遺伝子番号iの核酸に対応する遺伝子の標準化された発現量を示す。〕
により標準化され正規化された発現量を示し、Σ_iは各核酸に渡る総和を示す。｝
で表される判別式が挙げられる。かかる判別式を用いて発現量を解析する場合、前記判別式のｘ_i（ｉ＝１，２，・・・，９５）に、順に、当該検体における発現量の値を代入し、解Ｄを算出する。この場合、前記工程（Ｅ）においては、解Ｄが正の値のとき、予後が不良であると予測することができ、解Ｄがゼロまたは負の値のとき、予後が良好であると予測することができる。

【0049】

前記階層的クラスター分析は、例えば、被験者から採取した検体における前記発現量のデータ（あるいは前記発現量に関する蛍光強度のデータ）と、予後良好であることが既知の検体の群における前記発現量のデータ（あるいは前記発現量に関する蛍光強度のデータ）と、予後不良であることが既知の検体の群における前記発現量のデータ（あるいは前記発現量に関する蛍光強度のデータ）とを用い、前記発現量（あるいは前記発現量に関する蛍光強度）に基づいて検体間の類似度を示す距離を計算し、この距離に基づいて種々のクラスターを形成し、クラスターを統合して、樹形図を作成することにより行なうことができる。ここで、前記距離としては、例えば、Ｓｐｅａｒｍａｎ順位相関係数、ユークリッド距離などが挙げられる。また、クラスターの統合は、例えば、ｗａｒｄ法、最遠隣法、重心間距離法などにより行なうことができる。これらのなかでは、Ｓｐｅａｒｍａｎ順位相関係数およびｗａｒｄ法を用いることにより、予後良好であることが予測されるものと、予後不良であることが予測されるものとを良好に分けることができる。この場合、前記工程（Ｅ）においては、かかる階層的クラスター分析の結果によって、乳癌の予後を適切に予測することができる。

【0050】

前記スコア化手法として、公知の方法を用いることができる。かかるスコア化手法としては、例えば、主成分分析、重回帰分析、ロジスティック回帰分析、Ｐａｒｔｉａｌ Ｌｅａｓｔ Ｓｑｕａｒｅなどが挙げられる。これらのなかでは、予後良好であることが予測されるものと、予後不良であることが予測されるものとを良好に分けることができる観点から、主成分分析が好ましい。かかるスコア化手法を用いて前記発現量を解析する場合、前記発現量に基づいて、予後良好であることが予測される検体のスコアと、予後不良であることが予測される検体のスコアとが分けられるようスコア化がされる。したがって、この場合、前記工程（Ｅ）においては、かかるスコア化の結果によって、乳癌の予後を適切に予測することができる。

【0051】

以上のように、本発明の乳癌の予後の検査方法によれば、前記表１−１および表１−２に記載の遺伝子群の各遺伝子の発現量を解析するので、乳癌の予後を適切に検査することができる。

【0052】

なお、前記表１−１および表１−２に記載の遺伝子群の各遺伝子の発現量を用いることにより、乳癌の予後を判定するための指標を取得し、提供することができる。乳癌の予後を判定するための指標の取得方法は、
（ａ）乳癌に罹患した被験者から採取された検体からＲＮＡを抽出する工程、
（ｂ）前記工程（ａ）で抽出されたＲＮＡを用いて測定用試料を調製する工程、
（ｃ）前記工程（ｂ）で得られた測定用試料を用いて、表１−１および表１−２に記載の遺伝子群の各遺伝子の発現量を測定する工程、
（ｄ）前記工程（ｃ）で測定された前記発現量を解析する工程、および
（ｅ）前記工程（ｄ）で得られた解析結果に基づいて、被験者の乳癌の予後不良の可能性の指標または予後良好の可能性の指標を得る工程、
を含む。前記工程（ａ）〜工程（ｄ）は、前記乳癌の予後の検査方法における工程（Ａ）〜（Ｄ）と同様に行なうことができる。また、前記工程（ｅ）において、被験者の乳癌の予後不良の可能性の指標または予後良好の可能性の指標は、前記乳癌の予後の検査方法における工程（Ｅ）で用いられる手法により得ることができる。

【実施例】

【0053】

以下、実施例を挙げて本発明を詳細に説明するが、本発明はこれらに限定されるものではない。

【0054】

［実施例１］
マイクロアレイ実験の遺伝子発現情報データベース〔ＮＣＢＩ Ｇｅｎｅ Ｅｘｐｒｅｓｓｉｏｎ Ｏｍｎｉｂｕｓ （ｈｔｔｐ：／／ｗｗｗ．ｎｃｂｉ．ｎｌｍ．ｎｉｈ．ｇｏｖ／ｇｅｏ／）〕のアクセッション番号：ＧＳＥ２０３４、ＧＳＥ２９９０、ＧＳＥ４９２２、ＧＳＥ６５３２、ＧＳＥ７３９０およびＧＳＥ９１９５の６つのデータセットからリンパ節転移陰性・ＥＲ陽性の症例５４９症例のデータを抽出した。

【0055】

５４９症例の各データにおいて、前記データの取得に用いられたマイクロアレイ上の各核酸（プローブセット）の発現量について、解析ソフトウェア〔アフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社製、商品名：Ａｆｆｙｍｅｔｒｉｘ Ｅｘｐｒｅｓｓｉｏｎ Ｃｏｎｓｏｌｅソフトウェア〕のＲＭＡアルゴリズムを用いて前記データセット毎に正規化した。つぎに、前記データセット毎に、データの取得に用いられたアレイ上の核酸（プローブセット）の発現量の値から当該データセット内での前記核酸の平均発現量の値を減じることで、前記核酸の発現量を標準化した。

【0056】

その後、統計解析用ソフトウェア「Ｒ」で用いられる追加パッケージ集「ＢｉｏＣｏｎｄｕｃｔｏｒ」ｖｅｒ．２．４に含まれるパッケージ「ＧｅｎｅＭｅｔａ ｖ１．１６．０」（ｈｔｔｐ：／／ｗｗｗ．ｂｉｏｃｏｎｄｕｃｔｏｒ．ｏｒｇ／ｐａｃｋａｇｅｓ／２．４／ｂｉｏｃ／ｈｔｍｌ／ＧｅｎｅＭｅｔａ．ｈｔｍｌ）を用い、ジュン・キョン・チェ（Ｊｕｎｇ Ｋｙｏｏｎ Ｃｈｏｉ）らの文献〔「複数のマイクロアレイ研究の統合および研究間バリデーションのモデリング（Ｃｏｍｂｉｎｉｎｇ ｍｕｌｔｉｐｌｅ ｍｉｃｒｏａｒｒａｙ ｓｔｕｄｉｅｓ ａｎｄ ｍｏｄｅｌｉｎｇ ｉｎｔｅｒｓｔｕｄｙ ｖａｒｉａｔｉｏｎ）」バイオインフォマティックス（Ｂｉｏｉｎｆｏｒｍａｔｉｃｓ）、第１９巻、補遺１、２００３年、ｐ．ｉ８４−ｉ９０〕にしたがって、前記アレイ上の核酸（プローブセット）毎にｚ値を算出した。そして、前記核酸（プローブセット）をｚ値の絶対値の順に並べた。

【0057】

つぎに、Ｂｅｔｗｅｅｎ-ｇｒｏｕｐ ａｎａｌｙｓｉｓにより、判別式を構築した。最適な精度が得られるように、Ｓｅｑｕｅｎｔｉａｌ Ｆｏｒｗａｒｄ Ｆｉｌｔｅｒｉｎｇ ｍｅｔｈｏｄによって乳癌の予後の検査に最適なプローブセット数を求めた。ここでは、前記核酸（プローブセット）から、ｚ値の絶対値が大きい順に、３００個に達するまで、５個ずつ選択する核酸（プローブセット）の個数を増加させながら、当該核酸（プローブセット）を選択し、判別式の構築を行なった。

【0058】

得られた各判別式と、前記５４９症例の各データとを用い、Ｌｅａｖｅ-Ｏｎｅ-Ｏｕｔ Ｃｒｏｓｓ-Ｖａｌｉｄａｔｉｏｎにより、各判別式による検査精度を検証した。前記検査精度は、各判別式の感度および特異度を求め、前記感度および特異度からＲＯＣ曲線を作成し、このＲＯＣ曲線の曲線下面積を算出することにより評価した。そして、ＲＯＣ曲線の曲線下面積と、核酸（プローブセット）の個数との関係を調べて、ＲＯＣ曲線の曲線下面積が最大となる核酸（プローブセット）の個数（最適な個数）を調べた。

【0059】

なお、感度は、観察結果により術後１０年間に「再発あり」（予後不良）と判断され、かつ判別式により「再発あり」（予後不良）であると予測された検体の数を、観察結果により術後１０年間に「再発あり」（予後不良）と判断された検体の数で除算し、１００を乗ずることによって求めた。また、特異度は、観察結果により術後１０年間に「再発なし」（予後良好）と判断され、かつ判別式により「再発なし」（予後良好）と予測された検体の数を、観察結果により「再発なし」（予後良好）と判断された検体の数で除算し、１００を乗ずることによって求めた。

【0060】

実施例１において、プローブ数とＲＯＣ曲線の曲線下面積との関係を調べた結果を示すグラフを図１に示す。

【0061】

図１に示された結果から、核酸（プローブセット）の個数が９５個である場合、ＲＯＣ曲線の曲線下面積が最大となることがわかる。この結果から、核酸（プローブセット）の個数が９５個である場合、検査精度が最も高くなることがわかる。なお、これら９５個の核酸（プローブセット）は、前記表１−１および表１−２に示されるとおりである。

【0062】

また、これら９５個の核酸（プローブセット）のｚ値、術後１０年間に再発した再発群での発現傾向および判別式での重み係数を表２−１および表２−２に示す。

【0063】

【表3】

【0064】

【表4】

【0065】

前記結果に基づいて、最終的な判別式を構築した。構築された判別式は、下記式（１）：

【0066】

【数7】

【0067】

｛式（１）中、ｉは前記表１−１および表１−２に記載の遺伝子群の各遺伝子に付与された遺伝子番号を示し、ｗ_iは前記表１−１および表１−２に記載された遺伝子番号ｉの遺伝子に対応する重み係数を示し、Ｘ_iは下記式（２）：

【0068】

【数8】

【0069】

〔式（２）中、ｊは各検体に付与された検体番号を示し、ｙ_ijは遺伝子番号ｉの遺伝子の検体番号ｊの検体での標準化された発現量を示し、ｍｉｎは括弧内の値の最小値を示し、ｒｏｕｎｄは括弧内の値の小数点以下第一位を四捨五入した値を示し、ａｂｓは括弧内の値の絶対値を示し、ｙ_iは、下記式（３）：

【0070】

【数9】

【0071】

（式（３）中、ｘ_iは遺伝子番号ｉの核酸に対応する遺伝子の発現量を示し、ｕ_iは遺伝子番号ｉの核酸に対応する遺伝子の発現量の検体に渡る平均値を示す。）
に示される遺伝子番号iの核酸に対応する遺伝子の標準化された発現量を示す。〕
により標準化され正規化された発現量を示し、Σ_iは各核酸に渡る総和を示す。｝
で表される判別式である。

【0072】

ここで、前記判別式の解Ｄが正の値のとき、予後が不良であると予測され、解Ｄがゼロまたは負の値のとき、予後が良好であると予測される。

【0073】

［実施例２］
（１）核酸（プローブセット）の発現量のデータの取得
１０５人の乳癌患者それぞれの手術時に得られた腫瘍組織から、ＲＮＡ抽出キット〔キアゲン・サイエンシーズ（ＱＩＡＧＥＮ Ｓｃｉｅｎｃｅｓ）製、商品名：Ｑｉａｇｅｎ ＲＮｅａｓｙ ｍｉｎｉ ｋｉｔ〕を用いてＲＮＡを抽出した。

【0074】

なお、１０５人の乳癌患者は、リンパ節転移陰性・ＥＲ陽性であり、１９９６年〜２００５年の間に、乳房温存手術を受け、その後、放射線療法または乳腺切除術を受けた患者である。これらの患者の年齢の範囲は、３０〜８３歳であり、年齢の中央値は、５４歳である。前記１０５人の乳癌患者の臨床病理学的特徴を表３に示す。

【0075】

【表5】

【0076】

なお、前記腫瘍径Ｔは、マンモグラフィや超音波などの画像診断による測定結果に基づき、１〜４の４段階で表される。ここで、Ｔ＝１は腫瘍の最大径が２ｃｍ以下であること、Ｔ＝２は腫瘍の最大径２ｃｍを超えて５ｃｍ以下であること、Ｔ＝３は腫瘍最大径が５ｃｍを超えていること、Ｔ＝４は腫瘍の大きさに関わらず、胸壁または皮膚に浸潤していることを示す。
前記組織学的グレードは、核の異型度（１点：核多型性が低頻度、２点：核多型性が中頻度、３点：核多型性が高頻度）のスコア、組織構造の変化度（１点：＜１０％、２点：１０〜７５％、３点：＞７５％）のスコア、細胞分裂の頻度（１点：１０高倍率視野（ＨＰＦ）で細胞の有糸分裂数が０〜４、２点：１０ＨＰＦで細胞の有糸分裂数が５〜１０、３点：１０高倍率視野（ＨＰＦ）で細胞の有糸分裂数が１１≦）のスコアの合計によって、１〜３の３段階で表される。ここで、従来、ＨＧ＝１はスコア３〜５であり予後の良好な癌、ＨＧ＝２はスコア６〜７であり、ＨＧ＝３はスコア８〜９であり、最も予後の悪い癌を示すものとされていた。
前記ＥＲの有無は、免疫染色法による結果に基づき、陽性および陰性で表される。従来、ＥＲ陰性である場合、一般に、予後が不良であり、ＥＲ陽性である場合、予後が良好であるとされていた。
前記ＰＲの有無は、免疫染色法による結果に基づき、陽性および陰性で表される。従来、ＰＲ陰性である場合、一般に、予後が不良であり、ＰＲ陽性である場合、予後が良好であるとされていた。
前記ＨＥＲ２の有無は、免疫染色法による結果に基づき、陽性および陰性で表される。従来、ＨＥＲ２陽性である場合、一般に、予後が不良であり、ＨＥＲ２陰性である場合、予後が良好であるとされていた。
前記Ｋｉ６７は、免疫染色法による結果に基づき、陽性および陰性で表わされる。従来、Ｋｉ６７陽性である場合、一般に、予後が不良であり、Ｋｉ６７陰性である場合、予後が良好であるとされていた。

【0077】

つぎに、得られたＲＮＡ〔ＲＩＮ（ＲＮＡ Ｉｎｔｅｇｒｉｔｙ Ｎｕｍｂｅｒ））値＞６）１μｇと、発現解析用キット〔アフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社製、商品名：Ｏｎｅ−Ｃｙｃｌｅ Ｔａｒｇｅｔ Ｌａｂｅｌｉｎｇ ａｎｄ Ｃｏｎｔｒｏｌ Ｒｅａｇｅｎｔｓ〕とを用いて、ｃＲＮＡを増幅し、ビオチン化し、断片化した。

【0078】

得られた断片化ビオチン標識ｃＲＮＡを、ヒトゲノム発現解析用アレイ〔アフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社製、商品名：Ｈｕｍａｎ Ｇｅｎｏｍｅ Ｕ１３３ Ｐｌｕｓ ２．０ Ａｒｒａｙ〕上の核酸（プローブセット）と一晩ハイブリダイズさせた。なお、前記断片化ビオチン標識ｃＲＮＡと前記アレイ上の核酸（プローブセット）とのハイブリダイゼーションは、製造者〔アフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社〕による推奨条件に従って行なった。

【0079】

つぎに、ハイブリダイゼーション後のアレイを、マイクロアレイ洗浄・染色処理専用機器〔アフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社製、商品名：ＧｅｎｅＣｈｉｐ Ｆｌｕｉｄｉｃｓ Ｓｔａｔｉｏｎ ４５０〕に供して、前記アレイ上の核酸（プローブセット）にハイブリダイズしたｃＲＮＡを蛍光染色し、洗浄した。

【0080】

その後、前記アレイをレーザスキャナー〔アフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社製、商品名：ＧｅｎｅＣｈｉｐ（登録商標） Ｓｃａｎｎｅｒ ３０００〕に供して、前記アレイ上の核酸（プローブセット）にハイブリダイズしたｃＲＮＡの蛍光標識物質に基づくシグナルを読み取り、蛍光強度を定量化した。

【0081】

得られた蛍光強度のデータをソフトウェア〔アフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社製、商品名：ＧｅｎｅＣｈｉｐ（登録商標） Ｏｐｅｒａｔｉｎｇ Ｓｏｆｔｗａｒｅ〕によって処理して、ＣＥＬファイルを得た。

【0082】

得られた１０５例のＣＥＬファイルのデータと、解析ソフトウェア〔アフィメトリクス（Ａｆｆｙｍｅｔｒｉｘ）社製、商品名：Ａｆｆｙｍｅｔｒｉｘ Ｅｘｐｒｅｓｓｉｏｎ Ｃｏｎｓｏｌｅソフトウェア〕のＲＭＡアルゴリズムとを用いて、全１０５例の乳癌患者における核酸（プローブセット）の発現量のデータ（蛍光強度のデータ）を正規化した。

【0083】

（２）判別式の性能の検証
つぎに、前記（１）で得られた正規化後のデータと、前記判別式とを用いて、全１０５例の乳癌患者が再発を起こすかどうかを予測した。そして、病理学的な観察結果を真値として、当該病理学的な観察結果と前記判別式による予測結果とを比較することによって、当該判別式の性能を評価した。実施例２において、１０５例の乳癌患者について、判別式による予測結果と、観察結果との関係を調べた結果を図２に示す。

【0084】

図２に示された結果から、全１０５例の検体のうち、前記判別式によって、６１例の乳癌患者が「再発なし」と予測され、４４例の乳癌患者が「再発あり」と予測されることがわかる。すなわち、前記判別式を用いた場合には、前記６１例の乳癌患者については、予後が良好であることが予測され、一方、前記４４例の乳癌患者については、予後が不良であることが予測されることがわかる。

【0085】

また、病理学的な観察結果を真値として判別式の性能を評価した場合、感度８３．３％、特異度７０．４％、陰性適中率（ＮＰＶ）９３．４％、および陽性適中率（ＰＰＶ）４５．５％であることがわかる。

【0086】

したがって、これらの結果から、前記表１−１および表１−２に記載の遺伝子群の各遺伝子の発現量と前記判別式とを用いることによって、乳癌の予後を適切に予測することができることが示唆される。

【0087】

また、前記判別式により「再発なし」と予測された乳癌患者（低再発リスク群）および「再発あり」と予測された乳癌患者（高再発リスク群）それぞれについて、術後から観察し、Ｋａｐｌａｎ−Ｍｅｉｅｒプロットにより、非再発生存率を調べた。また、その結果をｌｏｇ−ｒａｎｋ検定により評価した。実施例２において、術後からの期間と非再発生存率との関係を図３に示す。

【0088】

図３に示された結果から、術後１０年間での非再発生存率は、高再発リスク群（図中、「高再発リスク」）では５３％であるのに対して、低再発リスク群（図中、「低再発リスク」）では９３％であることがわかる。また、ｌｏｇ−ｒａｎｋ検定によりｐ＝８．６ｅ−７であったため、低再発リスク群は、高再発リスク群と比べて有意に予後が良好であることがわかる。

【0089】

したがって、これらの結果から、前記表１−１および表１−２に記載の核酸（プローブセット）の発現量と前記判別式とを用いることによって、乳癌の予後を高い精度で予測することができ、乳癌の予後を適切に検査できることがわかる。

【0090】

［試験例１］
前記９５個の核酸を予後因子として用いる検査方法（実験番号１）、患者の閉経状態を予後因子として用いる検査方法（実験番号２）、腫瘍径を予後因子として用いる検査方法（実験番号３）、組織学的グレードを予後因子として用いる検査方法（実験番号４）、ＰＲの有無を予後因子として用いる検査方法（実験番号５）、ヒト上皮成長因子受容体２（ＨＥＲ２）の有無を予後因子として用いる検査方法（実験番号６）、全細胞中のＫｉ６７陽性細胞の割合が２０％以上であるかどうかを予後因子として用いる検査方法（実験番号７）または前記９５個の遺伝子とは異なる種類の９７個の遺伝子を予後因子として用いるＧｅｎｏｍｉｃ Ｇｒａｄｅ Ｉｎｄｅｘ（ＧＧＩ）による検査方法（実験番号８）のそれぞれについて、統計解析用ソフトウェア「Ｒ」にて使用する追加パッケージ「ｓｕｒｖｉｖａｌ ｖ２．３５−４」（ｈｔｔｐ：／／ｃｒａｎ．ｒ−ｐｒｏｊｅｃｔ．ｏｒｇ／ｗｅｂ／ｐａｃｋａｇｅｓ／ｓｕｒｖｉｖａｌ／ｉｎｄｅｘ．ｈｔｍｌ）を用い、多変量ＣＯＸ回帰ハザード解析を行なった。なお、前記ＧＧＩは、クリストス・ソテリウ（Ｓｏｔｉｒｉｏｕ Ｃｈｒｉｓｔｏｓ）らの文献〔ジャーナル・オブ・ザ・ナショナル・キャンサー・インスティチュート（Ｊｏｕｒｎａｌ ｏｆ ｔｈｅ Ｎａｔｉｏｎａｌ Ｃａｎｃｅｒ Ｉｎｓｔｉｔｕｔｅ）、２００６年、第９８巻、第４号、ｐ．２６２−２７２〕にしたがって求めた。

【0091】

その結果を表４に示す。なお、表中、「９５ｇｅｎｅｓ」は前記９５個の核酸を予後因子として用いる検査方法（実験番号１）、「Ｍｅｎｓ」は患者の閉経状態を予後因子として用いる検査方法（実験番号２）、「Ｔ」は腫瘍径を予後因子として用いる検査方法（実験番号３）、「ＨＧ」は組織学的グレードを予後因子として用いる検査方法（実験番号４）、「ＰｇＲ」はＰＲの有無を予後因子として用いる検査方法（実験番号５）、「ＨＥＲ２」はＨＥＲ２の有無を予後因子として用いる検査方法（実験番号６）、「Ｋｉ６７」は全細胞中のＫｉ６７陽性細胞の割合が２０％以上であるかどうかを予後因子として用いる検査方法（実験番号７）、「ｓｉｇｎ．ＧＧＩ」は前記９５個の遺伝子とは異なる種類の９７個の遺伝子を予後因子として用いるＧＧＩによる検査方法（実験番号８）をそれぞれ示す。各ハザード比は、表中、「参照」に該当する場合のハザードを１．０として算出した値である。前記閉経状態は、閉経前および閉経後で表される。また、従来、ｓｉｇｎ．ＧＧＩが低い場合、予後は良好であり、高い場合、予後は不良であるとされていた。

【0092】

【表6】

【0093】

表４に示される結果から、前記９５個の核酸を予後因子として用いた判別式を用いる検査方法（実験番号１；本発明の検査方法）は、ハザード比が７．７０であり、ｐ値が９．６Ｅ−０４であり、他の検査方法と比べて、精度が高い検査方法であることがわかる。

【0094】

［実施例３］
実施例２の（１）で得られた正規化後のデータについて、Ｓｐｅａｒｍａｎ順位相関係数およびｗａｒｄ法を用いて、階層的クラスター分析を行ない、樹形図を作成した。実施例３において、１０５例の乳癌患者それぞれにおける核酸（プローブセット）の発現量のデータの階層的クラスター分析を行なった結果を図４に示す。図４において、左側に、核酸（プローブセット）の発現量を表わすヒートマップを示し、右側に、観察結果による術後１０年間に再発あり（図中、「再発」）および術後１０年間に再発なし（図中、「非再発」）の判断結果を示す。

【0095】

図４に示された結果から、樹形図を分割するように付された太線を境界として、観察結果により術後１０年間に再発あり（予後不良）と判断された乳癌患者（上部に多い）と、観察結果により術後１０年間に再発なし（予後良好）と判断された乳癌患者（下部に多い）とを分けることができることがわかる。したがって、これらの結果から、前記表１−１および表１−２に記載の遺伝子群の各遺伝子の発現量を用いて、階層的クラスター分析を行なうことによって、乳癌の予後を高い精度で予測することができ、乳癌の予後を適切に検査できることがわかる。

【0096】

［実施例４］
実施例２の（１）で得られた正規化後のデータについて、表１−１および表１−２に記載の遺伝子を用いて主成分分析を行い、各遺伝子の変換係数を算出し、第一および第二主成分スコアを算出した。実施例４において算出された変換係数を表５−１および表５−２に示す。また、実施例４において、１０５例の乳癌患者の発現量のデータに基づき算出した第一主成分スコアおよび第二主成分スコアの散布図を図５に示す。図５中、ＰＣＡ１は、第一主成分スコアを示し、ＰＣＡ２は、第二主成分スコアを示す。図中、白丸は、観察結果により術後１０年間に「再発あり」と判断された乳癌患者、クロスは、観察結果により術後１０年間に「再発なし」と判断された乳癌患者である。

【0097】

【表7】

【0098】

【表8】

【0099】

図５に示された結果から、横軸の第一主成分スコアがゼロとなる所を境界として、観察結果により術後１０年間に「再発なし」（予後良好）と判断される乳癌患者と、観察結果により術後１０年間に「再発あり」（予後不良）と判断される乳癌患者とを分けることができることがわかる。

【0100】

したがって、これらの結果から、前記表１−１および表１−２に記載の遺伝子群の各遺伝子の発現量を用いて、主成分分析を行なうことによって、乳癌の予後を高い精度で予測することができ、乳癌の予後を適切に検査できることがわかる。

【図1】

【図2】

【図3】

【図5】

【図4】