(57)【要約】      （修正有）【課題】乳癌における遺伝子発現をゲノムワイドにかつ網羅的に解析した結果に基づき、乳癌患者の術後予後を遺伝子発現の観点から予測するシステムを提供する。【解決手段】ヒト遺伝子の遺伝子発現をＤＮＡマイクロアレイにより網羅的に解析し、様々な状態にある乳癌の遺伝子発現機能を比較することにより、乳癌術後予後予測システムを確立した。【選択図】なし【特許請求の範囲】
【請求項１】
乳癌の術後予後予測に関与する以下の定義の少なくとも１よりなる遺伝子；
１）エストロゲンレセプター陰性の乳癌において、外科手術後5年以内に死亡した乳癌患者からの遺伝子（5y-Dグループ）と数年以上無病（disease-free）生存した患者からの遺伝子（5y-Sグループ）とが、その発現機能によって区別することができたマーカー遺伝子群。
２）手術時にリンパ節への転移がなかった（node陰性）（n0）乳癌において、手術後5年以内に再発したn0乳癌患者からの遺伝子（5Y-Rグループ）と5年以上の間無病生存した患者からの遺伝子（5Y-Fグループ）とが、その発現機能によって区別することができたマーカー遺伝子群。
３）原発性乳癌において、外科手術後5年以内に死亡した乳癌患者からの遺伝子（5Dグループ）と数年以上無病生存した患者からの遺伝子（5Sグループ）とが、その発現機能によって区別することができたマーカー遺伝子群。
【請求項２】
原発性乳癌の術後予後予測に関与する以下の配列より選ばれる遺伝子；
pro-alpha-1 type 3 collagen (PIIIP)、
complement component C1r、
dihydropyrimidinase-like 3 (DPYSL3)、
protein tyrosine kinase 9-like (PTK9L)、
carboxypeptidase E (CPE)、
alpha-tubulin、
beta-tubulin、
heat shock protein HSP 90-alpha gene、
malate dehydrogenase、
NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 3 (NDUFB3)。
【請求項３】
原発性乳癌の術後予後予測に関与する予後の良い群で高発現する以下から選ばれる遺伝子；
pro-alpha-1 type 3 collagen (PIIIP)、
complement component C1r、
dihydropyrimidinase-like 3 (DPYSL3)、
protein tyrosine kinase 9-like (PTK9L)、
carboxypeptidase E (CPE)、
alpha-tubulin、
beta-tubulin。
【請求項４】
原発性乳癌の術後予後予測に関与する予後の悪い群で高発現する以下から選ばれる遺伝子；
heat shock protein HSP 90-alpha gene、
malate dehydrogenase、
NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 3 (NDUFB3)。
【請求項５】
手術時にリンパ節への転移がなかった（node陰性）（n0）乳癌において、術後予後予測に関与する以下の配列より選ばれる遺伝子；
AF058701/  DNA polymerase zeta catalytic subunit (REV3) 、
AI066764/  lectin, galactoside-binding, soluble, 1 (galectin 1)、
x15940/  ribosomal protein L31.、
Hs.94653/  neurochondrin(KIAA0607)、
M13436/  ovarian beta-A-inhibin、
Hs.5002/  copper chaperone for superoxide dismutase; CCS、
D67025/  proteasome (prosome, macropain) 26S subunit, non-ATPase, 3、
M80469/  MHC class I HLA-J gene、
Hs.4864/  ESTs、
Hs.106326/  ESTs。
【請求項６】
手術時にリンパ節への転移がなかった（node陰性）（n0）乳癌において、術後予後予測に関与する予後の悪い群で高発現する以下から選ばれる遺伝子；
AF058701/  DNA polymerase zeta catalytic subunit (REV3) 、
AI066764/  lectin, galactoside-binding, soluble, 1 (galectin 1)、
x15940/  ribosomal protein L31.。
【請求項７】
手術時にリンパ節への転移がなかった（node陰性）（n0）乳癌において、術後予後予測に関与する予後の良い群で高発現する以下から選ばれる遺伝子；
Hs.94653/  neurochondrin(KIAA0607)、
M13436/  ovarian beta-A-inhibin、
Hs.5002/  copper chaperone for superoxide dismutase; CCS、
D67025/  proteasome (prosome, macropain) 26S subunit, non-ATPase, 3、
M80469/  MHC class I HLA-J gene、
Hs.4864/  ESTs、
Hs.106326/  ESTs。
【請求項８】
エストロゲンレセプター陰性の乳癌において、術後予後予測に関与する以下の配列より選ばれる遺伝子；
Hs.108504/  FLJ20113/  ubiquitin-specific protease otubain 1
Hs.146550/  MYH9/  myosin, heavy polypeptide 9, non-muscle
Hs.194691/  RAI3/  retinoic acid induced 3
Hs.1975/  TDRD3/  tudor domain containing 3
Hs.203952/  TRRAP/  transformation/transcription domain-associated protein
Hs.278607/  GSA7/  ubiquitin activating enzyme E1-like protein
Hs.429/  ATP5G3/  
ATP synthase, H+transporting,mitochondrialF0complex,subunitc(subunit9)isoform3
Hs.75305/  AIP/  aryl hydrocarbon receptor interacting protein
Hs.81170/  PIM1/  pim-1 oncogene
Hs.99987/  ERCC2/
excision repaircross-complementingrodentrepairdeficiency,complementationgroup2
Y12781/  Transducin (beta) like 1 protein
Hs.104417/  KIAA1205 protein
cl.21783/  Hypothetical protein
Hs.112628/  Hypothetical protein: MGC43581
Hs.170345/  Hypothetical protein FLJ13710
Hs.53996/  weakly similar to zinc finger protein 135
Hs.55422/  Hypothetical protein
Hs.112718/  EST
Hs.115880/  EST
Hs.126495/  EST
【請求項９】
予後の悪い群で高発現する遺伝子である請求項８から選ばれる遺伝子。
【請求項１０】
請求項１〜９のいずれか一に記載の遺伝子が搭載されたＤＮＡマイクロアレイ。
【請求項１１】
請求項１〜９のいずれか一に記載の遺伝子に特異的なプローブが搭載されたＤＮＡマイクロアレイ。
【請求項１２】
ＤＮＡマイクロアレイが繊維型マイクロアレイである請求項１０または１１に記載のマイクロアレイ。
【請求項１３】
請求項１〜９のいずれか一に記載の遺伝子をマーカーにする乳癌の術後予後の検査方法。
【請求項１４】
請求項１０〜１２のいずれか一に記載のマイクロアレイを使用する乳癌の術後予後の検査方法。
【請求項１５】
請求項１〜９のいずれか一に記載の遺伝子をマーカーにする乳癌の術後予後を制御する癌治療薬のスクリーニング方法。
【請求項１６】
請求項１０〜１２のいずれか一に記載のマイクロアレイを使用する乳癌の術後予後を制御する癌治療薬のスクリーニング方法。
【請求項１７】
請求項１〜９のいずれか一に記載の遺伝子をマーカーにする試薬を含む乳癌の術後予後の診断キット。
【請求項１８】
請求項１６記載の診断キットがマイクロアレイを包含するものである診断キット。
【請求項１９】
マイクロアレイが繊維型マイクロアレイである請求項１８記載の診断キット。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、乳癌の術後予後予測に関与する遺伝子に関する。さらに、この遺伝子を使った乳癌の術後予後の検査方法、乳癌の術後予後を制御する癌治療薬のスクリーニング方法、及び乳癌の術後予後の診断キットに関する。
【背景技術】
【０００２】
乳癌は女性の癌死亡率の上位原因に位置する疾患であるが、いまだに生物学的見地からの悪性度、生存予後を規定する有力な因子は見出されていない。
【０００３】
エストロゲンレセプター（ＥＲ）の状態はヒト乳癌についての、臨床及び生物学的な症状の1つの決定要素である。 アジュバントホルモン治療は、年齢、閉経期の状態、 腋窩リンパ節（axillary  node）の関与、あるいは腫瘍径にかかわらずＥＲ陽性の乳癌患者において通常有効である；しかしながら、ＥＲ陰性乳癌は、この治療方法に対して抵抗性がある（非特許文献１−２）。ＥＲ陰性腫瘍を持っている患者が、化学療法に対して同じ応答をいつも示すものではない。そして、現存する指標では、臨床の症状によりこのタイプの乳癌を分別することができないことから、手術後の予後は、多様であると言える（非特許文献３−４）。
【０００４】
また、リンパ節転移の無い乳癌患者（node陰性乳癌；n0）の予後は、転移乳癌患者におけるよりは良い。しかし、日本において、本発明者らは、node陰性乳癌患者の16％が最初の手術後5年以内に再発することを見出している（非特許文献５）。
【０００５】
乳癌患者の術後予後予測は、現在利用できるアジュバント（adjuvant）治療の観点から重要性が増している。術後に再発しそうな患者を同定することに役立つ遺伝子マーカーは、ハイリスクな患者に適当な術前アジュバント療法を行い得る利益をもたらし、不要かつ煩雑で不快な副作用をもたらすことが阻止可能となる。
【０００６】
従来、個々の患者のための術後の処置決定は腫瘍径及びステージ、リンパ節への転移、臨床病理学的因子による診断、及びホルモンレセプターの検索などにより行われてきたが、決定的な方法ではなかった（非特許文献６−１１）。
【０００７】
近年、術後の乳癌患者の予後マーカーとして遺伝子の変異の重要性を決定しようとしたものがある。これらの遺伝子の変異にはp53の変異（非特許文献１２）、いくつかの対立遺伝子でのヘテロ接合性の欠失（非特許文献１３）、BRCA2遺伝子（非特許文献１４）、WT1遺伝子（非特許文献１５）、HER2/neu遺伝子（非特許文献１６）、及びKi-67遺伝子（非特許文献１７）の異常な発現を含む。しかしながら、癌が多遺伝子の異常の集積による疾患であることを考えると、有効な予後予測手段とは言いがたい。
【０００８】
さらに、近年各種生物におけるゲノムプロジェクトが進められており、ヒト遺伝子をはじめとして、多数の遺伝子とその塩基配列が急速に明らかにされつつある。配列の明らかにされた遺伝子の機能については、各種の方法で調べることができる。その有力な方法の一つとして、明らかにされ塩基配列情報を利用した遺伝子発現解析が知られている。例えば、ノーザンハイブリダイゼーションに代表されるような、各種の核酸−核酸間ハイブリダイゼーション反応や各種のＰＣＲ反応を利用した方法が開発され、当該方法により各種遺伝子とその生体機能発現との関係を調べることができる。これらの方法では適用し得る遺伝子の数に制限があるが、今日のゲノムプロジェクトを通して明らかにされつつあるような、一個体レベルという極めて多数の遺伝子の総合的・系統的解析を行うために、多数遺伝子の一括発現解析を可能とするＤＮＡマイクロアレイ法（ＤＮＡチップ法）と呼ばれる新しい分析法、及び方法論が開発されてきた。
【０００９】
ＤＮＡマイクロアレイとしては、リソグラフィー技術を応用して区画化された多数のセル上にＤＮＡ合成を行ったもの（特許文献１）、基盤上に溝又は穴で区画を形成し、該区画の内壁にプローブを固定化したもの（特許文献２−３）、チップ上に固定するプローブ量を多くするために、アクリルアミド等のゲルにプローブを固定したマイクロアレイ（特許文献４−５）等、多数の形状物が知られている。
【００１０】
本発明者らの一部は、先に新規なマイクロアレイ及びその製造法を開発した（特許文献６−７）。該発明は、核酸固定化ゲルを保持する核酸固定化ゲル保持繊維配列体を作製し、この配列体を配列体の繊維軸と交差する方向に切断することにより薄片をえるものである。この薄片は固定化核酸二次元高密度配列体、すなわちＤＮＡマイクロアレイとして使用される。
【００１１】
最近の研究により、癌診断のための新規な遺伝子マーカーの同定にｃＤＮＡマイクロアレイ技術が有効であることがわかった。現在までに、数人の研究者が乳癌のマイクロアレイ分析を行っているが、乳癌の術後予後予測の可能な乳癌遺伝子発現特性のデータを記述したものはない（非特許文献１８−２４）。一つの例外として、リンパ節の転移陰性腫瘍の特定のプロフィールが、遠隔転移への進行前の短い間隔を予測することを示す。（非特許文献２５）。
【特許文献１】USP 5445934
【特許文献２】特開平11-108928号
【特許文献３】特開2000-78998号
【特許文献４】USP5770721
【特許文献５】特開2000-60554号
【特許文献６】特開2000-270878号
【特許文献７】特開2000-270879号
【非特許文献１】J Clin Oncology (2001) 19, 3817-1827.
【非特許文献２】Breast Cancer (2001) 8,298-304.
【非特許文献３】J Natl Cancer Inst (1991) 83, 154-155.
【非特許文献４】J Natl Cancer Inst (2000) 93, 979-989.
【非特許文献５】Clin Cancer Res (2000) 6,3193-3198.
【非特許文献６】Cancer (1982) 50, 2131-2138.
【非特許文献７】Histopathology (1991) 19, 403-410.
【非特許文献８】Int JCancer (1996) 69, 135-141.
【非特許文献９】Am J Clin Oncol (1997) 20,546-551.
【非特許文献１０】Eur J Cancer (2002) 38, 1329-1334.
【非特許文献１１】Jpn JCancerRes (2000) 91,293-300.
【非特許文献１２】Breast Cancer Res Treat (2001) 69, 65-68.
【非特許文献１３】Int J Clin Oncol (2001) 6, 6-12.
【非特許文献１４】Int J Cancer (2002) I98, 879-882.
【非特許文献１５】Clin Cancer Res (2002) 8,1167-1171
【非特許文献１６】Arch Surg (2000) 135, 1469-1474.
【非特許文献１７】J Pathol (1999) 187, 207-216.
【非特許文献１８】Proc Natl Acad Sci U SA (1999) 96, 9212-9217.
【非特許文献１９】Nature (2000) 406, 747-752.
【非特許文献２０】Proc Natl Acad Sci U S A (2001) 98,11462-11467.
【非特許文献２１】Cancer Res (2001) 61, 5979-5984.
【非特許文献２２】Cancer Res (2000) 60, 2232-2238.
【非特許文献２３】Cancer Res (2001) 61, 5168-5178.
【非特許文献２４】Proc Natl Acad Sci U S A (2001) 98, 10869-10874.
【非特許文献２５】N Engl J Med (2002) 347, 1999-2009.
【発明の開示】
【発明が解決しようとする課題】
【００１２】
本発明の目的は、乳癌における遺伝子発現をゲノムワイドにかつ網羅的に解析した結果に基づき、乳癌患者の術後予後を遺伝子発現の観点から予測する画期的な手段を提供することを課題とする。
【課題を解決するための手段】
【００１３】
本発明は、ヒト遺伝子の遺伝子発現をＤＮＡマイクロアレイにより網羅的に解析し、様々な状態にある乳癌の遺伝子発現機能を比較することにより、乳癌術後予後予測システムを確立した。
【００１４】
すなわち、本発明は、以下の（１）〜（８）の遺伝子（群）である。
（１）  乳癌の術後予後予測に関与する以下の定義の少なくとも１よりなる遺伝子；
１）エストロゲンレセプター陰性の乳癌において、外科手術後5年以内に死亡した乳癌患者からの遺伝子（5y-Dグループ）と数年以上無病（disease-free）生存した患者からの遺伝子（5y-Sグループ）とが、その発現機能によって区別することができたマーカー遺伝子群。
２）手術時にリンパ節への転移がなかった（node陰性）（n0）乳癌において、手術後5年以内に再発したn0乳癌患者からの遺伝子（5Y-Rグループ）と5年以上の間無病生存した患者からの遺伝子（5Y-Fグループ）とが、その発現機能によって区別することができたマーカー遺伝子群。
３）原発性乳癌において、外科手術後5年以内に死亡した乳癌患者からの遺伝子（5Dグループ）と数年以上無病生存した患者からの遺伝子（5Sグループ）とが、その発現機能によって区別することができたマーカー遺伝子群。
（２）  原発性乳癌の術後予後予測に関与する以下の配列より選ばれる遺伝子；
pro-alpha-1 type 3 collagen (PIIIP)、
complement component C1r、
dihydropyrimidinase-like 3 (DPYSL3)、
protein tyrosine kinase 9-like (PTK9L)、
carboxypeptidase E (CPE)、
alpha-tubulin、
beta-tubulin、
heat shock protein HSP 90-alpha gene、
malate dehydrogenase、
NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 3 (NDUFB3)。
（３）  原発性乳癌の術後予後予測に関与する予後の良い群で高発現する以下から選ばれる遺伝子；
pro-alpha-1 type 3 collagen (PIIIP)、
complement component C1r、
dihydropyrimidinase-like 3 (DPYSL3)、
protein tyrosine kinase 9-like (PTK9L)、
carboxypeptidase E (CPE)、
alpha-tubulin、
beta-tubulin。
（４）  原発性乳癌の術後予後予測に関与する予後の悪い群で高発現する以下から選ばれる遺伝子；
heat shock protein HSP 90-alpha gene、
malate dehydrogenase、
NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 3 (NDUFB3)。
（５）  手術時にリンパ節への転移がなかった（node陰性）（n0）乳癌において、術後予後予測に関与する以下の配列より選ばれる遺伝子；
AF058701/  DNA polymerase zeta catalytic subunit (REV3) 、
AI066764/  lectin, galactoside-binding, soluble, 1 (galectin 1)、
x15940/  ribosomal protein L31.、
Hs.94653/  neurochondrin(KIAA0607)、
M13436/  ovarian beta-A-inhibin、
Hs.5002/  copper chaperone for superoxide dismutase; CCS、
D67025/  proteasome (prosome, macropain) 26S subunit, non-ATPase, 3、
M80469/  MHC class I HLA-J gene、
Hs.4864/  ESTs、
Hs.106326/  ESTs。
（６）  手術時にリンパ節への転移がなかった（node陰性）（n0）乳癌において、術後予後予測に関与する予後の悪い群で高発現する以下から選ばれる遺伝子；
AF058701/  DNA polymerase zeta catalytic subunit (REV3) 、
AI066764/  lectin, galactoside-binding, soluble, 1 (galectin 1)、
x15940/  ribosomal protein L31.。
（７）  手術時にリンパ節への転移がなかった（node陰性）（n0）乳癌において、術後予後予測に関与する予後の良い群で高発現する以下から選ばれる遺伝子；
Hs.94653/  neurochondrin(KIAA0607)、
M13436/  ovarian beta-A-inhibin、
Hs.5002/  copper chaperone for superoxide dismutase; CCS、
D67025/  proteasome (prosome, macropain) 26S subunit, non-ATPase, 3、
M80469/  MHC class I HLA-J gene、
Hs.4864/  ESTs、
Hs.106326/  ESTs。
（８）  エストロゲンレセプター陰性の乳癌において、術後予後予測に関与する以下の配列より選ばれる遺伝子；
Hs.108504/  FLJ20113/  ubiquitin-specific protease otubain 1
Hs.146550/  MYH9/  myosin, heavy polypeptide 9, non-muscle
Hs.194691/  RAI3/  retinoic acid induced 3
Hs.1975/  TDRD3/  tudor domain containing 3
Hs.203952/  TRRAP/  transformation/transcription domain-associated protein
Hs.278607/  GSA7/  ubiquitin activating enzyme E1-like protein
Hs.429/  ATP5G3/  
ATP synthase, H+transporting,mitochondrialF0complex,subunitc(subunit9)isoform3
Hs.75305/  AIP/  aryl hydrocarbon receptor interacting protein
Hs.81170/  PIM1/  pim-1 oncogene
Hs.99987/  ERCC2/
excision repaircross-complementingrodentrepairdeficiency,complementationgroup2
Y12781/  Transducin (beta) like 1 protein
Hs.104417/  KIAA1205 protein
cl.21783/  Hypothetical protein
Hs.112628/  Hypothetical protein: MGC43581
Hs.170345/  Hypothetical protein FLJ13710
Hs.53996/  weakly similar to zinc finger protein 135
Hs.55422/  Hypothetical protein
Hs.112718/  EST
Hs.115880/  EST
Hs.126495/  EST
【００１５】
また本発明は、上記（８）の内、予後の悪い群で高発現する遺伝子である。さらに本発明は、上記（１）〜（９）のいずれかの遺伝子が搭載されたＤＮＡマイクロアレイ、上記（１）〜（９）のいずれかの遺伝子に特異的なプローブが搭載されたＤＮＡマイクロアレイであり、好ましくは、ＤＮＡマイクロアレイが繊維型マイクロアレイである。
【００１６】
前記遺伝子は乳癌の術後予後のマーカー遺伝子として使用することができる。さらには、乳癌の術後予後を制御する癌治療薬のマーカー遺伝子としても使用することができる。
【００１７】
また、当該マーカー遺伝子は試薬として包含することができ、診断キットとして使用することができる。当該試薬キットにはマーカー遺伝子が搭載されたＤＮＡマイクロアレイ、好ましくは繊維型マイクロアレイを包含する。
【発明の効果】
【００１８】
本発明の手法により、完全に新規の乳癌関連遺伝子を見出すと同時に、それらの遺伝子が乳癌の悪性化に深く関与し、最終的には乳癌患者の予後に影響を及ぼしていることを発見した。さらには見出された遺伝子の発現状態を評価する数式を構築することにより、まったく新規かつ有効な乳癌術後予後予測システムを開発した。癌が遺伝子の異常による疾患であることを考慮すると、遺伝子発現の観点から見た本発明のシステムは従来の予後評価法とはまったく異なる、癌の生物学的本質を捉えた画期的な予後予測システムであると考えられる。
【発明を実施するための最良の形態】
【００１９】
本発明の一つの態様である乳癌の術後予後予測に関与するマーカー遺伝子群は、エストロゲンレセプター陰性乳癌、node陰性乳癌及び原発性乳癌において、外科手術後5年以内に死亡あるいは再発した患者、及び5年以上の間生存した患者からの遺伝子の発現機能をｃＤＮＡマイクロアレイで分析することにより得られたものである。
【００２０】
具体的には、本発明の態様の一つは、乳癌の術後予後予測に関与する既知配列から選ばれる以下の定義の少なくとも１よりなる遺伝子である；
１）エストロゲンレセプター陰性の乳癌において、外科手術後5年以内に死亡した乳癌患者からの遺伝子（5y-Dグループ）と数年以上無病（disease-free）生存した患者からの遺伝子（5y-Sグループ）とが、その発現機能によって区別することができたマーカー遺伝子群。
２）手術時にリンパ節への転移がなかった（node陰性）（n0）乳癌において、手術後5年以内に再発したn0乳癌患者からの遺伝子（5Y-Rグループ）と5年以上の間無病生存した患者からの遺伝子（5Y-Fグループ）とが、その発現機能によって区別することができたマーカー遺伝子群。
３）原発性乳癌において、外科手術後5年以内に死亡した乳癌患者からの遺伝子（5Dグループ）と数年以上無病生存した患者からの遺伝子（5Sグループ）とが、その発現機能によって区別することができたマーカー遺伝子群。
【００２１】
本発明の乳癌の術後予後予測に関与する遺伝子は、Random-permutationテスト、Mann-Whitneyテストを用いてｃＤＮＡマイクロアレイのデータを評価することにより得られたものである。本発明は遺伝子発現機能をｃＤＮＡマイクロアレイと半定量PCR実験とを組み合わせて評価することにより、臨床レベルにおいてより役立つアプローチを提示する。
【００２２】
本発明では、乳癌患者の遺伝子発現機能を評価することにより、原発性乳癌の術後予後予測に関与する遺伝子を同定した。
具体的には、本発明の態様の一つは原発性乳癌の術後予後予測に関与する既知配列から選ばれる以下の配列より選ばれる遺伝子である；
pro-alpha-1 type 3 collagen (PIIIP)、
complement component C1r、
dihydropyrimidinase-like 3 (DPYSL3)、
protein tyrosine kinase 9-like (PTK9L)、
carboxypeptidase E (CPE)、
alpha-tubulin、
beta-tubulin、
heat shock protein HSP 90-alpha gene、
malate dehydrogenase、
NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 3 (NDUFB3)。
【００２３】
上記の遺伝子のいくつかは腫瘍細胞の増殖又は遠隔転移に関連していると考えられる。例えば、heat shock protein HSP 90-alphaは、多くのキナーゼのシャペロンであり、癌細胞の成長を促進する可能性がある（Neckers,L.(2002)TrendsMolMed 8, S55-61.）。malate dehydrogenaseは、好気性及び嫌気性代謝に伴うエネルギーに関連する重要な酵素であり、malatedehydrogenaseの活性は扁平上細胞癌の腫瘍マーカーに関連する（Ross,C.D.,etal. (2000) Otolaryngol HeadNeck Surg 122, 195-200.）。NADHdehydrogenase(ubiquinone)1 beta subcomplex,3(NDUFB3)はミトコンドリアの電子伝達系に属しており、乳癌細胞株であるMDA-MB-231においてNDUFB3を含む領域の染色体異常が顕著である（Xie,D.,etal.(2002)IntJ Oncol 21, 499-507.）。
【００２４】
上記の原発性乳癌の術後予後予測に関与する10遺伝子は予後の良い群（5Sグループ）と予後の悪い群（5Yグループ）とで異なる発現を示し、10遺伝子のうち7遺伝子は予後の良い群（5Sグループ）で高発現する遺伝子である。
すなわち、本発明の態様の一つは、原発性乳癌の術後予後予測に関与する既知配列から選ばれる予後の良い群で高発現する以下から選ばれる遺伝子である；
pro-alpha-1 type 3 collagen (PIIIP)、
complement component C1r、
dihydropyrimidinase-like 3 (DPYSL3)、
protein tyrosine kinase 9-like (PTK9L)、
carboxypeptidase E (CPE)、
alpha-tubulin、
beta-tubulin。
【００２５】
上記の原発性乳癌の術後予後予測に関与する10遺伝子のうち3遺伝子は予後の悪い群（5Yグループ）で高発現する遺伝子である。すなわち、本発明の態様の一つは、原発性乳癌の術後予後予測に関与する既知配列から選ばれる予後の悪い群で高発現する以下から選ばれる遺伝子である；
heat shock protein HSP 90-alpha gene、
malate dehydrogenase、
NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 3 (NDUFB3)。
【００２６】
ここで、原発性乳癌の予測指標(PI)を以下のように定義し、乳癌の術後予後予測に用いることができる；
予測指標(PI)＝（乳癌組織における上記の予後の良い群で高発現する7遺伝子の正規化した発現比率の合計）−（乳癌組織における上記の予後の悪い群で高発現する3遺伝子の正規化した発現比率の合計）。
【００２７】
本発明では、乳癌患者の遺伝子発現機能を評価し、node陰性乳癌の術後予後予測に関与する10遺伝子を同定した。
具体的には、本発明の一つの態様は、手術時にリンパ節への転移がなかった（node陰性）（n0）乳癌において、術後予後予測に関与する既知配列から選ばれる以下の配列より選ばれる遺伝子である；
AF058701/  DNA polymerase zeta catalytic subunit (REV3) 、
AI066764/  lectin, galactoside-binding, soluble, 1 (galectin 1)、
x15940/  ribosomal protein L31.、
Hs.94653/  neurochondrin(KIAA0607)、
M13436/  ovarian beta-A-inhibin、
Hs.5002/  copper chaperone for superoxide dismutase; CCS、
D67025/  proteasome (prosome, macropain) 26S subunit, non-ATPase, 3、
M80469/  MHC class I HLA-J gene、
Hs.4864/  ESTs、
Hs.106326/  ESTs。
【００２８】
上記のnode陰性乳癌の術後予後予測に関与する遺伝子には腫瘍細胞の増殖、遠隔転移に関与する遺伝子が含まれる。例えば、galectin 1は細胞分化を制御するオートクライン型癌抑制因子である（AxelH,etal. (2003) Int. J. Cancer, 103: 370-379.）。また、癌転移を活性化する遺伝子が含まれる。
【００２９】
上記のnode陰性乳癌の術後予後予測に関与する10遺伝子は予後の良い群（5Y-Fグループ）と予後の悪い群（5Y-Rグループ）とで異なる発現を示し、10遺伝子のうち3遺伝子は予後の悪い群（5Y-Rグループ）で高発現する遺伝子である。すなわち、本発明の一つの態様は手術時にnode陰性乳癌において、術後予後予測に関与する既知配列から選ばれる予後の悪い群で高発現する以下から選ばれる遺伝子である；
AF058701/  DNA polymerase zeta catalytic subunit (REV3) 、
AI066764/  lectin, galactoside-binding, soluble, 1 (galectin 1)、
x15940/  ribosomal protein L31。
【００３０】
上記のnode陰性乳癌の術後予後予測に関与する10遺伝子のうち7遺伝子は予後の良い群（5Y-Fグループ）で高発現する遺伝子である。すなわち、本発明の一つの態様はnode陰性乳癌において、術後予後予測に関与する既知配列から選ばれる予後の良い群で高発現する以下から選ばれる遺伝子である；
Hs.94653/  neurochondrin(KIAA0607)、
M13436/  ovarian beta-A-inhibin、
Hs.5002/  copper chaperone for superoxide dismutase; CCS、
D67025/  proteasome (prosome, macropain) 26S subunit, non-ATPase, 3、
M80469/  MHC class I HLA-J gene、
Hs.4864/  ESTs、
Hs.106326/  ESTs。
【００３１】
ここで、node陰性乳癌の予後スコア(PS)を以下のように定義し、乳癌の術後予後予測に用いることができる；
予後スコア(PS)＝（乳癌組織における上記の予後の悪い群で高発現する3遺伝子の正規化した発現比率の合計）−（乳癌組織における上記の予後の良い群で高発現する7遺伝子の正規化した発現比率の合計）。
【００３２】
本発明では、乳癌患者の遺伝子発現機能を評価することにより、エストロゲンレセプター陰性の乳癌の術後予後予測に関与する20遺伝子を同定した。
具体的には、本発明の一つの態様は、エストロゲンレセプター陰性の乳癌において、術後予後予測に関与する既知配列から選ばれる以下の配列より選ばれる遺伝子である；
Hs.108504/  FLJ20113/  ubiquitin-specific protease otubain 1
Hs.146550/  MYH9/  myosin, heavy polypeptide 9, non-muscle
Hs.194691/  RAI3/  retinoic acid induced 3
Hs.1975/  TDRD3/  tudor domain containing 3
Hs.203952/  TRRAP/  transformation/transcription domain-associated protein
Hs.278607/  GSA7/  ubiquitin activating enzyme E1-like protein
Hs.429/  ATP5G3/    ATPsynthase,H+transporting,mitochondrialF0complex,subunitc(subunit9)isoform 3
Hs.75305/  AIP/  aryl hydrocarbon receptor interacting protein
Hs.81170/  PIM1/  pim-1 oncogene
Hs.99987/  ERCC2/  excisionrepaircross-complementingrodentrepairdeficiency,complementationgroup2
Y12781/  Transducin (beta) like 1 protein
Hs.104417/  KIAA1205 protein
cl.21783/  Hypothetical protein
Hs.112628/  Hypothetical protein: MGC43581
Hs.170345/  Hypothetical protein FLJ13710
Hs.53996/  weakly similar to zinc finger protein 135
Hs.55422/  Hypothetical protein
Hs.112718/  EST
Hs.115880/  EST
Hs.126495/  EST。
【００３３】
上記のエストロゲンレセプター陰性乳癌の術後予後予測に関与する遺伝子には腫瘍細胞の増殖又は遠隔転移に関連する遺伝子が含まれる。例えば、PIM1はセリン/スレオニンキナーゼであり、前立腺がんの臨床結果とその発現には相関性がある（Oesterreich,S.,etal.(1996) Clin Cancer Res, 2, 1199-1206.）。また、TRRAP蛋白質は哺乳類HAT複合体のサブユニットであり、TRRAPのアンチセンスRNAは乳癌細胞のエストロゲン依存性の成長を阻害する。
【００３４】
上記のエストロゲンレセプター陰性乳癌の術後予後予測に関与する20遺伝子は、予後の悪い群（5y-Dグループ）で高発現する。すなわち、本発明の一つの態様は、予後の悪い群で高発現する、前記のエストロゲンレセプターに対して陰性の乳癌において、術後予後予測に関与する既知配列から選ばれる遺伝子である。
【００３５】
ここで、上記のエストロゲンレセプター陰性乳癌の術後予後予測に関与する遺伝子の発現に基づいて、以下のように乳癌の術後予後を予測することができる；
（１）上記のエストロゲンレセプター陰性乳癌の術後予後予測に関与する20遺伝子の乳癌組織における発現量を母集団における平均値と比較し、各遺伝子発現量が母集団の平均値の２倍以上であるならば、1点を付与、
（２）（１）の操作を各20遺伝子について行い、合計点数が8点以上ならば、予後不良とする。
【００３６】
上記の乳癌の術後予後予測に関連する遺伝子は、マーカーとして乳癌の術後予後の検査に用いることができる。すなわち、本発明の一つの態様は、前記遺伝子をマーカーにする乳癌の術後予後の検査方法である。
【００３７】
上記の乳癌の術後予後予測に関連する遺伝子は、マーカーとして乳癌の術後予後を制御する癌治療薬のスクリーニングに用いることができる。すなわち、本発明の一つの態様は、前記の遺伝子をマーカーにする乳癌の術後予後を制御する癌治療薬のスクリーニング方法である。
【００３８】
上記の乳癌の術後予後予測に関連する遺伝子は、マーカーとして乳癌の術後予後の診断に用いることができる。また、上記遺伝子に対して特異的なプローブを設計し、それらプローブをマーカーとして使用することもできる。それらプローブは例えば、ダイナコム社製ProbeQuest（登録商標）により設計することが可能である。すなわち、本発明の一つの態様は、前記の遺伝子をマーカーにする試薬を含む乳癌の術後予後の診断キットである。
【００３９】
上記の診断キットはマイクロアレイを包含することができる。すなわち、本発明の一つの態様は、前記診断キットがマイクロアレイを包含するものである診断キットである。
【００４０】
上記マイクロアレイを包含するものである診断キットのマイクロアレイには繊維型マイクロアレイが含まれる。ここで、繊維型マイクロアレイの調製方法については、前記特許文献６−７を引用する。すなわち、本発明の一つの態様は、マイクロアレイが繊維型マイクロアレイである前記診断キットである。
【００４１】
次に、実施例により本発明の実施の態様を具体的に説明するが、本発明は実施例に限定されるものではない。
【実施例１】
【００４２】
エストロゲンレセプター陰性乳癌における術後予後予測の遺伝子発現機能の評価
【００４３】
（組織サンプル）
日本医科大学並びに癌研究会の倫理委員会より承認されたガイドラインに従ってインフォームドコンセントを得た後、癌研究会付属病院（東京）において1995-1997年に手術を受けた乳癌患者から原発性乳癌及び隣接する正常乳腺からの組織を採取した。組織は速やかに凍結され、-80°Cで保存された。954人の患者について、その全員を5年以上の間あるいは死亡時まで臨床的に追跡し、手術後5年以内に死亡したエストロゲンレセプター陰性乳癌の患者10人（5y-D）、及び5年以上の間無病生存した患者10人（5y-S）から試料を選んだ。両方の患者グループの臨床のバックグラウンドは、年齢、リンパ節転移、腫瘍径と組織型において一致させた（表１）。  
【００４４】
（乳癌の20症例の臨床的特徴）
【００４５】
【表１】
2005270093.tif 000002
【００４６】
すべての患者は癌研究会付属病院の「乳癌のための手術後の臨床のプロトコル」に従って手術後のアジュバント治療を受けた。 各症例での アジュバント治療の選択は、外科手術のタイプ、リンパ節関与の状態、及び局所あるいは遠隔転移の存在に基づいて厳密に決定した。本発明の研究においては、患者のいずれもがアジュバント化学療法の前に遠隔転移を持っておらず、外科手術の前に放射線療法あるいは化学療法を受けていなかった。
【００４７】
（臨床病理（Clinicopathological）パラメータ）
次のパラメータを調べた：組織型、腫瘍径及び浸潤（t因子）、リンパ節関与（lymph node involvement）、及びエストロゲンレセプター（ER）とプロゲステロンレセプター（PgR）の状態。腫瘍をTNM分類と日本の乳癌学会（1989）の組織分類によって、次のタイプに分類した；noninvasivetubular（1a）、invasivepapillotubular（a1）、invasivesolid-tubular（a2）、invasivescirrhouscarcinoma（a3）、及び他の特別なタイプ（b）。分類は基本的に世界保健機構の乳癌組織分類と同じである。ｔ因子は、組織学的TNM分類に従って次のタイプに分類した；最大寸法が2ｃｍ以下の腫瘍（t1）、皮膚又は胸筋への浸潤のない、最大寸法が2ｃｍ以上の腫瘍（t2）、皮膚あるいは胸筋への浸潤のあるもの（t3）。
【００４８】
（cDNAマイクロアレイのデザインと構成）
UniGene データベースから選択した25,344の cDNAsにより、” ゲノムワイド cDNA マイクロアレイ“を構築した。当該cDNAsは種々のヒト器官から分離されたpoly(A)+RNAを使いRT-PCRで作成された。PCR産物を、ArraySpotterGenerationIII（Amersham Biosciences）を使ってタイプ7のスライド・ガラス（AmershamBiosciences UKLimited、Buckinghamshire、UK）にスポットした。各スライドは、384のハウスキーピング遺伝子を含む。
【００４９】
（RNAの調製及び増幅）
腫瘍原料を採取後直ちに-80℃で急速凍結した。RNAを、TRIzol（Invitrogen Inc.、Carlsbad、CA、 USA）を使って抽出し、さらにRNeasykits(QuiagenInc.、Valencia、CA)を使って精製した。各RNAの純度を、分光測光法及び1.2％変性ホルムアミドゲル上で電気泳動することにより評価した。高純度RNAは、1.8-2.0の吸光度比（260nm/280nm）を持ち、ホルムアミドゲル電気泳動上で28S/18Sリボゾーマルバンドが1.8以上の比率を有するサンプルとして定義した。1ユニットのDNaseI（EpicentreTechnologies、Madison、WI）（1unit/μl）で処理後、出発原料として各試料からRNAの2μgを用いてT7RNAポリメラーゼによるRNA増幅を行った。増幅を2回行い、RNeasykits（QuiagenInc.、Valencia、CA）で、増幅されたRNA（aRNA）を精製した。各aRNAの量を分光光度計により測定し、その品質をホルムアミドゲル電気泳動によりチェックした。
【００５０】
（aRNAの標識、ハイブリダイゼーション及びスキャニング）
マイクロアレイ分析のcDNAを、aRNAから調製した。乳癌及び正常乳腺組織からのaRNA（5〜10μg）を、aminoallyl-cDNAlabelingkits（Ambion、Austin、TX）を使いCy5（癌試料）とCy3（正常試料）で標識した。Cy3-とCy5-標識cDNAプローブを混合し、95℃で5分加熱した後、30秒氷で急冷し、マイクロアレイにハイブリダイゼーションさせた。混合されたプローブを、microarrayhybridizationsolutionversion 2（Amersham Biosciences UK Limited、Buckinghamshire、UK）を50％終濃度のホルムアミド(Sigma-AldrichCorp.、St.Louis、MO、USA)に添加した。15時間40℃でのハイブリダイゼーションの後に、マイクロアレイスライドを、最初に1xSSCと0.2％SDSによって、10分間55℃で洗浄し、次いで2回0.1xSSC/0.2％SDSで各1分間室温で洗浄した。全ての処理は、AutomatedSlideProcessorSystem（Amersham）で行った。各ハイブリダイゼーションのシグナル強度を、Gene Pix 4000A（AxonInstruments,Inc.、FosterCity、CA、USA）でスキャニングし、Gene Pix 3.0 (Axon Instruments)で分光測光法により評価した。スキャンされたシグナルを、以下の文献記載の方法（thetotalgenenormalization method）で正規化した（Yang YH, Dudoit S, Luu P, et al. (2002)NucleicAcidsRes 30,e15; Manos EJ, Jones DA.(2001) Cancer Res 61: 433-438）。
【００５１】
（シグナル分析及び異なる発現を示す遺伝子の選択）
各ハイブリダイゼーションのシグナル強度を、Gene Pix Pro 3.0（Axon Instruments, Inc.、Foster City、 CA、USA）により測光法で評価した。癌とコントロール間のmRNA発現量を正規化するために、各遺伝子発現のCy5：Cy3比率を調整した。その結果、ハウスキーピング遺伝子の平均化されたLog（Cy5：Cy3比率）はゼロであった。27個のハウスキーピング遺伝子は、Webサイトhttp://www.nhgri.nih.gov/DIR/LCG/ARRAY/expn.htmlのハウスキーピングパネルから援用した。各マイクロアレイスライドについて、（Ｓ／Ｎ）比のカットオフ値を3.0に設定した。そして、Cy3とCy5のシグナル強度が、カットオフ値より低い遺伝子は検討から除外した。
【００５２】
（Mann-Whitney テスト）
5y-Dと5y-S腫瘍間で、明らかに異なる発現を示した遺伝子を検討するために、Mann-Whitneyテストを、一連のサンプルXに適用した。Xは、各遺伝子及び各サンプルのCy5/Cy3シグナル強度比率である（OnoK,TanakaT, Tsunoda T, et al. (2000) Cancer Res 2000; 60: 5007-5011）。 U値を、両グループの少なくとも5サンプルで有意なシグナルを与えた各遺伝子について計算した。Ｕ値が23より低いか、あるいは77より大きい遺伝子を選択した。U値は、各Ｘ値に基づく各遺伝子の、5y-Dグループに対する5y-Sグループを計算しているので、23より低いU値は、5y-Dグループに比べ5y-Sグループにおいて高発現するものと評価した。しかしながら、77より大きなU値を持つ遺伝子は、5y-Sグループに比べ5y-Dグループにおいて高発現するものと評価した。この基準によると、183遺伝子が5y-Sグループで高発現しており、31遺伝子が5y-Sグループで高発現していた。したがって、2つのグループ間での中間発現値が2倍以上の差異を示す遺伝子のみを（μXD/μXS  <  0.5あるいは  >  2.0、μXDとμXSがそれぞれ5y-Dあるいは5y-Sグループの平均Ｘ値を示す）予後関連の遺伝子と定義した。その結果、全部で110遺伝子が選ばれた。そのうち90遺伝子が5y-D腫瘍群で高いレベルで発現しており、20遺伝子が5y-S腫瘍群で高いレベルで発現していた。
【００５３】
（Random-permutationテスト）
さらにMann-Whitneyテストによって選択された110遺伝子の価値を評価するために、permutationテストを行った。そしてグループ差異に関連する遺伝子の可能性、Ps、もまた想定した。各遺伝子を、発現ベクターｖ（ｇ）＝（X1、X2、---、X20）（Xiは最初のサンプル群におけるi番目のサンプルの遺伝子発現レベルを示す）で表すとき、理想的発現パターンはc＝（c1、c2 、---、c20）（i番目のサンプルがS又はDグループに属するかによりci＝＋1又は0となる）で表現される。
遺伝子と グループ差異Pgc間の相関は、次のように定義された：すなわち、 Pgc  ＝ （μＳ−μD） / （ δS  ＋δD  ）；μＳ（μD）とδS（δD）は、新規に定義されたS（又はD）グループにおいて、各サンプルの該遺伝子“ｇ”のlog2Xの標準偏差を示す。
permutationテストは、ｃの座標を入れ換えることにより行った。すべてのpermutation間で相関値、 Pgcを計算した。これらの手順は10000回にわたり繰り返し行った。偶然に、2つのグループを分類する遺伝子の可能性を示すp値が、選択された110遺伝子の各々について評価された。最終的に、5y-Dケースで高発現する71遺伝子と、5y-Sケースで低発現する15遺伝子を選別した。
【００５４】
（半定量(Semi-quantitative)RT-PCR）
RNA（2μｇ）を、DNase I（Epicentre Technologies、Madison、WI、USA）で処理し、ReverscriptIIreversetranscriptase（和光純薬(株)、大阪、日本）とオリゴ（dT）12-18プライマーを使って単鎖ｃDNAを逆転写させた。単鎖cDNAsは、量的なコントロールとしてGAPD（glyceraldehyde-3-phosphatedehydrogenase）の発現をモニタリングすることにより次のPCR増幅のために濃度調整を行った。各PCRは、GeneAmpPCRシステム9700（AppliedBiosystems、Foster City、CA、USA）を用い1xPCRバッファー30μｌ量で、以下の反応条件で行った；
94℃5分、
（94℃30秒、60℃30秒、及び72℃30分）を25-35サイクル。
RT-PCRに使ったプライマー配列は以下である:  
配列番号1   GAPD (control) forward, 5'-GGA AGGTGA AGG TCG GAG T-3' 
配列番号2   reverse, 5'-TGG GTG GAA TCA TAT TGGAA-3'; 
配列番号3   Hs.108504F, 5’-ACA CTT CAT CTG CTCCCT CAT AG-3’;
配列番号4   Hs.108504R, 5’CTG CCT AGA CCT GAGGAC TGT AG-3’;
配列番号5   Hs.146550F, 5’ACT GAG GCC TTT TGGTAG TCG-3’; 
配列番号6   Hs.146550R, 5’TCT CTT TAT TGT GATGCT CAG TGG-3’; 
配列番号7   Hs.194691F, 5’AAA TCC TTC TCG TGT GTTGAC TG-3’; 
配列番号8   Hs.194691R, 5’CAG TCA TGA GGG CTA AAAACT GA-3’; 
配列番号9   Hs.1975F, 5'GAA GAC AAC AAG TTT TAC CGG G-3';
配列番号10   Hs.1975R, 5’ATG GTT TTA TTG ACG GCAGAA G-3’; 
配列番号11   Hs.203952F, 5’AGG ACA CGT CCT CTCCTC TCT C-3’; 
配列番号12   Hs.203952R, 5’TAA AGC TAG CGA AGGAAC GTA CA-3’; 
配列番号13   Hs.278607F, 5'TCC CTT CTG TTT CCT CAG TGT T-3'; 
配列番号14   Hs.278607R, 5'CCT GCC CCG ATA AAA ATA TCT AC -3'; 
配列番号15   Hs.429F, 5’TTG ACC TTA AGC CTC TTTTCC TC-3’; 
配列番号16   Hs.429R, 5’ATA ACG TAC ATT CCC ATGACA CC-3’; 
配列番号17   Hs.75305F, 5’ACT TTC AAG ATG GGACCA AGG-3’; 
配列番号18   Hs.75305R, 5’ATA TAC ACA GAA GCATGA CGC AG-3’; 
配列番号19   Hs.81170F, 5’TTG CTG GAC TCT GAAATA TCC C-3’; 
配列番号20   Hs.81170R, 5’TTC CCC TGT ACA GTATTT CAC TCA-3’; 
配列番号21   Hs.99987F, 5’GTT CCA GAT TCG TGA GAA TGA CT-3’; 
配列番号22   Hs.99987R, 5’CTG AGC AAT CTG CTC TAT TCC C-3’; 
配列番号23   Y12781F, 5’ACC AGT AAC AAC TGT GGGATG G-3’; 
配列番号24   Y12781R, 5’CAA ATG AGC TAC AAC ACACAA GG-3’; 
配列番号25   Hs.104417F, 5’CCC CCT CCA CCT TGTACA TAA T-3’; 
配列番号26   Hs.104417R, 5’GTT TTC GTT TGG CTGGTT GTG-3’; 
配列番号27   cl.21783F, 5’GTC TGA GAT TTT ACTGCA CCG-3’; 
配列番号28   cl.21783R, 5’GGA TGG AGC TGG AGGATA TTA-3’; 
配列番号29   Hs.112628F, 5’ATT GCT AAG GAT AAGTGC TGC TC-3’; 
配列番号30   Hs.112628R, 5’TGT CAG TAT AGA AGCCTG TGG GT-3’; 
配列番号31   Hs.170345F, 5’GCA TCT GAA TGT CTT TCT CCC TA-3’; 
配列番号32   Hs.170345R, 5’TTC TTA GGC CAT CCC TTT TCT AC-3’; 
配列番号33   Hs.53996F, 5’CCA TAG GAT CTT GACTCC AAC AG-3’; 
配列番号34   Hs.53996R, 5’ACT GGG AGT GGA GGAAAT TAG AG-3’; 
配列番号35   Hs.55422F, 5’CTA ATG TAA GCT CCATTG GGA TG-3’; 
配列番号36   Hs.55422R, 5’CAA ACT GCA AAC TAGCTC CCT AA-3’; 
配列番号37   Hs.112718F, 5'AAG ACT AAG AGG GAA AAT GTG GG-3'; 
配列番号38   Hs.112718R, 5'AGG TAA CCC AAA GTG ACA AAC CT-3'; 
配列番号39   Hs.115880F, 5'TTA AGT GAG TCT CCT TGG CTG AG-3'; 
配列番号40   Hs.115880R, 5'AGG GCC CCT ATA TCC AAT ACC TA-3'; 
配列番号41   Hs.126495F, 5’GAT CTT TCA AGA TGAGCC AAG GT-3’; 
配列番号42   Hs.126495R, 5’AGT CAT TCA GAA GCCATT GAG AC-3’
【００５５】
（RT-PCR産物のシグナル強度測定と予後スコアの計算）
PCR産物を、2％のアガロースゲル電気泳動とエチジウムブロマイド染色によって検出した。ゲルをSpot Density法でデジタル画像処理システム（AlphaImager3300；AlphaInnotech、San Leandro、CA、USA）でスキャンした。各バンドの2次元領域を構築し、その密度をIDV（IntegratedDensityValue）と定義したピクセル強度(遺伝子発現)を得た。各グループにおけるIDVの差異の重要性を、Student'st-テストで評価した。その結果、t-テストで0.05以下のp値を示した20遺伝子を候補として選択した（表２）；すなわち、当該20遺伝子の発現レベルは、5y-Dグループにおいて5y-Sグループに対して有意に高値であった。この情報を用いて、本発明者らは術後予後を予測するスコアリングシステムの確立を試みた。その際、各遺伝子は、各サンプルの発現レベルが20サンプルの平均発現レベルより高いか否かにより決定した。もしサンプルの発現レベルが、平均の2倍以上であるならば、付加的に＋1点を与えた。次に各サンプルのために全票（予後スコア）を得るためにすべての20遺伝子の点を合計した。その結果、8点以上のサンプルの場合は、悪い予後の表示と評価した。他方、8点以下の場合は、好ましい予後を表示していると評価した。
【００５６】
（予後スコアリングシステムの20候補遺伝子）
【００５７】
【表２】
2005270093.tif 000003
【００５８】
（結果）
エストロゲンレセプター陰性乳癌組織で、有意に高発現している257遺伝子を明らかにし、同様に低発現している378遺伝子を明らかにした。5y-Dと5y-Sグループ間で異なる発現を示す遺伝子を同定するために、マイクロアレイのデータを、Mann-WhitneyテストとRandom-permutationテストにより分析した。
その結果、5y-D腫瘍において全71遺伝子（10 EST と仮想タンパク質をコードしている9遺伝子を含む）が、共通して高発現のグループに分類された。これに対して、15遺伝子（3ESTを含む）が共通して低発現のグループに分類された（図１）。
【００５９】
5y-Dグループで高発現する遺伝子は、癌細胞の増殖及び転移に関与する以下の遺伝子を含む；matrix metalloproteinase 2（MMP2）、heatshockprotein27 (HSPB 1、Pim-1oncogene （PIM1）及びtransformation/transcriptiondomain-associatedprotein(TRRAP)。
5y-Dグループで低発現する遺伝子は、 HLA-C (major histocompatibility complex, class I, C)及び特異的なキナーゼの遺伝子を含む。ＤＮＡ修復、転写、シグナル伝達、細胞骨格及び接着性と関係がある多くの遺伝子が、2つのグループ間で異なる発現を示した。
【００６０】
マイクロアレイのデータの信頼性を確かめるために、5y-Dグループで高発現する20遺伝子を選び（ Hs.108504 、Hs.146550 、Hs.194691、Hs.1975、Hs.203952、Hs.278607、Hs.429、Hs.75305 、Hs.81170 、Hs.99987 、Y12781、Hs.104417、 cl.21783、Hs.112628 、Hs.170345、Hs.53996、Hs.55422、Hs.112718、Hs.115880及びHs.126495 ）、半定量RT-PCR によって該遺伝子の発現レベルを調べた。その結果はマイクロアレイのデータと一致し、5y-Dと5y-Sグループを区別するための統計学的な意義を有していた（代表的なデータを図２で示す）。
【００６１】
マーカー遺伝子の発現プロフィールを使って、手術後の予後を予測するスコアリングシステムを構築するために、前述の方法で予後スコアを計算した。簡略には、マーカー遺伝子が次の基準に従って選択された；
（１）調べた症例の少なくとも60％でカットオフレベルより高いシグナル強度を示す；
（２）│μD − μS │が1.0以下である。ここでμD（μS）は、5y-D（5y-S）のケースでの対数変換相対発現比率から導かれる平均値示す。
【００６２】
次に、発現機能により 5y-Dグループと 5y-Sグループとを区別することができたマーカー遺伝子を同定するためにMann-Whitney テストとRandom-permutationテストを行った。関連するマイクロアレイの結果を、半定量RT-PCR実験により確認した。Student'st-テストにより、20遺伝子を予後マーカーとして選別した（表２）。
【００６３】
本発明の予後スコア（PS）により、20人の患者を、予後不良と予測される10人（PSが11以上）と、予後良好と予測される10人（PSが11未満）に分けた。その結果、手術後の経過との対比により、本発明のスコアリングシステムは5y-Dケースでは80％、5y-Sケースでは100％の正確さにおいて信頼性があることを示した（図３Ａ）。
【００６４】
本発明の予後スコアリングシステムを使い、追加5症例について調べた（図３Ｂ ）。 当該システムは、2ケース（PS ＞ 11；患者TD-1、及びTD-2）で不良予後を、3ケース（PS＜11；患者TD-3、TS-1、及びTS-2）で良好予後を予測した。その結果、当該スコアリングシステムは、これら5症例の実際の臨床結果に関して80％の正確性であった。
【実施例２】
【００６５】
node陰性の乳癌における術後予後予測の遺伝子発現機能の評価
【００６６】
（組織サンプル）
組織サンプルは実施例１で記述した手法と同様に採取した。手術後5年以内に再発したnode陰性（ｎ0）癌の患者（5Y-R）12人、及び5年以上の間無病生存した患者（5Y-F）12人からの腫瘍について、遺伝子発現を検討した。両方の患者グループの臨床バックグラウンドは、年齢、リンパ節転移、腫瘍径、ホルモン受容体の状態、及び病理組織において一致させた（表３）。フォローアップの中間期間は、7.8年、最初の手術と再発の間の平均期間は、5Y-Rグループで2.7年であった。尚、すべての患者は実施例１で記載したアジュバント治療を受けた。
【００６７】
（臨床病理データ）
【００６８】
【表３】
2005270093.tif 000004
a  a1：invasive papillotubular carcinoma、a2：invasivesolid-tubularcarcinoma、a3：invasiveschirrhouscarcinoma 
b  TNM分類：日本乳癌学会のTNM分類に従って臨床学的に分類された
c  D.F.I：発病しない期間（disease free interval）
【００６９】
（臨床病理（Clinicopathological）パラメータ）
実施例１で記述した手法で臨床病理パラメータを調べた。組織学的グレードは、Elastonand Ellis（Abrams JS.BreastCancer2001; 8: 298-304.）の方法により評価した。リンパ管浸潤は、欠如か陽性で評価した（例えば、癌周辺のリンパ管に癌細胞が一つ又はそれ以上存在する場合、陽性と評価した）。脂肪浸潤（Fatinvasion）は、欠如か陽性かで評価した（例えば癌が、間質組織にまで浸潤している場合、陽性と評価した）。
【００７０】
（cDNAマイクロアレイの調製）
25,344の cDNAsを有する、”ゲノムワイド cDNAマイクロアレイキット（Amersham Biosciences UK Limited、Buckinghamshire、UK）“を使用した。PCR産物は、ArraySpotterGenerationIII（Amersham Biosciences）を使い、タイプ7ガラススライド（AmershamBiosciences）上にスポットした。
【００７１】
（RNAの調製と増幅）
実施例１で記述した手法と同様にRNAの調製と増幅を行った。
【００７２】
（aRNAの標識、ハイブリダイゼーション及びスキャニング）
実施例１で記述した手法と同様にaRNAのラベリング、ハイブリダイゼーション及びスキャニングを行った。
【００７３】
（Mann-Whitney テスト）
無病及び再発グループ間で異なる発現を示す遺伝子を同定するために、正規化したシグナルを、一連のXに適用されたMann-Whitneyテストにより分析した。ここでXは、各遺伝子と各試料についてのCy5/Cy3signal強度比である。2つのグループ間で、発現強度において2倍以上の差異を示す遺伝子を選んだ。3.0以下のシグナル−ノイズ比の遺伝子は分析から除外した。
U値を、両グループの少なくとも5試料で有意のシグナルを与えた各遺伝子について計算した。37より低い或いは107より大きいU値を持つ遺伝子を選択した。U値は、各X値に基づき各遺伝子について5Y-Rグループに対する5Y-Fグループについて計算したので、37より低いU値を持つ遺伝子は、5Y-Rグループに比べて5Y-Fグループで、高発現すると判断した（第一カテゴリー）。一方、107より大きいU値を持つ遺伝子は、5Y-Fグループに比べて5Y-Rグループで、高発現すると判断された（第二カテゴリー）。
この方法で、第一カテゴリーで78遺伝子、第二カテゴリーで55遺伝子を同定した。そこで、2つのグループ間の中間発現値の2倍以上の差異を示したもののみを予後関連遺伝子と定義した（μXR/μXF  <0.5あるいは>2.0、ここでμXRとμXFは、各5Y-R又は5Y-Fグループの平均的X値を示す）。全部で、98遺伝子が選ばれ、そのうち64遺伝子が5Y-F腫瘍で高い発現レベルを示し、34遺伝子が5Y-R腫瘍で有意に高い発現レベルを示した。
【００７４】
（Random-permutationテスト）
Mann-Whitney テストで選ばれた遺伝子の価値を評価するために、permutation テストを行い、そして各選択された遺伝子のグループ差異（Ps）への相関があることを評価した。各遺伝子が、発現ベクターｖ（ｇ）＝（X1、X2、----、X24）（Xiがサンプルの最初のセットでiサンプルの遺伝子の発現レベルを示す）で表現されるとき、理想化した発現パターンがｃ＝（c1、c2、----、c24）（iサンプルがFかRグループに属するかでci＝＋1又は0となる）で表現される。
遺伝子と グループ差異Pgc間の相関は、次のように定義された：すなわち、 Pgc ＝ (μF+μR)/(sF+sR)；μF（μR）と sF（sR）は、新規に定義された”F”又は“R”グループにおいて、各サンプルの該遺伝子”ｇ”のlog2Xの標準偏差を示す。
permutationテストは、ｃの座標を入れ換えることによって行われた。すべてのpermutationの間に、相関値、 Pgcsを計算した。これらの手順は10000回、繰り返して行われた。偶然に、2つのグループを分類する遺伝子の可能性を暗示するp値が、選択された58遺伝子の各々のために評価された。
【００７５】
（半定量RT-PCR）
RNA（5μｇ）を、DNase I （Epicentre Technologies、Madison、WI、USA）（1unit/μl ）で処理後、ReverscriptIIreversetranscriptase（和光純薬(株)、大阪、日本）と0.5μg/μloligo（dT）12-18プライマーを使って単鎖cDNAを逆転写させた。単鎖cDNAsの各調製物を、量的なコントロールとしてGAPDHをモニタリングすることにより、次のPCR増幅のために希釈した。全てのPCRは、GeneAmpPCRシステム9700（AppliedBiosystems、Foster City、CA、USA）を用い、1xPCRバッファー30μl量で以下の反応条件により行った；
94℃2分、
（94℃30秒、58-62℃30秒、及び72℃30秒）を27-35サイクル、
72℃5分。
GAPDHのRT-PCRのためのプライマー配列は以下である:
配列番号43   (forward) 5'-GGA AGG TGA AGG TCG GAG T-3' 
配列番号44   (reverse) 5'-TGG GTG GAA TCA TAT TGG AA -3'
【００７６】
（半定量PCRのプライマー（無病グループで高発現する遺伝子））
【００７７】
【表４Ａ】
2005270093.tif 000005
【００７８】
（半定量PCRのプライマー（再発グループで高発現する遺伝子））   
【００７９】
【表４Ｂ】
2005270093.tif 000006
【００８０】
（RT-PCR産物のシグナル強度測定と予後スコアの計算）
RT-PCR産物のシグナル強度を実施例１で記述した手法と同様に測定・評価し、t-テストで0.05以下のp値を示した10遺伝子を候補として選んだ；そのうち3遺伝子の発現レベルは、5y-Fグループに対して5y-Rグループで高かった。そして、7遺伝子の発現レベルは、5y-Rグループに対して5y-Fグループでより高かった。このインフォメーションを使って、node陰性乳癌の術後予後を予測するスコアリングシステムの確立を試みた。
【００８１】
各遺伝子の対象とする発現レベルを得るために、GAPDH発現に対するその発現比率（ER）を次の式によって計算した；
遺伝子AのER = 癌サンプルXの遺伝子Aの半定量PCR(エチジウムブロマイド染色したバンドの強度) の16-bit のイメージングスコア/癌サンプルXの遺伝子AのGAPDHの16-bitのイメージングスコア
【００８２】
（node陰性乳癌の術後予後を予測するスコアリングシステムの定義）
node陰性乳癌の術後の遺伝子予後指標を構築するため、予後スコア（PS）を定義した；（5Y-Fグループに比べて5Y-Rグループで高発現する遺伝子の正規化した発現比率の合計）−（5Y-Rグループに比べて5Y-Fグループで高発現する遺伝子の正規化した発現比率の合計）
【００８３】
2つのグループ間の発現比の有意さはStudent's t-テストで評価した。全ての統計手法は、Statview version 5.0（SASInstitute、Cary、NC）によった。
【００８４】
（結果）
ゲノムワイドな遺伝子発現を検討した24の乳がん患者の臨床病理所見を表３に要約した。本発明者らは、手術後5年以上の間無病生存した12例のnode陰性乳癌患者（5Y-F）、及び外科手術の後に5年の内に乳癌が再発した12例のnode陰性乳癌患者（5Y-R）からの腫瘍について、25,344のヒト遺伝子から構成されるcDNAマイクロアレイによる遺伝子発現を検討した。臨床のバックグラウンドは、2つのグループの間で、年齢、腫瘍径、エストロゲン受容体とプロゲステロン受容体、及び病理学的に合致させた。
【００８５】
cDNA マイクロアレイ のデータを、Mann-Whitney テストとRandom-permutationテストにより分析し、5Y-R及び5Y-F腫瘍間で異なる発現を示す遺伝子を同定した。このフィルターで全58遺伝子を選び、そのうち21遺伝子が5Y-R腫瘍で有意に強く発現した。そして37遺伝子が5Y-F腫瘍で高い発現レベルを示した。
【００８６】
5Y-R 腫瘍に比べて5Y-F腫瘍で高発現する37遺伝子には、6つのESTs と1つの仮想たんぱく質があった（表５Ａ、各グループ間での発現の差異を“foldchange”として示す）。
【００８７】
（5Y-R腫瘍に比べて5Y-F腫瘍で有意に高発現する遺伝子）
【００８８】
【表５Ａ】
2005270093.tif 000007
【００８９】
表５Ｂ は5Y-Rグループ において高発現する21遺伝子をリストにした。その5つがESTs で、一つが仮想たんぱく質をコードする。58遺伝子のこのパネルから、次のような基準で術後予後のマーカーを選んだ；（１）症例の少なくとも60％に位置するカットオフレベルより高いシグナル強度をもつ；（２）|μR-μF|>1.0、ここでμR(μF)は、5Y-R又は5Y-F症例における対数変換発現比率から導かれる平均値示す。
【００９０】
（5Y-F腫瘍に比べて5Y-R腫瘍で有意に高発現する遺伝子）
【００９１】
【表５Ｂ】
2005270093.tif 000008
【００９２】
5Y-Rに比べ5Y-F腫瘍で高発現する7遺伝子 (Hs.94653、M13436、Hs.5002、D67025、M80469、Hs.4864、及びHs.106326；p=0.0018、0.0011、0.001、0.008、0.0081、0.0018及び0.001；各Student'st-テストによる)及び5Y-R腫瘍で比較的高発現する3遺伝子（AF058701、AI066764、及びx15940；p=0.0351、0.00161及び0.0001；各Student'st-テストによる）が基準に合致し、予後マーカーとして選択した（表６）。
【００９３】
（node陰性乳癌の予後マーカーとして選別した遺伝子）
【００９４】
【表６】
2005270093.tif 000009
【００９５】
GAPDH発現に対する正規化後半定量RT-PCR実験でこれらのマーカーの発現を確認した。図３は、乳癌再発（5Y-Rグループ）の12患者からのサンプルで高発現する3つのマーカー遺伝子のRT-PCRの結果を示す。図４は、5Y-F（5年生存）グループで高発現する7つのマーカー遺伝子の結果を示す。これら10遺伝子の発現比率を予後指標の定義に用いた。
【００９６】
予後スコア (PS)を次のように定義した；
PS  =  （5Y-R腫瘍で高発現する3遺伝子の正規化した発現比率の合計）− （5Y-F腫瘍で高発現する7遺伝子の正規化した発現比率の合計）
【００９７】
検討した24ケースの予後スコアを、各マーカー遺伝子の発現比率と共に、表７に要約した。PSシステムは、3より高い予後スコアを有する症例R1からR12については、不良予後を予測した。一方、-16より低いスコアの症例F1からF12については、良好予後を予測した。この予測は、実際のこれらの臨床結果と、100％の正確性で一致した（図５）。5Y-Rグループの平均PSは9.44、そして5Y-Fグループの平均PSは-28.92であった。
【００９８】
（node陰性乳癌再発の予後スコア）
【００９９】
【表７】
2005270093.tif 000010
【実施例３】
【０１００】
原発性乳癌における術後予後予測の遺伝子発現機能の評価
【０１０１】
（組織サンプル）
組織サンプルは実施例１記載の手法と同様に採取した。1995-1997年の期間において、乳癌のための手術を受けて、5年以上の間または死亡時まで臨床的に追跡調査された954人の患者の中から、手術後5年以内に死亡した10例と手術後5年以上の間無病生存した10例を標本として選んだ。2つの患者グループの臨床背景は年齢、リンパ節への転移、腫瘍径と組織型に関して、可能な限り厳密に一致させた（表８）。最終的な予知システムをテストするのに用いた追加20ケースの臨床背景を表９に要約した。
【０１０２】
（マイクロアレイ分析に用いた患者の臨床的プロフィール）
【０１０３】
【表８】
2005270093.tif 000011
【０１０４】
（RT-PCR分析に用いた患者の臨床的プロフィール）
【０１０５】
【表９】
2005270093.tif 000012
【０１０６】
（臨床病理（Clinicopathological）パラメータ）
実施例１で記術した手法で臨床病理パラメータを調べた。
【０１０７】
（cDNAマイクロアレイの調製）
実施例２で記述した手法でcDNAマイクロアレイの調製を行った。
【０１０８】
（RNA抽出とRNA増幅）
RNAをTRIzol（Invitrogen、Carlsbad、CA、USA）で抽出した。変性RNAを除去するため、各々の抽出されたRNA（1μｇ）を、3.0%ホルムアルデヒド変性ゲル上で電気泳動した。DNA混入を除去するため、RNeasyキット（QIAGEN、Valencia、CA）を用いて精製した。MessageAmpaRNAキット（Ambion、Austin、TX）によるT7ＲＮＡポリメラーゼベースの増幅を行い、マイクロアレイ分析に用いるRNAを調製した。最初の増幅では、RNA（5μg）を鋳型として用いた。その後、最初に増幅されたRNA（aRNA）（2μg）を、二回目の増幅のための鋳型とした。増幅されたaRNAsはRNeasy精製キットで精製し、各aRNAの量を分光光度計により測定した。
【０１０９】
（aRNAの標識、ハイブリダイゼーション及びデータ分析）
Amino Allyl cDNA 標識キット（Ambion、Austin、TX）により、二回目の増幅による蛍光プローブ作成用aRNA（5μｇ）を用いて、ハイブリダイゼーションプローブを作成した。癌RNA及び正常なコントロールRNAに由来するプローブを、Cy5またはCy3のMono-Reactive'Dye（AmershamBioscienceUK Limited、Buckinghamshire、UK）で各々標識した。
【０１１０】
非結合色素を除去するために、標識プローブを、QIA quick PCR精製キット（QIAGEN、Valencia、CA）で精製した。腫瘍及び正常RNAからの蛍光標識プローブの各々10pmolを、4xマイクロアレイハイブリダイゼーション・バッファー（Amersham（UK））及び脱イオン化したホルムアミドと混合した。プローブ混合物を、40℃で15時間cDNAアレイにハイブリダイズさせた。その後、5分間で1回、その後10分間で2回にわたり0.2%SDSを含む0.1xSSCで洗浄した。全ての手順は、AutomatedSlideProcessorSystem（Amersham）で実行した。それぞれのハイブリダイゼーションのシグナル強度は、Gene Pix 4000（Amersham）で読み取り、GenePixPro3.0（Axon Instruments, Inc.、Foster City、CA、USA）によって評価した。読み取ったシグナルは、totalgenenormalizationmethod（Yang,Y.H.,Dudoit,S., Luu, P., Lin, D.M., Peng, V.,Ngai,J., and Speed,T.P. (2002). Nucleic Acids Res 30, e15.；Manos, E.J., andJones,D.A.(2001).Cancer Res 61, 433-438. ）によって正規化した。
【０１１１】
生存者と死亡者のグループ間で異なる発現を示す遺伝子を確認するために、正規化したシグナルを、Mann-Whitneyテストで分析した；正規化したシグナルを一連のXに適用した。Xとはそれぞれの遺伝子及びサンプルのためのCy5/Cy3シグナル強度比率である（Ono,K.,etal.(2000).Cancer Res 60, 5007-5011.）。Mann-Whitneyテストで0のＵ値を示した遺伝子と2つのグループ間で発現強度が2.0倍以上の違いを示したものを選んだ。S/N比が3.0未満の遺伝子は、検討から除外した。
【０１１２】
（半定量RT-PCR及び遺伝子発現比率）
マイクロアレイのデータを検証するため、本発明者らはRNA（10μg）を逆転写することにより、半定量RT-PCR実験をした。転写されたcDNAの濃度を調整するため、GAPDHを内部コントロールとして選び、半定量RT-PCRを実行した（Ono,K.,etal.(2000).CancerRes 60,5007-5011.）。GAPDHのプライマーは、5'-ggaaggtgaaggtcggagt-3（Foward）、及び5-tgggtggaatcatattggaa-3（Reverse）であった。プライマーの濃度を調整した後に、半定量RT-PCRを生存者と死亡者のグループからのサンプルで、選別した遺伝子について実行した。当該遺伝子（表１０）の各々のためのプライマーは、NCBIGenBank（http://www.ncbi.nlm.nih.gov/）のシーケンス情報とウェブサイト（http://www-genome.wi.mit.edu/cgi-bin/primer/primer3_www.cgi）上のプライマー3に基づいて設計した。各半定量PCR実験は、テンプレートとして濃度を調整したcDNA（1μl）、5UのTakaraEXTaq（Takara、大津、日本）、1xPCRバッファー（10mMのTris-HCl、50mMのKCl、1.5mMのMgCl2）、10nMのdNTPsと10pmolの前方及び後方プライマーで30μl総量にて行った。
配列番号160  ggaaggtgaaggtcggagt
配列番号161  tgggtggaatcatattggaa
【０１１３】
（半定量ＰＣＲのプライマー）
【０１１４】
【表１０】
2005270093.tif 000013
【０１１５】
生存者と死亡者のグループの間で遺伝子発現の強度を評価するために、各半定量PCR産物（8μl）を、2.5%アガロースゲル上で電気泳動し、エチジウムブロマイドで染色した。各々の染色されたサンプルの強度は、バックグラウンド校正を用いてAlphaImager3300（Alpha Inonotech、San Leandro、CA）により測定した。各遺伝子の発現レベルを得るために、その発現比率を、GAPDHの発現レベルで正規化した。
【０１１６】
発現比率は、以下の公式によって定義された：遺伝子Aの発現比率  =  癌サンプルX 中の遺伝子Aの半定量PCRの16ビットのイメージングスコア（エチジウムブロマイドで染色したバンドの強度）/癌サンプルX中のGAPDHの16ビットのイメージングスコア
【０１１７】
（原発性乳癌の予後指標（PI）の定義）
本発明者らは、5Ｓグループにおいて高発現する遺伝子の正規化した発現比率の合計から、5Dのグループにおいて高発現する遺伝子の正規化した発現比率の合計を差し引くことにより、原発性乳癌の予後指標（PI）を定義した。2つのグループ間の発現比率の有意性はStudent'st-テストで評価した。5S及び5Dのグループ間のPIの比較はMann-Whitneyテストにより行った。全ての統計はStatviewversion5.0を使って保管した（SASInstitute Inc.、Cary、NC）。
【０１１８】
（結果）
18,432のヒト遺伝子からなるcDNAマイクロアレイ上で、8人の乳癌患者からの腫瘍のゲノムワイドな遺伝子発現機能を調べた。患者のうちの4人は手術後5年以上の間無病生存し（5S）、4人は5年以内で、乳癌で死亡した（5D）。臨床病理学的背景は、2つのグループ間で年齢、腫瘍径、リンパ節転移、ホルモンレセプター状態と組織型について可能な限り厳密に一致させた（表８）。
【０１１９】
5Dと5Sのグループ間で異なる発現を示す遺伝子を同定するために、本発明者らはMann-Whitney テストによってcDNA マイクロアレイのデータを分析した。これらの遺伝子のうちの6つのESTs/仮想タンパク質である合計23遺伝子は、Mann-Whitneyテストで0のU値を示し、5Sグループにおいて高発現する遺伝子である（表１１）。
【０１２０】
（マイクロアレイ分析の解析による生存グループで高発現する遺伝子群）
【０１２１】
【表１１】
2005270093.tif 000014
【０１２２】
表１２は、5Dの腫瘍において一般に高発現しており、6つのESTs/仮想タンパク質を含み、Mann-Whitney テストでゼロのU値を持つ21遺伝子を記載する。表において、“foldchange”として、2つのグループ間の遺伝子発現の違いを示す。
【０１２３】
（マイクロアレイ分析の解析による死亡グループで高発現する遺伝子群）
【０１２４】
【表１２】
2005270093.tif 000015
【０１２５】
5Sグループにおいて高発現する23遺伝子及び5Dグループにおいて高発現する21遺伝子から、以下の基準に従って術後予後の予測マーカーを選んだ；（１）マイクロアレイ分析では、全ての症例の5S及び5D間のシグナル強度の違いが2.0倍より大きい；（２）半定量PCRの5S及び5D間のシグナル強度が有意に異なる（Student'st-テストによるp値<0.05）；（３）半定量PCRの結果は、独立した3回の実験によって再確認した。5S腫瘍において高発現する7遺伝子、及び5D腫瘍において高発現する3遺伝子は予後マーカーを選択するこれらの基準を満たした。
【０１２６】
5Sグループにおいて高発現する7遺伝子は、pro-alpha-1 type 3 collagen (PIIIP)、complement  component  C1r、dihydropyrimidinase-like  3（DPYSL3）、proteintyrosinekinase9-like（PTK9L）、carboxypeptidase E（CPE）、α-tubulinとβ-tubulinをコードしている遺伝子から構成されている。これらのマーカー遺伝子のStudent'st-テストのp値は、それぞれ0.00039、0.0012、0.0042、0.036、0.039、0.034と0.00069であった。
【０１２７】
5Dグループにおいて高発現する3つのマーカー遺伝子は、heat shock protein HSP 90-alpha gene、malatedehydrogenase、及びNADHdehydrogenase(ubiquinone) 1  beta  subcomplex, 3 (NDUFB3)をコードしていた。当該遺伝子のStudent'st-テストのp値は、それぞれ0.05、0.0055及び0.011であった。
【０１２８】
本発明者らは、半定量RT-PCRの実験結果を、内部コントロールとしてGAPDHで正規化し評価することによりマーカー遺伝子の選択を検証した。
【０１２９】
本発明者らは無作為に選んだ追加20症例を調べるために半定量PCRを行った。これらの患者のうちの10人は手術後5年以内に乳癌で死亡し、そして、残りの10人は5年以上の間無病生存した。図７は、5S腫瘍において高発現する7つのマーカー遺伝子のRT-PCRの結果を示す。図８は、5D腫瘍において高発現する3つのマーカー遺伝子のRT-PCRの結果示す。
【０１３０】
本発明者らは、次のように予後指標（PI）を定義した：（5Sグループにおいて高発現する遺伝子の正規化した発現比率の合計）−（5Dグループにおいて高発現する遺伝子の正規化した発現比率の合計）。選ばれたマーカー遺伝子の発現比率と共に、更なる試験例のための予後指標を、表１３にまとめた。
【０１３１】
（遺伝子の発現比率と予後指標）
【０１３２】
【表１３】
2005270093.tif 000016
【０１３３】
PIは、5Sグループの全10症例（S1からS10）の高い予後指標（>7）及び5Ｄグループの全10症例（D1からD10）の予後指標（<7）（図9）の実際の臨床結果を正確に予測した。5SグループのPIは21.2であった、そして、5DグループのPIは-0.7であった。ここでPI値7は、明らかに5S腫瘍と5D腫瘍を区別していた（p=0.0002）。
【産業上の利用可能性】
【０１３４】
本発明の術後予後予測システムは、乳癌患者の術後リスクの予測に有効である。さらに、本発明の乳癌関連遺伝子の広範囲に及び遺伝子発現リストは乳癌の進行についての様々な情報を提供するとともに、乳癌治療の潜在的ターゲット分子の予測が可能となる。
【図面の簡単な説明】
【０１３５】
【図１】5y-Sグループに比べて5y-Dグループで発現が上昇した遺伝子群（A）と減少した遺伝子群（B）を示す。
【図２】5y-Sグループと5y-Dグループ由来のRNAの半定量RT-PCRの分析結果を示す。
【図３】個々の患者における予後スコア
【図４】5Y-Fグループと5Y-Rグループ由来のRNAについての半定量RT-PCRの分析結果を示す。
【図５】5Y-Fグループと5Y-Rグループ由来のRNAについての半定量RT-PCRの分析結果を示す。
【図６】個々の患者の予後スコアを示す。
【図７】5S腫瘍において高発現する７遺伝子の半定量PCRの分析結果を示す。M：マーカーラダー（Marker ladder）S1−S10：手術後5年以上の間無病生存した患者の新たに検査した組織である。D1−D10：手術後5年以内に乳癌で死亡した患者を新たに検査したケースである。発現強度の違いはStudent's t-テストで評価した；p値が0.05以下の場合、統計的に有意であると考えられる。
【図８】5Dグループで高発現する3遺伝子の半定量PCRの分析結果を示す。符号の説明は図７の説明を参照。
【図９】新たに調べた20ケースの予後指標（PI）を描図した結果である。5年以上の間無病生存した全10人の患者の指標は7より高かった。一方、手術後5年以内に乳癌で死亡した患者の指標は7より低かった。2つのグループの分配は統計学的に有意である（ｐ＝0.0002）
【図３】
2005270093.tif 000017
【図４】
2005270093.tif 000018
【図５】
2005270093.tif 000019
【図６】
2005270093.tif 000020
【図７】
2005270093.tif 000021
【図８】
2005270093.tif 000022
【図９】
2005270093.tif 000023
【図１】
2005270093.tif 000024
【図２】
2005270093.tif 000025
【配列表】
2005270093000001.app