(57)【要約】【課題】ウィルス側や患者側の遺伝子に応じた治療の有効性を予測する。【解決手段】有効性判定手段２２は、ヒト遺伝子型と治療の有効性との相関度を示すヒト重み係数と、ウィルス遺伝子型と治療の有効性との相関度を示すウィルス重み係数とを読み出し、被検体のヒト遺伝子型情報とヒト重み係数を遺伝子型毎に乗算してヒト有効性予測個別値を算出し、被検体のウィルス遺伝子型情報とウィルス重み係数を遺伝子型毎に乗算してウィルス有効性予測個別値を算出し、ヒト有効性予測個別値の各々及びウィルス有効性予測個別値の各々を加算して有効性予測加算値を算出し、有効性予測加算値を出力する。【選択図】    図１
2008272510.tif 000002
【特許請求の範囲】
【請求項１】
コンピュータに、治療の有効性を判定するための有効性判定データベースを生成させる処理を実行させる有効性判定データベース生成プログラムであって、この有効性判定データベース生成プログラムは前記コンピュータに、
検体毎に、検体の感染宿主の遺伝子型を示す感染宿主遺伝子型情報と、検体が感染した感染病原体の遺伝子型を示す感染病原体遺伝子型情報と、前記感染病原体に対する治療の有効性を示す有効性情報が与えられた検体データを複数の検体について読み出させ、
前記複数の検体についての検体データに基づき、前記感染宿主遺伝子型情報及び前記感染病原体遺伝子型情報の少なくとも一方と前記治療の有効性との相関度を示す重み係数を、遺伝子型毎に算出させ、
前記重み係数を、遺伝子型毎に記憶装置に格納させる
ことを特徴とする有効性判定データベース生成プログラム。
【請求項２】
１あるいは複数の感染宿主遺伝子型と、治療の有効性との相関度を示す感染宿主重み係数とを記録した第１のフィールドと、
１あるいは複数の感染病原体遺伝子型と、前記治療の有効性との相関度を示す感染病原体重み係数とを記録した第２のフィールドと
を有することを特徴とする有効性判定データのデータ構造。
【請求項３】
１あるいは複数の感染宿主遺伝子型と、治療の有効性との相関度を示す感染宿主重み係数と、１あるいは複数の感染病原体遺伝子型と、前記治療の有効性との相関度を示す感染病原体重み係数とを記録した第１のフィールドと、
前記治療の有効性を示す有効性情報を記録した第２のフィールドと
を有することを特徴とする有効性判定データのデータ構造。
【請求項４】
コンピュータに、治療の有効性を判定する処理を実行させる有効性判定プログラムであって、この有効性判定プログラムは前記コンピュータに、
１あるいは複数の感染宿主遺伝子型と治療の有効性との相関度を示す感染宿主重み係数と、１あるいは複数の感染病原体遺伝子型と前記治療の有効性との相関度を示す感染病原体重み係数とを読み出させ、
被検体の感染宿主遺伝子型情報と、前記感染宿主重み係数を遺伝子型毎に乗算して感染宿主有効性予測個別値を算出させ、
被検体の感染病原体遺伝子型情報と、前記感染病原体重み係数を遺伝子型毎に乗算して感染病原体有効性予測個別値を算出させ、
前記感染宿主有効性予測個別値の各々及び感染病原体有効性予測個別値の各々を加算して有効性予測加算値を算出させ、
前記有効性予測加算値あるいは前記有効性予測加算値に所定の数値を加算した値に基づき治療の有効性を判定し、
前記判定結果を出力させる
ことを特徴とする有効性判定プログラム。
【請求項５】
コンピュータを用いて治療の有効性を判定するための有効性判定データベースを生成する有効性判定データベース生成システムであって、
検体毎に、検体の感染宿主遺伝子型を示す感染宿主遺伝子型情報と、検体が感染した感染病原体の遺伝子型を示す感染病原体遺伝子型情報と、前記感染病原体に対する治療の有効性を示す有効性情報が与えられた検体データを複数の検体について読み出し、前記複数の検体についての検体データに基づき、前記感染宿主遺伝子型情報及び前記感染病原体遺伝子型情報の少なくとも一方と前記治療の有効性との相関度を示す重み係数を、遺伝子型毎に算出するデータベース生成手段と、
前記データベース生成手段により算出された前記重み係数を、遺伝子型毎に記憶する記憶手段と
を具備してなることを特徴とする有効性判定データベース生成システム。
【請求項６】
コンピュータを用いて治療の有効性を判定する治療法有効性判定システムであって、
１あるいは複数の感染宿主遺伝子型と治療の有効性との相関度を示す感染宿主重み係数と、１あるいは複数の感染病原体遺伝子型と前記治療の有効性との相関度を示す感染病原体重み係数とを記憶した有効性判定データベースと、
被検体の感染宿主遺伝子型情報と、前記感染宿主重み係数を遺伝子型毎に乗算して感染宿主有効性予測個別値を算出させ、被検体の感染病原体遺伝子型情報と、前記感染病原体重み係数を遺伝子型毎に乗算して感染病原体有効性予測個別値を算出させ、前記感染宿主有効性予測個別値の各々及び感染病原体有効性予測個別値の各々を加算して有効性予測加算値を算出させ、前記有効性予測加算値あるいは前記有効性予測加算値に所定の数値を加算した値に基づき治療の有効性を判定させる有効性判定手段と、
前記有効性予測加算値を出力する出力装置と
を具備してなることを特徴とする治療法有効性判定システム。
【請求項７】
コンピュータを用いて治療の有効性を判定するための有効性判定データベースを生成する有効性判定データベース生成方法であって、
検体毎に、検体の感染宿主遺伝子型を示す感染宿主遺伝子型情報と、検体が感染した感染病原体の遺伝子型を示す感染病原体遺伝子型情報と、前記感染病原体に対する治療の有効性を示す有効性情報が与えられた検体データを複数の検体について第１の記憶装置から読み出し、
前記複数の検体についての検体データに基づき、前記感染宿主遺伝子型情報及び前記感染病原体遺伝子型情報の少なくとも一方と前記治療の有効性との相関度を示す重み係数を、遺伝子型毎に算出し、
前記重み係数を、遺伝子型毎に第２の記憶装置に記憶する
有効性判定データベース生成方法。
【請求項８】
コンピュータを用いて治療の有効性を判定する治療法有効性判定方法であって、
１あるいは複数の感染宿主遺伝子型と治療の有効性との相関度を示す感染宿主重み係数と、１あるいは複数の感染病原体遺伝子型と前記治療の有効性との相関度を示す感染病原体重み係数とを読み出し、
被検体の感染宿主遺伝子型情報と、前記感染宿主重み係数を遺伝子型毎に乗算して感染宿主有効性予測個別値を算出し、
被検体の感染病原体遺伝子型情報と、前記感染病原体重み係数を遺伝子型毎に乗算して感染病原体有効性予測個別値を算出し、
前記感染宿主有効性予測個別値の各々及び感染病原体有効性予測個別値の各々を加算して有効性予測加算値を算出し、
前記有効性予測加算値又は前記有効性予測加算値に所定の数値を加算した値に基づき治療の有効性を判定し、
前記判定結果を出力する
治療法有効性判定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、遺伝子の存在状態に基づき治療法の有効性を予測するためのプログラム、データベース、システム及び方法に関する。
【背景技術】
【０００２】
患者への薬剤投与などの治療法の有効性は、患者の状態によって変わってくる場合が多い。そのため、患者の臨床データに基づき、これを解析して患者側にフィードバックする仕組みが考案されている（例えば特許文献１参照）。これによれば、患者毎に蓄積された臨床データを解析及び評価し、臨床的危険度をスコア化し、医療者に適切な医療処置を表示する。
【特許文献１】特開２００２−９５６５０号公報（第１頁、左上欄１行〜１８行）
【発明の開示】
【発明が解決しようとする課題】
【０００３】
前述した通り、臨床データをスコア化し、患者の臨床的危険度を表示することにより、医療者に適切な医療処置を施すことが高まる。
【０００４】
しかしながら、例えば患者がウィルスに感染している場合、そのウィルスに対する治療の有効性は臨床データのみにより定まるものではない。この場合、患者側の遺伝子や、ウィルス側の遺伝子などにより治療の有効性が変わるものと考えられている。しかしながら、これらウィルス側や患者側の遺伝子に基づき、治療の有効性を決定する仕組みは存在しなかった。
【０００５】
本発明は上記課題を解決するためになされたもので、その目的とするところは、ウィルス側や患者側の遺伝子に応じた治療の有効性を予測するためのプログラム、データベース、システム及び方法を提供することにある。
【課題を解決するための手段】
【０００６】
この発明の一の観点によれば、コンピュータに、治療の有効性を判定するための有効性判定データベースを生成させる処理を実行させる有効性判定データベース生成プログラムであって、この有効性判定データベース生成プログラムは前記コンピュータに、検体毎に、検体の感染宿主遺伝子型を示す感染宿主遺伝子型情報と、検体が感染した感染病原体の遺伝子型を示す感染病原体遺伝子型情報と、前記感染病原体に対する治療の有効性を示す有効性情報が与えられた検体データを複数の検体について読み出させ、前記複数の検体についての検体データに基づき、前記感染宿主遺伝子型情報及び前記感染病原体遺伝子型情報の少なくとも一方と前記治療の有効性との相関度を示す重み係数を、遺伝子型毎に算出させ、前記重み係数を、遺伝子型毎に記憶装置に格納させることを特徴とする有効性判定データベース生成プログラムが提供される。
【０００７】
また、本発明の別の観点によれば、１あるいは複数の感染宿主遺伝子型と、治療の有効性との相関度を示す感染宿主重み係数とを記録した第１のフィールドと、１あるいは複数の感染病原体遺伝子型と、前記治療の有効性との相関度を示す感染病原体重み係数とを記録した第２のフィールドとを有することを特徴とする有効性判定データのデータ構造が提供される。
【０００８】
また、本発明のさらに別の観点によれば、コンピュータに、治療の有効性を判定する処理を実行させる有効性判定プログラムであって、この有効性判定プログラムは前記コンピュータに、１あるいは複数の感染宿主遺伝子型と治療の有効性との相関度を示す感染宿主重み係数と、１あるいは複数の感染病原体遺伝子型と前記治療の有効性との相関度を示す感染病原体重み係数とを読み出させ、被検体の感染宿主遺伝子型情報と、前記感染宿主重み係数を遺伝子型毎に乗算して感染宿主有効性予測個別値を算出させ、被検体の感染病原体遺伝子型情報と、前記感染病原体重み係数を遺伝子型毎に乗算して感染病原体有効性予測個別値を算出させ、前記感染宿主有効性予測個別値の各々及び感染病原体有効性予測個別値の各々を加算して有効性予測加算値を算出させ、前記有効性予測加算値あるいは前記有効性予測加算値に所定の数値を加算した値に基づき治療の有効性を判定し、前記判定結果を出力させることを特徴とする有効性判定プログラムが提供される。
【０００９】
また、プログラムに係る本発明は、そのプログラムを実行するためのコンピュータにより構成される装置、そのプログラムによりコンピュータで実行される手順からなる方法、そのプログラムを記録した記録媒体の発明としても成立する。
【発明の効果】
【００１０】
本発明によれば、ウィルス側や患者側の遺伝子に応じた治療の有効性を予測することができる。
【発明の実施の形態】
【００１１】
以下、図面を参照しながら本発明の実施形態を説明する。
【００１２】
なお、本実施形態では、本発明をＣ型肝炎ウィルスに対するインターフェロンを用いた治療法の有効性の予測に適用する例として説明する。
【００１３】
（第１実施形態）
図１は本発明の第１実施形態に係る治療法有効性予測装置の全体構成を示す図である。
【００１４】
図１に示すように、治療法有効性予測装置１００は、コンピュータ１０と、このコンピュータ１０と通信ネットワーク１１を介して接続された検体統計データベース１２から構成される。
【００１５】
コンピュータ１０は、入力装置１と、この入力装置１に接続された処理装置２と、この処理装置２に接続された有効性判定データベース３と、処理装置２に接続された出力装置４及び記憶媒体読取装置５から構成される。このコンピュータ１０は例えばパーソナルコンピュータにより実現される。コンピュータ１０は通信インタフェース（不図示）を介して通信インタフェースに接続される。コンピュータ１０は、この通信インタフェースを介して通信ネットワーク１１との間でデータを送受信することができる。
【００１６】
入力装置１は、例えばキーボードやマウスなどにより実現される。処理装置２は、ＣＰＵなど、一般的なコンピュータの演算処理を実現するハードウェアにより実現される。有効性判定データベース３は、磁気ディスク、光学式ディスクなどにより実現される。出力装置４は、例えばディスプレイやプリンタなどにより実現される。記憶媒体読取装置５は、例えば磁気ディスク読取装置、ＣＤ−ＲＯＭ読取装置、ＤＶＤ読取装置などにより実現される。
【００１７】
入力装置１は、治療が施される個体である患者の遺伝子型に関する情報と、その患者が感染しているウィルスの遺伝子型に関する情報や、その他処理装置２における処理に必要な各種データを入力するための装置である。
【００１８】
処理装置２は、検体統計データベース１２に基づく有効性判定データベース３の生成処理や、入力装置１から入力されたデータ及び有効性判定データベース３に基づく有効性判定処理を実行する。
【００１９】
処理装置２は、例えば処理装置に接続された記憶媒体読取装置５からＣＤ−ＲＯＭやＤＶＤ、磁気ディスクなどの記憶媒体６に記憶されたデータベース生成プログラムや有効性判定プログラムを読み取り、それらプログラムを実行する。これにより、処理装置２が前述した生成処理を行うデータベース生成手段２１や有効性判定処理を行う有効性判定手段２２として機能する。もちろん、処理装置２内蔵のメモリや記憶装置などからこれらプログラムを読み取り実行されるようにしてもよい。また、処理装置２は、これらデータベース生成手段２１や有効性判定手段２２による処理以外の他の処理も、前述のプログラムに基づき実行する。
【００２０】
データベース生成手段２１は、著効率算出処理、オッズ算出処理、χ２値算出処理、重み係数算出処理などを実行する。有効性判定手段２２は、遺伝子型情報入力処理、有効性判定処理などを実行する。
【００２１】
出力装置４は、処理装置２で処理された結果を表示し、あるいは処理結果を印字出力する。
【００２２】
検体統計データベース（ＤＢ）１２の構成の一例を図２に示す。図２に示すように、検体番号、複数のヒト遺伝子型識別情報及び有効性情報の各々がデータベースの各フィールドに検体毎に関連づけて格納されている。図２に示す検体統計データベース１２は、ウィルス遺伝子型識別情報で識別されるウィルス遺伝子型毎に与えられている。具体的には、図２の例では、インターフェロン治療に関連するＣ型肝炎ウィルスの遺伝子型を識別するウィルス遺伝子型として、ＨＣＶ−１ｂ、ＨＣＶ−２ａ、ＨＣＶ−２ｂなどが挙げられる。
【００２３】
検体番号は、検体としての患者を識別するための検体識別情報として用いられる番号である。ヒト遺伝子型識別情報は、患者の遺伝子型を識別する情報である。有効性情報は、その患者への治療の有効性（薬剤の投与であれば薬剤の有効性）を示す情報である。以下の実施形態では、既にヒト遺伝子型情報、ウィルス遺伝子型情報、薬剤投与などの治療を施した結果としての著効／非著効の別が分かっている対象（患者）を検体と呼び、以下の実施形態の発明による有効性判定の対象（患者）を被検体と呼ぶ。
【００２４】
図２には、Ｃ型肝炎のインターフェロンを用いた治療に関連するヒト遺伝子型識別情報として、ＭｘＡ−８８、ＭｘＡ−１２３、ＭＢＬ、ＬＭＰ７、ＩＦＮＡＲ１（ＧＴ）ｎ、ＩＦＮＡＲ１  Ｃ／ＴなどのＳＮＰｓに関する情報が挙げられている。ＩＦＮＡＲ１（ＧＴ）ｎ及びＩＦＮＡＲ１  Ｃ／Ｔは、インターフェロン受容体遺伝子の遺伝子型を識別する情報である。有効性情報は、「著効」あるいは「非著効」のいずれかで示される。より具体的には、図２の例では、インターフェロン療法、すなわちＣ型肝炎の患者にインターフェロンを投与する治療の有効性が「著効」あるいは「非著効」で示されている。
【００２５】
ヒト遺伝子型情報は、それぞれ２つあるいはそれ以上の情報で区別される。ウィルス遺伝子型情報も同様に、２つあるいはそれ以上の情報で区別される。例えばヒト遺伝子型のＭｘＡ−８８は、Ｇ／Ｔ、Ｔ／Ｔ及びＧ／Ｇの３種類の対立遺伝子型のいずれかに分類される。ＭｘＡ−１２３及びＬＭＰ７は、Ｃ／Ａ、Ａ／Ａ及びＣ／Ｃの３種類の対立遺伝子型のいずれかに分類される。ＭＢＬは、ＹＡ及びＸＢの２種類の対立遺伝子型のいずれかに分類される。ＩＦＮＡＲ１（ＧＴ）ｎは、５／５、５／１４及びそれ以外（others）の３種類の対立遺伝子型のいずれかに分類される。ＩＦＮＡＲ１  Ｃ／Ｔは、Ｃ／Ｔ、Ｔ／Ｔ及びＣ／Ｃの３種類の対立遺伝子型のいずれかに分類される。なお、対立遺伝子型による分類の他に対立ハプロタイプにより分類してもよい。
【００２６】
「遺伝子型」とは、全対立遺伝子、あるいは、注目している遺伝子座の遺伝子の存在状態を指す。「遺伝子座」とは、染色体、もしくは遺伝子地図上での遺伝子の位置を指す。例えば、互いに対立遺伝子である遺伝子の場合は同一の遺伝子座にある。「対立遺伝子」とは、相同染色体の相同な場所に位置し、機能的にも相同な遺伝子をいう。
【００２７】
有効性判定データベース３の構成の一例を図３に示す。図３に示すように、複数のヒト遺伝子型重み係数の各々がデータベースの各フィールドに格納されている。図３に示すヒト遺伝子型重み係数は、ウィルス遺伝子型としてＨＣＶ−１ｂに対して与えられた例を示している。このようなヒト遺伝子型重み係数がＨＣＶ−２ａ／ｂなどの他のウィルス遺伝子型に対しても同様に与えられている。以下の実施形態で、重み係数とは、治療の有効性に対してある遺伝子型が与える影響、すなわち治療の有効性に対する相関度を数値化したものである。また、図３の例では、重み係数に符号が乗算された値として示されている。
【００２８】
次に、データベース生成手段２１による有効性判定データベース３の作成方法について図４のフローチャートに沿って説明する。
【００２９】
まず、検体統計データベース１２のヒト遺伝子型の各々に対する治療の著効率を算出する著効率算出処理を実行する（ｓ１）。具体的には、あるヒト遺伝子型に着目し、その遺伝子型の対立遺伝子型あるいは対立ハプロタイプの各々について、著効率を算出する。
【００３０】
例えばＭｘＡ−８８という遺伝子型に着目すると、その対立遺伝子型はＧ／Ｔ、Ｔ／Ｔ及びＧ／Ｇが存在する。この３種類の対立遺伝子型をＧ／Ｔ及びＴ／Ｔの第１の型と、Ｇ／Ｇの第２の型に類型化した場合、第１の型に属する全検体データのうち、著効であった検体データの割合を算出する。第１の型に属する検体のうちＸａ個が著効（ＳＲ）でＹａ個が非著効（ＮＲ）であった場合には、第１の型の著効率（ＳＲrate１）は、
ＳＲrate１＝１００×Ｘａ／（Ｘａ＋Ｙａ）（％）
となる。同様に、第２の型に属する検体のうちＸｂ個が著効（ＳＲ）でＹｂ個が非著効（ＮＲ）であった場合には、第２の型の著効率（ＳＲrate２）は、
ＳＲrate２＝１００×Ｘｂ／（Ｘｂ＋Ｙｂ）（％）
により求められる。このように、検体統計データベース１２の対立遺伝子型が２つ以上の型に分類されている場合、２つの型１，２に分類し、その型毎に著効率ＳＲrate１及びＳＲrate２を算出する。２つの型への分類は、治療の有効性に対する相関度が互いに差が生じるように分類するのが望ましい。より具体的には、２つの型の各々に対する著効率の差が大きくなるように分類するのが望ましい。以下では、この２つの型をそれぞれ「第１の型」及び「第２の型」と呼ぶ。
【００３１】
なお、この２つの型への分類を自動化してもよい。自動化する場合には、遺伝子型の各々を２つの型に分類したと仮定した場合に得られる著効率をすべての分類の組合せについて算出し、互いの著効率が最も大きくなるような分類を選択する。例えばＡ，Ｂ，Ｃ，Ｄ４つの遺伝子型を２つの型に分類する場合、型（Ａ，Ｂ）と型（Ｃ，Ｄ）の分類、型（Ａ，Ｃ）と型（Ｂ，Ｄ）の分類、型（Ａ，Ｄ）と型（Ｂ，Ｃ）の分類、型（Ａ）と型（Ｂ，Ｃ，Ｄ）の分類、型（Ｂ）と型（Ａ，Ｃ，Ｄ）の分類、型（Ｃ）と型（Ａ，Ｂ，Ｄ）の分類、型（Ｄ）と型（Ａ，Ｂ，Ｃ）の分類という７種類の分類が考えられる。この７種類の分類のそれぞれについて、２つの型の著効率及びその差を算出する。その著効率の差が最大となるような分類を選択する。
【００３２】
以上はＭｘＡ−８８を例に説明したが、他の遺伝子型、すなわちＭｘＡ−１２３、ＭＢＬ、ＬＭＰ７、ＩＦＮＡＲ１（ＧＴ）ｎ、ＩＦＮＡＲ１  Ｃ／Ｔなどについても同様の２つの型への分類、各型に対する著効率の算出を実行する。
【００３３】
ＭｘＡ−１２３の場合、Ｃ／Ａ及びＡ／Ａからなる型と、Ｃ／Ｃからなる型の２つの型に分類される。ＭＢＬの場合、ＹＡの型とＸＢの型の２つの型に分類される。ＬＭＰ７は、Ｃ／Ａ及びＡ／Ａからなる型と、Ｃ／Ｃからなる型の２つの型に分類される。ＩＦＮＡＲ１（ＧＴ）ｎは、５／５及び５／１４からなる型と、それ以外の型の２つの型に分類される。ＩＦＮＡＲ１  Ｃ／Ｔは、Ｃ／Ｔ及びＴ／Ｔからなる型と、Ｃ／Ｃからなる型の２つの型に分類される。
【００３４】
得られた著効率は、その型毎に図示しないデータベース（有効性判定データベース３でもよい）に格納される。
【００３５】
次に、オッズ算出処理を実行する（ｓ２）。オッズＯｋとは、２つの型の著効率のうち、大きい方の著効率を小さな方の著効率で除した数値であり、以下の式で示される。なお、ｋは、遺伝子型毎に付与される識別番号を示す。
【００３６】
Ｏｋ＝ＳＲrate１／ＳＲrate２      （ＳＲrate１＞ＳＲrate２の場合）
Ｏｋ＝ＳＲrate２／ＳＲrate１      （ＳＲrate１＜ＳＲrate２の場合）
ＳＲrate１＝ＳＲrate２の場合は、上記２式のいずれで定義してもよい。
【００３７】
以上の式により、すべての遺伝子型についてオッズＯｋを算出し、得られたオッズＯｋは図示しないデータベース（有効性判定データベース３でもよい）に格納される。
【００３８】
次に、χ２値算出処理を実行する（ｓ３）。以下、ｃｋ＝χ２と定義する。χ２値ｃｋは以下の式により算出される。
【００３９】
ｃｋ＝χ２＝（ｎ−１）ｓ２／σ２
ここで、確率変数はＸａ，Ｙａ，Ｘｂ，Ｙｂ、ｎは標本の大きさ、σ２は分散で、ｓ２は標本分散である。より具体的には、Ｘａ，Ｙａ，Ｘｂ，Ｙｂの各数値をカイ二乗検定の２×２分割表に代入することによりカイ二乗値χ２を得ることができる。
【００４０】
次に、オッズＯｋ及びχ２値に基づき、重み係数Ｓｋを算出する重み係数算出処理を実行する（ｓ４）。重み係数Ｓｋは、以下の式により算出される。
【００４１】
Ｓｋ＝（Ｏｋ−１）ｃｋ／２０
次に、得られた重み係数Ｓｋに対する２つの型への符号割り当て処理を実行する（ｓ５）。重み係数Ｓｋは正の値をとる。正の符号＋は、著効率（ＳＲrate）の高い型に割り当てられ、負の符号は、著効率（ＳＲrate）の低い型に割り当てられる。具体的には、各遺伝子型における２つの型の著効率を比較し、大きな著効率を有する型に正の符号を割り当て、小さな著効率を有する型に負の符号を割り当てる。
【００４２】
このように符号割り当て処理を各遺伝子型について実行し、得られた重み係数Ｓｋ、−Ｓｋを有効性判定データベース３に格納する。
【００４３】
以上により、有効性判定データベース３の作成が終了する。
【００４４】
図５は、有効性判定データベース３作成に際し算出された著効率、オッズ、χ２値ｃｋ、重み係数Ｓｋなどを表にしたものである。また、図５には、著効とされた検体と非著効とされた検体の各々の標本数も参考のために示している。重み係数Ｓｋの欄には、各型に対して割り当てられる符号に重み係数Ｓｋを乗算した値が示されている。
【００４５】
次に、有効性判定手段２２による有効性判定データベース３を用いた有効性判定処理について図６のフローチャートに沿って説明する。
【００４６】
まず、被検体としての患者から、ヒト遺伝子型情報及びウィルス遺伝子型情報を得る（ｓ６１）。具体的には、個体に由来するサンプルを取得し、得られたサンプルをＤＮＡチップのセル内に注入し、セル内に固定化されたＤＮＡプローブとのハイブリダイゼーション反応を行わせる。ハイブリダイゼーション反応の後、バッファや挿入剤をセル内に充填し、セル内に固定された電極からの電気化学信号を検出する。この電気化学信号を信号処理することにより、ヒト遺伝子型情報及びウィルス遺伝子型情報を得ることができる。もちろん、他の手法によりこの被検体からの情報を取得してもよい。
【００４７】
次に、得られた被検体のヒト遺伝子型情報及びウィルス遺伝子型情報からなる被検体データを入力装置１を用いて入力する（ｓ６２）。なお、入力装置１の代わりに、ＤＮＡチップを信号線により接続し、ＤＮＡチップから信号線を介してヒト遺伝子型情報及びウィルス遺伝子型情報を入力してもよい。
【００４８】
処理装置２は、入力された遺伝子型情報及びウィルス遺伝子型情報に基づき有効性の判定を行う（ｓ６３）。
【００４９】
具体的には、有効性判定プログラムが記憶媒体６から記憶媒体読取装置５により読み取られ、これにより処理装置２は有効性判定手段２２として機能する。有効性判定手段２２は、まず入力されたウィルス遺伝子型情報に対応するデータテーブルを有効性判定データベース３から読み出す（ｓ６３ａ）。そして、入力された遺伝子型に対応付けられた符号と重み係数を乗算することにより各々の遺伝子型についてヒト有効性予測個別値を得る（ｓ６３ｂ）。そして、得られた各々の遺伝子型についての有効性予測個別値を加算し、有効性予測加算値を算出する（ｓ６３ｃ）。
【００５０】
この有効性予測加算値の算出概念図を図７に示す。図７に示すように、入力された被検体データがＭｘＡ−８８、Ｍｘａ−１２３、ＭＢＬ、ＬＭＰ７、ＩＦＮＡＲ１（ＧＴ）ｎ、ＩＦＮＡＲ１  Ｃ／Ｔの順にＧ／Ｔ、Ｃ／Ｃ、ＹＡ、Ａ／Ａ、５／５、Ｔ／Ｔである場合、各々についてＭｘａ−８８、ＭＢＬ、ＬＭＰ７、ＩＦＮＡＲ１（ＧＴ）ｎ、ＩＦＮＡＲ１  Ｃ／Ｔについては重み係数Ｓｋに＋１を乗算した値が割り当てられ、Ｍｘａ−１２３については重み係数Ｓｋに−１が乗算した値が割り当てられる。その結果、重み係数は図７の通りとなる。この重み係数（有効性予測個別値）を加算すると、有効性個別加算値は（＋０．４１）＋（−０．４０）＋（＋０．２８）＋（＋０．０２）＋（＋０．２８）＋（＋０．００）＝０．５９となる。
【００５１】
次に、処理装置２は、出力装置４に有効性判定結果として有効性予測加算値を出力する（ｓ６４）。主治医あるいは患者は、出力された有効性予測加算値に基づき、治療の有効性を予測し、その治療方法を選択するか否かを判断する際の参考に供することができる。一般に有効性予測加算値が高いほど、治療の有効性が高い。
【００５２】
なお、有効性の判定は、有効性予測加算値の算出に限定されない。例えば、有効性予測加算値に対応する著効率を算出してもよい。
【００５３】
図８は、有効性判定結果の出力例を示す図である。図８に示すように、検体である患者の識別ＩＤ、主治医、その患者の遺伝子型に関する情報、有効性予測加算値、ウィルス遺伝子型と、ウィルス遺伝子型毎の有効性予測加算値に対する著効率が示されている。また、有効性予測加算値を所定の範囲毎に分類し、各類型毎に著効率が算出されている。そして、その著効率を所定のしきい値に基づき治療が有効か否かの目安を示す有効性予測境界を定め（例えば図８では５０％以上と５０％未満）、その境界線がいずれかの類型の間に線引きされて示されている。
【００５４】
有効性予測加算値の範囲毎の著効率の算出は簡単に行える。具体的には、すべての検体データについて有効性予測加算値を算出する。そして、得られた有効性予測加算値を、所定の範囲毎に分類する。図８の例では、０〜１．００、１．００〜２．００、２．００〜３．００…というように、分類している。これは、検体統計データベース１２に格納された検体データの有効性情報に基づき容易に分類可能である。そして、その分類された検体データにおけるすべての標本数に対する著効の標本数の割合を算出することにより、範囲別の著効率が算出される。この範囲別の著効率を有効性予測加算値とともに出力することにより、ウィルスや宿主である患者のヒトの遺伝子に応じた治療の有効性を確率により簡単に把握することができる。また、宿主はヒトに限らず、他の生物でもよい。
【００５５】
以上により、有効性判定データベース３の生成及びそれを用いた有効性判定は終了する。
【００５６】
検体統計データベース１２が更新される毎に上述した（ｓ１）〜（ｓ５）に示す有効性判定データベース３の生成処理を行うことにより、有効性判定データベース３を更新することができる。
【００５７】
また、検体統計データベース１２以外の検体データをさらに入力して有効性判定データベース３を更新することもできる。
【００５８】
このように、検体の増加に伴い有効性判定データベース３をその都度更新することにより、最新の統計情報に基づく治療有効性判定が可能となり、また母集団の増加に伴い判定の制度が向上する。
【００５９】
以上説明したように本実施形態によれば、治療の有効性をその治療対象のウィルスの遺伝子型やヒトの遺伝子型に応じて確率により簡便に判断できる。
【００６０】
（第２実施形態）
図９は本発明の第２実施形態に係る治療法有効性予測装置の全体構成を示す図である。
【００６１】
図９に示すように、本実施形態の治療法有効性予測装置２００の基本的な構成は、図１に示した治療法有効性予測装置１００の構成と共通する。共通する構成には同一符号を付し、詳細な説明は省略する。
【００６２】
治療法有効性予測装置２００では、データベース生成手段２１０、有効性判定手段２２０の処理内容、記憶媒体６０に格納されたプログラム、有効性判定データベース３００に格納される有効性判定データベースの構成、検体統計データベース１２０の構成が図１の治療法有効性予測装置１００のそれと異なる。記憶媒体６０には、処理装置２をデータベース生成手段２１０及び有効性判定手段２２０として機能させるためのデータベース生成プログラム及び有効性判定プログラムなどが格納されている。有効性判定データベース３００には、データベース生成手段２１０により生成された有効性判定データベースが格納される。
【００６３】
データベース生成手段２１０は、データ分類処理、分散共分散行列算出処理、平均差行列算出処理、判別係数算出処理、判別点算出処理、判別関数算出処理、著効率、誤判別率算出処理などを実行する。有効性判定手段２２０は、被検体遺伝子情報入力処理、有効性判定データ選択処理、判別関数値算出処理、判別関数値正負判定処理、判定結果出力処理などを実行する。
【００６４】
検体統計データベース１２０の構成の一例を図１０に示す。図１０に示すように、検体番号、複数のヒト遺伝子型識別情報、複数のウィルス遺伝子型識別情報及び有効性情報の各々がデータベースの各フィールドに検体毎に関連づけて格納されている。
【００６５】
検体番号は、検体としての患者を識別するための検体識別番号として用いられる番号である。ヒト遺伝子型識別情報は、患者の遺伝子型を識別する情報である。ウィルス遺伝子型識別情報は、患者が感染したウィルスの遺伝子型を識別する情報である。有効性情報は、その患者への治療の有効性を示す情報である。
【００６６】
ウィルスやヒトの遺伝子は、通常Ａ，Ｔ，Ｇ，Ｃと略記される４つの塩基により定まるもので、検体データに関しても、この塩基の配列で与えられる。これら遺伝子の配列は、部分的に個々人により差があり、この差が生物学的な特徴付けに関わっており、医療応用の立場からはこのような遺伝子配列のわずかな差がその人の薬剤応答性にかかわっていることが分かっている。図１０では、このような配列の差を数量化して示している。
【００６７】
図１０の例では、Ｃ型肝炎のインターフェロンを用いた治療に関連するヒト遺伝子型識別情報として、ＭｘＡ−８８、ＭｘＡ−１２３、ＭＢＬ、ＬＭＰ７などのＳＮＰｓに関する情報が挙げられる。有効性情報は、著効が「２」、非著効が「１」で示されている。より具体的には、図１０の例では、インターフェロン療法、すなわちＣ型肝炎の患者にインターフェロンを投与する治療の有効性が著効あるいは非著効の２値で示されている。
【００６８】
ヒト遺伝子型情報は、それぞれ２以上の情報で区別される。例えばヒト遺伝子型のＭｘＡ−８８は、Ｇ／Ｇ、Ｇ／Ｔ及びＴ／Ｔの３種類の対立遺伝子型の各々に数値「１」、「２」及び「３」が対応付けられている。なお、医学的な知見に基づき、この３つの類型を２つに分類することも可能である。例えば、Ｇ／Ｇ及びＧ／Ｔを「１」に、Ｔ／Ｔを「２」に対応付ける等である。ＭｘＡ−１２３は、Ｃ／Ｃ、Ｃ／Ａ及びＡ／Ａの３種類の対立遺伝子型の各々に「１」、「２」及び「３」が対応付けられている。ＭＢＬは、ＸＢ及びＹＡの２種類の対立遺伝子型の各々に「１」及び「２」が対応付けられている。ＬＭＰ７は、Ｃ／Ｃ、Ｃ／Ａ及びＡ／Ａの３種類の対立遺伝子型の各々に「１」、「２」及び「３」が対応付けられている。
【００６９】
図１０の例では、Ｃ型肝炎のウィルスの遺伝子型情報として、遺伝子型、ＩＳＤＲが挙げられる。遺伝子型は、「１」及び「２」で示される。それぞれ、ＨＣＶ−１ｂ及びＨＣＶ−２ａ／２ｂが対応している。ＩＳＤＲは、単純な塩基の相違で示されるものではないため、Ｃ型肝炎のＩＳＤＲでのアミノ酸の変異の数を用いている。もちろん、ＭｘＡ−８８等と同様に、医学的な見地に基づき２以上の類型に分類してもよい。例えば、変異数３以下に「１」を割り当て、４以上に「２」を割り当てもよい。
【００７０】
なお、この図１０に示した数量化の手法はほんの一例にすぎず、遺伝子型の各類型に他の数値を割り当てたり、遺伝子の存在状態を示す他の類型や数値を用いたりする等、他の数量化の手法を適用してもよい。望ましくは、誤判別率が最小となるような数値の割り当てが選択される。また、検体データを格納した統計データベースとしては、検体統計データベース１２０に示される形式で格納されていない場合もある。この場合、その統計データベースから、上述した類型化、数量化手法により、検体統計データベース１２０を生成すればよい。
【００７１】
有効性判定データベース３００の構成の一例を図１１に示す。図１１に示すように、線形判別関数と、変数と、有効群著効率、有効群誤判別率、無効群著効率、無効群誤判別率の各々がデータベースの各フィールドに対応付けられて格納されている。
【００７２】
線形判別関数は、入力としての検体の遺伝子情報に対する治療法の有効性を判別するための関数である。この線形判別関数は、複数の変数ａ，ｂ，ｃ，ｄの関数として表現されている。この判別関数を一般化した式（１）を以下に示す。
【００７３】
ｆ（ａ，ｂ，ｃ，ｄ）＝ｐ×ａ＋ｑ×ｂ＋ｒ×ｃ＋ｓ×ｄ＋ｔ  …（１）
式（１）で、ｐ，ｑ，ｒ，ｓは、有効性判定データベース３００の生成処理により決定される係数、ｔは定数である。なお、関数ｆは、ａ，ｂ，ｃ，ｄの４つの変数すべてで示される必要はなく、少なくともそのうちの１つの変数で示されていればよい。定数ｔは、治療法が有効か無効かを関数値の正負により判別するために与えられた数値である。判別関数の値が正であると、その治療法が有効であり、その値が負であると、その治療法が有効でないことを示している。これら判別関数の係数ｐ，ｑ，ｒ，ｓ及び定数ｔは、通常の統計学の手法に従って決定される（例えば栗原考次著“データの科学”放送大学教育振興会参照）。
【００７４】
有効群著効率は、線形判別関数により有効と判別される領域に属する検体データに対する著効が対応付けられた検体データの割合を示している。
【００７５】
無効群著効率は、線形判別関数により無効と判別される領域に属する検体データに対する著効が対応付けられた検体データの割合を示している。
【００７６】
有効群誤判別率は、本当は著効なのに、判別関数により誤って無効と判別される確率、無効群誤判別率は、本当は非著効なのに、判別関数により誤って有効と判別される確率を示す数値である。
【００７７】
図１１（ａ）〜（ｃ）は、いずれもＣ型肝炎ウィルスに対するインターフェロン治療の有効性を判別するための有効性判定データを示している。図１１（ｂ）は、Ｃ型肝炎ウィルスのうちの遺伝子型がＨＣＶ−１ｂに該当する検体データに対して生成された有効性判定データを示している。図１１（ｃ）は、ＨＣＶ−２に該当する検体データに対して生成された有効性判定データを示している。図１１（ａ）は、ＨＣＶ−１ｂに該当する検体データとＨＣＶ−２に該当する検体データの全データに対して生成された有効性判定データを示している。
【００７８】
判別関数は、ウィルスの遺伝子型やヒトの遺伝子型の１つの組合せに対して１つ与えられる。従って、遺伝子型の組合せが複数ある場合には、その組合せ毎に複数の判別関数が与えられる。図１１（ａ）に示すように、遺伝子型としてgeno-type（ウィルス遺伝子型）のみを考える場合、ＭｘＡ−１２３のみを考える場合、ＭｘＡ−８８のみを考える場合に対して１つの判別関数が与えられ、geno-type、ＭｘＡ−１２３、ＭｘＡ−８８の組合せに対して１つの判別関数が与えられ、ＭｘＡ−１２３、Ｍｘａ−８８の組合せに対して１つの判別関数が与えられる。同様に、図１１（ｂ）では、遺伝子型の８つの組合せに対してそれぞれ判別関数が与えられる。図１１（ｃ）では、遺伝子型の５つの組合せに対してそれぞれ判別関数が与えられる。
【００７９】
次に、有効性判定データベース３００の作成方法について図１２のフローチャートに沿って説明する。
【００８０】
まず、検体統計データベース１２から検体データを読み出し、各検体を、著効及び非著効の別で２つのデータ群に分類する（ｓ１１１）。分類された各データ群は図１３に示される。図１３（ａ）は著効である検体データのデータ群［Ａ］、図１３（ｂ）は非著効である検体データのデータ群［Ｂ］を示している。遺伝子型情報は、検体のヒト遺伝子型情報及びウィルス遺伝子型情報の図１０に従った数値を示している。なお、被検体のウィルス遺伝子型が予め分かっている場合などには、その被検体のウィルス遺伝子型に関する検体データのみを抽出し、それに基づきデータ群［Ａ］及び［Ｂ］を算出してもよい。
【００８１】
次に、分類されたデータ群［Ａ］について分散共分散行列ＶＡを算出し、分類されたデータ群［Ｂ］について分散共分散行列ＶＢを算出する（ｓ１１２）。
【００８２】
分散共分散行列ＶＡは、データ群［Ａ］におけるヒト遺伝子型情報、あるいはウィルス遺伝子型情報を変数とした場合に、そのうちの２つの変数の間の直線的傾向の度合いを示す共分散を、すべての変数の組合せについて求めた行列である。また、分散共分散行列ＶＢは、データ群［Ｂ］におけるヒト遺伝子型情報、あるいはウィルス遺伝子型情報を変数とした場合に、そのうちの２つの変数の間の直線的傾向の度合いを示す共分散を、すべての変数の組合せについて求めた行列である。
【００８３】
分散共分散行列ＶＡを以下の式（２）で、分散共分散行列ＶＢを以下の式（３）で定義する。
【００８４】
【数１】
2008272510.tif 000003
【００８５】
この場合、分散共分散行列ＶｉｊＡ及びＶｉｊＢは、以下の式（４）及び（５）で算出される。
【００８６】
【数２】
2008272510.tif 000004
【００８７】
次に、データ群［Ａ］及び［Ｂ］全体の分散共分散行列Ｖｔｏｔを算出する（ｓ１１３）。全体の分散共分散行列の行列要素Ｖｉｊｔｏｔは以下の式（６）で示される。
【００８８】
【数３】
2008272510.tif 000005
【００８９】
次に、平均差行列ｄを算出する（ｓ１１４）。平均差行列ｄは以下の式（７）で示される。
【００９０】
【数４】
2008272510.tif 000006
【００９１】
次に、判別係数Ａを算出する（ｓ１１５）。判別係数Ａは、以下の式（８）で示され、式（９）で示される行列により表現される。
【００９２】
【数５】
2008272510.tif 000007
【００９３】
Ａ＝Ｖｔｏｔ−１×ｄ  …（９）
次に、判別点Ｙ０を算出する（ｓ１１６）。判別点Ｙ０は、以下の式（１０）で示される。
【００９４】
【数６】
2008272510.tif 000008
【００９５】
（ｓ１１５）及び（ｓ１１６）で得られた判別係数Ａ及び判別点Ｙ０により、判別関数Ｚが決定される（ｓ１１７）。判別関数Ｚは以下の式（１１）で示される。
【００９６】
Ｚ＝ａ１Ｘ１＋ａ２Ｘ２＋…＋ａｐＸｐ−Ｙ０  …（１１）
次に、判別関数Ｚに対して検体データを適用し、有効群著効率、有効群誤判別率、無効群著効率及び無効群誤判別率を算出する（ｓ１１８）。
【００９７】
具体的には、判別関数Ｚにより有効、すなわちＺ＞０と判別された検体のうち、実際に著効であった検体数をＳＲ＋、実際には非著効であった検体数をＮＲ＋とすると、有効群著効率はＳＲ＋×１００／（ＳＲ＋＋ＮＲ＋)で計算できる。また、判別関数Ｚにより無効、すなわちＺ＜０と判別された検体のうち、実際に非著効だった検体数をＮＲ−、予想に反して著効だった検体数をＳＲ−とすると、無効群著効率はＳＲ−×１００／（ＳＲ−＋ＮＲ−)で計算できる。理想的には、有効群著効率は１００で、無効群著効率は０になるべきである。
【００９８】
有効群誤判別率はＳＲ−×１００／（ＳＲ＋＋ＳＲ−)で算出される。無効群誤判別率は、ＮＲ＋×１００／（ＮＲ＋＋ＮＲ−)で算出される。
【００９９】
得られた判別関数Ｚ、有効群著効率、有効群誤判別率、無効群著効率及び無効群誤判別率は、それぞれ関連づけられて有効性判定データとして有効性判定データベース３００に格納される（ｓ１１９）。
【０１００】
なお、以上（ｓ１１２）〜（ｓ１１８）に示される処理は、遺伝子型の組合せに応じて繰り返し実行される。すなわち、（ｓ１１２）で生成する分散共分散行列ＶＡ及びＶＢを、複数の遺伝子型の組合せについて設定する。そして、その設定された分散共分散行列ＶＡ及びＶＢに対して（ｓ１１３）〜（ｓ１１８）で示される処理を実行することにより、遺伝子型の組合せの各々について有効性判定データが得られる。
【０１０１】
次に、有効性判定データベース３００を用いた有効性判定手段２２０による有効性判定処理について図１４のフローチャートに沿って説明する。
【０１０２】
まず、検体としての患者から、ヒト遺伝子型情報及びウィルス遺伝子型情報を得る（ｓ６１）。取得の手法は第１実施形態で述べた通りであるので詳細な説明は省略する。
【０１０３】
次に、得られたヒト遺伝子型情報及びウィルス遺伝子型情報を入力装置１を用いて入力する（ｓ６２）。
【０１０４】
処理装置２は、入力された遺伝子型情報及びウィルス遺伝子型情報に基づき有効性の判定を行う（ｓ１３１）。
【０１０５】
具体的には、有効性判定プログラムが記憶媒体６０から記憶媒体読取装置５により読み取られ、これにより処理装置２は有効性判定手段２２０として機能する。有効性判定手段２２０は、入力されたヒト遺伝子型情報及びウィルス遺伝子型情報に対応する変数の組合せを有する有効性判定データを検索し、例えば複数の有効性判定データを選択する（ｓ１３１ａ）。
【０１０６】
なお、この（ｓ１３１ａ）では、被検体データが入力される毎に、前述の（ｓ１１１）〜（ｓ１１９）に示す有効性判定データベース３００生成処理を実行してもよい。その場合、被検体データのウィルス遺伝子型に応じて生成する有効性判定データベース３００を特定のウィルス遺伝子型のものに限定することもできる。
【０１０７】
次に、選択された有効性判定データの変数として定められる遺伝子型に対応するヒト遺伝子型情報を、線形判別関数の変数にそれぞれ代入し、判別関数値Ｚ１を得る（ｓ１３１ｂ）。
【０１０８】
次に、判別関数値Ｚ１の正負の判定を行う（ｓ１３１ｃ）。
【０１０９】
判別関数値Ｚ１＞０であれば、その判別関数Ｚが対応付けられた有効著効率及び有効誤判別率を読み出し、判別結果「有効」とともに、出力装置４に出力させる（ｓ１３２）。
【０１１０】
判別関数値Ｚ１＜０であれば、その判別関数Ｚが対応付けられた無効著効率及び無効誤判別率を読み出し、判別結果「無効」とともに、出力装置４に出力する（ｓ１３２）。
【０１１１】
ＨＣＶ−１ｂに感染した患者の被検体データとして、ＩＳＤＲ＝６、ＭＢＬ＝２、ＭｘＡ−１２３＝３及びＭｘａ−８８＝２が与えられている場合の判定手法を具体例として説明する。図１１の９番の判別関数Ｚ＝０．７６ａ＋１．７１ｂ−０．０４ｃ＋０．９２ｄ−６．３４に着目する。（ｓ１３１ｂ）で、被検体データを判別関数Ｚに代入すると、有効性予測加算値として、判別関数値Ｚ１＝（０．７６×６）＋（１．７１×２）−（０．０４×３）＋（０．９２×２）−６．３４＝３．３６となる。０．７６×６はウィルス有効性予測個別値、１．７１×２、−（０．０４×３）、０．９２×２はヒト有効性予測個別値である。Ｚ１＝３．３６＞０であるため、判定結果「有効」とともに、この９番の判別関数に関連づけられた有効群著効率５２％、有効群誤判別率１８％、無効群著効率４％、無効群誤判別率１６％が読み出されて出力装置４に出力される。
【０１１２】
この判定結果は、被検体に対する治療が有効と判定されたことを示している。さらに、過去にこの判別式で有効と判別されたデータ群に着目すると、判別どおり実際に著効であった割合は５２％であり、実際には著効であるのに誤って無効と判別された割合が１８％であることを示している。また、過去に判別式で無効と判別されたデータ群に着目すると、判別結果に反して著効だった割合は４％であり、実際には非著効であるのに誤って有効と判別された割合は１６％であることを示している。この判定結果により、医者は適切な治癒手段を選択することができる。
【０１１３】
なお、判別結果の出力の前に、読み出された有効誤判別率あるいは無効誤判別率が所定のしきい値以下（例えば２０％以下）か否かを判定し、所定のしきい値以下の判別結果とそれに関連するデータを出力するようにしてもよい。また、読み出された有効誤判別率あるいは無効誤判別率の低い順にデータを出力するようにしてもよい。これにより、誤判別率が高く、信頼性のあまり高く無いデータに基づき有効性を判断するおそれが低減し、誤判別率の低い信頼性の高いデータに基づく有効性予測が可能となる。
【０１１４】
有効性判定データベース３００の更新は、第１実施形態の有効性判定データベース３の更新処理と同様に行うことができる。
【０１１５】
以上説明したように本実施形態によれば、治療の有効性をその治療対象のウィルスの遺伝子型やヒトの遺伝子型に応じて有効／無効の別とその判定結果の誤判別率により、簡便に判断することができる。
【０１１６】
また、このような仕組みにより、薬剤有効性と遺伝子レベルでの関係がさらに明らかになり、かつ検体の母集団が増大することによってデータベースは更新され、より高精度かつ高信頼性のある予測システムを構築することが可能となる。
【０１１７】
この第２実施形態では、判別関数Ｚの係数は、通常の統計学の手法に従って決定され、ヒト遺伝子型情報同士、ウィルス遺伝子型情報同士、あるいはヒト遺伝子型情報及びウィルス遺伝子型情報同士に先験的な相関はないとした。すなわち、１のヒト遺伝子型情報と治療の有効性との相関、あるいは１のウィルス遺伝子型情報と治療の有効性との相関に基づき、係数を決定した。これにより、判別関数の算出に際して計算量の少ない線形判別関数が求まる。これは、事前にデータ群の相関が分かっていないときには、それぞれが正規分布することを前提としたものである。
【０１１８】
しかしながら、本発明はこれに限定されるものではない。例えば、あるヒト遺伝子型情報と別のヒト遺伝子型情報の間の相関、あるヒト遺伝子型情報とウィルス遺伝子型情報の間の相関、あるいはあるウィルス遺伝子型情報と別のウィルス遺伝子型情報の相関が事前に分かっている場合のように、明らかにデータ因子間に相関がある場合には、独立データ因子への縮約作業や非線形判別関数による判定作業を行えばよい。より具体的には、それらの相関に対する係数を改めて求め、その係数とその相関に対応する変数の乗じた項を判別関数Ｚに付加すればよい。従って、第２実施形態で示した線形判別関数のみならず、二次関数を用いた判別関数や、多次元解析の手法であるクラスタ分析法を取り入れてもよい。ここで採用するデータ解析的な手法がアルゴリズムとして既知のものであるならが、それは本発明の範囲内にあることはもちろんである。また、この有効性予測システムでは、判別関数はシステムの中に予め与えられているものである。本発明者らは、ここでは母集団の大きさをある程度変化させて、それが判別関数の係数にどの程度影響するかを調べたが、数％の影響しかないことを確認した。従って、ある程度の大きさの母集団を用いれば、汎用的な判別関数を構築することが可能である。しかし、検体データの母集団に結果は依存し、また人種や地域、生活習慣による相違は当然あり得る。本予測システムは、臨床用遺伝子検査装置への実装が可能であるが、このような装置に通信手段を装備させ、新たな検体データを母集団として決定した判別関数のセットを通信手段を介して遺伝し検査装置に送り込むということも可能である。
【０１１９】
本実施形態で述べた判別関数はあくまである検体データを用いた判別関数であって、その値に普遍性は存在しないが、ここで述べた方法によって、それぞれの診断に適した、また常に最新の情報に基づいた判別関数を提供することができる。また、このような判別関数及び誤判別率のデータベースを通信によらず、磁気媒体などの記録媒体によって配布することも可能である。
【０１２０】
本発明は上記第１，第２実施形態に限定されるものではない。
【０１２１】
検体統計データベース１２は、コンピュータ１０と通信ネットワーク１１を介して接続される形態を示したが、これに限定されず、コンピュータ１０に通信ネットワーク１１を介さず、コンピュータ１０に備えられた記憶装置に格納されていてもよい。
【０１２２】
また、上記実施形態では、Ｃ型肝炎に対してインターフェロンを投与して治療を行う方法が有効か否かの予測に本発明を適用する例を示したが、これに限定されるものではない。治療法は、インターフェロンの投与のみならず、他の治療法も該当する。また、治療法は、薬物の投与を伴わない治療も含まれる。また、治療により治癒すべきウィルスは、Ｃ型肝炎に限らず、他のいかなるウィルスも本発明の予測の対象となり得る。また、ウィルスのみならず、細菌や真菌など遺伝子を有するすべての感染病原体による感染治療に簡単に適用できる。また、感染宿主をヒトとして説明したが、これに限定されず、動物や他の生物などを感染宿主とした場合にも同様に適用可能である。
【０１２３】
また、第１実施形態では、検体統計データベース１２に基づき有効性判定データベース３を生成する例を示したが、これに限定されない。例えば、第２実施形態の検体統計データベース１２０のようなデータベースに基づき有効性判定データベース３を生成してもよい。この場合、検体統計データベース１２０に示されるＩＳＤＲやウィルス遺伝子型についても図５のＭｘａ−８８などのヒト遺伝子型と同様の著効率、オッズ、χ２値、重み係数などを算出する。この場合、第１実施形態と第２実施形態では同一の検体統計データベースから生成され得る。
【０１２４】
ウィルスを類型化する第１実施形態の変形例における検体統計データベース１２ａを図１５に示す。図１５に示すように、ヒト遺伝子型情報に関するＭｘＡ−８８、ＭｘＡ−１２３、ＭＢＬ、ＬＭＰ７、ＩＦＮＡＲ１（ＧＴ）ｎ、ＩＦＮＡＲ１  Ｃ／Ｔなどの遺伝子型のみならず、ウィルス遺伝子型情報に関するＨＣＶ−１ｂ、ＨＣＶ−２ａ、ＨＣＶ−２ｂの別に関する情報が検体毎に対応付けられている。このウィルス遺伝子型情報も例えばＨＣＶ−１ｂを第１の型に、ＨＣＶ−２ａ／２ｂを第２の型に類型化し、ウィルス遺伝子型も含め重み係数を算出した例を図１６に示す。この図１６に示す重み係数と、被検体から得られるヒト遺伝子型情報及びウィルス遺伝子型情報を乗算し、その乗算値をそれぞれ加算することにより、ウィルス遺伝子型も含めた有効性予測加算値を得ることが出来る。なお、図１５や図１６の例では、ウィルス遺伝子型を１つのみ取り上げて重み係数を算出する例を示したが、他のウィルス遺伝子型に関する情報についても同様に重み係数を算出してもよいことはもちろんである。
【０１２５】
また、上述した有効性判定データベース３や３００の生成は、検体統計データベース１２が更新される毎に繰り返し自動的に実行することにより、有効性判定の精度を向上させることができる。また、検体統計データベース１２が更新された場合のみならず、別の経路から検体データが得られた場合には、その検体データを検体統計データベース１２から得られた検体データに加算して有効性判定データベース３や３００の生成を行えばよい。
【０１２６】
また、本実施形態では有効性判定データベースを有するコンピュータ１０で有効性判定を実行する場合を示したが、これに限定されない。例えば、通信ネットワーク１１に接続された端末から被検体の遺伝子型データを送信してコンピュータ１０に判定要求を行い、この要求に応答してコンピュータ１０が上述の有効性判定処理を実行し、判定結果をその端末に送信するようにしてもよい。
【０１２７】
また、第１実施形態では、複数のヒト遺伝子型情報を２種類の型に分類し、その型の各々について分散共分散行列を算出する例を示したが、これに限定されるものではない。例えば３種類以上の型に分類し、その型の各々について第１実施形態と同様に著効率ＳＲrate、オッズＯｋ、重み係数Ｓｋを算出してもよい。第２実施形態でも同様に、（ｓ１１２）３種類以上のデータ群に対して分散共分散行列ＶＡ，ＶＢ，ＶＣ，…を算出し、これらに基づき全体の分散共分散行列Ｖｔｏｔを算出してもよい。
【０１２８】
また、上記第１、第２実施形態における有効／無効の結果の表示は、数値で与えてもよいし、あるいは視覚的に色分けして見分けられる仕組みにしてもよい。
【図面の簡単な説明】
【０１２９】
【図１】本発明の第１実施形態に係る治療法有効性予測装置の全体構成を示す図。
【図２】同実施形態に係る検体統計データベース（ＤＢ）の構成の一例を示す図。
【図３】同実施形態に係る有効性判定データベースの構成の一例を示す図。
【図４】同実施形態に係る有効性判定データベースの作成方法のフローチャートを示す図。
【図５】同実施形態に係る有効性判定データベース作成に際し算出された著効率、オッズ、χ２値、重み係数などを表にして示した図。
【図６】同実施形態に係る有効性判定処理のフローチャートを示す図。
【図７】同実施形態に係る有効性予測加算値の算出概念図。
【図８】同実施形態に係る有効性判定結果の出力例を示す図。
【図９】本発明の第２実施形態に係る治療法有効性予測装置の全体構成を示す図。
【図１０】同実施形態に係る検体統計データベースの構成の一例を示す図。
【図１１】同実施形態に係る有効性判定データベースの構成の一例を示す図。
【図１２】同実施形態に係る有効性判定データベースの生成方法のフローチャートを示す図。
【図１３】同実施形態に係る分類された各データ群のデータの一例を示す図。
【図１４】同実施形態に係る有効性判定データベースを用いた有効性判定処理のフローチャートを示す図。
【図１５】第１実施形態の変形例に関わる検体統計データベースの構成の一例を示す図。
【図１６】第１実施形態の変形例に関わる有効性判定データベースの構成の一例を示す図。
【符号の説明】
【０１３０】
１…入力装置
２…処理装置
３…有効性判定データベース
４…出力装置
５…記憶媒体読取装置
６…記憶媒体
１０…コンピュータ
１１…通信ネットワーク
１２…検体統計データベース
【図１】
2008272510.tif 000009
【図２】
2008272510.tif 000010
【図３】
2008272510.tif 000011
【図４】
2008272510.tif 000012
【図５】
2008272510.tif 000013
【図６】
2008272510.tif 000014
【図７】
2008272510.tif 000015
【図８】
2008272510.tif 000016
【図９】
2008272510.tif 000017
【図１０】
2008272510.tif 000018
【図１１】
2008272510.tif 000019
【図１２】
2008272510.tif 000020
【図１３】
2008272510.tif 000021
【図１４】
2008272510.tif 000022
【図１５】
2008272510.tif 000023
【図１６】
2008272510.tif 000024
【手続補正書】
【提出日】平成20年8月13日(2008.8.13)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
治療の有効性を判定するための有効性判定データベースの作成方法であって、
（１）検体番号、複数のヒト遺伝子型識別情報、複数のウイルス遺伝子型識別情報および有効性情報を関連付けて格納する検体統計データベースから、検体データを読み出すことと、
（２）各検体を著効及び非著効の別で分類して、著効群をデータ群［Ａ］とし、非著効群をデータ群［Ｂ］とすることと、
（３）データ群［Ａ］について分散共分散行列ＶＡを求め、データ群［Ｂ］について分散共分散行列ＶＢを求めることと、
（４）前記（３）から得られた値から、データ群［Ａ］および［Ｂ］の全体の分散共分散行列Ｖｔｏｔを算出することと、
（５）平行差行列を算出することと、
（６）判別係数Ａを算出することと、
（７）判別点Ｙ０を算出することと、
（８）判別関数Ｚを算出することと、
（９）有効群著効率、有効群誤判別率、無効群著効率及び無効群誤判別率を算出することと、
（１０）前記（９）で得られた判別関数Ｚ、有効群著効率、有効群誤判別率、無効群著効率及び無効群誤判別率を、それぞれ関連付けて有効性判定データとして有効性判定データベースに格納することと、
を具備する方法。
【請求項２】
請求項１に記載の方法であって、
（ａ）前記（３）の分散共分散行列ＶＡが以下の式により求められ、
【数１】
2008272510.tif 000025
分散共分散行列ＶＢが以下の式により求められることと、
【数２】
2008272510.tif 000026
（ｂ）前記（４）の全体の分散共分散行列Ｖｔｏｔが以下の式により求められることと、
【数３】
2008272510.tif 000027
（ｃ）前記（５）の平行差行列ｄが以下の式により算出されることと、
【数４】
2008272510.tif 000028
（ｄ）前記（６）の判別係数Ａが以下の式により算出されることと、
【数５】
2008272510.tif 000029
（ｅ）前記（７）の判別点Ｙ０が以下の式により算出されることと、
【数６】
2008272510.tif 000030
（ｆ）前記（８）の判別関数Ｚが以下の式により算出されることと、
Ｚ＝ａ１Ｘ１＋ａ２Ｘ２＋…＋ａｐＸｐ−Ｙ０
（ｇ）前記（９）算出がそれぞれ以下の式により行われることと、
有効群著効率＝ＳＲ＋×１００／（ＳＲ＋＋ＮＲ＋）
（ここで、ＳＲ＋は、Ｚ＞０と判別された検体のうち、実際に著効であった検体数であり、ＮＲ＋は、Ｚ＞０と判定された検体のうち、実際には非著効であった検体数である）
無効群著効率＝ＳＲ−×１００／（ＳＲ−＋ＮＲ−）
（ここで、ＮＲ−は、Ｚ＜０と判別された検体のうち、実際に非著効だった検体数であり、ＳＲ−は、Ｚ＜０と判別された検体のうち、実際には著効であった検体数である）
有効群誤判定率＝ＳＲ−×１００／（ＳＲ＋＋ＳＲ−）
無効群誤判定率＝ＮＲ＋×１００／（ＮＲ＋＋ＮＲ−）
を特徴とする方法。
【請求項３】
コンピュータに、治療の有効性を判定するための有効性判定データベースを生成させる処理を実行させる有効性データベース生成プログラムであって、この有効性判定データベース生成プログラムは、前記コンピュータに、
（１）検体番号、複数のヒト遺伝子型識別情報、複数のウイルス遺伝子型識別情報および有効性情報を関連付けて格納される検体統計データベースから、検体データを読み出させ、
（２）各検体を著効及び非著効の別で分類し、著効群をデータ群［Ａ］とし、非著効群をデータ群［Ｂ］と設定させ、
（３）データ群［Ａ］について分散共分散行列ＶＡを求め、データ群［Ｂ］について分散共分散行列ＶＢを求めさせ、
（４）前記（３）から得られた値から、データ群［Ａ］および［Ｂ］の全体の分散共分散行列Ｖｔｏｔを算出させ、
（５）平行差行列を算出させ、
（６）判別係数Ａを算出させ、
（７）判別点Ｙ０を算出させ、
（８）判別関数Ｚを算出させ、
（９）有効群著効率、有効群誤判別率、無効群著効率及び無効群誤判別率を算出させ、
（１０）前記（９）で得られた判別関数Ｚ、有効群著効率、有効群誤判別率、無効群著効率及び無効群誤判別率を、それぞれ関連付けて有効性判定データとして有効性判定データベースに格納させる、
ことを特徴とする有効性判定データベース生成プログラム。
【請求項４】
請求項３に記載のプログラムであって、
（ａ）前記（３）の分散共分散行列ＶＡが以下の式により求められ、
【数７】
2008272510.tif 000031
分散共分散行列ＶＢが以下の式により求められることと、
【数８】
2008272510.tif 000032
（ｂ）前記（４）の全体の分散共分散行列Ｖｔｏｔが以下の式により求められることと、
【数９】
2008272510.tif 000033
（ｃ）前記（５）の平行差行列ｄが以下の式により算出されることと、
【数１０】
2008272510.tif 000034
（ｄ）前記（６）の判別係数Ａが以下の式により算出されることと、
【数１１】
2008272510.tif 000035
（ｅ）前記（７）の判別点Ｙ０が以下の式により算出されることと、
【数１２】
2008272510.tif 000036
（ｆ）前記（８）の判別関数Ｚが以下の式により算出されることと、
Ｚ＝ａ１Ｘ１＋ａ２Ｘ２＋…＋ａｐＸｐ−Ｙ０
（ｇ）前記（９）算出がそれぞれ以下の式により行われることと、
有効群著効率＝ＳＲ＋×１００／（ＳＲ＋＋ＮＲ＋）
（ここで、ＳＲ＋は、Ｚ＞０と判別された検体のうち、実際に著効であった検体数であり、ＮＲ＋は、Ｚ＞０と判定された検体のうち、実際には非著効であった検体数である）
無効群著効率＝ＳＲ−×１００／（ＳＲ−＋ＮＲ−）
（ここで、ＮＲ−は、Ｚ＜０と判別された検体のうち、実際に非著効だった検体数であり、ＳＲ−は、Ｚ＜０と判別された検体のうち、実際には著効であった検体数である）
有効群誤判定率＝ＳＲ−×１００／（ＳＲ＋＋ＳＲ−）
無効群誤判定率＝ＮＲ＋×１００／（ＮＲ＋＋ＮＲ−）
を特徴とするプログラム。
【請求項５】
治療の有効性を判定するための治療法有効性判定方法であって、
（１）入力されたヒト遺伝子情報及びウイルス遺伝子情報に対応する変数の組合せを有する有効性判定データを、
（Ａ）検体番号、複数のヒト遺伝子型識別情報、複数のウイルス遺伝子型識別情報および有効性情報を関連付けて格納する検体統計データベースから、検体データを読み出すことと、
（Ｂ）各検体を著効及び非著効の別で分類して、著効群をデータ群［Ａ］とし、非著効群をデータ群［Ｂ］とすることと、
（Ｃ）データ群［Ａ］について分散共分散行列ＶＡを求め、データ群「Ｂ」について分散共分散行列ＶＢを求めることと、
（Ｄ）前記（Ｃ）から得られた値から、データ群［Ａ］および［Ｂ］の全体の分散共分散行列Ｖｔｏｔを算出することと、
（Ｅ）平行差行列を算出することと、
（Ｆ）判別係数Ａを算出することと、
（Ｇ）判別点Ｙ０を算出することと、
（Ｈ）判別関数Ｚを算出することと、
（Ｉ）有効群著効率、有効群誤判別率、無効群著効率及び無効群誤判別率を算出することと、
（Ｊ）前記（Ｉ）で得られた判別関数Ｚ、有効群著効率、有効群誤判別率、無効群著効率及び無効群誤判別率を、それぞれ関連付けて有効性判定データとして有効性判定データベースに格納することと、
により得られたデータベースから抽出することと、
（２）抽出された変数を線形判別関数の変数に代入し、判別関数値Ｚ１を算出することと、
（３）判別関数値Ｚ１の正負の判定をすることと、
（４）判別関数Ｚ、有効群著効率、有効群誤判別率、無効群著効率及び無効群誤判別率を求めることと、
（５）前記（４）の結果から治療の有効性を判定することと、
ことを特徴とする方法。
【請求項６】
コンピュータを用いて治療の有効性を判定するための治療法有効性判定方法であって、当該コンピュータに
（１）入力されたヒト遺伝子情報及びウイルス遺伝子情報に対応する変数の組合せを有する有効性判定データを、
（Ａ）検体番号、複数のヒト遺伝子型識別情報、複数のウイルス遺伝子型識別情報および有効性情報を関連付けて格納する検体統計データベースから、検体データを読み出すことと、
（Ｂ）各検体を著効及び非著効の別で分類して、著効群をデータ群［Ａ］とし、非著効群をデータ群［Ｂ］とすることと、
（Ｃ）データ群［Ａ］について分散共分散行列ＶＡを求め、データ群「Ｂ」について分散共分散行列ＶＢを求めることと、
（Ｄ）前記（Ｃ）から得られた値から、データ群［Ａ］および［Ｂ］の全体の分散共分散行列Ｖｔｏｔを算出することと、
（Ｅ）平行差行列を算出することと、
（Ｆ）判別係数Ａを算出することと、
（Ｇ）判別点Ｙ０を算出することと、
（Ｈ）判別関数Ｚを算出することと、
（Ｉ）有効群著効率、有効群誤判別率、無効群著効率及び無効群誤判別率を算出することと、
（Ｊ）前記（Ｉ）で得られた判別関数Ｚ、有効群著効率、有効群誤判別率、無効群著効率及び無効群誤判別率を、それぞれ関連付けて有効性判定データとして有効性判定データベースに格納することと、
により得られたデータベースから抽出させ、
（２）抽出された変数を線形判別関数の変数に代入し、判別関数値Ｚ１を算出させ、
（３）判別関数値Ｚ１の正負の判定をさせ、
（４）判別関数Ｚ、有効群著効率、有効群誤判別率、無効群著効率及び無効群誤判別率を求めさせ、
（５）前記（４）の結果から治療の有効性を判定させ、
（６）判定結果を出力させる、
ことを特徴とする方法。