疾患リスクの提示方法およびそのプログラム
【課題】個人の疾患リスクを提示する方法及びプログラムを提供すること。
【解決手段】疾患リスクを提示する方法は、対象者の検査項目の検査値、対象者が有する遺伝子多型セットの遺伝子情報、並びに、検査項目及び遺伝子多型セットを変数とし、検査項目と少なくとも遺伝子多型セットを含む多項式との積を含む判定式を、記録部から読み出す第１ステップ（Ｓ２、Ｓ４）と、読み出された遺伝子情報を用いて多項式の値を計算する第２ステップ（Ｓ６）と、読み出された検査値及び遺伝子情報を用いて、検査項目を含む全ての項の値を計算する第３ステップ（Ｓ７）と、第２及び第３ステップで計算された値の各々に対応する第１及び第２図形を提示するステップ（Ｓ９）とを含み、判定式が重回帰分析または多重ロジスティック回帰分析によって得られた式である。
【解決手段】疾患リスクを提示する方法は、対象者の検査項目の検査値、対象者が有する遺伝子多型セットの遺伝子情報、並びに、検査項目及び遺伝子多型セットを変数とし、検査項目と少なくとも遺伝子多型セットを含む多項式との積を含む判定式を、記録部から読み出す第１ステップ（Ｓ２、Ｓ４）と、読み出された遺伝子情報を用いて多項式の値を計算する第２ステップ（Ｓ６）と、読み出された検査値及び遺伝子情報を用いて、検査項目を含む全ての項の値を計算する第３ステップ（Ｓ７）と、第２及び第３ステップで計算された値の各々に対応する第１及び第２図形を提示するステップ（Ｓ９）とを含み、判定式が重回帰分析または多重ロジスティック回帰分析によって得られた式である。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、個人の遺伝子情報およびその他の生体検査結果に応じて、その個人が疾患に罹患するリスクを判定して分かり易く提示する方法およびそのプログラムに関する。
【背景技術】
【０００２】
人の遺伝子に関する解析が進み、個人の有する遺伝子と疾患との関連性の解析が精力的に行なわれている。例えば、動脈硬化の発症や進展には生活習慣の違いと遺伝因子の個人差が関連することが知られている。従って、個人の遺伝子検査の結果や生活習慣の情報を利用して、各個人の身体状態の維持・管理（健康増進や罹患の防止、罹患者に関してはその進展の管理など）に役立つ情報を提供することが要望されている。これは、動脈硬化に限らず、心筋梗塞、脳梗塞などの、生活習慣の違いによって発症する可能性や発症後の進展の程度が異なる疾患に関しても同様に要望されている。なお、生活習慣の情報は、血液検査などの検査結果（検査値）や、問診による喫煙や飲酒に関する情報によって得られる情報を意味する。
【０００３】
例えば、下記特許文献１には、糖尿病に起因する動脈硬化、心筋梗塞、腎症及び網膜症について、疾患危険度（疾患の発症しやすさ、進行しやすさを表す指標）を視覚的に分かり易く提示する方法が開示されている。具体的には、被験者の遺伝子情報と環境因子との交互作用を含む疾患危険度の判定式を用いて、レーダーチャート、発症予測グラフ、及びバブルチャートで提示する方法が開示されている。また、特許文献１には、環境因子が変化した場合に疾患危険度の判定式の変化の程度を予測し、その予測値を基準値と比較して、疾患への対処の必要性を判定する方法が開示されている。
【特許文献１】国際公開ＷＯ２００６／１２６６１８
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかし、上記特許文献１に開示された視覚的な提示方法では、疾患毎の危険度の判定式を用いて求めた現在（評価時）の危険度を提示することしかできていなかった。また、疾患への対処の必要性の判定方法では、１つの判定式の変化のみで判断しているので、各環境因子の変化による影響の程度、特に管理可能な環境因子の変化による影響の程度を、分かり易く提示することができなかった。即ち、生活習慣の変化などによって検査値（血圧、総コレステロールなど）が変化した場合の危険度を予測するのに役立ち得る情報を分かり易く提示することはできなかった。
【０００５】
本発明は、上記した問題を解決するためのものであって、その目的は、個人の将来の疾患リスクをわかり易く提示することができる方法およびそのプログラムを提供することにある。
【課題を解決するための手段】
【０００６】
前記課題を解決するための手段は以下のとおりである。なお、ここでは後述する実施の形態に関する図面、符号を引用しているが、これは本発明の理解を容易にするためであり、本発明の限定を意図するものではない。
【０００７】
（１）記録部、演算部および出力部を備えた装置を用いて対象者の疾患リスクを求めて提示する方法であって、
前記演算部が、
前記対象者の検査項目（SBP、Tchなど）の値である検査値、
前記対象者が有する遺伝子多型セット（CSNPi）を表す遺伝子情報、並びに、
前記検査項目および前記遺伝子多型セットを変数とし、前記検査項目（SBPなど）と少なくとも前記遺伝子多型セットを含む多項式（Σ(ai×CSNPi)など）との積（SBP×[a0＋Σ(ai×CSNPi)]など）を含む判定式（例えば式１で表されるｙ）
を、前記記録部から読み出す第１ステップと、
前記演算部が、読み出された前記遺伝子情報を用いて前記多項式の値を計算する第２ステップと、
前記演算部が、前記第２ステップで計算された値に対応する第１図形データ（例えば図３の遺伝子リスクの欄に示された図形データ）を決定する第３ステップと、
前記演算部が、前記検査項目に対応させて、前記検査値および前記第１図形データを、前記出力部から出力する第４ステップとを含み、
前記判定式が、前記検査項目および前記遺伝子多型セットを説明変数として、重回帰分析または多重ロジスティック回帰分析によって得られた式である疾患リスクの提示方法。
【０００８】
（２）前記第４ステップの前に、前記演算部が、読み出された前記検査値および前記遺伝子情報を用いて、前記検査項目（SBPなど）を含む全ての項（SBP×[a0＋Σ(ai×CSNPi)]など）の値を計算する第５ステップと、
前記演算部が、前記第５ステップで計算された値に対応する第２図形データ（顔マーク）を決定する第６ステップとをさらに含み、
前記第４ステップにおいて、前記演算部が、前記検査項目に対応させて前記第２図形データを前記出力部から出力することを特長とする疾患リスクの提示方法。
【０００９】
（３）前記第４ステップの前に、前記演算部が、複数の前記検査値のうち所定範囲を超えている検査値を決定し、決定された該検査値を、前記判定式（ｙ）の値が低下する方向に所定値だけ変化させて目標検査値を求め、該目標検査値を用いて前記判定式から改善度（Δｙ）を求める第７ステップをさらに含み、
前記第４ステップにおいて、前記演算部が、決定された前記検査値に該当する検査項目に対応させて、前記改善度を前記出力部から出力する疾患リスクの提示方法。
【００１０】
（４）前記判定式が、複数の被験者の検査値および遺伝子情報の集合を用いて多重ロジスティック回帰分析によって得られた式であり、
年代毎に分類された前記被験者に関する検査値の平均値（ｔav）が、前記記録部に記録されており、
前記演算部が、前記対象者に関して、前記判定式中の所定の検査項目を含む項の係数値（α）を計算し、該係数値（α）を指数として自然対数の底のべき乗値（ｅα）を計算して前記所定の検査項目に関するオッズ比を求め、前記対象者の年齢が該当する年代に対応する前記平均値を前記記録部から読み出し、該平均値を前記対象者の前記検査値（ｔ）から減算した値（ｔ−ｔav）を指数として、前記オッズ比のべき乗値（β＝（ｅα）(ｔ−ｔav)）を計算する第８ステップと、
前記演算部が、前記検査項目毎の前記べき乗値（β）を、全て乗算して評価値を計算する第９ステップと、
前記演算部が、前記評価値を前記出力部から出力する第１０ステップと
をさらに含む疾患リスクの提示方法。
【００１１】
（５）前記判定式が、複数の被験者の検査値および遺伝子情報の集合を用いて重回帰分析によって得られた式であり、
年代毎に分類された前記被験者に関する前記判定式の値の平均値が、前記記録部に記録されており、
前記演算部が、読み出された前記検査値および前記遺伝子情報を用いて前記判定式の値を計算して評価値とする第１１ステップと、
前記演算部が、前記対象者の年齢が該当する年代に対応する前記平均値を、前記記録部から読み出す第１２ステップと、
前記演算部が、前記評価値および前記平均値を前記出力部から出力する第１３ステップと
をさらに含む疾患リスクの提示方法。
【００１２】
（６）前記疾患が心筋梗塞であり、
前記判定式に含まれる各々の各遺伝子多型セットが、図７に記載された遺伝子多型セットの１つである疾患リスクの提示方法。
【００１３】
（７）前記疾患が脳梗塞であり、
前記判定式に含まれる各々の遺伝子多型セットが、図８に記載された遺伝子多型セットの１つである疾患リスクの提示方法。
【００１４】
（８）前記疾患が動脈硬化であり、
前記判定式に含まれる各々の遺伝子多型セットが、図６に記載された遺伝子多型セットの１つである疾患リスクの提示方法。
【００１５】
（９）記録部および出力部を備えたコンピュータに、対象者の疾患リスクを求めて提示させるプログラムであって、
前記コンピュータに、
前記対象者の検査項目の値である検査値、
前記対象者が有する遺伝子多型セットを表す遺伝子情報、並びに、
前記検査項目および前記遺伝子多型セットを変数とし、前記検査項目と少なくとも前記遺伝子多型セットを含む多項式との積を含む判定式を、前記記録部から読み出す第１機能と、
読み出された前記遺伝子情報を用いて前記多項式の値を計算する第２機能と、
前記第２ステップで計算された値に対応する第１図形データを決定する第３機能と、
前記検査項目に対応させて、前記検査値および前記第１図形データを、前記出力部から出力する第４機能とを実現させ、
前記判定式が、前記検査項目および前記遺伝子多型セットを説明変数として、重回帰分析または多重ロジスティック回帰分析によって得られた式である疾患リスクの提示プログラム。
【発明の効果】
【００１６】
本発明によれば、特定の個人について、遺伝因子および環境因子を総合的に考慮した疾患リスクを、同世代の平均値とともにグラフィカルに提示することができる。
【００１７】
また、特定の疾患について、関連する複数の検査項目の各々に対応させて、遺伝因子を考慮した遺伝因子リスクと、遺伝因子および環境因子を考慮した総合リスクとを提示することができる。
【００１８】
また、遺伝因子リスク、総合リスクの数値が該当するランクをマークで分かり易く提示することができる。
【００１９】
従って、これらの提示に基づいて医師などが、その個人（評価対象者）が将来どのように身体を管理すべきか（生活習慣の改善など）を、わかり易く且つ適切に指導することが容易となる。
【発明を実施するための最良の形態】
【００２０】
以下、まず用語について説明した後、図面を参照して本発明についてより詳細に説明する。
【００２１】
「遺伝子多型」とは、同一集団内において、一つの遺伝子座に２種類以上の対立遺伝子（アレル）が存在する遺伝子の多様性を意味する。具体的には、ある集団において一定の頻度以上で存在する遺伝子の変異を示す。ここでいう遺伝子の変異は、ＲＮＡとして転写される領域に限定されるものではなく、プロモーター、エンハンサー等の制御領域などを含むヒトゲノム上で特定しうるすべてのＤＮＡにおける変異を含むものである。ヒトゲノムＤＮＡの９９．９％は各個人間で共通しており、残る０．１％がこのような多様性の原因となり、特定の疾患に対する感受性、薬物や環境因子に対する反応性の個人差として関与し得る。遺伝子多型があっても表現型に差が出るとは限らない。なお、ＳＮＰ（一塩基多型）も遺伝子多型の一種であるが、本発明が対象とする遺伝子多型はこれに限定されない。
【００２２】
このように、「遺伝子多型」は、遺伝子型（genotype）を考慮したものであって、それを含むものであり、本明細書において「遺伝子多型」とは特定の遺伝子型（genotype）を有する遺伝子多型を意味する。
【００２３】
本明細書で示す遺伝子型（genotype）は、「１」は置換塩基のうち塩基のアルファベット順（Ａ、Ｃ、Ｇ、Ｔ）で前にくる塩基を有する多型のホモを、「２」はヘテロを、「３」は、置換塩基のうち塩基のアルファベット順で後にくる塩基を有する多型のホモを表す。例えば、対象とする遺伝子（遺伝子多型）がXXXX(C3413G)（この表記は、遺伝子XXXXが、その塩基配列の3413番目の塩基がＣまたはＧである遺伝子多型を有することを意味する。）であるとすると、その遺伝子型（genotype）がCCのホモである場合を遺伝子型１、CGまたはGCのヘテロである場合を遺伝子型２、GGのホモである場合を遺伝子型３という。なお、「１２」は前記１と２の遺伝子型（genotype）の両方の遺伝子型、「２３」は前記２と３の遺伝子型（genotype）の両方の遺伝子型を表す。
【００２４】
なお、Aldose redactase遺伝子多型（Aldose redactase(AC)nで表され、ｎはACの繰り返し回数を表し、ｎ＝１３〜３０である）に関しては、上記の規則と異なり、次の規則で遺伝子型を付与している。すなわち、遺伝子の塩基配列中の「AC」の繰り返し回数がｎ回のホモである場合を遺伝子型「１」、「AC」の繰り返し回数がｎ回とｎ回以外とのヘテロである場合を遺伝子型「２」、これら以外の場合、即ち「AC」の繰り返し回数がｎ回以外のホモまたはヘテロである場合を遺伝子型「３」とする。
【００２５】
「遺伝子多型セット」とは、複数の遺伝子多型の組合せをいう。ここで複数の遺伝子多型とは、異なる遺伝子座を有する２種以上の遺伝子多型を意味する。
【００２６】
本明細書における塩基配列（ヌクレオチド配列）、核酸などの略号による表示は、ＩＵＰＡＣ−ＩＵＢの規定〔IUPAc-IUB communication on Biological Nomenclature, Eur. J. Biochem., 138; 9 (1984)〕、「塩基配列又はアミノ酸配列を含む明細書等の作製のためのガイドライン」（特許庁編）及び当該分野における慣用記号に従うものとする。
【００２７】
本明細書中において「遺伝子」は、特に言及しない限り、ヒトゲノムＤＮＡを含む２本鎖ＤＮＡ、及び１本鎖ＤＮＡ（センス鎖）、並びに当該センス鎖と相補的な配列を有する１本鎖ＤＮＡ（アンチセンス鎖）、及びそれらの断片のいずれもが含まれる、また、本明細書で「遺伝子」とは、特に言及しない限り、調節領域、コード領域、エクソン、及びイントロンを区別することなく示すものとする。
【００２８】
また、「動脈硬化」とは「頚動脈肥厚」を意味する。
【００２９】
図１は、本発明の実施の形態に係る疾患リスクの提示方法を実施するための装置（以下、提示装置と記す）を示すブロック図である。本提示装置１は、演算部（以下、ＣＰＵと記す）２と、メモリ３と、コンピュータ用ハードディスクなどの記録部４と、コンピュータ用キーボードやマウスなどの操作部５と、液晶ディスプレイなどの表示部６と、プリンタなどの印刷部７と、入出力のインタフェースである入出部８と、各部間でデータ（制御データを含む）を交換するための内部バス９とを備えている。
【００３０】
記録部４には、疾患リスクの計算に使用する遺伝子多型セット（図６〜８参照）の情報が参照テーブルとして記録されている。ここで、図６〜８には、動脈硬化、心筋梗塞および脳梗塞のそれぞれについて、図９Ａ〜９Ｃに示した１３３個の遺伝子多型の中から、上記特許文献１に開示された方法と同様の方法で決定された遺伝子多型セットを示している。なお、上記特許文献１では、有効な遺伝子多型セットの決定に正および負のｔ値を用いているが、図６〜８に示した遺伝子多型セットは、正のｔ値のみを用いて決定されたものである。
【００３１】
図９Ａ〜９Ｃにおいて、SNP-No.、GENETIC POLYMORPHISM、GENE、SITEはそれぞれ、各遺伝子多型に付した連続番号、遺伝子多型、遺伝子多型の略称、塩基配列上のアレルの位置を表す。図６〜８において、SSNP、Gtype、t-valueはそれぞれ、１つのセットを構成する遺伝子多型の数、遺伝子型、ｔ値を表す。SNP-No.は図９Ａ〜９Ｃに示した番号であり、GENEは図９Ａ〜９Ｃに示したGENE(略称)とSITEとを組み合わせた表記である。
【００３２】
なお、図９Ｃに示したSNP-No.が124の遺伝子は、Aldose redactase遺伝子であり、その遺伝子型は上記したように「AC」の繰り返し回数に応じて付与されている。さらに、図６〜８において、Aldose redactase遺伝子をALR_Z_4、ALR_Z_2、ALR_Z、ALR_Z2、ALR_Z4と表記しており、これらの遺伝子型（Gtype）は、図１０に示した規則で表記されている。例えば、ALR_Z_4に関しては、「AC」の繰り返し回数が２２回のホモである場合に遺伝子型を「１」とし、「AC」の繰り返し回数が２２回と２２回以外（即ち、１３〜２１回及び２３〜３０回の何れか）とのヘテロである場合に遺伝子型を「２」とし、「AC」の繰り返し回数が２２回以外のホモまたはヘテロである場合に遺伝子型を「３」とする。また、例えば、ALR_Z_2に関しては、「AC」の繰り返し回数が２３回のホモである場合に遺伝子型を「１」とし、「AC」の繰り返し回数が２３回以外（即ち、１３〜２２回及び２４〜３０回の何れか）のホモである場合に遺伝子型を「３」とし、「AC」の繰り返し回数が２３回と２３回以外とのヘテロである場合に遺伝子型を「２」とする。ALR_Z、ALR_Z2、ALR_Z4に関しても、繰り返し回数が異なるが同様である。
【００３３】
また、図９Ｂに示したSNP-No.が71のResistin VNTR 遺伝子多型（Resistin(ATG)n repeat)で表わされ、ｎはATGの繰り返し回数を表し、ｎ＝６〜８である）に関しては、一塩基多型と異なり、Aldose redatase(AC)nと同じく、即ち、遺伝子の塩基配列の「ATG」の繰り返し回数がｎ回のホモである場合を遺伝子型「１」、「ATG」の繰り返し回数がｎ回とｎ回以外とのヘテロである場合を遺伝子型「２」、これら以外の場合、即ち「ATG」の繰り返し回数がｎ回以外のホモまたはヘテロである場合を遺伝子型「３」とする。図７では、ｎ＝６のResistin遺伝子多型をresistin_6と表している。
【００３４】
記録部４にはさらに、判定対象者の検査値（血液検査結果および問診結果の情報）及び遺伝子情報が、対象者毎に付した情報（対象者ＩＤ）と対応させて記録されている。また、記録部４には、疾患危険度の判定式ｙが記録されている。判定式ｙは、後述するように検査項目および遺伝子多型セットを変数とする多項式であり、使用するコンピュータプログラム言語に応じた形式で記録されている。
【００３５】
以下、図２に示したフローチャートを用いて提示装置１を用いた、疾患リスクの提示方法について説明する。ここでは、特に断らない限り、ＣＰＵ２が行なう処理として記載する。ＣＰＵ２は、操作部５を介して指示やデータ入力を受け、処理対象データを記録部４からメモリ３に読み出し、メモリ３をワーク領域として使用して各処理を実行し、処理の途中結果や最終結果を、必要に応じて記録部４の所定領域に記録する。ＣＰＵ２は、処理結果を記録部４から読み出して、表示部６の入力信号形式に応じた所定形式のデータ（アナログ又はディジタル）に変換し、入出力部８を介して表示部６に伝送し、表示させる。また、ＣＰＵ２は、処理結果を記録部４から読み出して、印刷部７の入力信号形式に応じた所定形式のデータに変換し、印刷部７に伝送し、例えば紙に印刷させる。
【００３６】
ステップＳ１において、疾患リスクを判定して提示する対象者を指定するために操作部５が操作されて入力される対象者ＩＤを、入出力部８を介してメモリ３に一時記憶する。
【００３７】
ステップＳ２において、ステップＳ１でメモリ３に一時記憶された対象者ＩＤに対応するデータを記録部４からメモリ３に読み出す。このとき読み出されるデータは、対象者の遺伝子情報（遺伝子多型の有無を表す情報）および検査値である。検査値は、例えば、血液検査の結果である血圧（mmHg）、血中総コレステロール値（mg/dl）、血中中性脂肪値（mg/dl）、血糖コントロール（空腹時血糖値（mg/dl）、及びヘモグロビンA1c値（HbA1c）（％））、身体検査によって得られるBMI（Body-mass-index）（＝体重（kg）/（身長（m））2）、問診によって得られる喫煙歴である。ここで、血圧には収縮期血圧（SBP）および拡張期血圧（DBP）が含まれる。また、喫煙歴はこの判定時の喫煙習慣を表す情報であり、判定時に喫煙習慣を有してれば１、有していなければ０が設定される。なお、判定時に喫煙習慣を有していないとしても、過去に喫煙習慣を有していたことがあれば、喫煙歴には１が設定される。
【００３８】
ステップＳ３において、予め指定された複数の疾患のうち、処理対象とする最初の疾患を指定する。ここでは、動脈硬化が指定されたとする。また、その他に心筋梗塞および脳梗塞を対象とし、疾患リスクが判定されるとする。
【００３９】
ステップＳ４において、ステップＳ３で指定された疾患、すなわち動脈硬化に対応する危険度の判定式（回帰式）ｙを記録部４からメモリ３に読み出す。判定式ｙは予め決定されており、上記特許文献１に開示されているように、ｙ＝ΣａiＧi＋ΣｂjＥj＋ΣｄmnＥmＧn＋ｃ で表される。各係数ａi、ｂj、ｄmn、ｃは、所定の母集団を用いて重回帰分析によって疾患毎に決定された値である。例えば、動脈硬化の場合、ＩＭＴ（Intima Media Thickness）を目的変数とする。
【００４０】
上記した検査項目などの環境因子を説明変数として、判定式ｙを書き直せば、次式のように記述することができる。
ｙ＝SBP×[a0＋Σ(ai×CSNPi)]＋Tch×[b0＋Σ(bi×CSNPi)]
＋TG×[c0＋Σ(ci×CSNPi)]＋HbA1c×[d0＋Σ(di×CSNPi)]
＋BMI×[e0＋Σ(ei×CSNPi)]＋SM×[f0＋Σ(fi×CSNPi)]
−HDL×[g0＋Σ(gi×CSNPi)]＋const ・・・・・・・・・（式１）
ここで、SBP、Tch、TG、HbA1c、BMI、SMおよびHDLはそれぞれ、収縮期血圧（mmHg）、血中総コレステロール値（mg/dl）、中性脂肪（mg/dl）、HbA1c値（％）、BMI（kg/m2）、喫煙歴（１又は０）およびHDLコレステロール値（mg/dl）である。CSNPiは、i番目の遺伝子多型セットCSNPの有無に応じて１又は０を取る変数である。ai、bi、ci、di、ei、ｆi、gi（iは０以上の整数）は回帰係数であり、constは定数である。また、Σはiに関する和を求める演算子である。HDLを含む項の符号がマイナスになっているのは、HDLの疾患リスクへの影響が他の検査値とは逆であること、即ちHDLが大きくなると疾患リスクが低減することを明示的に示すためである。なお、各項の係数は、疾患および判定式ｙを求めるときに使用した被験者の集合に応じて異なる。また、疾患によっては、判定式ｙに一部の項が含まれていない場合もある。
【００４１】
上記の記述では、本発明の疾患リスクの判定に用いる７種類の環境因子（SBP、Tch、TG、HbA1c、BMI、SMおよびHDL）を含む項を明示的に記述しているが、判定式ｙには、上記の７種類以外の環境因子（例えば、性別（SEX）、年齢（AGE）、罹病期間など）を含む項や、環境因子を含まない遺伝子多型セットのみの項が含まれていてもよい。その場合、上記の７種類以外の環境因子を含む項（遺伝因子を含まない環境因子のみの項、及び、環境因子と遺伝因子との積である交互作用項を含む）はconstに含める。上記では、constは定数であるとしたが、「定数」とはこのような意味である。本発明では、個人が管理可能な情報を提示することを目的としているので、性別、年齢などの管理不可能な環境因子を含む項は、本発明の処理対象とはしない。
【００４２】
ステップＳ５において、ステップＳ２で読み出された対象者のデータとステップＳ４で読み出された判定式ｙとを用いて、ステップＳ１で指定された疾患に関する評価値を計算し、対象者ＩＤと対応させて記録部４に記録する。即ち、対象者の遺伝子多型セットのデータ（CSNPi）と検査値（SBP、Tch、TG、HbA1c、BMI、SM、HDL）とを判定式ｙに代入し、評価値を計算する。
【００４３】
ステップＳ６において、ステップＳ２で読み出された対象者のデータとステップＳ４で読み出された判定式ｙとを用いて、遺伝因子リスクを計算し、検査項目を特定する情報（検査項目ＩＤ）及び対象者ＩＤと対応させて記録部４に記録する。即ち、検査項目（BP、Tch、TG、HbA1c、BMI、SM）毎に、検査項目を含む項のうち検査項目を除いた部分の値を計算する。従って、この遺伝因子リスクは、対象者が持っている遺伝子に応じて決定される検査項目毎のリスクといえる。
【００４４】
例えば、収縮期血圧に関しては、SBPを含む項であるSBP×[a0＋Σ(ai×CSNPi)]のうち、a0＋Σ(ai×CSNPi)の値を計算し、血中総コレステロール値に関しては、Tchを含む項であるTch×[b0＋Σ(bi×CSNPi)]のうち、b0＋Σ(bi×CSNPi)を計算する。中性脂肪、HbA1c値、BMIおよび喫煙歴に関しても、同様に検査項目を除いた部分の値を計算する。
【００４５】
ステップＳ７において、ステップＳ２で読み出された対象者のデータとステップＳ４で読み出された判定式ｙとを用いて、総合リスクを計算し、検査項目ＩＤ及び対象者ＩＤと対応させて記録部４に記録する。即ち、検査項目（BP、Tch、TG、HbA1c、BMI、SM）毎に、特定の検査項目を含む全ての項、即ち、特定の検査項目と遺伝子多型セットの多項式（この多項式は特定の検査項目にとって係数に相当する。）との積を計算する。この総合リスクは、対象者が持っている遺伝子および検査値の両方を考慮して決定される検査項目毎のリスクといえる。
【００４６】
例えば、収縮期血圧に関しては、SBPを含む全ての項であるSBP×[a0＋Σ(ai×CSNPi)]の値を計算し、血中総コレステロール値に関しては、Tchを含む全ての項であるTch×[b0＋Σ(bi×CSNPi)]を計算する。中性脂肪、HbA1c値、BMIおよび喫煙歴に関しても、同様に計算する。
【００４７】
以上、ステップＳ５〜Ｓ７での処理によって、ステップＳ１で指定された対象者について、動脈硬化に関する評価値、遺伝因子リスクおよび総合リスクが求められる。次に、これらの値を提示する方法を説明する。動脈硬化に関する提示方法の一例を図３に示す。以下、図３を参照して具体的に説明する。
【００４８】
ステップＳ８において、対象者ＩＤに対応する評価値を記録部４から読み出し、表示部６の画面に、例えば図３の上段（左端に「評価」と表示された枠内）に示されているように提示する。具体的には、横長の細い長方形（以下、スケールと記す）を、例えばグラデーションを付して描画する。スケールは、基準尺度とするために、所定範囲のＩＭＴ値を表示したものである。スケールの左端に数値０を表示し、右側に向かって順に数値０．８、１．１、１．４、２．８を表示する。そして、対象者と同世代の人のＩＭＴの平均値に対応するスケール上の位置に矢印を表示する。さらに、評価値に対応するスケール上の位置に、人型の図形を表示し、その下に評価値、例えば１．２を表示する。
【００４９】
ここで、スケールに付す数値は、動脈硬化の程度を表す基準値であるＩＭＴの値であり、これらの値は予め設定されて、記録部４に記録されている。なお、図３に示した値（０．８、１．１、１．４、２．８）は、動脈硬化の程度を表す基準値（ＩＭＴ値）として一般に使用されている値である。
【００５０】
また、図３の例では、評価値と平均値とを比較して、評価を文章で自動的に表示している。図３の上段に表示された文章の下線部は条件によって変化する部分であり、それ以外の部分は固定である。即ち、左側の下線部分には疾患名称を表示し、右側の下線部分には、評価値と平均値との大小関係に応じた表示をする。例えば、評価値が平均値よりも大きい場合、「すすんでいる」と表示する。
【００５１】
ステップＳ９において、対象者ＩＤに対応する検査値、遺伝因子リスク及び総合リスクを記録部４から読み出し、表示部６の画面に、例えば図３の下段（左端に「分析表」と表示された枠内）に示したように提示する。即ち、テーブル形式で表示し、最上段に検査項目を、２段目に検査値を、３段目に遺伝因子リスクを、最下段に総合リスクをそれぞれ表示する。
【００５２】
図３のテーブルの２段目には、記録部４から読み出された検査値を検査項目毎に提示している。このとき、検査値が、検査項目毎に予め設定された所定範囲内にあるか否かの情報を併せて提示する。図３では、所定範囲から外れている検査値を白抜き文字で表示しているが、別の方法として、異なる色で表示してもよい。なお、血圧に関しては収縮期血圧（SBP）値および拡張期血圧（DBP）値を、血糖コントロールに関しては空腹時血糖値およびHbA1c値を表示し、喫煙歴に関しては、読み出された値が０であれば「無」、１であれば「有」を表示する。
【００５３】
図３のテーブルの３段目には、遺伝因子リスクの数値を直接表示せずに、リスクの程度をランク付けして表示するために、数値に応じたマーク（図形）を表示する。即ち、ステップＳ６で検査項目（ｉ）毎に求められた遺伝因子リスクＢiと、予め検査項目（ｉ）毎に決定された各ランク（ｊ）に対応する基準値（Ｒij）とを比較して、Ｒi,j≦Ｂi＜Ｒi,j+1である場合、ランク（ｊ）に対応するマークを記録部４から読み出して表示する。具体的には、判定式ｙを重回帰分析によって求めるときに使用した被験者の集合に関して、各被験者の検査項目毎の遺伝因子リスクを計算し（ステップＳ６参照）、その度数分布において、例えば２５、５０、７５パーセンタイルに該当する遺伝因子リスクの値を、検査項目毎の基準値（Ｒij）とする。例えば、血圧に関しては、判定式ｙ中の収縮期血圧SBPを含む項であるSBP×[a0＋Σ(ai×CSNPi)]のうち、a0＋Σ(ai×CSNPi)の値を、集合の各被験者について計算する。これによって、対象者の検査項目毎の遺伝因子リスクをマークで表示することができる。なお、図３では、一部の欄（中性脂肪、ＢＭＩに対応する欄）が空白になっているが、これは、使用した判定式ｙ中に該当する項が存在しないこと、若しくは、判定式ｙ中に該当する項は存在するが、その対象者が該当する遺伝子多型セットを有していないためにその項が０になったことを意味する。
【００５４】
図３のテーブルの最下段には、総合リスクの数値を直接表示せずに、リスクの程度をランク付けして表示するために、数値に応じたマークを表示する。遺伝因子リスクの提示に関して上記したのと同様に、ステップＳ７で検査項目（ｉ）毎に求められた総合リスクＣiと、予め検査項目（ｉ）毎に決定された各ランク（ｊ）に対応する基準値（Ｔij）とを比較して、Ｔi,j≦Ｃi＜Ｔi,j+1である場合、ランク（ｊ）に対応するマークを記録部４から読み出して表示する。具体的には、判定式ｙを重回帰分析によって求めるときに使用した被験者の集合に関して、各被験者の検査項目毎の遺伝因子リスクを計算し（ステップＳ７参照）、その度数分布において、例えば２５、５０、７５パーセンタイルに該当する遺伝因子リスクの値を、検査項目毎の基準値（Ｔij）とする。例えば、血圧に関しては、判定式ｙ中の収縮期血圧SBPを含む項であるSBP×[a0＋Σ(ai×CSNPi)]の値を、集合の各被験者について計算する。これによって、対象者の検査項目毎の総合リスクをマークで表示することができる。遺伝因子リスクが判定式y中に該当する項が存在しないために空白の欄（図３参照）があれば、それに対応する総合リスクとして表示するマークは、検査値のみから決定する。
【００５５】
以上、ステップＳ９およびＳ１０での処理によって、例えば図３に示したように、動脈硬化に関して対象者の遺伝因子リスク及び総合リスクを、ランクに対応するマークを用いてわかり易く提示することができる。
【００５６】
ステップＳ１０において、予め指定された複数の疾患（動脈硬化、心筋梗塞および脳梗塞）の全てについて処理が完了したか否かを判断し、完了していなければステップＳ１１に移行して別の疾患を指定した後、ステップＳ４に移行する。予め指定された複数の疾患の全てについて処理が完了していればステップＳ１２に移行し、終了の指示があったか否かを判断し、指示が無かった場合、ステップＳ１に移行して新たに対象者の指定を受け付ける。
【００５７】
以上によって、特定の対象者について、動脈硬化、心筋梗塞および脳梗塞の各疾患に関するリスクを計算し、図３と同様に提示することができる。なお、心筋梗塞および脳梗塞に関しても、上記した動脈硬化に関する処理と同様の処理を行えばよいが、動脈硬化の場合と異なる点もあるので、以下にその点について説明する。
【００５８】
心筋梗塞、脳梗塞に関しては、説明変数としては動脈硬化の場合と同様の環境因子（検査項目など）を用いるが、動脈硬化の場合と異なり多重ロジスティック回帰分析によって判定式ｙを求める。最終的に求められた判定式ｙは、上記の式１と同様の多項式の形ではあるが、左辺ｙの値はロジットである。疾患の発症確率をｐとすれば、ロジットはlog[ｐ／(１−ｐ)]で表される。ここで、logは自然対数を意味する。また、ｐ／(１−ｐ)はオッズ（Odds）である。
【００５９】
即ち、所定の被験者の集合に関して、環境因子（検査項目など）と脳梗塞の発症確率ｐから求められるオッズとを用いて、多重ロジスティック回帰分析によって脳梗塞の判定式ｙを求める。また、所定の被験者の集合に関して、環境因子（検査項目など）と心筋梗塞の発症確率ｐから求められるオッズとを用いて、多重ロジスティック回帰分析によって心筋梗塞の判定式ｙを求める。例えば、心筋梗塞に関して図７に示した遺伝子多型セットを使用して、中性脂肪およびＢＭＩの各々と遺伝子多型セットとの積の項を含まない判定式ｙが得られた。また、脳梗塞に関して図８に示した遺伝子多型セットを使用して、中性脂肪と遺伝子多型セットとの積の項を含まない判定式ｙが得られた。
【００６０】
従って、脳梗塞、心筋梗塞に関しては、判定式ｙの左辺がロジットである点で動脈硬化の場合と異なるが、遺伝因子リスクおよび総合リスクを求める処理はステップＳ６およびＳ７の処理と同様であり、遺伝因子リスクおよび総合リスクをマークで提示する処理はステップＳ９の処理と同様である。
【００６１】
しかし、脳梗塞、心筋梗塞に関して、評価値を求める処理はステップＳ４とは異なる。脳梗塞、心筋梗塞に関して、評価値を求める処理の一例を図５に示す。図５のフローチャートは、図２のフローチャートにおけるステップＳ５での処理を、脳梗塞、心筋梗塞に関する処理に修正したものである。
【００６２】
ステップＳ５１において、処理対象とする最初の検査項目を指定する。具体的には、収縮期血圧、血中総コレステロール値、中性脂肪、HbA1c値、BMIおよび喫煙歴のうちの１つが指定される。
【００６３】
ステップＳ５２において、ステップＳ４で読み出された判定式ｙ（多重ロジスティック回帰分析によって予め決定された判定式）中の、ステップＳ５１で指定された検査項目を含む項の係数値αを計算する。例えば、収縮期血圧に関しては、SBPを含む項であるSBP×[a0＋Σ(ai×CSNPi)]のうち、a0＋Σ(ai×CSNPi)の値がαである。血中総コレステロール値に関しては、Tchを含む項であるTch×[b0＋Σ(bi×CSNPi)]のうち、b0＋Σ(bi×CSNPi)の値がαである。
【００６４】
ステップＳ５３において、記録部４から、指定された検査項目に関する、対象者と同年代の平均値ｔavを読み出す。ここで、平均値ｔavは、判定式ｙを求めるのに使用した被験者の集合中の被験者を年代毎に分類し、同じ年代の全被験者の検査値を平均して求めた値である。平均値ｔavは、年代および検査項目毎に予め計算されて、記録部４に記録されている。
【００６５】
ステップＳ５４において、ステップＳ５２で求めた値αを用いてｅαを計算し、さらに、ステップＳ５３で読み出された値ｔavおよびステップＳ２で読み出された対象者の検査値ｔを用いてβ＝（ｅα）(ｔ−ｔav)を計算する。
【００６６】
ステップＳ５５において、全ての検査項目についてステップＳ５２〜Ｓ５４の処理を終了したか否かを判断し、終了していなければステップＳ５６に移行し、別の検査項目を指定し、ステップＳ５２〜Ｓ５４の処理を繰り返す。
【００６７】
全ての検査項目についてステップＳ５２〜Ｓ５４の処理を終了した後、ステップＳ５７において、ステップＳ５４で計算した各検査項目に関する指数関数値βを全て乗算して評価値を求め、記録部４に記録し、ステップＳ６に移行する。
【００６８】
そして、ステップＳ５７で求められた評価値を提示する方法は、ステップＳ８での処理と同様である。心筋梗塞の場合の提示例を図４に示す。脳梗塞に関しても、心筋梗塞と同様に提示することができる。
【００６９】
図４に示したスケールは、判定式ｙを求めるのに使用した集合における、対象者と同じ同年代の平均値と比較した疾患になりやすさ（疾患危険度）を表している。スケールの左端に付された値１は、同年代の平均値を表している。このことは、次のことから理解できる。
【００７０】
ある事象を多重ロジスティック回帰分析した場合、説明変数ｘの回帰係数ａxのべき乗exp(ａx)は、説明変数ｘが１単位変化した場合の、その事象に関するオッズ比（ステップＳ５４で求めたｅαに対応）を表す。従って、説明変数である各検査項目に関して求めた値β（ステップＳ５４参照）は、所定の検査項目に関して、同年代の人（平均値ｔav）と比較した場合の、対象者（検査値ｔを有する）の疾患危険度のオッズ比である。よって、全検査項目の値βを乗算して最終的に得られた評価値は、同年代の人と比較した、全検査項目の影響を考慮した、対象者の疾患になりやすさを表しているといえる。
【００７１】
以上のように、本発明によれば、特定の対象者について、遺伝因子および環境因子を総合的に考慮した評価値を、同世代の平均値とともに、疾患リスクとしてグラフィカルに提示することができる。また、特定の疾患について、発症に関連する複数の検査項目の各々に対応させて、遺伝子を考慮した遺伝因子リスクと、遺伝因子および環境因子を考慮した総合リスクとを提示することができる。また、遺伝因子リスク、総合リスクの数値に対応するランクを表すマークとして提示することができる。従って、これらの提示に基づいて、医師は、その対象者が将来どのように身体を管理すべきか（生活習慣の改善など）を容易に知ることができ、対象者にわかり易く且つ適切に指導することが可能となる。
【００７２】
例えば、図３から、次のことが容易に分かる。即ち、評価値から、対象者は動脈硬化の状態が、同年代の人に比べて進んでいると推測できる。また、テーブルの３段目の遺伝因子リスクにマークが提示されていることから、対象者は血圧、総コレステロール、ヘモグロビンA1cおよび喫煙が高く（多く）なった場合に、動脈硬化が進展する遺伝因子を保有していることが容易にわかる。また、テーブルの２段目の検査値から、判定時点で、総コレステロール、中性脂肪および血糖コントロールが正常範囲を超えていることが容易に分かる。そして、テーブルの最下段の総合リスクに表示されたマークから、正常範囲を超えている総コレステロール、中性脂肪および血糖コントロールのうち、総コレステロールおよびヘモグロビンA1cが遺伝因子と相乗的に作用して、動脈硬化を進展させることが予想でき、中性脂肪を重点的に管理する必要があることが容易に分かる。また、血圧および喫煙は、判定時点で正常範囲内にあるが、これらに対応する遺伝因子リスクを有するため、注意が必要であることが容易にわかる。
【００７３】
以上、実施の形態を用いて本発明を説明したが、本発明は上記の実施の形態に限定されず、種々変更して実施することが可能である。
【００７４】
例えば、各処理の順序は、図２、図５に示したフローチャートの順序に限定されない。評価値、遺伝因子リスク、総合リスクの各値を計算する毎に、直ちに提示するようにしてもよい。
【００７５】
また、ステップＳ６での遺伝因子リスクの計算において、遺伝子多型セットを含まない検査値のみの項、即ち係数a0、b0、c0、d0、e0、ｆ0、ｇ0を含まない項の値を計算してもよい。例えば、収縮期血圧に関しては、SBPを含む項であるSBP×[a0＋Σ(ai×CSNPi)]のうち、少なくとも遺伝子多型セットを含むΣ(ai×CSNPi)の値を計算してもよい。これは、係数a0、b0、c0、d0、e0、ｆ0、ｇ0は定数であり、計算結果が一定値だけシフトされるだけであるので、所定の集合に関する平均値との比較や、度数分布からランクの境界値をパーセンタイル値で指定する場合には影響が無いからである。
【００７６】
また、多重ロジスティック回帰分析を使用する心筋梗塞、脳梗塞の場合、ステップＳ６の処理はステップＳ５２の処理と同じ処理であるので、ステップＳ５２の処理結果を記録部に記録しておき、ステップＳ６の処理を省略してもよい。
【００７７】
また、動脈硬化に関しては、ＩＭＴ測定を行なって得られた実測値があれば、ステップＳ５の処理を行わずに、その実測値を評価値として提示してもよい。
【００７８】
また、上記では、図３に示したような提示によって、医師が、対象者が将来どのように身体を管理すべきか（生活習慣の改善など）を容易に知ることができ、対象者にわかり易く且つ適切に指導することができることを示したが、さらに、検査値の改善による疾患リスクの軽減を定量的に提示することもできる。例えば、正常範囲を超えている総コレステロールの検査値が、所定値だけ減少した場合の疾患リスクの改善度Δｙ（％）を次の式で求めることができる。
Δｙ＝１００×[ｙ(Tch0)−ｙ(Tch0−ΔTch)]／ｙ(Tch0)
ここで、ｙ(Tch0)は、総コレステロール値がTch0である評価時の疾患リスクである。ｙ(Tch1)は、総コレステロール値として、評価時の検査値からΔTchだけ減少させた値（目標検査値）を用い、総コレステロール値以外の検査値（SBP、TG、HbA1c、BMI、SM、HDL）には評価時と同じ値を用いて得られた疾患リスクである。
【００７９】
このように、特定の対象者の検査値のうち、所定範囲を超えている検査値について、その検査値から所定の値だけ変化させた値である目標検査値を用いて疾患リスクを求め、上記の式によって、疾患リスク改善度（％）を求めて、検査項目毎に提示すれば、医師および検査者にとって有効である。
【００８０】
また、改善度を提示する場合、対応する検査項目を管理する必要がある旨と共に、目標検査値（変更値）および改善度を提示すればより有効である。例えば、「総コレステロールを管理する必要があります。総コレステロールが現状の値から××まで下がれば（△△だけ減少すれば）、疾患リスクが○○％改善します。」のように提示する。なお、改善度を求める式は上式に限定されず、判定式の値が評価時の値からどの程度変化するかを算出できればよい。また、検査項目によっては（例えばHDLコレステロール値）値が増大すれば疾患リスクが改善される検査値もあり、その場合にはその検査値を所定値だけ増大させて、改善度を求めればよい。
【００８１】
また、図３および図４に示した各疾患に関する評価値のスケール値は、一例であり、判定式ｙを求めるときの被験者の集合に応じて設定した値や表示であればよい。また、スケール値やスケール表示は、国や地域が異なれば、違った値や表示になり得る。
【００８２】
また、上記では、図３を用いて具体的に説明したが、図３はあくまで一例に過ぎない。判定式ｙを決定するために用いる被験者の集合が異なれば、決定された判定式ｙが異なる場合もある。例えば、その集合を構成する被験者の国籍や人種が異なれば、また日本においても被験者の居住地域が異なれば、判定式ｙが異なる可能性がある。従って、動脈硬化に関して、例えば、血糖コントロールおよびＢＭＩの各々と遺伝因子との積の項を含む判定式ｙを使用してもよい。同様に、心筋梗塞に関して、例えば、中性脂肪およびＢＭＩの各々と遺伝子多型セットとの積の項を含む判定式ｙを使用してもよく、脳梗塞に関して、例えば、中性脂肪と遺伝子多型セットとの積の項を含む判定式ｙを使用してもよい。
【００８３】
また、図３では、６種類の検査値を提示する場合を説明したが、これらに限定されない。これらに加えて、若しくはこれらの一部に代えてHDLコレステロール値を提示してもよく、その他の管理可能な検査値を提示してもよい。
【００８４】
また、疾患として動脈硬化、心筋梗塞および脳梗塞の３種類を対象とする場合を説明したが、これらに限定されず、その他の疾患に関しても本発明を適用することが可能である。その場合、疾患に応じて遺伝子多型セットおよび判定式ｙを決定すればよい。
【００８５】
また、提示するマークは任意の図形であってもよい。
【００８６】
また、提示装置の構成も図１に示した構成に限定されない。処理に必要なデータを記録した記録部が、ネットワークを介して別の装置（コンピュータ）に装備された記録装置であっても、コンピュータ読み取り可能な可搬性の記録媒体（光学式ディスク、半導体メモリなど）であってもよい。
【００８７】
また、提示装置を構成する表示部や印刷部を用いて、疾患リスクを表示し、紙に印刷する場合に限らず、計算結果のデータ（評価値、遺伝因子リスク、総合リスク）を所定構造のファイルデータに変換して、コンピュータ読み取り可能な記録媒体や通信回線を介して別の装置に移し、提示してもよい。
【図面の簡単な説明】
【００８８】
【図１】本発明の実施の形態に係る疾患リスクの提示装置を示すブロック図である。
【図２】本発明の実施の形態に係る疾患リスクの提示方法を示すフローチャートである。
【図３】本発明の実施の形態に係る疾患リスクの提示例を示す図である。
【図４】本発明の実施の形態に係る疾患リスク（評価値）の提示例を示す図であり、心筋梗塞の場合を示す。
【図５】本発明の実施の形態に係る、心筋梗塞または脳梗塞に関する評価値の計算方法を示すフローチャートである。
【図６】動脈硬化リスクの判定に使用する遺伝子多型セットの一例を示す表である。
【図７】心筋梗塞リスクの判定に使用する遺伝子多型セットの一例を示す表である。
【図８】脳梗塞リスクの判定に使用する遺伝子多型セットの一例を示す表である。
【図９Ａ】図６〜８の遺伝子多型セットを決定する際に用いた遺伝子多型を示す表である。
【図９Ｂ】図６〜８の遺伝子多型セットの決定する際に用いた遺伝子多型を示す表である。
【図９Ｃ】図６〜８の遺伝子多型セットの決定する際に用いた遺伝子多型を示す表である。
【図１０】Aldose redactase遺伝子の遺伝子型を決める規則を示す表である。
【符号の説明】
【００８９】
１ 提示装置
２ 演算部（ＣＰＵ）
３ メモリ
４ 記録部
５ 操作部
６ 表示部
７ 印刷部
８ 入出力部
９ 内部バス
【技術分野】
【０００１】
本発明は、個人の遺伝子情報およびその他の生体検査結果に応じて、その個人が疾患に罹患するリスクを判定して分かり易く提示する方法およびそのプログラムに関する。
【背景技術】
【０００２】
人の遺伝子に関する解析が進み、個人の有する遺伝子と疾患との関連性の解析が精力的に行なわれている。例えば、動脈硬化の発症や進展には生活習慣の違いと遺伝因子の個人差が関連することが知られている。従って、個人の遺伝子検査の結果や生活習慣の情報を利用して、各個人の身体状態の維持・管理（健康増進や罹患の防止、罹患者に関してはその進展の管理など）に役立つ情報を提供することが要望されている。これは、動脈硬化に限らず、心筋梗塞、脳梗塞などの、生活習慣の違いによって発症する可能性や発症後の進展の程度が異なる疾患に関しても同様に要望されている。なお、生活習慣の情報は、血液検査などの検査結果（検査値）や、問診による喫煙や飲酒に関する情報によって得られる情報を意味する。
【０００３】
例えば、下記特許文献１には、糖尿病に起因する動脈硬化、心筋梗塞、腎症及び網膜症について、疾患危険度（疾患の発症しやすさ、進行しやすさを表す指標）を視覚的に分かり易く提示する方法が開示されている。具体的には、被験者の遺伝子情報と環境因子との交互作用を含む疾患危険度の判定式を用いて、レーダーチャート、発症予測グラフ、及びバブルチャートで提示する方法が開示されている。また、特許文献１には、環境因子が変化した場合に疾患危険度の判定式の変化の程度を予測し、その予測値を基準値と比較して、疾患への対処の必要性を判定する方法が開示されている。
【特許文献１】国際公開ＷＯ２００６／１２６６１８
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかし、上記特許文献１に開示された視覚的な提示方法では、疾患毎の危険度の判定式を用いて求めた現在（評価時）の危険度を提示することしかできていなかった。また、疾患への対処の必要性の判定方法では、１つの判定式の変化のみで判断しているので、各環境因子の変化による影響の程度、特に管理可能な環境因子の変化による影響の程度を、分かり易く提示することができなかった。即ち、生活習慣の変化などによって検査値（血圧、総コレステロールなど）が変化した場合の危険度を予測するのに役立ち得る情報を分かり易く提示することはできなかった。
【０００５】
本発明は、上記した問題を解決するためのものであって、その目的は、個人の将来の疾患リスクをわかり易く提示することができる方法およびそのプログラムを提供することにある。
【課題を解決するための手段】
【０００６】
前記課題を解決するための手段は以下のとおりである。なお、ここでは後述する実施の形態に関する図面、符号を引用しているが、これは本発明の理解を容易にするためであり、本発明の限定を意図するものではない。
【０００７】
（１）記録部、演算部および出力部を備えた装置を用いて対象者の疾患リスクを求めて提示する方法であって、
前記演算部が、
前記対象者の検査項目（SBP、Tchなど）の値である検査値、
前記対象者が有する遺伝子多型セット（CSNPi）を表す遺伝子情報、並びに、
前記検査項目および前記遺伝子多型セットを変数とし、前記検査項目（SBPなど）と少なくとも前記遺伝子多型セットを含む多項式（Σ(ai×CSNPi)など）との積（SBP×[a0＋Σ(ai×CSNPi)]など）を含む判定式（例えば式１で表されるｙ）
を、前記記録部から読み出す第１ステップと、
前記演算部が、読み出された前記遺伝子情報を用いて前記多項式の値を計算する第２ステップと、
前記演算部が、前記第２ステップで計算された値に対応する第１図形データ（例えば図３の遺伝子リスクの欄に示された図形データ）を決定する第３ステップと、
前記演算部が、前記検査項目に対応させて、前記検査値および前記第１図形データを、前記出力部から出力する第４ステップとを含み、
前記判定式が、前記検査項目および前記遺伝子多型セットを説明変数として、重回帰分析または多重ロジスティック回帰分析によって得られた式である疾患リスクの提示方法。
【０００８】
（２）前記第４ステップの前に、前記演算部が、読み出された前記検査値および前記遺伝子情報を用いて、前記検査項目（SBPなど）を含む全ての項（SBP×[a0＋Σ(ai×CSNPi)]など）の値を計算する第５ステップと、
前記演算部が、前記第５ステップで計算された値に対応する第２図形データ（顔マーク）を決定する第６ステップとをさらに含み、
前記第４ステップにおいて、前記演算部が、前記検査項目に対応させて前記第２図形データを前記出力部から出力することを特長とする疾患リスクの提示方法。
【０００９】
（３）前記第４ステップの前に、前記演算部が、複数の前記検査値のうち所定範囲を超えている検査値を決定し、決定された該検査値を、前記判定式（ｙ）の値が低下する方向に所定値だけ変化させて目標検査値を求め、該目標検査値を用いて前記判定式から改善度（Δｙ）を求める第７ステップをさらに含み、
前記第４ステップにおいて、前記演算部が、決定された前記検査値に該当する検査項目に対応させて、前記改善度を前記出力部から出力する疾患リスクの提示方法。
【００１０】
（４）前記判定式が、複数の被験者の検査値および遺伝子情報の集合を用いて多重ロジスティック回帰分析によって得られた式であり、
年代毎に分類された前記被験者に関する検査値の平均値（ｔav）が、前記記録部に記録されており、
前記演算部が、前記対象者に関して、前記判定式中の所定の検査項目を含む項の係数値（α）を計算し、該係数値（α）を指数として自然対数の底のべき乗値（ｅα）を計算して前記所定の検査項目に関するオッズ比を求め、前記対象者の年齢が該当する年代に対応する前記平均値を前記記録部から読み出し、該平均値を前記対象者の前記検査値（ｔ）から減算した値（ｔ−ｔav）を指数として、前記オッズ比のべき乗値（β＝（ｅα）(ｔ−ｔav)）を計算する第８ステップと、
前記演算部が、前記検査項目毎の前記べき乗値（β）を、全て乗算して評価値を計算する第９ステップと、
前記演算部が、前記評価値を前記出力部から出力する第１０ステップと
をさらに含む疾患リスクの提示方法。
【００１１】
（５）前記判定式が、複数の被験者の検査値および遺伝子情報の集合を用いて重回帰分析によって得られた式であり、
年代毎に分類された前記被験者に関する前記判定式の値の平均値が、前記記録部に記録されており、
前記演算部が、読み出された前記検査値および前記遺伝子情報を用いて前記判定式の値を計算して評価値とする第１１ステップと、
前記演算部が、前記対象者の年齢が該当する年代に対応する前記平均値を、前記記録部から読み出す第１２ステップと、
前記演算部が、前記評価値および前記平均値を前記出力部から出力する第１３ステップと
をさらに含む疾患リスクの提示方法。
【００１２】
（６）前記疾患が心筋梗塞であり、
前記判定式に含まれる各々の各遺伝子多型セットが、図７に記載された遺伝子多型セットの１つである疾患リスクの提示方法。
【００１３】
（７）前記疾患が脳梗塞であり、
前記判定式に含まれる各々の遺伝子多型セットが、図８に記載された遺伝子多型セットの１つである疾患リスクの提示方法。
【００１４】
（８）前記疾患が動脈硬化であり、
前記判定式に含まれる各々の遺伝子多型セットが、図６に記載された遺伝子多型セットの１つである疾患リスクの提示方法。
【００１５】
（９）記録部および出力部を備えたコンピュータに、対象者の疾患リスクを求めて提示させるプログラムであって、
前記コンピュータに、
前記対象者の検査項目の値である検査値、
前記対象者が有する遺伝子多型セットを表す遺伝子情報、並びに、
前記検査項目および前記遺伝子多型セットを変数とし、前記検査項目と少なくとも前記遺伝子多型セットを含む多項式との積を含む判定式を、前記記録部から読み出す第１機能と、
読み出された前記遺伝子情報を用いて前記多項式の値を計算する第２機能と、
前記第２ステップで計算された値に対応する第１図形データを決定する第３機能と、
前記検査項目に対応させて、前記検査値および前記第１図形データを、前記出力部から出力する第４機能とを実現させ、
前記判定式が、前記検査項目および前記遺伝子多型セットを説明変数として、重回帰分析または多重ロジスティック回帰分析によって得られた式である疾患リスクの提示プログラム。
【発明の効果】
【００１６】
本発明によれば、特定の個人について、遺伝因子および環境因子を総合的に考慮した疾患リスクを、同世代の平均値とともにグラフィカルに提示することができる。
【００１７】
また、特定の疾患について、関連する複数の検査項目の各々に対応させて、遺伝因子を考慮した遺伝因子リスクと、遺伝因子および環境因子を考慮した総合リスクとを提示することができる。
【００１８】
また、遺伝因子リスク、総合リスクの数値が該当するランクをマークで分かり易く提示することができる。
【００１９】
従って、これらの提示に基づいて医師などが、その個人（評価対象者）が将来どのように身体を管理すべきか（生活習慣の改善など）を、わかり易く且つ適切に指導することが容易となる。
【発明を実施するための最良の形態】
【００２０】
以下、まず用語について説明した後、図面を参照して本発明についてより詳細に説明する。
【００２１】
「遺伝子多型」とは、同一集団内において、一つの遺伝子座に２種類以上の対立遺伝子（アレル）が存在する遺伝子の多様性を意味する。具体的には、ある集団において一定の頻度以上で存在する遺伝子の変異を示す。ここでいう遺伝子の変異は、ＲＮＡとして転写される領域に限定されるものではなく、プロモーター、エンハンサー等の制御領域などを含むヒトゲノム上で特定しうるすべてのＤＮＡにおける変異を含むものである。ヒトゲノムＤＮＡの９９．９％は各個人間で共通しており、残る０．１％がこのような多様性の原因となり、特定の疾患に対する感受性、薬物や環境因子に対する反応性の個人差として関与し得る。遺伝子多型があっても表現型に差が出るとは限らない。なお、ＳＮＰ（一塩基多型）も遺伝子多型の一種であるが、本発明が対象とする遺伝子多型はこれに限定されない。
【００２２】
このように、「遺伝子多型」は、遺伝子型（genotype）を考慮したものであって、それを含むものであり、本明細書において「遺伝子多型」とは特定の遺伝子型（genotype）を有する遺伝子多型を意味する。
【００２３】
本明細書で示す遺伝子型（genotype）は、「１」は置換塩基のうち塩基のアルファベット順（Ａ、Ｃ、Ｇ、Ｔ）で前にくる塩基を有する多型のホモを、「２」はヘテロを、「３」は、置換塩基のうち塩基のアルファベット順で後にくる塩基を有する多型のホモを表す。例えば、対象とする遺伝子（遺伝子多型）がXXXX(C3413G)（この表記は、遺伝子XXXXが、その塩基配列の3413番目の塩基がＣまたはＧである遺伝子多型を有することを意味する。）であるとすると、その遺伝子型（genotype）がCCのホモである場合を遺伝子型１、CGまたはGCのヘテロである場合を遺伝子型２、GGのホモである場合を遺伝子型３という。なお、「１２」は前記１と２の遺伝子型（genotype）の両方の遺伝子型、「２３」は前記２と３の遺伝子型（genotype）の両方の遺伝子型を表す。
【００２４】
なお、Aldose redactase遺伝子多型（Aldose redactase(AC)nで表され、ｎはACの繰り返し回数を表し、ｎ＝１３〜３０である）に関しては、上記の規則と異なり、次の規則で遺伝子型を付与している。すなわち、遺伝子の塩基配列中の「AC」の繰り返し回数がｎ回のホモである場合を遺伝子型「１」、「AC」の繰り返し回数がｎ回とｎ回以外とのヘテロである場合を遺伝子型「２」、これら以外の場合、即ち「AC」の繰り返し回数がｎ回以外のホモまたはヘテロである場合を遺伝子型「３」とする。
【００２５】
「遺伝子多型セット」とは、複数の遺伝子多型の組合せをいう。ここで複数の遺伝子多型とは、異なる遺伝子座を有する２種以上の遺伝子多型を意味する。
【００２６】
本明細書における塩基配列（ヌクレオチド配列）、核酸などの略号による表示は、ＩＵＰＡＣ−ＩＵＢの規定〔IUPAc-IUB communication on Biological Nomenclature, Eur. J. Biochem., 138; 9 (1984)〕、「塩基配列又はアミノ酸配列を含む明細書等の作製のためのガイドライン」（特許庁編）及び当該分野における慣用記号に従うものとする。
【００２７】
本明細書中において「遺伝子」は、特に言及しない限り、ヒトゲノムＤＮＡを含む２本鎖ＤＮＡ、及び１本鎖ＤＮＡ（センス鎖）、並びに当該センス鎖と相補的な配列を有する１本鎖ＤＮＡ（アンチセンス鎖）、及びそれらの断片のいずれもが含まれる、また、本明細書で「遺伝子」とは、特に言及しない限り、調節領域、コード領域、エクソン、及びイントロンを区別することなく示すものとする。
【００２８】
また、「動脈硬化」とは「頚動脈肥厚」を意味する。
【００２９】
図１は、本発明の実施の形態に係る疾患リスクの提示方法を実施するための装置（以下、提示装置と記す）を示すブロック図である。本提示装置１は、演算部（以下、ＣＰＵと記す）２と、メモリ３と、コンピュータ用ハードディスクなどの記録部４と、コンピュータ用キーボードやマウスなどの操作部５と、液晶ディスプレイなどの表示部６と、プリンタなどの印刷部７と、入出力のインタフェースである入出部８と、各部間でデータ（制御データを含む）を交換するための内部バス９とを備えている。
【００３０】
記録部４には、疾患リスクの計算に使用する遺伝子多型セット（図６〜８参照）の情報が参照テーブルとして記録されている。ここで、図６〜８には、動脈硬化、心筋梗塞および脳梗塞のそれぞれについて、図９Ａ〜９Ｃに示した１３３個の遺伝子多型の中から、上記特許文献１に開示された方法と同様の方法で決定された遺伝子多型セットを示している。なお、上記特許文献１では、有効な遺伝子多型セットの決定に正および負のｔ値を用いているが、図６〜８に示した遺伝子多型セットは、正のｔ値のみを用いて決定されたものである。
【００３１】
図９Ａ〜９Ｃにおいて、SNP-No.、GENETIC POLYMORPHISM、GENE、SITEはそれぞれ、各遺伝子多型に付した連続番号、遺伝子多型、遺伝子多型の略称、塩基配列上のアレルの位置を表す。図６〜８において、SSNP、Gtype、t-valueはそれぞれ、１つのセットを構成する遺伝子多型の数、遺伝子型、ｔ値を表す。SNP-No.は図９Ａ〜９Ｃに示した番号であり、GENEは図９Ａ〜９Ｃに示したGENE(略称)とSITEとを組み合わせた表記である。
【００３２】
なお、図９Ｃに示したSNP-No.が124の遺伝子は、Aldose redactase遺伝子であり、その遺伝子型は上記したように「AC」の繰り返し回数に応じて付与されている。さらに、図６〜８において、Aldose redactase遺伝子をALR_Z_4、ALR_Z_2、ALR_Z、ALR_Z2、ALR_Z4と表記しており、これらの遺伝子型（Gtype）は、図１０に示した規則で表記されている。例えば、ALR_Z_4に関しては、「AC」の繰り返し回数が２２回のホモである場合に遺伝子型を「１」とし、「AC」の繰り返し回数が２２回と２２回以外（即ち、１３〜２１回及び２３〜３０回の何れか）とのヘテロである場合に遺伝子型を「２」とし、「AC」の繰り返し回数が２２回以外のホモまたはヘテロである場合に遺伝子型を「３」とする。また、例えば、ALR_Z_2に関しては、「AC」の繰り返し回数が２３回のホモである場合に遺伝子型を「１」とし、「AC」の繰り返し回数が２３回以外（即ち、１３〜２２回及び２４〜３０回の何れか）のホモである場合に遺伝子型を「３」とし、「AC」の繰り返し回数が２３回と２３回以外とのヘテロである場合に遺伝子型を「２」とする。ALR_Z、ALR_Z2、ALR_Z4に関しても、繰り返し回数が異なるが同様である。
【００３３】
また、図９Ｂに示したSNP-No.が71のResistin VNTR 遺伝子多型（Resistin(ATG)n repeat)で表わされ、ｎはATGの繰り返し回数を表し、ｎ＝６〜８である）に関しては、一塩基多型と異なり、Aldose redatase(AC)nと同じく、即ち、遺伝子の塩基配列の「ATG」の繰り返し回数がｎ回のホモである場合を遺伝子型「１」、「ATG」の繰り返し回数がｎ回とｎ回以外とのヘテロである場合を遺伝子型「２」、これら以外の場合、即ち「ATG」の繰り返し回数がｎ回以外のホモまたはヘテロである場合を遺伝子型「３」とする。図７では、ｎ＝６のResistin遺伝子多型をresistin_6と表している。
【００３４】
記録部４にはさらに、判定対象者の検査値（血液検査結果および問診結果の情報）及び遺伝子情報が、対象者毎に付した情報（対象者ＩＤ）と対応させて記録されている。また、記録部４には、疾患危険度の判定式ｙが記録されている。判定式ｙは、後述するように検査項目および遺伝子多型セットを変数とする多項式であり、使用するコンピュータプログラム言語に応じた形式で記録されている。
【００３５】
以下、図２に示したフローチャートを用いて提示装置１を用いた、疾患リスクの提示方法について説明する。ここでは、特に断らない限り、ＣＰＵ２が行なう処理として記載する。ＣＰＵ２は、操作部５を介して指示やデータ入力を受け、処理対象データを記録部４からメモリ３に読み出し、メモリ３をワーク領域として使用して各処理を実行し、処理の途中結果や最終結果を、必要に応じて記録部４の所定領域に記録する。ＣＰＵ２は、処理結果を記録部４から読み出して、表示部６の入力信号形式に応じた所定形式のデータ（アナログ又はディジタル）に変換し、入出力部８を介して表示部６に伝送し、表示させる。また、ＣＰＵ２は、処理結果を記録部４から読み出して、印刷部７の入力信号形式に応じた所定形式のデータに変換し、印刷部７に伝送し、例えば紙に印刷させる。
【００３６】
ステップＳ１において、疾患リスクを判定して提示する対象者を指定するために操作部５が操作されて入力される対象者ＩＤを、入出力部８を介してメモリ３に一時記憶する。
【００３７】
ステップＳ２において、ステップＳ１でメモリ３に一時記憶された対象者ＩＤに対応するデータを記録部４からメモリ３に読み出す。このとき読み出されるデータは、対象者の遺伝子情報（遺伝子多型の有無を表す情報）および検査値である。検査値は、例えば、血液検査の結果である血圧（mmHg）、血中総コレステロール値（mg/dl）、血中中性脂肪値（mg/dl）、血糖コントロール（空腹時血糖値（mg/dl）、及びヘモグロビンA1c値（HbA1c）（％））、身体検査によって得られるBMI（Body-mass-index）（＝体重（kg）/（身長（m））2）、問診によって得られる喫煙歴である。ここで、血圧には収縮期血圧（SBP）および拡張期血圧（DBP）が含まれる。また、喫煙歴はこの判定時の喫煙習慣を表す情報であり、判定時に喫煙習慣を有してれば１、有していなければ０が設定される。なお、判定時に喫煙習慣を有していないとしても、過去に喫煙習慣を有していたことがあれば、喫煙歴には１が設定される。
【００３８】
ステップＳ３において、予め指定された複数の疾患のうち、処理対象とする最初の疾患を指定する。ここでは、動脈硬化が指定されたとする。また、その他に心筋梗塞および脳梗塞を対象とし、疾患リスクが判定されるとする。
【００３９】
ステップＳ４において、ステップＳ３で指定された疾患、すなわち動脈硬化に対応する危険度の判定式（回帰式）ｙを記録部４からメモリ３に読み出す。判定式ｙは予め決定されており、上記特許文献１に開示されているように、ｙ＝ΣａiＧi＋ΣｂjＥj＋ΣｄmnＥmＧn＋ｃ で表される。各係数ａi、ｂj、ｄmn、ｃは、所定の母集団を用いて重回帰分析によって疾患毎に決定された値である。例えば、動脈硬化の場合、ＩＭＴ（Intima Media Thickness）を目的変数とする。
【００４０】
上記した検査項目などの環境因子を説明変数として、判定式ｙを書き直せば、次式のように記述することができる。
ｙ＝SBP×[a0＋Σ(ai×CSNPi)]＋Tch×[b0＋Σ(bi×CSNPi)]
＋TG×[c0＋Σ(ci×CSNPi)]＋HbA1c×[d0＋Σ(di×CSNPi)]
＋BMI×[e0＋Σ(ei×CSNPi)]＋SM×[f0＋Σ(fi×CSNPi)]
−HDL×[g0＋Σ(gi×CSNPi)]＋const ・・・・・・・・・（式１）
ここで、SBP、Tch、TG、HbA1c、BMI、SMおよびHDLはそれぞれ、収縮期血圧（mmHg）、血中総コレステロール値（mg/dl）、中性脂肪（mg/dl）、HbA1c値（％）、BMI（kg/m2）、喫煙歴（１又は０）およびHDLコレステロール値（mg/dl）である。CSNPiは、i番目の遺伝子多型セットCSNPの有無に応じて１又は０を取る変数である。ai、bi、ci、di、ei、ｆi、gi（iは０以上の整数）は回帰係数であり、constは定数である。また、Σはiに関する和を求める演算子である。HDLを含む項の符号がマイナスになっているのは、HDLの疾患リスクへの影響が他の検査値とは逆であること、即ちHDLが大きくなると疾患リスクが低減することを明示的に示すためである。なお、各項の係数は、疾患および判定式ｙを求めるときに使用した被験者の集合に応じて異なる。また、疾患によっては、判定式ｙに一部の項が含まれていない場合もある。
【００４１】
上記の記述では、本発明の疾患リスクの判定に用いる７種類の環境因子（SBP、Tch、TG、HbA1c、BMI、SMおよびHDL）を含む項を明示的に記述しているが、判定式ｙには、上記の７種類以外の環境因子（例えば、性別（SEX）、年齢（AGE）、罹病期間など）を含む項や、環境因子を含まない遺伝子多型セットのみの項が含まれていてもよい。その場合、上記の７種類以外の環境因子を含む項（遺伝因子を含まない環境因子のみの項、及び、環境因子と遺伝因子との積である交互作用項を含む）はconstに含める。上記では、constは定数であるとしたが、「定数」とはこのような意味である。本発明では、個人が管理可能な情報を提示することを目的としているので、性別、年齢などの管理不可能な環境因子を含む項は、本発明の処理対象とはしない。
【００４２】
ステップＳ５において、ステップＳ２で読み出された対象者のデータとステップＳ４で読み出された判定式ｙとを用いて、ステップＳ１で指定された疾患に関する評価値を計算し、対象者ＩＤと対応させて記録部４に記録する。即ち、対象者の遺伝子多型セットのデータ（CSNPi）と検査値（SBP、Tch、TG、HbA1c、BMI、SM、HDL）とを判定式ｙに代入し、評価値を計算する。
【００４３】
ステップＳ６において、ステップＳ２で読み出された対象者のデータとステップＳ４で読み出された判定式ｙとを用いて、遺伝因子リスクを計算し、検査項目を特定する情報（検査項目ＩＤ）及び対象者ＩＤと対応させて記録部４に記録する。即ち、検査項目（BP、Tch、TG、HbA1c、BMI、SM）毎に、検査項目を含む項のうち検査項目を除いた部分の値を計算する。従って、この遺伝因子リスクは、対象者が持っている遺伝子に応じて決定される検査項目毎のリスクといえる。
【００４４】
例えば、収縮期血圧に関しては、SBPを含む項であるSBP×[a0＋Σ(ai×CSNPi)]のうち、a0＋Σ(ai×CSNPi)の値を計算し、血中総コレステロール値に関しては、Tchを含む項であるTch×[b0＋Σ(bi×CSNPi)]のうち、b0＋Σ(bi×CSNPi)を計算する。中性脂肪、HbA1c値、BMIおよび喫煙歴に関しても、同様に検査項目を除いた部分の値を計算する。
【００４５】
ステップＳ７において、ステップＳ２で読み出された対象者のデータとステップＳ４で読み出された判定式ｙとを用いて、総合リスクを計算し、検査項目ＩＤ及び対象者ＩＤと対応させて記録部４に記録する。即ち、検査項目（BP、Tch、TG、HbA1c、BMI、SM）毎に、特定の検査項目を含む全ての項、即ち、特定の検査項目と遺伝子多型セットの多項式（この多項式は特定の検査項目にとって係数に相当する。）との積を計算する。この総合リスクは、対象者が持っている遺伝子および検査値の両方を考慮して決定される検査項目毎のリスクといえる。
【００４６】
例えば、収縮期血圧に関しては、SBPを含む全ての項であるSBP×[a0＋Σ(ai×CSNPi)]の値を計算し、血中総コレステロール値に関しては、Tchを含む全ての項であるTch×[b0＋Σ(bi×CSNPi)]を計算する。中性脂肪、HbA1c値、BMIおよび喫煙歴に関しても、同様に計算する。
【００４７】
以上、ステップＳ５〜Ｓ７での処理によって、ステップＳ１で指定された対象者について、動脈硬化に関する評価値、遺伝因子リスクおよび総合リスクが求められる。次に、これらの値を提示する方法を説明する。動脈硬化に関する提示方法の一例を図３に示す。以下、図３を参照して具体的に説明する。
【００４８】
ステップＳ８において、対象者ＩＤに対応する評価値を記録部４から読み出し、表示部６の画面に、例えば図３の上段（左端に「評価」と表示された枠内）に示されているように提示する。具体的には、横長の細い長方形（以下、スケールと記す）を、例えばグラデーションを付して描画する。スケールは、基準尺度とするために、所定範囲のＩＭＴ値を表示したものである。スケールの左端に数値０を表示し、右側に向かって順に数値０．８、１．１、１．４、２．８を表示する。そして、対象者と同世代の人のＩＭＴの平均値に対応するスケール上の位置に矢印を表示する。さらに、評価値に対応するスケール上の位置に、人型の図形を表示し、その下に評価値、例えば１．２を表示する。
【００４９】
ここで、スケールに付す数値は、動脈硬化の程度を表す基準値であるＩＭＴの値であり、これらの値は予め設定されて、記録部４に記録されている。なお、図３に示した値（０．８、１．１、１．４、２．８）は、動脈硬化の程度を表す基準値（ＩＭＴ値）として一般に使用されている値である。
【００５０】
また、図３の例では、評価値と平均値とを比較して、評価を文章で自動的に表示している。図３の上段に表示された文章の下線部は条件によって変化する部分であり、それ以外の部分は固定である。即ち、左側の下線部分には疾患名称を表示し、右側の下線部分には、評価値と平均値との大小関係に応じた表示をする。例えば、評価値が平均値よりも大きい場合、「すすんでいる」と表示する。
【００５１】
ステップＳ９において、対象者ＩＤに対応する検査値、遺伝因子リスク及び総合リスクを記録部４から読み出し、表示部６の画面に、例えば図３の下段（左端に「分析表」と表示された枠内）に示したように提示する。即ち、テーブル形式で表示し、最上段に検査項目を、２段目に検査値を、３段目に遺伝因子リスクを、最下段に総合リスクをそれぞれ表示する。
【００５２】
図３のテーブルの２段目には、記録部４から読み出された検査値を検査項目毎に提示している。このとき、検査値が、検査項目毎に予め設定された所定範囲内にあるか否かの情報を併せて提示する。図３では、所定範囲から外れている検査値を白抜き文字で表示しているが、別の方法として、異なる色で表示してもよい。なお、血圧に関しては収縮期血圧（SBP）値および拡張期血圧（DBP）値を、血糖コントロールに関しては空腹時血糖値およびHbA1c値を表示し、喫煙歴に関しては、読み出された値が０であれば「無」、１であれば「有」を表示する。
【００５３】
図３のテーブルの３段目には、遺伝因子リスクの数値を直接表示せずに、リスクの程度をランク付けして表示するために、数値に応じたマーク（図形）を表示する。即ち、ステップＳ６で検査項目（ｉ）毎に求められた遺伝因子リスクＢiと、予め検査項目（ｉ）毎に決定された各ランク（ｊ）に対応する基準値（Ｒij）とを比較して、Ｒi,j≦Ｂi＜Ｒi,j+1である場合、ランク（ｊ）に対応するマークを記録部４から読み出して表示する。具体的には、判定式ｙを重回帰分析によって求めるときに使用した被験者の集合に関して、各被験者の検査項目毎の遺伝因子リスクを計算し（ステップＳ６参照）、その度数分布において、例えば２５、５０、７５パーセンタイルに該当する遺伝因子リスクの値を、検査項目毎の基準値（Ｒij）とする。例えば、血圧に関しては、判定式ｙ中の収縮期血圧SBPを含む項であるSBP×[a0＋Σ(ai×CSNPi)]のうち、a0＋Σ(ai×CSNPi)の値を、集合の各被験者について計算する。これによって、対象者の検査項目毎の遺伝因子リスクをマークで表示することができる。なお、図３では、一部の欄（中性脂肪、ＢＭＩに対応する欄）が空白になっているが、これは、使用した判定式ｙ中に該当する項が存在しないこと、若しくは、判定式ｙ中に該当する項は存在するが、その対象者が該当する遺伝子多型セットを有していないためにその項が０になったことを意味する。
【００５４】
図３のテーブルの最下段には、総合リスクの数値を直接表示せずに、リスクの程度をランク付けして表示するために、数値に応じたマークを表示する。遺伝因子リスクの提示に関して上記したのと同様に、ステップＳ７で検査項目（ｉ）毎に求められた総合リスクＣiと、予め検査項目（ｉ）毎に決定された各ランク（ｊ）に対応する基準値（Ｔij）とを比較して、Ｔi,j≦Ｃi＜Ｔi,j+1である場合、ランク（ｊ）に対応するマークを記録部４から読み出して表示する。具体的には、判定式ｙを重回帰分析によって求めるときに使用した被験者の集合に関して、各被験者の検査項目毎の遺伝因子リスクを計算し（ステップＳ７参照）、その度数分布において、例えば２５、５０、７５パーセンタイルに該当する遺伝因子リスクの値を、検査項目毎の基準値（Ｔij）とする。例えば、血圧に関しては、判定式ｙ中の収縮期血圧SBPを含む項であるSBP×[a0＋Σ(ai×CSNPi)]の値を、集合の各被験者について計算する。これによって、対象者の検査項目毎の総合リスクをマークで表示することができる。遺伝因子リスクが判定式y中に該当する項が存在しないために空白の欄（図３参照）があれば、それに対応する総合リスクとして表示するマークは、検査値のみから決定する。
【００５５】
以上、ステップＳ９およびＳ１０での処理によって、例えば図３に示したように、動脈硬化に関して対象者の遺伝因子リスク及び総合リスクを、ランクに対応するマークを用いてわかり易く提示することができる。
【００５６】
ステップＳ１０において、予め指定された複数の疾患（動脈硬化、心筋梗塞および脳梗塞）の全てについて処理が完了したか否かを判断し、完了していなければステップＳ１１に移行して別の疾患を指定した後、ステップＳ４に移行する。予め指定された複数の疾患の全てについて処理が完了していればステップＳ１２に移行し、終了の指示があったか否かを判断し、指示が無かった場合、ステップＳ１に移行して新たに対象者の指定を受け付ける。
【００５７】
以上によって、特定の対象者について、動脈硬化、心筋梗塞および脳梗塞の各疾患に関するリスクを計算し、図３と同様に提示することができる。なお、心筋梗塞および脳梗塞に関しても、上記した動脈硬化に関する処理と同様の処理を行えばよいが、動脈硬化の場合と異なる点もあるので、以下にその点について説明する。
【００５８】
心筋梗塞、脳梗塞に関しては、説明変数としては動脈硬化の場合と同様の環境因子（検査項目など）を用いるが、動脈硬化の場合と異なり多重ロジスティック回帰分析によって判定式ｙを求める。最終的に求められた判定式ｙは、上記の式１と同様の多項式の形ではあるが、左辺ｙの値はロジットである。疾患の発症確率をｐとすれば、ロジットはlog[ｐ／(１−ｐ)]で表される。ここで、logは自然対数を意味する。また、ｐ／(１−ｐ)はオッズ（Odds）である。
【００５９】
即ち、所定の被験者の集合に関して、環境因子（検査項目など）と脳梗塞の発症確率ｐから求められるオッズとを用いて、多重ロジスティック回帰分析によって脳梗塞の判定式ｙを求める。また、所定の被験者の集合に関して、環境因子（検査項目など）と心筋梗塞の発症確率ｐから求められるオッズとを用いて、多重ロジスティック回帰分析によって心筋梗塞の判定式ｙを求める。例えば、心筋梗塞に関して図７に示した遺伝子多型セットを使用して、中性脂肪およびＢＭＩの各々と遺伝子多型セットとの積の項を含まない判定式ｙが得られた。また、脳梗塞に関して図８に示した遺伝子多型セットを使用して、中性脂肪と遺伝子多型セットとの積の項を含まない判定式ｙが得られた。
【００６０】
従って、脳梗塞、心筋梗塞に関しては、判定式ｙの左辺がロジットである点で動脈硬化の場合と異なるが、遺伝因子リスクおよび総合リスクを求める処理はステップＳ６およびＳ７の処理と同様であり、遺伝因子リスクおよび総合リスクをマークで提示する処理はステップＳ９の処理と同様である。
【００６１】
しかし、脳梗塞、心筋梗塞に関して、評価値を求める処理はステップＳ４とは異なる。脳梗塞、心筋梗塞に関して、評価値を求める処理の一例を図５に示す。図５のフローチャートは、図２のフローチャートにおけるステップＳ５での処理を、脳梗塞、心筋梗塞に関する処理に修正したものである。
【００６２】
ステップＳ５１において、処理対象とする最初の検査項目を指定する。具体的には、収縮期血圧、血中総コレステロール値、中性脂肪、HbA1c値、BMIおよび喫煙歴のうちの１つが指定される。
【００６３】
ステップＳ５２において、ステップＳ４で読み出された判定式ｙ（多重ロジスティック回帰分析によって予め決定された判定式）中の、ステップＳ５１で指定された検査項目を含む項の係数値αを計算する。例えば、収縮期血圧に関しては、SBPを含む項であるSBP×[a0＋Σ(ai×CSNPi)]のうち、a0＋Σ(ai×CSNPi)の値がαである。血中総コレステロール値に関しては、Tchを含む項であるTch×[b0＋Σ(bi×CSNPi)]のうち、b0＋Σ(bi×CSNPi)の値がαである。
【００６４】
ステップＳ５３において、記録部４から、指定された検査項目に関する、対象者と同年代の平均値ｔavを読み出す。ここで、平均値ｔavは、判定式ｙを求めるのに使用した被験者の集合中の被験者を年代毎に分類し、同じ年代の全被験者の検査値を平均して求めた値である。平均値ｔavは、年代および検査項目毎に予め計算されて、記録部４に記録されている。
【００６５】
ステップＳ５４において、ステップＳ５２で求めた値αを用いてｅαを計算し、さらに、ステップＳ５３で読み出された値ｔavおよびステップＳ２で読み出された対象者の検査値ｔを用いてβ＝（ｅα）(ｔ−ｔav)を計算する。
【００６６】
ステップＳ５５において、全ての検査項目についてステップＳ５２〜Ｓ５４の処理を終了したか否かを判断し、終了していなければステップＳ５６に移行し、別の検査項目を指定し、ステップＳ５２〜Ｓ５４の処理を繰り返す。
【００６７】
全ての検査項目についてステップＳ５２〜Ｓ５４の処理を終了した後、ステップＳ５７において、ステップＳ５４で計算した各検査項目に関する指数関数値βを全て乗算して評価値を求め、記録部４に記録し、ステップＳ６に移行する。
【００６８】
そして、ステップＳ５７で求められた評価値を提示する方法は、ステップＳ８での処理と同様である。心筋梗塞の場合の提示例を図４に示す。脳梗塞に関しても、心筋梗塞と同様に提示することができる。
【００６９】
図４に示したスケールは、判定式ｙを求めるのに使用した集合における、対象者と同じ同年代の平均値と比較した疾患になりやすさ（疾患危険度）を表している。スケールの左端に付された値１は、同年代の平均値を表している。このことは、次のことから理解できる。
【００７０】
ある事象を多重ロジスティック回帰分析した場合、説明変数ｘの回帰係数ａxのべき乗exp(ａx)は、説明変数ｘが１単位変化した場合の、その事象に関するオッズ比（ステップＳ５４で求めたｅαに対応）を表す。従って、説明変数である各検査項目に関して求めた値β（ステップＳ５４参照）は、所定の検査項目に関して、同年代の人（平均値ｔav）と比較した場合の、対象者（検査値ｔを有する）の疾患危険度のオッズ比である。よって、全検査項目の値βを乗算して最終的に得られた評価値は、同年代の人と比較した、全検査項目の影響を考慮した、対象者の疾患になりやすさを表しているといえる。
【００７１】
以上のように、本発明によれば、特定の対象者について、遺伝因子および環境因子を総合的に考慮した評価値を、同世代の平均値とともに、疾患リスクとしてグラフィカルに提示することができる。また、特定の疾患について、発症に関連する複数の検査項目の各々に対応させて、遺伝子を考慮した遺伝因子リスクと、遺伝因子および環境因子を考慮した総合リスクとを提示することができる。また、遺伝因子リスク、総合リスクの数値に対応するランクを表すマークとして提示することができる。従って、これらの提示に基づいて、医師は、その対象者が将来どのように身体を管理すべきか（生活習慣の改善など）を容易に知ることができ、対象者にわかり易く且つ適切に指導することが可能となる。
【００７２】
例えば、図３から、次のことが容易に分かる。即ち、評価値から、対象者は動脈硬化の状態が、同年代の人に比べて進んでいると推測できる。また、テーブルの３段目の遺伝因子リスクにマークが提示されていることから、対象者は血圧、総コレステロール、ヘモグロビンA1cおよび喫煙が高く（多く）なった場合に、動脈硬化が進展する遺伝因子を保有していることが容易にわかる。また、テーブルの２段目の検査値から、判定時点で、総コレステロール、中性脂肪および血糖コントロールが正常範囲を超えていることが容易に分かる。そして、テーブルの最下段の総合リスクに表示されたマークから、正常範囲を超えている総コレステロール、中性脂肪および血糖コントロールのうち、総コレステロールおよびヘモグロビンA1cが遺伝因子と相乗的に作用して、動脈硬化を進展させることが予想でき、中性脂肪を重点的に管理する必要があることが容易に分かる。また、血圧および喫煙は、判定時点で正常範囲内にあるが、これらに対応する遺伝因子リスクを有するため、注意が必要であることが容易にわかる。
【００７３】
以上、実施の形態を用いて本発明を説明したが、本発明は上記の実施の形態に限定されず、種々変更して実施することが可能である。
【００７４】
例えば、各処理の順序は、図２、図５に示したフローチャートの順序に限定されない。評価値、遺伝因子リスク、総合リスクの各値を計算する毎に、直ちに提示するようにしてもよい。
【００７５】
また、ステップＳ６での遺伝因子リスクの計算において、遺伝子多型セットを含まない検査値のみの項、即ち係数a0、b0、c0、d0、e0、ｆ0、ｇ0を含まない項の値を計算してもよい。例えば、収縮期血圧に関しては、SBPを含む項であるSBP×[a0＋Σ(ai×CSNPi)]のうち、少なくとも遺伝子多型セットを含むΣ(ai×CSNPi)の値を計算してもよい。これは、係数a0、b0、c0、d0、e0、ｆ0、ｇ0は定数であり、計算結果が一定値だけシフトされるだけであるので、所定の集合に関する平均値との比較や、度数分布からランクの境界値をパーセンタイル値で指定する場合には影響が無いからである。
【００７６】
また、多重ロジスティック回帰分析を使用する心筋梗塞、脳梗塞の場合、ステップＳ６の処理はステップＳ５２の処理と同じ処理であるので、ステップＳ５２の処理結果を記録部に記録しておき、ステップＳ６の処理を省略してもよい。
【００７７】
また、動脈硬化に関しては、ＩＭＴ測定を行なって得られた実測値があれば、ステップＳ５の処理を行わずに、その実測値を評価値として提示してもよい。
【００７８】
また、上記では、図３に示したような提示によって、医師が、対象者が将来どのように身体を管理すべきか（生活習慣の改善など）を容易に知ることができ、対象者にわかり易く且つ適切に指導することができることを示したが、さらに、検査値の改善による疾患リスクの軽減を定量的に提示することもできる。例えば、正常範囲を超えている総コレステロールの検査値が、所定値だけ減少した場合の疾患リスクの改善度Δｙ（％）を次の式で求めることができる。
Δｙ＝１００×[ｙ(Tch0)−ｙ(Tch0−ΔTch)]／ｙ(Tch0)
ここで、ｙ(Tch0)は、総コレステロール値がTch0である評価時の疾患リスクである。ｙ(Tch1)は、総コレステロール値として、評価時の検査値からΔTchだけ減少させた値（目標検査値）を用い、総コレステロール値以外の検査値（SBP、TG、HbA1c、BMI、SM、HDL）には評価時と同じ値を用いて得られた疾患リスクである。
【００７９】
このように、特定の対象者の検査値のうち、所定範囲を超えている検査値について、その検査値から所定の値だけ変化させた値である目標検査値を用いて疾患リスクを求め、上記の式によって、疾患リスク改善度（％）を求めて、検査項目毎に提示すれば、医師および検査者にとって有効である。
【００８０】
また、改善度を提示する場合、対応する検査項目を管理する必要がある旨と共に、目標検査値（変更値）および改善度を提示すればより有効である。例えば、「総コレステロールを管理する必要があります。総コレステロールが現状の値から××まで下がれば（△△だけ減少すれば）、疾患リスクが○○％改善します。」のように提示する。なお、改善度を求める式は上式に限定されず、判定式の値が評価時の値からどの程度変化するかを算出できればよい。また、検査項目によっては（例えばHDLコレステロール値）値が増大すれば疾患リスクが改善される検査値もあり、その場合にはその検査値を所定値だけ増大させて、改善度を求めればよい。
【００８１】
また、図３および図４に示した各疾患に関する評価値のスケール値は、一例であり、判定式ｙを求めるときの被験者の集合に応じて設定した値や表示であればよい。また、スケール値やスケール表示は、国や地域が異なれば、違った値や表示になり得る。
【００８２】
また、上記では、図３を用いて具体的に説明したが、図３はあくまで一例に過ぎない。判定式ｙを決定するために用いる被験者の集合が異なれば、決定された判定式ｙが異なる場合もある。例えば、その集合を構成する被験者の国籍や人種が異なれば、また日本においても被験者の居住地域が異なれば、判定式ｙが異なる可能性がある。従って、動脈硬化に関して、例えば、血糖コントロールおよびＢＭＩの各々と遺伝因子との積の項を含む判定式ｙを使用してもよい。同様に、心筋梗塞に関して、例えば、中性脂肪およびＢＭＩの各々と遺伝子多型セットとの積の項を含む判定式ｙを使用してもよく、脳梗塞に関して、例えば、中性脂肪と遺伝子多型セットとの積の項を含む判定式ｙを使用してもよい。
【００８３】
また、図３では、６種類の検査値を提示する場合を説明したが、これらに限定されない。これらに加えて、若しくはこれらの一部に代えてHDLコレステロール値を提示してもよく、その他の管理可能な検査値を提示してもよい。
【００８４】
また、疾患として動脈硬化、心筋梗塞および脳梗塞の３種類を対象とする場合を説明したが、これらに限定されず、その他の疾患に関しても本発明を適用することが可能である。その場合、疾患に応じて遺伝子多型セットおよび判定式ｙを決定すればよい。
【００８５】
また、提示するマークは任意の図形であってもよい。
【００８６】
また、提示装置の構成も図１に示した構成に限定されない。処理に必要なデータを記録した記録部が、ネットワークを介して別の装置（コンピュータ）に装備された記録装置であっても、コンピュータ読み取り可能な可搬性の記録媒体（光学式ディスク、半導体メモリなど）であってもよい。
【００８７】
また、提示装置を構成する表示部や印刷部を用いて、疾患リスクを表示し、紙に印刷する場合に限らず、計算結果のデータ（評価値、遺伝因子リスク、総合リスク）を所定構造のファイルデータに変換して、コンピュータ読み取り可能な記録媒体や通信回線を介して別の装置に移し、提示してもよい。
【図面の簡単な説明】
【００８８】
【図１】本発明の実施の形態に係る疾患リスクの提示装置を示すブロック図である。
【図２】本発明の実施の形態に係る疾患リスクの提示方法を示すフローチャートである。
【図３】本発明の実施の形態に係る疾患リスクの提示例を示す図である。
【図４】本発明の実施の形態に係る疾患リスク（評価値）の提示例を示す図であり、心筋梗塞の場合を示す。
【図５】本発明の実施の形態に係る、心筋梗塞または脳梗塞に関する評価値の計算方法を示すフローチャートである。
【図６】動脈硬化リスクの判定に使用する遺伝子多型セットの一例を示す表である。
【図７】心筋梗塞リスクの判定に使用する遺伝子多型セットの一例を示す表である。
【図８】脳梗塞リスクの判定に使用する遺伝子多型セットの一例を示す表である。
【図９Ａ】図６〜８の遺伝子多型セットを決定する際に用いた遺伝子多型を示す表である。
【図９Ｂ】図６〜８の遺伝子多型セットの決定する際に用いた遺伝子多型を示す表である。
【図９Ｃ】図６〜８の遺伝子多型セットの決定する際に用いた遺伝子多型を示す表である。
【図１０】Aldose redactase遺伝子の遺伝子型を決める規則を示す表である。
【符号の説明】
【００８９】
１ 提示装置
２ 演算部（ＣＰＵ）
３ メモリ
４ 記録部
５ 操作部
６ 表示部
７ 印刷部
８ 入出力部
９ 内部バス
【特許請求の範囲】
【請求項１】
記録部、演算部および出力部を備えた装置を用いて対象者の疾患リスクを求めて提示する方法であって、
前記演算部が、
前記対象者の検査項目の値である検査値、
前記対象者が有する遺伝子多型セットを表す遺伝子情報、並びに、
前記検査項目および前記遺伝子多型セットを変数とし、前記検査項目と少なくとも前記遺伝子多型セットを含む多項式との積を含む判定式
を、前記記録部から読み出す第１ステップと、
前記演算部が、読み出された前記遺伝子情報を用いて前記多項式の値を計算する第２ステップと、
前記演算部が、前記第２ステップで計算された値に対応する第１図形データを決定する第３ステップと、
前記演算部が、前記検査項目に対応させて、前記検査値および前記第１図形データを、前記出力部から出力する第４ステップとを含み、
前記判定式が、前記検査項目および前記遺伝子多型セットを説明変数として、重回帰分析または多重ロジスティック回帰分析によって得られた式であることを特長とする疾患リスクの提示方法。
【請求項２】
前記第４ステップの前に、前記演算部が、読み出された前記検査値および前記遺伝子情報を用いて、前記検査項目を含む全ての項の値を計算する第５ステップと、
前記演算部が、前記第５ステップで計算された値に対応する第２図形データを決定する第６ステップとをさらに含み、
前記第４ステップにおいて、前記演算部が、前記検査項目に対応させて前記第２図形データを前記出力部から出力することを特長とする請求項１に記載の疾患リスクの提示方法。
【請求項３】
前記第４ステップの前に、前記演算部が、複数の前記検査値のうち所定範囲を超えている検査値を決定し、決定された該検査値を、前記判定式の値が低下する方向に所定値だけ変化させて目標検査値を求め、該目標検査値を用いて前記判定式から改善度を求める第７ステップをさらに含み、
前記第４ステップにおいて、前記演算部が、決定された前記検査値に該当する検査項目に対応させて、前記改善度を前記出力部から出力することを特徴とする請求項１または２に記載の疾患リスクの提示方法。
【請求項４】
前記判定式が、複数の被験者の検査値および遺伝子情報の集合を用いて多重ロジスティック回帰分析によって得られた式であり、
年代毎に分類された前記被験者に関する検査値の平均値が、前記記録部に記録されており、
前記演算部が、前記対象者に関して、前記判定式中の所定の検査項目を含む項の係数値を計算し、該係数値を指数として自然対数の底のべき乗値を計算して前記所定の検査項目に関するオッズ比を求め、前記対象者の年齢が該当する年代に対応する前記平均値を前記記録部から読み出し、該平均値を前記対象者の前記検査値から減算した値を指数として、前記オッズ比のべき乗値を計算する第８ステップと、
前記演算部が、前記検査項目毎の前記べき乗値を、全て乗算して評価値を計算する第９ステップと、
前記演算部が、前記評価値を前記出力部から出力する第１０ステップと
をさらに含むことを特長とする請求項１〜３の何れか１項に記載の疾患リスクの提示方法。
【請求項５】
前記判定式が、複数の被験者の検査値および遺伝子情報の集合を用いて重回帰分析によって得られた式であり、
年代毎に分類された前記被験者に関する前記判定式の値の平均値が、前記記録部に記録されており、
前記演算部が、読み出された前記検査値および前記遺伝子情報を用いて前記判定式の値を計算して評価値とする第１１ステップと、
前記演算部が、前記対象者の年齢が該当する年代に対応する前記平均値を、前記記録部から読み出す第１２ステップと、
前記演算部が、前記評価値および前記平均値を前記出力部から出力する第１３ステップと
をさらに含むことを特長とする請求項１〜３の何れか１項に記載の疾患リスクの提示方法。
【請求項６】
前記疾患が心筋梗塞であり、
前記判定式に含まれる各々の各遺伝子多型セットが、図７に記載された遺伝子多型セットの１つであることを特徴とする請求項４に記載の疾患リスクの提示方法。
【請求項７】
前記疾患が脳梗塞であり、
前記判定式に含まれる各々の遺伝子多型セットが、図８に記載された遺伝子多型セットの１つであることを特徴とする請求項４に記載の疾患リスクの提示方法。
【請求項８】
前記疾患が動脈硬化であり、
前記判定式に含まれる各々の遺伝子多型セットが、図６に記載された遺伝子多型セットの１つであることを特徴とする請求項５に記載の疾患リスクの提示方法。
【請求項９】
記録部および出力部を備えたコンピュータに、対象者の疾患リスクを求めて提示させるプログラムであって、
前記コンピュータに、
前記対象者の検査項目の値である検査値、
前記対象者が有する遺伝子多型セットを表す遺伝子情報、並びに、
前記検査項目および前記遺伝子多型セットを変数とし、前記検査項目と少なくとも前記遺伝子多型セットを含む多項式との積を含む判定式を、前記記録部から読み出す第１機能と、
読み出された前記遺伝子情報を用いて前記多項式の値を計算する第２機能と、
前記第２ステップで計算された値に対応する第１図形データを決定する第３機能と、
前記検査項目に対応させて、前記検査値および前記第１図形データを、前記出力部から出力する第４機能とを実現させ、
前記判定式が、前記検査項目および前記遺伝子多型セットを説明変数として、重回帰分析または多重ロジスティック回帰分析によって得られた式であることを特長とする疾患リスクの提示プログラム。
【請求項１】
記録部、演算部および出力部を備えた装置を用いて対象者の疾患リスクを求めて提示する方法であって、
前記演算部が、
前記対象者の検査項目の値である検査値、
前記対象者が有する遺伝子多型セットを表す遺伝子情報、並びに、
前記検査項目および前記遺伝子多型セットを変数とし、前記検査項目と少なくとも前記遺伝子多型セットを含む多項式との積を含む判定式
を、前記記録部から読み出す第１ステップと、
前記演算部が、読み出された前記遺伝子情報を用いて前記多項式の値を計算する第２ステップと、
前記演算部が、前記第２ステップで計算された値に対応する第１図形データを決定する第３ステップと、
前記演算部が、前記検査項目に対応させて、前記検査値および前記第１図形データを、前記出力部から出力する第４ステップとを含み、
前記判定式が、前記検査項目および前記遺伝子多型セットを説明変数として、重回帰分析または多重ロジスティック回帰分析によって得られた式であることを特長とする疾患リスクの提示方法。
【請求項２】
前記第４ステップの前に、前記演算部が、読み出された前記検査値および前記遺伝子情報を用いて、前記検査項目を含む全ての項の値を計算する第５ステップと、
前記演算部が、前記第５ステップで計算された値に対応する第２図形データを決定する第６ステップとをさらに含み、
前記第４ステップにおいて、前記演算部が、前記検査項目に対応させて前記第２図形データを前記出力部から出力することを特長とする請求項１に記載の疾患リスクの提示方法。
【請求項３】
前記第４ステップの前に、前記演算部が、複数の前記検査値のうち所定範囲を超えている検査値を決定し、決定された該検査値を、前記判定式の値が低下する方向に所定値だけ変化させて目標検査値を求め、該目標検査値を用いて前記判定式から改善度を求める第７ステップをさらに含み、
前記第４ステップにおいて、前記演算部が、決定された前記検査値に該当する検査項目に対応させて、前記改善度を前記出力部から出力することを特徴とする請求項１または２に記載の疾患リスクの提示方法。
【請求項４】
前記判定式が、複数の被験者の検査値および遺伝子情報の集合を用いて多重ロジスティック回帰分析によって得られた式であり、
年代毎に分類された前記被験者に関する検査値の平均値が、前記記録部に記録されており、
前記演算部が、前記対象者に関して、前記判定式中の所定の検査項目を含む項の係数値を計算し、該係数値を指数として自然対数の底のべき乗値を計算して前記所定の検査項目に関するオッズ比を求め、前記対象者の年齢が該当する年代に対応する前記平均値を前記記録部から読み出し、該平均値を前記対象者の前記検査値から減算した値を指数として、前記オッズ比のべき乗値を計算する第８ステップと、
前記演算部が、前記検査項目毎の前記べき乗値を、全て乗算して評価値を計算する第９ステップと、
前記演算部が、前記評価値を前記出力部から出力する第１０ステップと
をさらに含むことを特長とする請求項１〜３の何れか１項に記載の疾患リスクの提示方法。
【請求項５】
前記判定式が、複数の被験者の検査値および遺伝子情報の集合を用いて重回帰分析によって得られた式であり、
年代毎に分類された前記被験者に関する前記判定式の値の平均値が、前記記録部に記録されており、
前記演算部が、読み出された前記検査値および前記遺伝子情報を用いて前記判定式の値を計算して評価値とする第１１ステップと、
前記演算部が、前記対象者の年齢が該当する年代に対応する前記平均値を、前記記録部から読み出す第１２ステップと、
前記演算部が、前記評価値および前記平均値を前記出力部から出力する第１３ステップと
をさらに含むことを特長とする請求項１〜３の何れか１項に記載の疾患リスクの提示方法。
【請求項６】
前記疾患が心筋梗塞であり、
前記判定式に含まれる各々の各遺伝子多型セットが、図７に記載された遺伝子多型セットの１つであることを特徴とする請求項４に記載の疾患リスクの提示方法。
【請求項７】
前記疾患が脳梗塞であり、
前記判定式に含まれる各々の遺伝子多型セットが、図８に記載された遺伝子多型セットの１つであることを特徴とする請求項４に記載の疾患リスクの提示方法。
【請求項８】
前記疾患が動脈硬化であり、
前記判定式に含まれる各々の遺伝子多型セットが、図６に記載された遺伝子多型セットの１つであることを特徴とする請求項５に記載の疾患リスクの提示方法。
【請求項９】
記録部および出力部を備えたコンピュータに、対象者の疾患リスクを求めて提示させるプログラムであって、
前記コンピュータに、
前記対象者の検査項目の値である検査値、
前記対象者が有する遺伝子多型セットを表す遺伝子情報、並びに、
前記検査項目および前記遺伝子多型セットを変数とし、前記検査項目と少なくとも前記遺伝子多型セットを含む多項式との積を含む判定式を、前記記録部から読み出す第１機能と、
読み出された前記遺伝子情報を用いて前記多項式の値を計算する第２機能と、
前記第２ステップで計算された値に対応する第１図形データを決定する第３機能と、
前記検査項目に対応させて、前記検査値および前記第１図形データを、前記出力部から出力する第４機能とを実現させ、
前記判定式が、前記検査項目および前記遺伝子多型セットを説明変数として、重回帰分析または多重ロジスティック回帰分析によって得られた式であることを特長とする疾患リスクの提示プログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９Ａ】
【図９Ｂ】
【図９Ｃ】
【図１０】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９Ａ】
【図９Ｂ】
【図９Ｃ】
【図１０】
【公開番号】特開２００８−１９１７１６（Ｐ２００８−１９１７１６Ａ）
【公開日】平成２０年８月２１日（２００８．８．２１）
【国際特許分類】
【出願番号】特願２００７−２２３１５（Ｐ２００７−２２３１５）
【出願日】平成１９年１月３１日（２００７．１．３１）
【出願人】（５０５１９５７６４）株式会社サインポスト (6)
【Ｆターム（参考）】
【公開日】平成２０年８月２１日（２００８．８．２１）
【国際特許分類】
【出願日】平成１９年１月３１日（２００７．１．３１）
【出願人】（５０５１９５７６４）株式会社サインポスト (6)
【Ｆターム（参考）】
[ Back to top ]