(57)【特許請求の範囲】
【請求項１】
Ａｍｅｓ試験陽性の少なくとも１の発がん物質（Ｃ＋Ｍ＋）及びＡｍｅｓ試験陰性の少なくとも１の発がん物質（Ｃ＋Ｍ−)を試験動物（ヒトを除く）に投与し、所定時間経過後に試験動物からｍＲＮＡを採取してその発現パターンをそれぞれ予め準備しておき、発がん性未知の被検物質の発がん性を予測するに際しては、
（１）発がん性を予測する被検物質にＡｍｅｓ試験を行い、該被検物質がＡｍｅｓ試験陽性又は陰性の何れに属するかを決定する第１工程、
（２）発がん性を予測する被検物質を試験動物に投与して所定期間経過後、試験動物からｍＲＮＡを採取してその発現パターンを得る第２工程、
（３）第１工程で決定される被検物質のＡｍｅｓ属性に応じて、被検物質が前記第１工程でＡｍｅｓ試験陽性に属する場合は、予め準備したＡｍｅｓ試験陽性の発がん物質（Ｃ＋Ｍ＋）のｍＲＮＡの発現パターンと前記第２工程で得られたｍＲＮＡの発現パターンとの一致度を算出し、被検物質が前記第１工程でＡｍｅｓ試験陰性に属する場合は、予め準備したＡｍｅｓ試験陰性の発がん物質（Ｃ＋Ｍ−）のｍＲＮＡの発現パターンと前記第２工程で得られたｍＲＮＡの発現パターンとの一致度を算出し、算出した一致度から被検物質の発がん性を予測する第３工程、を行い、
ｍＲＮＡの発現パターンが、Ａｍｅｓ試験陽性の発がん物質（Ｃ＋Ｍ＋）を投与した試験動物とＡｍｅｓ試験陽性の非発がん物質（Ｃ−Ｍ＋）を投与した試験動物の間又はＡｍｅｓ試験陰性の発がん物質（Ｃ＋Ｍ−）を投与した試験動物とＡｍｅｓ試験陰性の非発がん物質（Ｃ−Ｍ−）を投与した試験動物の間で発現量の差が有意なｍＲＮＡを用いて得られる、配列番号１〜１４３５又は配列番号１４３６〜２３５７に記載する塩基配列の遺伝子から選定された３つ以上の遺伝子を搭載したＤＮＡマイクロアレイに、前記第２工程で試験動物から採取したｍＲＮＡから調製した蛍光標識ｃＤＮＡ又はｃＲＮＡをハイブリダイゼーションさせて得られるＤＮＡマイクロアレイの蛍光パターンを使用するものである被検物質の発がん性予測方法。
【請求項２】
ｍＲＮＡの発現パターンが、Ａｍｅｓ試験陽性の発がん物質（Ｃ＋Ｍ＋）を投与した試験動物とＡｍｅｓ試験陽性の非発がん物質（Ｃ−Ｍ＋）を投与した試験動物の間又はＡｍｅｓ試験陰性の発がん物質（Ｃ＋Ｍ−）を投与した試験動物とＡｍｅｓ試験陰性の非発がん物質（Ｃ−Ｍ−）を投与した試験動物の間で発現量の差が有意なｍＲＮＡを逆転写し、前記逆転写反応により合成したｃＤＮＡ配列の全部又は一部を、ＰＣＲにより、少なくともフォワードプライマー又はリバースプライマーのいずれかを用いて増幅して得たＤＮＡを搭載したＤＮＡマイクロアレイ、又は、前記逆転写反応により合成したｃＤＮＡ配列の一部を化学合成して得たＤＮＡを搭載したオリゴマイクロアレイに、前記第２工程で試験動物から採取したｍＲＮＡから調製した蛍光標識ｃＤＮＡ又はｃＲＮＡをハイブリダイゼーションさせて得られるＤＮＡマイクロアレイ、又はオリゴマイクロアレイの蛍光パターンを使用するものである、請求項１に記載の被検物質の発がん性予測方法。
【請求項３】
Ａｍｅｓ試験陽性の発がん物質（Ｃ＋Ｍ＋）又はＡｍｅｓ試験陰性の発がん物質（Ｃ＋Ｍ−)を複数試験動物に投与し、Ａｍｅｓ試験陽性の発がん物質（Ｃ＋Ｍ＋）又はＡｍｅｓ試験陰性の発がん物質（Ｃ＋Ｍ−）について得られた複数のｍＲＮＡの発現パターンをその類似性により２以上のグループに分類し、グループ毎にｍＲＮＡの発現パターンを準備する請求項１又は請求項２に記載の被検物質の発がん性予測方法。
【請求項４】
前記第２工程において、被検物質を１〜９０日間試験動物に連続投与した後、試験動物からｍＲＮＡを採取する請求項１乃至３のいずれかに記載の被検物質の発がん性予測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検物質の投与により試験動物に発現するmRNAの発現パターンを、発がん性が既知の化学物質の投与により発現するmRNAの発現パターンと比較することにより、被検物質の発がん性を予測する発がん性予測方法に関する。
【背景技術】
【０００２】
化学物質の有害性（ハザード）評価においては、発がん性など長期毒性を評価するには多額の費用と長期の試験期間を要する動物実験の実施が必要とされる。また、動物実験の結果をヒトへ外挿する際は、種差や作用機序などを考慮しなければならないなど様々な課題がある。
【０００３】
近年ではゲノム情報に関する技術の著しい発展により、遺伝子レベルで化学物質の有害性評価が行われるようになっている。例えば、特許文献１には、化学物質を曝露した組織や細胞の遺伝子の発現の差異を検出することにより化学物質の毒性作用を予測する方法が記載されている。
【０００４】
化学物質の発がん性の評価においても発がんメカニズムに関与している遺伝子の存在が予想されることからこれらの遺伝子の発現の差異を検出することにより遺伝子レベルでの評価が可能であると考えられる。しかしながら、現段階では、化学物質が引き起こす遺伝子レベルでの発がんメカニズムはほとんど解明されておらず、遺伝子の発現の差異から化学物質の発がん性予測を行うことは非常に困難である。
【０００５】
最近では、遺伝子の発現プロファイルデータを網羅的に収集し解析する手法としてＤＮＡマイクロアレイが開発され、様々な分野で使用されるようになっている。ＤＮＡマイクロアレイ（ＤＮＡチップ）はスライドガラス等の基板上に固定化した数百から数万種類の微量のＤＮＡに、核酸分子同士の相補性を利用してｍＲＮＡから作成したｃＤＮＡ或いはｃＲＮＡを結合させ、ｍＲＮＡの発現量に関する情報を得る技術である。
【０００６】
既知の遺伝子を搭載したＤＮＡマイクロアレイが既に商品化され一般に使用できる環境にあるが、市販のアレイは化学物質の毒性評価を目的に設計・製作されたものではない。また、通常の毒性試験規模で使用するにはデータの解析を含めたシステム全体を構築するコストが多大であることから、毒性試験分野で広く活用されていないのが現状である。
【０００７】
発がん性物質は遺伝毒性あるいは突然変異誘起性を有していると考えられおり、従来、化学物質のがん原性の第一次スクリーニング試験には復帰変異試験であるAmes試験が広く使用されている。
【特許文献１】特開２００３−３０４８８８号公報（特許請求の範囲）
【発明の開示】
【発明が解決しようとする課題】
【０００８】
一般に、化学物質の発がん性を評価する動物実験では、試験動物に癌が出現するまで、あるいは試験動物が死亡するまで化学物質の連続投与を行うが、癌は長い潜伏期間を経て出現するため長期にわたる動物実験が必要となる。
【０００９】
しかしながら、癌が出現する前段階においても、発がん物質の刺激により遺伝子レベルでは何らかの変化が生じていることが予想される。
【００１０】
本発明の目的は、がん発生の初期段階に発がんメカニズムに関与している可能性が高い遺伝子群について発現の差異を検出してデータ解析することにより、長期にわたる動物実験を行うことなく被検物質の発がん性を予測する方法を提供することにある。
【課題を解決するための手段】
【００１１】
本発明者らは発がん物質の投与により発現量が変動するmRNAに着目し、その変動から被検物質の発がん性を予測することを試みた。当初は、これらmRNAの変動パターンが全ての発がん物質で類似性を持つものと考えていた。しかし、実際には異なる挙動を示し、１つの類似性を共有することはなかった。
【００１２】
本発明者らは、化学物質の復帰変異試験として広く採用されているAmes試験において陰性を示す発がん物質が存在することに着目し、Ames試験陽性の発がん物質（C+M+）とAmes試験陰性の発がん物質（C+M-）との間で発がんに到るメカニズムが異なることに想到した。そこで、被検物質を試験動物に投与したときのmRNAの発現パターンを、被検物質のAmes属性に応じてAmes試験陽性の発がん物質（C+M+）の発現パターン又はAmes試験陰性の発がん物質（C+M-）の発現パターンと比較することにより、動物実験の初期段階においても被検物質の発がん性が高い精度で予測できることを見出し本発明を完成するに到った。
【００１３】
即ち、上記課題を解決する本発明は以下に記載するものである。
【００１４】
〔１〕  Ames試験陽性の少なくとも１の発がん物質（C+M+）及びAmes試験陰性の少なくとも１の発がん物質（C+M-)を試験動物に投与し、所定時間経過後に試験動物からmRNAを採取してその発現パターンをそれぞれ予め準備しておき、発がん性未知の被検物質の発がん性を予測するに際しては、
（１）発がん性を予測する被検物質にAmes試験を行い、該被検物質がAmes試験陽性又は陰性の何れに属するかを決定する第１工程、
（２）発がん性を予測する被検物質を試験動物に投与して所定期間経過後、試験動物からmRNAを採取してその発現パターンを得る第２工程、
（３）被検物質が前記第１工程でAmes試験陽性に属する場合は、予め準備したAmes試験陽性の発がん物質（C+M+）のmRNAの発現パターンと前記第２工程で得られたmRNAの発現パターンとの一致度を算出し、被検物質が前記第１工程でAmes試験陰性に属する場合は、予め準備したAmes試験陰性の発がん物質（C+M-）のmRNAの発現パターンと前記第２工程で得られたmRNAの発現パターンとの一致度を算出し、算出した一致度から被検物質の発がん性を予測する第３工程、
を行うことを特徴とする被検物質の発がん性予測方法。
【００１５】
〔２〕  Ames試験陽性の発がん物質（C+M+）又はAmes試験陰性の発がん物質（C+M-)を複数試験動物に投与し、Ames試験陽性の発がん物質（C+M+）又はAmes試験陰性の発がん物質（C+M-)について得られた複数のmRNAの発現パターンをその類似性により２以上のグループに分類し、グループ毎にmRNAの発現パターンを準備する〔１〕に記載の被検物質の発がん性予測方法。
【００１６】
〔３〕  mRNAの発現パターンが、Ames試験陽性の発がん物質（C+M+）を投与した試験動物とAmes試験陽性の非発がん物質（C-M+）を投与した試験動物の間又はAmes試験陰性の発がん物質（C+M-）を投与した試験動物とAmes試験陰性の非発がん物質（C-M-）を投与した試験動物の間で発現量の差が有意なmRNAを逆転写したcDNA又はcRNAを搭載したDNAマイクロアレイに、前記第２工程で試験動物から採取したmRNAから調製した蛍光標識cDNA又はcRNAをハイブリダイゼーションさせて得られるDNAマイクロアレイの蛍光パターンを使用するものである〔１〕に記載の発がん性予測方法。
【００１７】
〔４〕  DNAマイクロアレイが、Ames試験陽性の発がん物質（C+M+）を投与した試験動物とAmes試験陽性の非発がん物質（C-M+）を投与した試験動物の間及びAmes試験陰性の発がん物質（C+M-）を投与した試験動物とAmes試験陰性の非発がん物質（C-M-）を投与した試験動物の間で発現量の差が有意なmRNAから逆転写反応により合成したcDNA配列の全部又は一部をPCRにより増幅又は化学合成したDNAを搭載したものである〔３〕に記載の被検物質の発がん性予測方法。
【００１８】
〔５〕  DNAマイクロアレイが、配列番号１〜１４３５又は配列番号１４３６〜２３５７に記載する塩基配列から選定された３つ以上を搭載したものである〔３〕又は〔４〕に記載の被検物質の発がん性予測方法。
【００１９】
〔６〕前記第２工程において、被検物質を１〜９０日間試験動物に連続投与した後、試験動物からmRNAを採取する〔１〕乃至〔５〕のいずれかに記載の被検物質の発がん性予測方法。
【発明の効果】
【００２０】
本発明においては、被検物質のAmes試験における属性によりAmes試験陽性の発がん物質（C+M+)又はAmes試験陰性の発がん物質のmRNAの発現パターンと、被検物質のmRNAの発現パターンと比較を行うので、従来予測が難しかったmRNAの発現量からの被検物質の発がん性予測が可能となる。
【００２１】
本発明においては、mRNAの発現量から被検物質の発がん性を予測するので、試験動物に癌が出現するまで化合物の連続投与を行うような長期にわたる動物実験を必要とせず、１〜９０日程度の動物実験で容易に化学物質の発がん性を予測することができる。
【発明を実施するための最良の形態】
【００２２】
被検物質の発がん性を予測するに際しては、予め以下の手順に従ってAmes試験陽性の発がん物質（C+M+）とAmes試験陰性の発がん物質（C+M-)の発現パターンを準備する。
【００２３】
まず最初に、Ames試験陽性の少なくとも１の発がん物質（C+M+）及びAmes試験陰性少なくとも１の発がん物質（C+M-)を各々試験動物に投与して所定時間経過後に試験動物の組織からmRNAを採取する。
【００２４】
Ames試験陽性の発がん物質（C+M+）又はAmes試験陰性の発がん物質（C+M-)を投与後、試験動物の組織を採取するまでの期間としては、１〜９０日程度とするが、より迅速に試験する観点から１〜２８日程度とすることが好ましい。
【００２５】
発がん物質の試験動物への投与は、１日１〜数回（好ましくは１日１回）、１〜９０日程度の所定期間反復して連続投与する。
【００２６】
発がん物質の投与方法は特に制限されず、経口投与、腹腔内投与、静脈内投与等の汎用的な方法を使用できる。
【００２７】
試験動物には、ラット、マウス、イヌ、サル、モルモット、ウサギ等が使用できる。また、mRNAの発現量を測定するために採取する試験動物の組織としては、肝臓、腸、肺、腎臓、胃、脾臓、脳、血液等が使用できる。
【００２８】
試験動物に投与するAmes試験陽性の発がん物質（C+M+）、Ames試験陰性の発がん物質（C+M-)は少なくとも１以上とするが、予測の精度を高めるため２物質以上とすることが好ましく、３物質以上とすることがより好ましい。使用する発がん物質の数は多いほど被検物質の発がん性予測の精度が高くなるので投与する物質の数に上限はないが、いずれの発がん物質も３０物質程度を投与すれば充分な精度で被検物質の発がん性が予測できる。
【００２９】
試験動物に投与する発がん物質は、Ames試験の属性が既知の公知の発がん物質を使用することが可能である。
【００３０】
Ames試験陽性の発がん物質（C+M+）としては、例えば、2,4-ジアミノトルエン、キノリン、ジエチルニトロサミン、2-ニトロプロパン、N-ニトロソモルホリン、サフロール、N-ニトロソジメチラミン、N-ニトロソピペリジン、2-アミノ-3,8-ジメチルイミダゾ[4,5-f]キノキサリン、2-アミノ-1-メチル-6-フェニルイミダゾ[4,5-b]-ピリジン、ベンツ(a)アントラセン、7,12-ジメチルベンズアントラセン、3-メチルコラントレン、4-ニトロキノリン-1-オキサイド、N-エチル-N-ニトロソ尿素、タンニン酸、ベンゾ（a）ピレン、N-メチル-N'-ニトロ-N-ニトロソグアニジン、4-ジメチルアミノアゾベンゼン等を挙げることができる。
【００３１】
Ames試験陰性の発がん物質（C+M-）としては、クロフィブレート、ジ（2-エチルヘキシル）フタレート、四塩化炭素、フェノバルビタール、アルドリン、アジピン酸ジ(2-エチルヘキシル)、エチニルエストラジオール、ヘキサクロロベンゼン、α-ヘキサクロロシクロヘキサン、トリクロロエチレン、ブチル化ヒドロキシアニソール、リモネン、トリクロロ酢酸、ウレタン、ペンタクロルエタン、クロロホルム、テトラクロロエチレン、アセタミド、ジエチルスチルベストロール、フェニトインナトリウム塩、D,L-エチオニン等を挙げることができる。
【００３２】
試験動物の組織から公知の方法によりmRNAを抽出、精製した後、mRNAの発現量を測定し、Ames試験陽性の発がん物質（C+M+）、Ames試験陰性の発がん物質（C+M-）について、それぞれmRNAの発現パターンを得る。発現パターンは、Ames試験陽性の発がん物質（C+M+）を投与した試験動物とAmes試験陽性の非発がん物質（C-M+）を投与した試験動物の間、又はAmes試験陰性の発がん物質（C+M-）を投与した試験動物とAmes試験陰性の非発がん物質（C-M-）を投与した試験動物の間で、発現量が有意に増加又は減少するmRNAの種類とその発現量により構成されるデータの集まりである。
【００３３】
mRNAの発現パターンには、少なくとも１以上のmRNAを選択して使用するが、使用するmRNAの数は、好ましくは３以上、より好ましくは５以上である。
【００３４】
試験動物にラットを使用する場合、Ames試験陽性の発がん物質（C+M+）とAmes試験陽性の非発がん物質（C-M+）との間、又はAmes試験陰性の発がん物質（C+M-)とAmes試験陰性の非発がん物質（C-M-)との間で、mRNAの発現量が変動する傾向がある遺伝子の塩基配列の一部を１〜１４３５又は１４３６〜２３５７に示す。
【００３５】
mRNAの発現量の測定にDNAマイクロアレイを使用する場合には、多くのmRNAの発現量に関する情報を同時に取得できる。そのため、アレイの蛍光をスキャナで取り込んで取得した画像は、必要に応じてスポットの選択、画像処理等を行った後、そのまま発現パターンとして使用することができる。
【００３６】
DNAマイクロアレイには、試験動物から採取したmRNAを逆転写したcDNA又はDNAが固定してあるDNAマイクロアレイを使用する。mRNAの発現量の測定及び発現パターンの取得には、試験動物から採取したｍＲＮＡから蛍光標識ｃＤＮＡ又はｃＲＮＡを調製し、これらをアレイ上のcDNA又はDNAと結合させることにより行う。なお、マイクロアレイに代えてマイクロプレートを用いることも可能である。
【００３７】
DNAマイクロアレイに搭載するcDNA又はDNAとしては、試験動物の組織から抽出したmRNAから逆転写反応により合成したcDNAの配列の全部又は一部をPCR等により増幅したDNAのほかに、cDNAの配列の一部を化学合成して得たDNAも使用できる。
【００３８】
DNAマイクロアレイには、試験動物の種類に応じて選択した市販のアレイを使用することが可能である。例えば、試験動物にラットを使用する場合には、市販のアレイとして、Gene Chip（商品名、アフィメトリクス社製）、Rat Oligo MicroArray Kit（商品名、アジレント社製）等が例示できる。
【００３９】
本発明においては、DNAマイクロアレイを用いる方法等により各mRNAの発現量を一度に測定してmRNAの発現パターンを取得してもよいし、各mRNAの発現量を別々に測定することによりmRNAの発現パターンを取得してもよい。mRNAの発現量の測定方法には、例えばノーザンブロッティング、定量的RT-PCR、RNaseプロテクションアッセイ等の公知の方法を制限することなく使用できる。
【００４０】
次に、発がん性未知の被検物質の発がん性を予測する際の手順について説明する。被検物質の発がん性を予測するに際しては、まず最初に被検物質にAmes試験を行い、該被検物質がAmes試験陽性又は陰性の何れに属するかを決定する。
【００４１】
Ames試験の方法自体は公知であり、OECDの試験法ガイドラインNo.471等の試験法が定められている。Ames試験は、医薬品の製造（輸入）承認申請に必要な遺伝毒性試験等として実施されている。
【００４２】
被検物質は上述したAmes試験陽性の発がん物質（C+M+）又はAmes試験陰性の発がん物質（C+M-）と同様の方法を用いて試験動物へ投与し、所定期間経過後に試験動物からmRNAの採取を行う。Ames試験陽性の発がん物質（C+M+）又はAmes試験陰性の発がん物質（C+M-）に特徴的な変動を示すmRNAについてその発現量を測定し、被検物質についてmRNAの発現パターンを得る。
【００４３】
その後、被検物質がAmes試験陽性に属する場合は、予め準備してあるAmes試験陽性の発がん物質（C+M+）のmRNAの発現パターンと、被検物質がAmes試験陰性に属する場合は、Ames試験陰性の発がん物質（C+M-）のmRNAの発現パターンと、被検物質のmRNAの発現パターンとの一致度を統計学的手法により算出し、算出した一致度から被検物質の発がん性を予測する。被検物質のmRNAの発現パターンが、Ames試験陽性の発がん物質（C+M+）のmRNAの発現パターン又はAmes試験陰性の発がん物質（C+M-）のmRNAの発現パターンと一致度が高い場合には被検物質は発がん性有りと判定し、いずれの発現パターンとも一致度が低い場合には発がん性無しと判定する。
【００４４】
発現パターンの一致度を評価する方法としては、例えば被検物質と発がん物質とを一緒にクラスタ分析や決定木による解析を行って得られた樹形図の階層から評価する方法、mRNAの発現量の相関を算出する方法等を挙げることができるが、サポートベクターマシーン（Support Vector Machine :SVM）等の公知の分類方法を用いて予測式を作成する方法が最も簡便である。
【００４５】
なお、Ames試験陽性の発がん物質（C+M+）、又はAmes試験陰性の発がん物質（C+M-)のmRNAの発現パターンを準備する際に、これらの発がん物質を複数使用してmRNAの発現パターンを取得する場合には、各発がん物質について得られたmRNAの発現パターンをその類似性により複数のグループに分類し、各グループ毎に発現パターンを準備してもよい。
【００４６】
Ames試験陽性の発がん物質（C+M+）又はAmes試験陰性の発がん物質（C+M-)のmRNAの発現量から、mRNAの発現パターンをグループ化する方法につき以下説明する。
【００４７】
まず、Ames試験陽性の発がん物質（C+M+）又はAmes試験陰性の発がん物質（C+M-)について発がん物質の投与により有意に発現量が増加又は減少（変動）したmRNAを選択する。mRNAの発現量が有意に変動したかどうかの判断は、Ames試験陽性の発がん物質（C+M+）についてはAmes試験陽性の非発がん物質（C-M+）との間、Ames試験陰性の発がん物質（C+M-)についてはAmes試験陰性の非発がん物質（C-M-)との間でmRNAの発現量を比較することにより行う。
【００４８】
Ames試験陽性の非発がん物質（C-M+）としては、例えば2,6-ジアミノトルエン、8-ヒドロキシキノリン、2-クロロエタノール、2-(クロロメチル)ピリジン塩酸塩、4-ニトロ-o-フェニレンジアミン、2-クロロ-p-フェニレンジアミン硫酸塩、p-フェニレンジアミン二塩酸塩、2,5-トルエンジアミン硫酸塩、4-（クロロアセチル）アセトアニリド、3-クロロ-p-トルイジン、グルタルアルデヒド、4-ニトロアントラニル酸、１−ニトロナフタレン等を挙げることができる。
【００４９】
Ames試験陰性の非発がん物質（C-M-)としては、D-マンニトール、L-アスコルビン酸、1-ニトロプロパン、リンダン、アスピリン、フタルアミド、カプロラクタム、安息香酸ナトリウム、インドメタシン等を挙げることができる。
【００５０】
発現パターンの分類に使用するmRNAは、発現量が発がん物質と非発がん物質の間で有意に変動したmRNAの中から選択した少なくとも3以上とするが、4〜200程度とすることがより好ましく、5〜50程度とすることが特に好ましい。
【００５１】
選択したmRNAについて、発がん物質間で、クラスタ解析、決定木等の解析を行い、発現パターンを分類する。
【実施例】
【００５２】
製造例１（cDNAマイクロアレイの作製）
（１）完全長cDNA クローン及びライブラリーの作製（理化学研究所法）
発がん剤投与成獣（6週齢）ラット由来の肝臓、腎臓、脾臓及び大腸組織から、株式会社ダナフォームが作製したラット完全長cDNAライブラリーを利用して完全長cDNAクローンを得た。なお当該ラット完全長ｃDNAライブラリーの作製手順の概略は以下の通りであった。
【００５３】
Ｆ344 ラット成獣（6 週齢）に作用機構が異なる４種の既知発がん物質（Diethylnitrosoamine：投与量100mg/kg/day：投与期間3 時間, 1 日, 3 日、2-amino-3,8-dimethyl-imadazo[4,5-f]quinoxaline：投与量20mg/kg/day：投与期間 3 時間, 1 日, 4 日、Clofibrate：投与量100mg/kg/day：投与期間1 日, 2 日, 4 日、Phenobarbital：投与量80mg/kg/day：投与期間1 日, 2 日, 4 日）を投与し、肝臓、腎臓、脾臓及び大腸を採取した。採取した組織からAGPC 法によりtotal RNA を抽出し、MACS mRNA Isolation Kit (Miltenyi Biotec 社製)を用いてmRNA を精製した。その後、Anchored oligo(dT)をプライマーとしてトレハロース存在下で耐熱化した逆転写酵素により逆転写反応を行った。このAnchored oligo(dT)プライマーには、6 塩基のタグ配列を用いてライブラリー及びcDNA の起源に関する情報を付した。次にキャップトラッパー法を用いて完全長cDNA を選択し、第2 鎖合成、制限酵素切断を行った後、方向を定めてクローニングし、完全長cDNA ライブラリーを作製した。この完全長cDNAライブラリーは、化合物未処理ラット臓器由来のcDNA ライブラリーをドライバーとすることでサブトラクション処理を行い、特異性を高めた。
【００５４】
（２）cDNAマイクロアレイに搭載する遺伝子の選択
発がん物質投与成獣由来完全長ｃDNA ライブラリー中のクローンの3'端の塩基配列からクローン間の相同性を調べた。発がん物質投与成獣由来遺伝子の合計数は15,762 で、このうちの肝臓由来の遺伝子数は4,139 であった。また、UniGene ID を元にした非重複塩基配列数は6,954 であった。
【００５５】
未処理ラット臓器由来クローンと、発がん物質処理ラット臓器由来クローンを搭載する遺伝子の選定対象の遺伝子プールとして、以下の基準から8,862 クローンを選定した。
【００５６】
a)  発がん物質の投与によりコントロール動物と比較して2倍以上あるいは1/2 以下に発現量が変動した遺伝子。
【００５７】
b)  発がん物質処理ラット臓器由来クローン。
【００５８】
c)  文献等で発がん性に関連していることが報告されている遺伝子。これらの遺伝子の中にはがん遺伝子やがん抑制遺伝子、転写因子、増殖因子などが含まれる。なお、文献等で発がん性に関連していることが報告されている遺伝子で、上記の遺伝子プールに入っていないが発がん性評価には重要と考えられる24 種の遺伝子については、個別にクローンを購入あるいは単離した。
【００５９】
d)  市販の毒性評価用マイクロアレイに共通に搭載されている遺伝子。市販アレイとしては、Clontech Rat Toxicology1.2 array、TAKARA Rat Toxicology CHIP ver1.0、NIEHS Rat Chip、MWG Rat Liver、Mergen Rat R01 を調査した。
【００６０】
e)  ポジティブコントロールとして3 遺伝子(ベータ-アクチン、GAPDH、ユビキチン)及びネガティブコントロールとしてのラムダファージ遺伝子。
【００６１】
a)に含まれる遺伝子は2,890 クローン、b)に含まれるのは5,759 クローン、c)に含まれる遺伝子は44 クローン、d)に含まれる遺伝子は185 クローンであった。そして、上記の基準で選定した遺伝子について、PCRにより増幅した遺伝子をマイクロアレイに搭載した。
【００６２】
（３）cDNAマイクロアレイの作製
選定した遺伝子の塩基配列からお互いに非相補的な部分配列を設計し、（２） で選択した8,862 クローンをテンプレートとして、ポリメラーゼ連鎖反応（PCR）法により遺伝子断片を増幅した。増幅遺伝子の長さは、一部の遺伝子を除いて約300bpに設計した。cDNAマイクロアレイに搭載する遺伝子断片をPCR法で増幅するためのプライマーはタカラバイオ株式会社の設計アルゴリズムを用いて設計した。フォワードプライマー、リバースプライマーを各35pmol、Ex-Taq ポリメラーゼ2.5 unit 及び酵素に添付されたバッファーを用いて、反応液量を100μL としてPCR 反応を行った。反応は、95℃ 30 秒間の後、95℃ 45 秒間、54℃ 30 秒間、72℃ 60 秒間を1 サイクルとしてこれを37 サイクル行い、さらに72℃ 3 分間行った。反応終了後、反応液の一部をアガロース電気泳動に供することで、PCR 産物を確認した。
【００６３】
PCR 反応後、産物が増幅されたものは8,261クローンであった。さらに産物を詳しく検討した結果、PCR 産物が単一の遺伝子断片でないものやスメアーなものを除き、最終的な解析対象としては8,051 クローンとした。なお2003/01/06 版のUniGeneデータベースでアノテーションし直した結果、非重複なUniGene ID 数としては6,353であった。
【００６４】
得られたPCR 産物はQIAGEN 社製  QIAquick PCR Purification Kit で精製したのち、ガラス基板上にスポットし、cDNAマイクロアレイを得た。１グラス当たり2 アレイをスポットした。
【００６５】
実施例１
（１）total RNAの取得
表１〜４に示す化学物質をそれぞれ媒体に溶解し、溶液を調製した。化学物質とその物質番号、発がん性、使用した媒体、試験動物に投与した用量を表１〜４に示す。
【００６６】
【表１】
0004841279.tif 000002
【００６７】
【表２】
0004841279.tif 000003
【００６８】
【表３】
0004841279.tif 000004
【００６９】
【表４】
0004841279.tif 000005
【００７０】
日本チャールス・リバー社から入手した５週令の雄性ラット（Ｆ３４４、ＳＰＦ系統）を４匹／群に分け、各群のラットに調製した溶液又は媒体を強制経口投与した。投与は１日１回、１、３、７、１４、２８日間行った。投与開始から２、４、８、１５、又は２９日目に各個体の肝臓を採取した。肝臓から断片を切り出し、total RNA 保存安定用試薬であるRNAlater（登録商標） RNA Stabilization Reagent に浸漬させた。これを室温で24 時間放置後、RNAlater RNA Stabilization Reagent に浸漬させたまま-20℃で保存した。
【００７１】
ラット組織からtotal RNAを抽出し、QIAGEN 社BioRobot 3000 を用いてQIAGEN 社RNeasy Midi Kit により精製した。精製法は製造会社の使用法に準じて行ったが、概略は以下のとおりである。
【００７２】
切り出した組織サンプル50-100mg と破砕用ビーズ（直径5 mm、ジルコニア製）を2 mL エッペンドルフチューブに入れた。RLT バッファーを添加した後、Mixer Mill (QIAGEN 社製)により25 Hz 4 分間で2 回組織破砕を行った。5,000×g 5分間遠心し、未破砕物を沈殿させた。上清に70%エタノールを添加、混合した後、RNeasy Midi column へ添加した。3,000-5,000×g 5 分間遠心した後、RW1 バッファーを添加した。3,000-5,000×g 5 分間遠心した後、DNase I 溶液(QIAGEN 社製) を添加し室温で15 分間放置しさらにRW1 バッファーを添加した。3,000-5,000×g 5 分間遠心した後、RPE バッファーを添加した。3,000-5,000×g 2 分間遠心した後、再度RPE バッファーを添加した。3,000-5,000×g 2 分間遠心した後、3,000-5,000×g でさらに5 分間遠心しカラムを乾燥させた。RNase-free 水を添加し5 分間放置した後、3,000-5,000×g 2 分間遠心することでRNA 溶液を溶出させた。濃度及び回収率を向上させるために溶出したRNA 溶液を再度カラムに添加し5 分間放置した後、3,000-5,000×g 2 分間遠心することで最終的な精製RNA 溶液を得た。
【００７３】
精製total RNA の純度を示す260 nm/280 nm 比は光学的測定装置を用いて測定した。さらに精製total RNA はAgilent 社BioAnalyser2100 を用いてその泳動パターンを検査した。泳動にはRNA6000nano チップ(Agilent 社製) を用い、操作法は製造会社の使用法に準じて行った。精製total RNA の濃度測定は光学的測定装置を用いた方法で行った。
【００７４】
光学的測定装置として、分光光度計による場合は光路長1cm の場合1{吸光度(260 nm)-吸光度(320 nm)} = 40 μg RNA/mL として計算した。
【００７５】
（２）蛍光標識cDNA の作製
精製total RNA 20 または10 μg を用いて、CyScript 逆転写酵素を用いて蛍光ラベル化cDNA を作製した。
【００７６】
total RNA またはmRNA(自家精製、クロンテック社製または北山ラベス株式会社製)を用い、以下の反応条件で蛍光ラベル化を行った。最終反応液量は50μL で行った。
【００７７】
反応溶液の終濃度はtotal RNA 10-20μg、Anchored Oligo(dT) (Amersham Bioscience 社製) 0.075μg、dCTP Nucleotide Mix 1μL、1×CyScript buffer (Amersham Bioscience 社製)、10 mM DTT、Cy3-dCTP またはCy5-dCTP (Amersham Bioscience 社製) 1 nmol である。
【００７８】
RNA とAnchored Oligo(dT)の混合液を70℃ 5 分静置した後、氷上に1 分静置した。その他の成分を添加した後CyScript 逆転写酵素(Amersham Bioscience 社製)を100 unit添加し42℃ 90 分遮光静置した。1 N NaOH 12.5μL を添加し、65℃ 10 分遮光静置し、RNA を加水分解した。1 N HCl 15μL を添加し中和した後QIAGEN 社製MinElute PCR Purification Kit を用いて蛍光ラベル体を精製した。未反応の蛍光標識核酸を除去するためにPE バッファーによる洗浄は2 回行った。
【００７９】
投与群（高用量、低用量）をCy3 で標識し、コントロール群をCy5 で標識した。コントロール群の4 匹から作製したCy5 標識cDNA は精製が終了した時点で等量づつ混合し、ハイブリダイゼーション溶液作製には混合液から分注した。
【００８０】
（３）cDNAマイクロアレイへのハイブリダイゼーション及び検出
最終的なハイブリダイゼーション溶液12μL 当たりにCy3 ラベル体4.2 μL、Cy5ラベル体4.2 μL、20×SSC 3 μL、10%SDS 0.6μL になるように混合した後、95℃ 2分遮光静置した。室温3 分遮光静置した後、BSA ブロッキング処理済みのcDNAマイクロアレイに12μL/アレイでスポットした。速やかに24 mm × 32 mm カバーグラス(松浪硝子製)でカバーした後、ハイブリカセット(Telechem 社製、日立ソフトウェア社製等)にセットし55℃で一晩遮光静置した。その後ハイブリカセットからcDNAマイクロアレイを取り出し，2×SSC/0.1%SDS 溶液に浸漬してカバーグラスを取り除いた。室温で2×SSC/0.1%SDS溶液に20 分遮光浸漬し、さらに0.2×SSC/0.1%SDS 溶液に20 分遮光浸漬した。42℃の0.2×SSC/0.1%SDS 溶液に2 回20 分遮光浸漬した後、0.2×SSC/0.1%SDS 溶液、0.05×SSC溶液で洗浄した。遠心器を用いて乾燥させた後、Agilent 社製MicroArray Scanner を用いて検出した。検出感度を示すPMT 値は100％を用いた。スキャン画像から各スポットの蛍光値の数値化はAxon Instruments 社製GenePix Pro ver.4.0.1.17 を用いて行った。
【００８１】
cDNAマイクロアレイのBSA ブロッキングは0.1 または0.22μm のフィルターで膜ろ過した1％BSA ブロッキング溶液(1％BSA、4×SSC、2％SDS)中で42℃ 45 分遮光静置することで行った。2 回水中で振とう洗浄し遠心器を用いて乾燥させた後、ハイブリダイゼーションに用いた。
【００８２】
（４）解析に有効な遺伝子の選定
cDNAマイクロアレイに搭載した各遺伝子の発現量データに対して、以下に示すデータクレンジングを行い、解析に使用する遺伝子の選定を行った。
【００８３】
a)  cDNAマイクロアレイの数値化発現データには、それぞれのアレイ上のスポット（遺伝子）の有効性を表すフラグ（flag）がついており、フラグが"0"（有効）を示すスポットを選定した。選定後の解析可能なスポット数が、3000未満のアレイについてはデータを用いないこととした。
【００８４】
b)  陰性対照遺伝子のlambda DNAの発現量に標準偏差の2倍を加えた値以下の発現量を示す遺伝子は、ノイズデータと判断して除外した。
【００８５】
c)  4個体から抽出した2個体間の相関計数を計算し、総当りで6通りの組み合わせのうち相関係数が0.5以下を2回以上示す場合は、相関の悪いデータとみなし解析対象から除外した。
【００８６】
d)  実験毎にばらつきの少ない遺伝子を選定するため、サンプル投与群のシグナル値（n=4）と媒体対照群（コントロール）のシグナル値(n=4)間でウエルチｔ検定（Welch's t-test）を行い、p≦0.05である遺伝子を選定した（帰無仮説：サンプル投与群（n=4）のシグナル値と媒体対照群（n=4）の遺伝子発現量（シグナル値）には差がない）。
【００８７】
e)  統計的に有意といえるn=4の平均発現比を示す遺伝子を選定するため、検出力を基準にしきい値を算出した。結果、LogRatio（化合物投与群の発現量と媒体対照群の発現量の比をLog変換した値）が0を基準に0.8以上離れていれば、有意な発現比であることが確認できた為、LogRatio=|0.8|をしきい値とし、LogRatioが0.8以上-0.8以下変動している遺伝子を選定した。
【００８８】
上記a)〜e)に従ってデータクレンジングを行ったのち、以降の解析で得られる発がん物質グループで変動していることが確認された遺伝子を選定した結果、1435であった。これらの遺伝子の塩基配列を配列番号1〜1435に示す。
【００８９】
更に、配列番号1〜1435に示す塩基配列に相補的な塩基配列を有する遺伝子のユニジーン（UniGene）番号と、それぞれの塩基配列に付与したクローンIDとを表5〜26に示す。
【００９０】
【表５】
0004841279.tif 000006
【００９１】
【表６】
0004841279.tif 000007
【００９２】
【表７】
0004841279.tif 000008
【００９３】
【表８】
0004841279.tif 000009
【００９４】
【表９】
0004841279.tif 000010
【００９５】
【表１０】
0004841279.tif 000011
【００９６】
【表１１】
0004841279.tif 000012
【００９７】
【表１２】
0004841279.tif 000013
【００９８】
【表１３】
0004841279.tif 000014
【００９９】
【表１４】
0004841279.tif 000015
【０１００】
【表１５】
0004841279.tif 000016
【０１０１】
【表１６】
0004841279.tif 000017
【０１０２】
【表１７】
0004841279.tif 000018
【０１０３】
【表１８】
0004841279.tif 000019
【０１０４】
【表１９】
0004841279.tif 000020
【０１０５】
【表２０】
0004841279.tif 000021
【０１０６】
【表２１】
0004841279.tif 000022
【０１０７】
【表２２】
0004841279.tif 000023
【０１０８】
【表２３】
0004841279.tif 000024
【０１０９】
【表２４】
0004841279.tif 000025
【０１１０】
【表２５】
0004841279.tif 000026
【０１１１】
【表２６】
0004841279.tif 000027
【０１１２】
（５）発がん物質の分類
発がん物質39物質（表１〜４に示した発がん物質からジ（2-エチルヘキシル）フタレートとフェノバルビタールを除外）及び非発がん物質20物質（表１〜４に示した非発がん物質）の投与期間28日間の投与群について、各物質のAmes試験結果をもとに化合物の分類を行った。その結果、Ames陽性発がん物質は19物質、Ames陽性非発がん物質は14物質であった。また、Ames陰性発がん物質は20物質、そしてAmes陰性非発がん物質は6物質であった。
【０１１３】
（６）予測式作成手順と検証方法
まず、Ames陽性発がん物質群とAmes陽性非発がん物質群の間で発現量差が有意である遺伝子セットを統計量の一つであるウェルチt値（Welch's t-value）とフィルタリングの手法をもちいて表5〜26に示す遺伝子群より選定し、それぞれの手法で選定した遺伝子セットを併せたものを次の予測式作成へと用いた。ウエルチｔ値については、発がん物質グループのLogRatioの平均値の絶対値が、非発がん物質グループのLogRatioの平均値の絶対値より大きい遺伝子を選定した。フィルタリングの手法とは、ウェルチt値では選定しにくい一部の発がん物質の予測に有効な遺伝子セットを探索する為の方法である。具体的には、発がん物質群でLogRatioの絶対値が0.8以上であった物質がX個以上で、非発がん物質群でLogRatioの絶対値が0.8以下であった物質がY個以上を示す遺伝子をXとYを変化させて選定し、ウェルチt値で選定した遺伝子セットを合わせて予測式を作成し、最も予測率のよい遺伝子セットを探索した。また、Ames陰性発がん物質群とAmes陰性非発がん物質群の間で発現量差が有意である遺伝子セットについても同様にウェルチt値とフィルタリングの手法で選定を行った。
【０１１４】
Support Vector Machine（SVM）という教師付き分類方法に基づいて作成されたフリーソフトSVMlight（URL http://svmlight.joachims.org/ から入手）を用い、以下の手順でAmes陽性物質の発がん性予測式を作成した。
【０１１５】
a)  発がん物質グループと、非発がん物質グループに属する各物質数が等しくない場合、物質数の少ないグループの物質数と同じ数だけ、物質数の多い方の物質グループからランダムサンプリングを行い、両グループの物質数を等しくした。
【０１１６】
b)  物質数を等しくした両物質グループと、前もって選定した発がん、非発がん物質群間で発現量差が有意である遺伝子セットを用いてSVMlightで学習させ、予測式を作成した。
【０１１７】
c)  予測式の有効性を検証する為に、予測式を作成する為に用いた両物質グループ（トレーニングセット）の発がん性予測と、検証手法の一つである物質のリーブ  ワン  アウト（Leave one out）によるクロスバリデーションを行なった。
【０１１８】
d)  予測式作成に用いていない残りの物質データ（テストセット）の予測を行なった。
【０１１９】
e)  最も高い予測率を示す遺伝子セットが得られるまで選定条件を調整し、a)〜e)までの手順を繰り返し行った。
【０１２０】
（７）各Ames分類毎の予測式の予測結果
Ames陽性物質の発がん性予測式は、（６）a)〜e)の手順に従い遺伝子を選定し、クローンIDが11648、14929、17343、22524、23126、26054、29775、31823、33592、33615、33649、35603、35953、36476、38445、38980、39147、40318、40694、41351、41608、41764、41875、42265、42759、42763、45500、45533、46291、46451、46526の発現量データを用いて作成した。そして、構築した予測式を用いて、トレーニングセットとテストセットの予測を行なった結果を表２７に示す。
【０１２１】
【表２７】
0004841279.tif 000028
【０１２２】
結果、トレーニングセットは、ほぼ予測できた。Ames陽性物質の発がん性予測式は、L00クロスバリデーションにより81.8％の確率で、未知のAmes陽性の発がん性を予測できることが確認できた。また、テストセットのAmes陽性発がん物質は100％予測できた。
【０１２３】
さらに同様の手順で、Ames陰性物質の発がん性予測式は、クローンIDが02620、05151、11648、15551、15562、15777、17912、18209、18327、18636、18810、20938、22117、22953、24206、24614、27762、29598、30778、32657、32667、32851、33469、33615、35953、36351、37092、37685、37936、38674、38764、38846、38980、39701、39766、40078、40337、40673、40981、41275、41433、42065、42510、43282、43512、44281、45238、45463、45718の遺伝子の発現量データを用いて作成した。その予測式を用いて各々のトレーニングセット、テストセットを予測した結果を表２８に示す。
【０１２４】
【表２８】
0004841279.tif 000029
【０１２５】
Ames陰性物質の発がん性予測結果についても、Ames陽性物質の結果と同じくトレーニングセットとテストセットの非発がん物質は比較的高い予測率を示した。L00クロスバリデーションの結果、76.9％予測できることが確認できた。
【０１２６】
（８）予測フローの構築
（７）で作成したAmes陽性物質グループと、Ames陰性物質グループの予測式を組み合わせ、発がん物質を予測する予測フローを構築した（図１）。予測式作成に用いた59物質のデータを、予測フローに適用した結果、Ames陽性物質の発がん性予測式では発がん物質が18物質予測できたが、非発がん物質4物質を発がん物質と誤って予測した。一方、Ames陰性物質の発がん性予測式では発がん物質が18物質予測できた。両予測式の結果を集約させた結果、36物質の発がん物質を予測することができた。
【０１２７】
実施例２
実施例１（１）〜（５）で得られたデータを使用してAmes陽性発がん物質、Ames陽性非発がん物質の発現パターンを更に下記の手順に従ってグループに分類し、化学物質の発がん性の予測を行った。
【０１２８】
（１）発がん物質の分類
発がん物質39物質（表１〜４に示した発がん物質からジ（2-エチルヘキシル）フタレートとフェノバルビタールを除外）及び非発がん物質20物質（表１〜４に示した非発がん物質）について、遺伝子の発現パターンが近い物質ごとに発がん物質を分類して予測式を作成するために、次のような処理を行った。
【０１２９】
まず、Ames陽性発がん物質群とAmes陽性非発がん物質群の間で発現量差が有意である遺伝子セットを、ウエルチｔ値により選定した。選定した遺伝子セットの発現パターンの分類をGene Maths（Applied Maths, Sint-Martens-Latem, Belgium）のクラスタ解析機能を用いて行い、樹形図の形状からクラスタを決定した。ウエルチｔ値の選定条件を変化させ、３パターンの遺伝子セットにより物質のクラスタ構成について確認を行い、遺伝子セットを変化させても安定して同じクラスタを形成する２つの発がん物質のグループ１、２を選定した。その遺伝子セットの選定条件（ウエルチｔ値の値）と、選定した遺伝子のクローンＩＤを表２９に示す。更に、それぞれのウエルチｔ値で選定した遺伝子セットを用いてクラスタ解析した樹形図を図２〜４に示す。図２〜４中、●はAmes陽性発がん物質、○はAmes陽性非発がん物質を表す。
【０１３０】
【表２９】
0004841279.tif 000030
【０１３１】
次に、Ames陰性発がん物質とAmes陰性非発がん物質を用いて発現パターンの近い物質グループを同じ方法で選定した。遺伝子セットの選定条件（ウエルチｔ値の値）と、選定した遺伝子のクローンＩＤを表３０に示す。更に、それぞれのウエルチｔ値で選定した遺伝子セットを用いてクラスタ解析した樹形図を図５〜７に示す。図５〜７中、●はAmes陰性発がん物質、○はAmes陰性非発がん物質を表す。
【０１３２】
【表３０】
0004841279.tif 000031
【０１３３】
各発がん物質グループに分類された化学物質は以下の通りである。
（Ames陽性発がん物質グループ１）
TS011,013,014,025,041,042,046,047,048,049,063,065
（Ames陽性発がん物質グループ２）
TS004,006,045,052,077
（Ames陰性発がん物質グループ１）
TS001,003,019,021,022,023,024,051,055,056,066,067,068,070,078
（２）予測式作成手順と検証方法
予測に有効な遺伝子は、Ames陽性発がん物質グループ１、２に対してはAmes陽性非発がん物質グループ、そしてAmes陰性発がん物質グループ１に対してはAmes陰性非発がん物質グループを組み合わせ、各々の発がん物質グループと非発がん物質グループの間で、発現量差が有意である遺伝子セットを統計量の一つであるウエルチｔ値とフィルタリングの手法によって、表5〜26に示す遺伝子群から選定した。そして、それぞれの手法で選定した遺伝子セットを併せたものを次の予測式作成へと用いた。ウエルチｔ値については、発がん物質グループのLogRatioの平均値の絶対値が、非発がん物質グループのLogRatioの平均値の絶対値より大きい遺伝子を選定した。また、フィルタリングの手法により発がん物質群でLogRatioの絶対値が0.8以上であった物質がX個以上で、非発がん物質群でLogRatioの絶対値が0.8以下であった物質がY個以上を示す遺伝子をXとYを変化させて選定し、ウェルチt値で選定した遺伝子セットを合わせて予測式を作成し、最も予測率のよい遺伝子セットを探索した。Ames陰性発がん物質群とAmes陰性非発がん物質群の間で発現量差が有意である遺伝子セットについても同様にウェルチt値とフィルタリングの手法で選定を行った。
【０１３４】
フリーソフトSVMlightを用い、実施例１（６）のa)〜e)と同じ手順でAmes陽性物質の発がん性予測式を作成した。
【０１３５】
（３）各Ames分類毎の予測式の予測結果
Ames陽性物質グループ１の発がん性予測式は、（２）の手順に従い遺伝子を選定し、クローンIDが01371、05376、17343、34871、40104、41764、46526の発現量データを用いて作成した。構築した予測式を用いて、トレーニングセットとテストセットの予測を行なった結果を表３１に示す。
【０１３６】
【表３１】
0004841279.tif 000032
【０１３７】
結果、トレーニングセットは、ほぼ予測できた。Ames陽性物質グループ１の発がん性予測式は、L00クロスバリデーションにより92.3％の確率で、未知のAmes陽性の発がん性グループ１を予測できることが確認できた。また、テストセットのAmes陽性発がん物質は28.6％予測できた。
【０１３８】
同様にAmes陽性物質グループ２の発がん性予測式は、クローンIDが09104、20106、29121、30146、31305、31578、31831、32851、36121、38445、39766、40668、40694、40805、41875、42265、45500、45533、46291の発現量データを用いて作成した。構築した予測式を用いて、トレーニングセットとテストセットの予測を行なった結果を表３２に示す。
【０１３９】
【表３２】
0004841279.tif 000033
【０１４０】
結果、トレーニングセットは、ほぼ予測できた。Ames陽性物質グループ２の発がん性予測式は、L00クロスバリデーションにより73.7％の確率で、未知のAmes陽性グループ２の発がん性を予測できることが確認できた。また、テストセットのAmes陽性発がん物質は57.1％予測できた。
【０１４１】
さらに同様の手順で、Ames陰性物質グループ１の発がん性予測式は、クローンＩＤが17917、18636、20938、33469、35953、36351、37685、38764、38846、39766、40337、41253、42510、45238、46254、46675の遺伝子の発現量データを用いて作成した。その予測式を用いて各々のトレーニングセット、テストセットを予測した結果を表３３に示す。
【０１４２】
【表３３】
0004841279.tif 000034
【０１４３】
Ames陰性物質グループ１の予測結果についても、Ames陽性物質の結果と同じくトレーニングセットとテストセットの非発がん物質は比較的高い予測率を示した。L00クロスバリデーションの結果、95.2％予測できることが確認できた。
【０１４４】
（４）予測フローの構築
（３）で作成したAmes陽性物質グループと、Ames陰性物質グループの予測式を組み合わせ、発がん物質を予測する予測フローを構築した（図８）。予測式作成に用いた59物質のデータを、予測フローに適用した結果、Ames陽性物質の発がん性予測式では発がん物質が19物質予測できたが、非発がん物質6物質を発がん物質と誤って予測した。一方、Ames陰性物質の発がん性予測式では発がん物質が18物質予測できたが、非発がん物質1物質を発がん物質と誤って予測した。両予測式の結果を集約させた結果、37物質の発がん物質を予測することができた。
【０１４５】
製造例２（オリゴマイクロアレイの作製）
（１）オリゴマイクロアレイに搭載する遺伝子の選択
オリゴマイクロアレイに搭載する遺伝子は、以下の基準から6,785遺伝子を選定した。
【０１４６】
a)  製造例１で得たcDNAマイクロアレイを用いた試験において、発がん物質の内2物質以上の投与によりコントロール動物と比較して2倍以上あるいは1/2 以下に発現量が変動した遺伝子。
【０１４７】
b)  Affymetrix社製 GeneChip U34A アレイを用いた試験において、発がん物質の内2物質以上の投与によりコントロール動物と比較して2倍以上あるいは1/2 以下に発現量が変動した遺伝子。
【０１４８】
c)  文献等で発がん性に関連していることが報告されている遺伝子。これらの遺伝子の中にはがん遺伝子やがん抑制遺伝子、転写因子、増殖因子などが含まれる。
【０１４９】
d)  ポジティブコントロールとして3 遺伝子(ベータ-アクチン、GAPDH、ユビキチン)及び製造例１で得たcDNAマイクロアレイを用いた試験において発現が変動しなかった遺伝子。ただしベータ-アクチン及びGAPDHについては3'末端、中央部、5'末端部分についてそれぞれプローブを作成した。
【０１５０】
e)  ネガティブコントロールとしてラムダファージ遺伝子、バクテリオファージ遺伝子（cre）及びバクテリア由来の8遺伝子（bioB、bioC、bioD、jojG、lys、pheA、thrB、trpF）。
【０１５１】
a)に含まれる遺伝子は4,461遺伝子、b)に含まれるのは4,386遺伝子、c)に含まれる遺伝子は1,714遺伝子、d)に含まれる遺伝子は76遺伝子であった。ただし、それぞれのグループに重複して含まる遺伝子が存在するため、単純合計は6,785遺伝子とはならない。
【０１５２】
（２）オリゴマイクロアレイの作製
選定した遺伝子の塩基配列からタカラバイオ株式会社の設計アルゴリズムを用いてお互いに非相補的な60塩基長の部分配列を設計した。（１） で選択した6,785遺伝子の内有効な配列を設計できなかった遺伝子を除く6,722遺伝子について6,726の60塩基長の部分配列を設計し、オリゴを合成した。60塩基オリゴの合成はオリゴの5'末端に基盤との共有結合に用いるアミノ基を加えてタカラバイオ株式会社に委託した。得られた60塩基オリゴを3次元表面構造のガラス基板（Amersham Bioscience社製）上にスポットし、オリゴマイクロアレイを得た。1ガラス当たり1アレイをスポットし、ハイブリダイゼーション操作の簡略化のために表面にチャンバーシールを貼り付けた。
【０１５３】
実施例３
（１）total RNAの取得
表１〜４に示す化学物質から1-ニトロプロパンを除く62物質と、表３４に示す化学物質11物質とをそれぞれ媒体に溶解し、溶液を調製した。
【０１５４】
【表３４】
0004841279.tif 000035
【０１５５】
日本チャールス・リバー社から入手した５週令の雄性ラット（Ｆ３４４、ＳＰＦ系統）を４匹／群に分け、各群のラットに調製した溶液又は媒体を強制経口投与した。投与は１日１回、１、３、７、１４、２８日間行った。投与開始から２、４、８、１５、又は２９日目に各個体の肝臓を採取した。肝臓から断片を切り出し、実施例１（１）と同様にしてtotal RNA の抽出、精製を行った。
【０１５６】
（２）蛍光標識cRNA の作製
（１）で得た精製total RNA 5μgを用いて、SuperScript Indirect RNA Amplification System（Invitrogen社製）を用いて蛍光cRNAを作製した。投与群では各検体毎に、コントロール群は等重量を混合した後に、以下の反応条件で蛍光ラベル化を行った。
【０１５７】
total RNA 5μgとT7-Oligo(dT) primer 1μLの混合液 10μLを70℃ 10分静置した後、氷上に2分静置した。5X First-Strand buffer 4μL、0.1 M DTT 2μL、10 mM dNTP Mix 1μL、RNaseOUT (40 U/μL) 1μL、SuperScript III逆転写酵素 (200 U/μL) 2μLを添加し、46℃ 2時間遮光静置した。70℃ 10分静置した後、氷上に2分静置した。DEPC処理水 91μL、5X Second-Strand Buffer 30μL、10 mM dNTP Mix 3μL、E. coli DNA Polymerase I (10 units/μL) 4μL、E. coli DNA Ligase (10 units/μL) 1μL、E. coli RNase H (2 units/μL) 1μLを添加し、16℃ 2時間遮光静置した。作成された2本鎖cDNAにcDNA Loading Buffer 500μL を添加した後、cDNA Purification Moduleを用いて精製した。20μLのDEPC処理水で溶出した精製cDNA溶液の内 14μLにDEPC処理水 7.75μL、100 mM ATP 1.5μL、100 mM CTP 1.5μL、100 mM GTP 1.5μL、100 mM UTP 0.75μL、50 mM aa-UTP 2μL、10X T7 Reaction Buffer 4μL、T7 Enzyme Mix 7μLを添加し、37℃ 3時間静置した。さらにDNase I (1 unit/μL) 2μLを添加し、37℃ 30分間静置した。得られた溶液にDEPC処理水 58μL、Buffer RLT (QIAGEN社) 350μLを添加した後、RNeasy MinElute Cleanup Kit (QIAGEN社)を用いてcRNAを精製した。20μLのDEPC処理水で溶出した精製cRNA溶液の内 6μLにCoupling Buffer 10μL及びDMSO 4μLに溶解したCy3 reactive dye for one labelling reactionまたはCy5 reactive dye for one labelling reaction (Amersham Bioscience社製)を添加し、室温で30分間遮光で静置した。投与群をCy3 で標識し、コントロール群をCy5 で標識した。4M Hydroxyamine 4.5μLを添加し、室温で15分間遮光で静置することで反応を停止した。得られた溶液にDEPC処理水 75.5μL、Buffer RLT (QIAGEN社) 350μLを添加した後、RNeasy MinElute Cleanup Kit (QIAGEN社)を用いて標識cRNAを精製した。14μLのDEPC処理水で溶出することで精製標識cRNAを得た。
【０１５８】
（３）オリゴマイクロアレイへのハイブリダイゼーション及び検出
分光光度計を用いて得られた精製標識cRNAの濃度を決定した後、Cy3標識cRNA 0.75μg、Cy5標識cRNA 0.75μg及び5X Fragmentation Buffer 8μLを混合しDEPC処理水で総量40μLにした。94℃ 15分間遮光で静置した後、氷上に2分静置した。得られた溶液から30μLを分取し、DEPC処理水 46.5μL、20XSSC 37.5μL、10％ SDS 7.5μL及び50Xデンハルト溶液 12μを添加した後、95℃ 2分間遮光で静置した。氷上に2分静置した後、Salmon Sperm DNA (Invitrogen社) 1.5μL及びホルムアミド 15μLを添加し、60℃ 3分間遮光で静置することで最終的なハイブリダイゼーション溶液を得た。オリゴマイクロアレイに120-140μLのハイブリダイゼーション溶液を注入した後に、注入口をシールし、60℃で一晩遮光静置した。その後オリゴマイクロアレイの表面シールを剥離し、素早く5×SSC/0.1%SDS 溶液に浸漬した。30℃の5×SSC/0.1%SDS溶液に10分遮光浸漬し、さらに30℃の0.5×SSC溶液に2回2分遮光浸漬した。0.5×SSC/0.01%Tween20溶液で洗浄した後、遠心器を用いて乾燥させた。Agilent Technologies社製Microarray Scanner を用いて検出し、検出感度を示すPMT 値は70-100％を用いた。スキャン画像から各スポットの蛍光値の数値化はAxon Instruments 社製GenePix Pro ver.4.0.1.17 を用いて行った。
【０１５９】
（４）解析に有効な遺伝子の選定
b)の陰性対照遺伝子にjojG、bioD、creを使用した以外は、実施例１（４）のa)〜e)と同様にしてオリゴマイクロアレイに搭載した各遺伝子の発現量データに対してデータクレンジングを行い、解析に使用する遺伝子の選定を行った。さらに、以降の解析で得られる発がん物質グループで変動していることが確認された遺伝子を選定した結果、解析に有効な遺伝子数は922であった。これらの遺伝子の塩基配列を配列番号1436〜2357に示す。
【０１６０】
更に、配列番号1436〜2357に示す塩基配列に相補的な塩基配列を有する遺伝子のユニジーン（UniGene）番号と、それぞれの塩基配列に付与したクローンIDとを表35〜48に示す。
【０１６１】
【表３５】
0004841279.tif 000036
【０１６２】
【表３６】
0004841279.tif 000037
【０１６３】
【表３７】
0004841279.tif 000038
【０１６４】
【表３８】
0004841279.tif 000039
【０１６５】
【表３９】
0004841279.tif 000040
【０１６６】
【表４０】
0004841279.tif 000041
【０１６７】
【表４１】
0004841279.tif 000042
【０１６８】
【表４２】
0004841279.tif 000043
【０１６９】
【表４３】
0004841279.tif 000044
【０１７０】
【表４４】
0004841279.tif 000045
【０１７１】
【表４５】
0004841279.tif 000046
【０１７２】
【表４６】
0004841279.tif 000047
【０１７３】
【表４７】
0004841279.tif 000048
【０１７４】
【表４８】
0004841279.tif 000049
【０１７５】
（５）発がん物質の分類
発がん物質47物質（表１〜４及び３４に示した発がん物質から四塩化炭素とグルタルアルデヒドを除外）及び非発がん物質24物質（表１〜４及び３４に示した非発がん物質）の投与期間28日間の投与群について、各物質のAmes試験結果をもとに化合物の分類を行った。その結果、Ames陽性発がん物質は23物質、Ames陽性非発がん物質は13物質であった。また、Ames陰性発がん物質は24物質、そしてAmes陰性非発がん物質は11物質であった。
【０１７６】
（６）予測式作成手順と検証方法
予測に有効な遺伝子は、Ames陽性発がん物質群とAmes陽性非発がん物質群の間で発現量差が有意である遺伝子セットを統計量の一つであるウェルチt値とフィルタリングの手法によって、表35〜48に示す遺伝子群から選定した。そして、それぞれの手法で選定した遺伝子セットを併せたものを次の予測式作成へと用いた。ウエルチｔ値については、発がん物質グループのLogRatioの平均値の絶対値が、非発がん物質グループのLogRatioの平均値の絶対値より大きい遺伝子を選定した。フィルタリングの手法とは、ウェルチt値では選定しにくい一部の発がん物質の予測に有効な遺伝子セットを探索する為の方法である。具体的には、発がん物質群でLogRatioの絶対値が0.8以上であった物質がX個以上で、非発がん物質群でLogRatioの絶対値が0.8以下であった物質がY個以上を示す遺伝子をXとYを変化させて選定し、ウェルチt値で選定した遺伝子セットを合わせて予測式を作成し、最も予測率のよい遺伝子セットを探索した。また、Ames陰性発がん物質群とAmes陰性非発がん物質群の間で発現量差が有意である遺伝子セットについても同様にウェルチt値とフィルタリングの手法により選定した。
【０１７７】
フリーソフトSVMlightを用い、実施例１（６）のa)〜e)と同じ手順でAmes陽性物質の発がん性予測式を作成した。
【０１７８】
（７）各Ames分類毎の予測式の予測結果
Ames陽性物質の発がん性予測式は、（６）の手順に従い遺伝子を選定し、クローンIDが0136、0297、0336、0518、0523、0580、0688、0719、0987、1091、1496、1512、1643、1662、1671、1754、1912、1940、2286、2372、2471、2654、2667、2754、2786、2788、2879、3211、3460、3469、3629、3942、4047、4096、4173、4476、4673、4688、4784、4830、4846、4975、5171、5231、5299、5500、5574、5600、5643、5670、5706、5791、6182、6196、6570、6670の発現量データを用いて作成した。そして、構築した予測式を用いて、トレーニングセットとテストセットの予測を行なった結果を表４９に示す。
【０１７９】
【表４９】
0004841279.tif 000050
【０１８０】
結果、トレーニングセットは、ほぼ予測できた。Ames陽性物質の発がん性予測式は、L00クロスバリデーションにより75％の確率で、未知のAmes陽性の発がん性を予測できることが確認できた。また、テストセットのAmes陽性発がん物質は70％予測できた。
【０１８１】
同様の手順で、Ames陰性物質の発がん性予測式は、クローンIDが0007、0187、0725、0745、0762、1294、1764、1845、2873、3201、3286、3305、3903、4273、4275、4293、4697、4770、4846、5224、5515、6109、6182、6343の遺伝子の発現量データを用いて作成した。その予測式を用いて各々のトレーニングセット、テストセットを予測した結果を表５０に示す。
【０１８２】
【表５０】
0004841279.tif 000051
【０１８３】
Ames陰性物質の発がん性予測結果についても、Ames陽性物質の結果と同じくトレーニングセットとテストセットの非発がん物質は比較的高い予測率を示した。L00クロスバリデーションの結果、82.9％予測できることが確認できた。
【０１８４】
（８）予測フローの構築
（７）で作成したAmes陽性物質グループと、Ames陰性物質グループの予測式を組み合わせ、発がん物質を予測する予測フローを構築した（図９）。予測式作成に用いた71物質のデータを、予測フローに適用した結果、Ames陽性物質の発がん性予測式では発がん物質が20物質予測できたが、非発がん物質1物質を発がん物質と誤って予測した。一方、Ames陰性物質の発がん性予測式では発がん物質が22物質予測できたが、非発がん物質1物質を発がん物質と誤って予測した。両予測式の結果を集約させた結果、42物質の発がん物質を予測することができた。
【０１８５】
実施例４
実施例３（１）〜（５）で得られたデータを使用してAmes陽性発がん物質、Ames陽性非発がん物質の発現パターンを更に下記の手順に従ってグループに分類し、化学物質の発がん性の予測を行った。
【０１８６】
（１）発がん物質の分類
発がん物質47物質（表１〜４及び３４に示した発がん物質から四塩化炭素とグルタルアルデヒドを除外）及び非発がん物質24物質（表１〜４及び３４に示した非発がん物質）について、遺伝子の発現パターンが近い物質ごとに発がん物質を分類して予測式を作成するために、次のような処理を行った。
【０１８７】
まず、Ames陽性発がん物質群とAmes陽性非発がん物質群の間で発現量差が有意である遺伝子セットを、ウエルチｔ値により選定した。選定した遺伝子セットの発現パターンの分類をGene Maths（Applied Maths, Sint-Martens-Latem, Belgium）のクラスタ解析機能を用いて行い、樹形図の形状からクラスタを決定した。ウエルチｔ値の選定条件を変化させ、数パターンの遺伝子セットにより物質のクラスタ構成について確認を行い、遺伝子セットを変化させても安定して同じクラスタを形成する発がん物質のグループ１を選定した。その遺伝子セットの選定条件（ウエルチｔ値の値）と、選定した遺伝子のクローンＩＤを表５１に示す。更に、それぞれのウエルチｔ値で選定した遺伝子セットを用いてクラスタ解析した樹形図を図１０〜１２に示す。図１０〜１２中、●はAmes陽性発がん物質、○はAmes陽性非発がん物質を表す。
【０１８８】
【表５１】
0004841279.tif 000052
【０１８９】
Ames陽性発がん物質グループ１を除いた残りのAmes陽性発がん物質と、全非発がん物質を用いて、残りの発がん物質中に存在する発現パターンの近い物質グループを同じ方法で選定した。即ち、残りのAmes陽性発がん物質投与群と、全非発がん物質投与群の間で特徴的に変動している遺伝子をウエルチｔ値で選定した後、その遺伝子セットを用いてクラスタ解析を行った。ウエルチｔ値の選定条件を変化させ、数パターンの遺伝子セットにより物質のクラスタ構成の確認を行い、常に安定してクラスタを形成する発がん物質グループ２、３を選定した。遺伝子セットの選定条件（ウエルチｔ値の値）と、選定した遺伝子のクローンＩＤを表５２に示す。更に、それぞれのウエルチｔ値で選定した遺伝子セットを用いてクラスタ解析した樹形図を図１３〜１５に示す。図１３〜１５中、●は発がん物質、○は非発がん物質を表す。
【０１９０】
【表５２】
0004841279.tif 000053
【０１９１】
次に、Ames陰性発がん物質とAmes陰性非発がん物質を用いて発現パターンの近い物質グループを同じ方法で選定した。遺伝子セットの選定条件（ウエルチｔ値の値）と、選定した遺伝子のクローンＩＤを表５３に示す。更に、それぞれのウエルチｔ値で選定した遺伝子セットを用いてクラスタ解析した樹形図を図１６〜１８に示す。図１６〜１８中、●はAmes陰性発がん物質、○はAmes陰性非発がん物質を表す。
【０１９２】
【表５３】
0004841279.tif 000054
【０１９３】
各発がん物質グループに分類された化学物質は以下の通りである。
（Ames陽性発がん物質グループ１）
TS004,006,011,013,014,042,043,076
（Ames陽性発がん物質グループ２）
TS041,047,048,049,065
（Ames陽性発がん物質グループ３）
TS028,044,050,053,077
（Ames陰性発がん物質グループ１）
TS008,018,020,021,023,024,055,056,070
（Ames陰性発がん物質グループ２）
TS027,029,030,054,067
（２）予測式作成手順と検証方法
予測に有効な遺伝子は、Ames陽性発がん物質群とAmes陽性非発がん物質群の間で発現量差が有意である遺伝子セットを統計量の一つであるウェルチt値とフィルタリングの手法によって、表35〜48に示す遺伝子群から選定した。そして、それぞれの手法で選定した遺伝子セットを併せたものを次の予測式作成へと用いた。ウエルチｔ値については、発がん物質グループのLogRatioの平均値の絶対値が、非発がん物質グループのLogRatioの平均値の絶対値より大きい遺伝子を選定した。フィルタリングの手法により、発がん物質群でLogRatioの絶対値が0.8以上であった物質がX個以上で、非発がん物質群でLogRatioの絶対値が0.8以下であった物質がY個以上を示す遺伝子をXとYを変化させて選定し、ウェルチt値で選定した遺伝子セットを合わせて予測式を作成し、最も予測率のよい遺伝子セットを探索した。また、Ames陰性発がん物質群とAmes陰性非発がん物質群の間で発現量差が有意である遺伝子セットについても同様にウェルチt値とフィルタリングの手法により選定した。
【０１９４】
フリーソフトSVMlightを用い、実施例１（６）a)〜e)と同じ手順でAmes陽性物質の発がん性予測式を作成した。
【０１９５】
（３）各Ames分類毎の予測式の予測結果
Ames陽性物質グループ１の発がん性予測式は、（２）の手順に従い遺伝子を選定し、クローンIDが0065、0688、1974、2286、2654、3460、4673、4771、5231、5569の発現量データを用いて作成した。そして、構築した予測式を用いて、トレーニングセットとテストセットの予測を行なった結果を表５４に示す。
【０１９６】
【表５４】
0004841279.tif 000055
【０１９７】
結果、トレーニングセットは、ほぼ予測できた。Ames陽性物質グループ１の発がん性予測式は、L00クロスバリデーションにより95.2％の確率で、未知のAmes陽性の発がん性グループ１を予測できることが確認できた。また、テストセットのAmes陽性非発がん物質は100％予測できた。しかし、残りのAmes陽性発がん物質は13.3％しか予測できなかった。この結果から、Ames陽性物質グループ１の予測式は、グループ１の変動パターンに特化した予測式であり、グループ２、３は、変動パターンが異なる為に予測できなかったと考えられる。
【０１９８】
同様にAmes陽性物質グループ２の発がん性予測式は、クローンIDが0477、0543、1875、1962、2007、2435、4578、5726、5731、6196、6275の発現量データを用いて作成した。構築した予測式を用いて、トレーニングセットとテストセットの予測を行なった結果を表５５に示す。
【０１９９】
【表５５】
0004841279.tif 000056
【０２００】
結果、トレーニングセットは、ほぼ予測できた。Ames陽性物質グループ２の発がん性予測式は、L00クロスバリデーションにより94.4％の確率で、未知のAmes陽性グループ２の発がん性を予測できることが確認できた。また、テストセットのAmes陽性非発がん物質は100％予測できた。しかし、グループ１の時と同様に、残りのAmes陽性発がん物質は16.7％しか予測できなかった。
【０２０１】
Ames陽性物質グループ３の発がん性予測式は、クローンIDが2171、2877、4047、6701の発現量データを用いて作成した。構築した予測式を用いて、トレーニングセットとテストセットの予測を行なった結果を表５６に示す。
【０２０２】
【表５６】
0004841279.tif 000057
【０２０３】
結果、トレーニングセットは、ほぼ予測できた。Ames陽性物質グループ３の発がん性予測式は、L00クロスバリデーションにより100％の確率で、未知のAmes陽性グループ３の発がん性を予測できることが確認できた。また、テストセットのAmes陽性非発がん物質は100％予測できた。しかし、グループ１、２の時と同様に、残りのAmes陽性発がん物質は38.9％しか予測できなかった。
【０２０４】
Ames陰性物質グループ１の発がん性予測式についても、クローンIDが0187、0725、0745、1431、3201、4275の遺伝子の発現量データを用いて作成した。その予測式を用いて各々のトレーニングセット、テストセットを予測した結果を表５７に示す。
【０２０５】
【表５７】
0004841279.tif 000058
【０２０６】
Ames陰性物質グループ１の予測結果は、これまでの結果と同じくトレーニングセットとテストセットの非発がん物質は高い予測率を示した。L00クロスバリデーションの結果、95％予測できることが確認できた。
【０２０７】
Ames陰性物質グループ２の発がん性予測式についても、クローンIDが1882、2281、4846、5515、5650の遺伝子の発現量データを用いて作成した。その予測式を用いて各々のトレーニングセット、テストセットを予測した結果を表５８に示す。
【０２０８】
【表５８】
0004841279.tif 000059
【０２０９】
Ames陰性物質グループ２の発がん性予測結果についても、これまでの結果と同じくトレーニングセットとテストセットの非発がん物質は高い予測率を示した。L00クロスバリデーションの結果、100％予測できることが確認できた。
【０２１０】
（４）予測フローの構築
（３）で作成したAmes陽性物質グループと、Ames陰性物質グループの予測式を組み合わせ、発がん物質を予測する予測フローを構築した（図１９）。予測式作成に用いた71物質のデータを、予測フローに適用した結果、Ames陽性物質の発がん性予測式では発がん物質が20物質予測できたが、非発がん物質1物質を発がん物質と誤って予測した。一方、Ames陰性物質の発がん性予測式では発がん物質が21物質予測でき、非発がん物質は全て正しく予測できた。両予測式の結果を集約させた結果、41物質の発がん物質を予測することができた。
【図面の簡単な説明】
【０２１１】
【図１】実施例１において構築した発がん性予測方法を示すフロー図である。
【図２】実施例２においてクラスタ解析により化学物質を分類した樹形図である。
【図３】実施例２においてクラスタ解析により化学物質を分類した樹形図である。
【図４】実施例２においてクラスタ解析により化学物質を分類した樹形図である。
【図５】実施例２においてクラスタ解析により化学物質を分類した樹形図である。
【図６】実施例２においてクラスタ解析により化学物質を分類した樹形図である。
【図７】実施例２においてクラスタ解析により化学物質を分類した樹形図である。
【図８】実施例２において構築した発がん性予測方法を示すフロー図である。
【図９】実施例３において構築した発がん性予測方法を示すフロー図である。
【図１０】実施例４においてクラスタ解析により化学物質を分類した樹形図である。
【図１１】実施例４においてクラスタ解析により化学物質を分類した樹形図である。
【図１２】実施例４においてクラスタ解析により化学物質を分類した樹形図である。
【図１３】実施例４においてクラスタ解析により化学物質を分類した樹形図である。
【図１４】実施例４においてクラスタ解析により化学物質を分類した樹形図である。
【図１５】実施例４においてクラスタ解析により化学物質を分類した樹形図である。
【図１６】実施例４においてクラスタ解析により化学物質を分類した樹形図である。
【図１７】実施例４においてクラスタ解析により化学物質を分類した樹形図である。
【図１８】実施例４においてクラスタ解析により化学物質を分類した樹形図である。
【図１９】実施例４において構築した発がん性予測方法を示すフロー図である。
【図１】
0004841279.tif 000060
【図２】
0004841279.tif 000061
【図３】
0004841279.tif 000062
【図４】
0004841279.tif 000063
【図５】
0004841279.tif 000064
【図６】
0004841279.tif 000065
【図７】
0004841279.tif 000066
【図８】
0004841279.tif 000067
【図９】
0004841279.tif 000068
【図１０】
0004841279.tif 000069
【図１１】
0004841279.tif 000070
【図１２】
0004841279.tif 000071
【図１３】
0004841279.tif 000072
【図１４】
0004841279.tif 000073
【図１５】
0004841279.tif 000074
【図１６】
0004841279.tif 000075
【図１７】
0004841279.tif 000076
【図１８】
0004841279.tif 000077
【図１９】
0004841279.tif 000078
【配列表】
0004841279000001.app