(57)【要約】本発明は、生物個体から採取した検体の遺伝子発現産物量測定データに基づき生物個体の体内時刻を推定する装置１１である。標準検体における遺伝子発現産物量測定データの経時的変化が、所定時間を周期とする余弦波に類似する概日振動遺伝子を選択する概日振動遺伝子選択手段と、特定時間を周期とし、位相の異なる複数の余弦波から、選択された概日振動遺伝子の発現産物量の経時的変化に類似する概日発現変動曲線を選択する概日発現変動曲線選択手段と、選択された前記概日発現変動曲線を特定する情報を登録する登録手段を備える。
2004012128.tif 000002
【特許請求の範囲】
【請求項１】
生物個体から採取した標準検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定するための分子時刻表を作成する装置であって、
所定生物種の複数個体の所定部位から採取した標準検体中の遺伝子発現産物量測定データを入力するデータ入力手段と、
前記標準検体中において発現する遺伝子の中から、その遺伝子発現産物量測定データの経時的変化が、所定時間を周期とする余弦波に類似する概日振動遺伝子を選択する概日振動遺伝子選択手段と、
特定時間を周期とし、位相の異なる複数の余弦波から、前記選択された概日振動遺伝子の発現産物量の経時的変化に類似する前記概日発現変動曲線を選択する概日発現変動曲線選択手段と、
選択された前記概日発現変動曲線を特定する情報を登録する登録手段と、
を備えることを特徴とする分子時刻表作成装置。
【請求項２】
前記登録手段は、前記概日発現変動曲線上の値が最大となる時刻を、前記概日振動遺伝子の標準分子時刻として前記体内時刻の推定のために用いられる分子時刻表に登録すると共に、
前記発現産物量の前記概日振動遺伝子毎の平均値及び標準偏差を、前記概日振動遺伝子の標準発現量及び標準変動量として前記分子時刻表に登録することを特徴とする請求項１記載の分子時刻表作成装置。
【請求項３】
生物個体から採取した検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定する装置であって、
所定生物種の所定部位の概日振動遺伝子の発現産物量の経時的変化を示す概日発現変動曲線を特定する分子時刻表を記憶する分子時刻表記憶手段と、
前記生物個体の前記所定部位から採取された検体に含まれる前記概日振動遺伝子の前記遺伝子発現産物量測定データを入力するデータ入力手段と、
入力された前記測定データを、前記分子時刻表で特定される前記概日発現変動曲線と照合することにより、前記生物個体が概日リズム障害であるかの判定結果及び概日リズム障害でない場合には前記生物個体の体内時刻推定結果を含む体内時刻情報を導出する体内時刻情報導出手段と、を備えることを特徴とする体内時刻推定装置。
【請求項４】
生物個体から採取した標準検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定するための分子時刻表を、情報処理装置を用いて作成する方法であって、
所定生物種の複数個体の所定部位から採取した標準検体中の遺伝子発現産物量測定データを入力するデータ入力手順と、
前記標準検体中において発現する遺伝子の中から、その遺伝子発現産物量測定データの経時的変化が、所定時間を周期とする余弦波に類似する概日振動遺伝子を選択する概日振動遺伝子選択手順と、
特定時間を周期とし、位相の異なる複数の余弦波から、前記選択された概日振動遺伝子の発現産物量の経時的変化に類似する前記概日発現変動曲線を選択する概日発現変動曲線選択手順と、
選択された前記概日発現変動曲線を特定する情報を登録する登録手順と、
を行うことを特徴とする分子時刻表作成方法。
【請求項５】
生物個体から採取した検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を、情報処理装置を用いて推定する方法であって、
所定生物種の所定部位の概日振動遺伝子の発現産物量の経時的変化を示す概日発現変動曲線を特定する分子時刻表を記憶する分子時刻表記憶手順と、
前記生物個体の前記所定部位から採取された検体に含まれる前記概日振動遺伝子の前記遺伝子発現産物量測定データを入力するデータ入力手順と、
入力された前記測定データを、前記分子時刻表で特定される前記概日発現変動曲線と照合することにより、前記生物個体が概日リズム障害であるかの判定結果及び概日リズム障害でない場合には前記生物個体の体内時刻推定結果を含む体内時刻情報を導出する体内時刻情報導出手順と、を行うことを特徴とする体内時刻推定方法。
【請求項６】
生物個体から採取した標準検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定するための分子時刻表を作成する情報処理装置に、
所定生物種の複数個体の所定部位から採取した標準検体中の遺伝子発現産物量測定データを入力するデータ入力手順と、
前記標準検体中において発現する遺伝子の中から、その遺伝子発現産物量測定データの経時的変化が、所定時間を周期とする余弦波に類似する概日振動遺伝子を選択する概日振動遺伝子選択手順と、
特定時間を周期とし、位相の異なる複数の余弦波から、前記選択された概日振動遺伝子の発現産物量の経時的変化に類似する前記概日発現変動曲線を選択する概日発現変動曲線選択手順と、
選択された前記概日発現変動曲線を特定する情報を登録する登録手順と、
を実行させるための分子時刻表作成プログラム。
【請求項７】
生物個体から採取した検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定する情報処理装置に、
所定生物種の所定部位の概日振動遺伝子の発現産物量の経時的変化を示す概日発現変動曲線を特定する分子時刻表を記憶する分子時刻表記憶手順と、
前記生物個体の前記所定部位から採取された検体に含まれる前記概日振動遺伝子の前記遺伝子発現産物量測定データを入力するデータ入力手順と、
入力された前記測定データを、前記分子時刻表で特定される前記概日発現変動曲線と照合することにより、前記生物個体が概日リズム障害であるかの判定結果及び概日リズム障害でない場合には前記生物個体の体内時刻推定結果を含む体内時刻情報を導出する体内時刻情報導出手順と、を実行させるための体内時刻推定プログラム。
【請求項８】
生物個体から採取した標準検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定するための分子時刻表を作成すると共に、前記生物個体から採取した検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定するシステムであって、
前記体内時刻情報を提供する情報センターに設置されるサーバコンピュータと、該サーバコンピュータに情報送受信可能に接続された端末コンピュータとを備え、
前記サーバコンピュータは、
所定生物種の複数個体の所定部位から採取した標準検体中の遺伝子発現産物量測定データを入力する標準データ入力手段と、
前記標準検体中において発現する遺伝子の中から、その遺伝子発現産物量測定データの経時的変化が、所定時間を周期とする余弦波に類似する概日振動遺伝子を選択する概日振動遺伝子選択手段と、
特定時間を周期とし、位相の異なる複数の余弦波から、前記選択された概日振動遺伝子の発現産物量の経時的変化に類似する前記概日発現変動曲線を選択する概日発現変動曲線選択手段と、
前記概日発現変動曲線上の値が最大となる時刻を、前記概日振動遺伝子の標準分子時刻として前記体内時刻の推定のために用いられる分子時刻表に登録すると共に、前記発現産物量の前記概日振動遺伝子毎の平均値及び標準偏差を、前記概日振動遺伝子の標準発現量及び標準変動量として前記分子時刻表に登録する登録手段と、
前記生物個体の前記所定部位から採取された検体に含まれる前記概日振動遺伝子の前記遺伝子発現産物量測定データを入力する測定データ入力手段と、
入力された前記測定データを、前記分子時刻表で特定される前記概日発現変動曲線と照合することにより、前記生物個体が概日リズム障害であるかの判定結果及び概日リズム障害でない場合には前記生物個体の体内時刻推定結果を含む体内時刻情報を導出する体内時刻情報導出手段と、
前記導出された体内時刻情報を、前記端末コンピュータに送信する体内時刻情報送信手段と、を備えることを特徴とする体内時刻推定システム。
【発明の詳細な説明】
【技術分野】  本発明は分子時刻表作成装置、体内時刻推定装置、分子時刻表作成方法、体内時刻推定方法、分子時刻表作成プログラム、体内時刻推定プログラム及び体内時刻推定システムに係り、特に１回の採取で得た検体の検査結果から、コンピュータを用いて簡易に体内時刻を推定可能な分子時刻表作成装置、体内時刻推定装置、分子時刻表作成方法、体内時刻推定方法、分子時刻表作成プログラム、体内時刻推定プログラム及び体内時刻推定システムに関するものである。
【背景技術】  多くの生物は体内時計を有しており、睡眠覚醒リズム、血圧リズム、体温リズム、及び一部のホルモン分泌リズムなど種々の生体リズムを制御している。
体内時計の発振周期は約２４時間であるが、生物種によって多少異なり、ヒトの場合は約２５時間といわれている。環境サイクルのない定常環境下における約２４時間のリズムを概日リズムといい、日常生活では太陽光が、日々概日リズムを２４時間サイクルに同調させる最も強力な同調因子として働く。
概日リズムの異常は、睡眠覚醒リズム障害（睡眠相後退症候群、睡眠相前進症候群、非２４時間睡眠・覚醒症候群など）、季節性うつ病、時差症候群（ＪＥＴ−ＬＡＧ）、昼夜交代勤務労働者における睡眠障害、痴呆老人にみられる夜間徘徊、及びせん妄等の疾患の原因や、不登校，不出社の原因の一つになることが知られている。
更に、投薬時刻によって薬効や副作用の発現の強弱、薬物動態が影響されることが知られている。深夜勤務者，交替制勤務者等の場合、体内時刻が一般的な人とは異なる場合があり、一般的な人の時刻に合わせて投薬しても、その患者の体内時刻が投薬に適した時間でない場合には、十分な薬効の発現が期待できないことがある。
従って、概日リズム障害患者の診断、個人毎の投薬時間の決定、痴呆老人の夜間徘徊防止等のため、個人の体内時刻を個別に測定し、個人の体内時刻の情報を医療機関の神経科、投薬のためのカルテを記載する医師、老人病院、老人保健施設に提供できるシステムの開発が望まれており、体内時刻を推定する方法の研究が進められている。
体内時刻推定方法としては、経時的に採血して血中のメラトニン量を測定し、ピークを示す時刻及び量から個人の体内時刻を推定する方法、同様に経時的に活動量や睡眠を測定して睡眠・覚醒リズムを観察し、個人の体内時刻を推定する方法、体温を測定して体温リズムを観察することにより体内時刻を推定する方法などが考えられる。
しかし、これらの方法は、２４時間以上にわたって経時的に採血又は体温測定等を行う必要があり、患者及び医療従事者側にとって負担が大きく、医療現場で用いるには現実的な方法とはいえないという問題点があった。
一方、オリゴＤＮＡアレイを用いたショウジョウバエ頭部の経時的遺伝子発現産物量の測定（ＭｃＤｏｎａＬＤ，Ｍ．ａｎｄ  Ｒｏｓｂａｓｈ，Ｍ．，Ｃｅｌｌ，１０７，５６７−５７８．（２００１），ＣｌａｒＩＤｇｅ−Ｃｈａｎｇ，Ａ．ら、Ｎｅｕｒｏｎ，３２，６５７−６７１（２００１））やディファレンシャルディスプレイ法を用いたマウスの肝臓における概日振動遺伝子の探索（Ｋｏｒｎｍａｎｎ，Ｂ．ら、Ｎｕｃｌｅｉｃ  ＡｃＩＤｓ  Ｒｅｓ．，２９（１１），ｅ５１（２００１））により該組織において概日振動遺伝子が多数見出され、それぞれの概日振動遺伝子が最大発現値を示す時刻についても報告されている。
しかし、これらの測定データを用いて体内時刻の測定をする方法は知られておらず、体内時刻の測定をする方法の確立が望まれている。
本発明の目的は、上記問題点を解決することにあり、複数回の検体採取をする必要なく、一回の検体採取手順のみで個人の体内時刻を推定可能な分子時刻表作成装置、体内時刻推定装置、分子時刻表作成方法、体内時刻推定方法、分子時刻表作成プログラム、体内時刻推定プログラム及び体内時刻推定システムを提供することにある。
本発明の他の目的は、個人の体内時刻の推定結果を医療機関、老人保健施設、スポーツクラブ等の顧客に提供可能な分子時刻表作成装置、体内時刻推定装置、分子時刻表作成方法、体内時刻推定方法、分子時刻表作成プログラム、体内時刻推定プログラム及び体内時刻推定システムを提供することを目的とする。
【発明の開示】
本発明における分子時刻表作成装置は、請求項１に係る発明によれば、生物個体から採取した標準検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定するための分子時刻表を作成する装置であって、所定生物種の複数個体の所定部位から採取した標準検体中の遺伝子発現産物量測定データを入力するデータ入力手段と、前記標準検体中において発現する遺伝子の中から、その遺伝子発現産物量測定データの経時的変化が、所定時間を周期とする余弦波に類似する概日振動遺伝子を選択する概日振動遺伝子選択手段と、特定時間を周期とし、位相の異なる複数の余弦波から、前記選択された概日振動遺伝子の発現産物量の経時的変化に類似する前記概日発現変動曲線を選択する概日発現変動曲線選択手段と、選択された前記概日発現変動曲線を特定する情報を登録する登録手段と、を備えることにより解決される。
上記課題は、請求項４に係る分子時刻表作成方法によれば、生物個体から採取した標準検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定するための分子時刻表を、情報処理装置を用いて作成する方法であって、所定生物種の複数個体の所定部位から採取した標準検体中の遺伝子発現産物量測定データを入力するデータ入力手順と、前記標準検体中において発現する遺伝子の中から、その遺伝子発現産物量測定データの経時的変化が、所定時間を周期とする余弦波に類似する概日振動遺伝子を選択する概日振動遺伝子選択手順と、特定時間を周期とし、位相の異なる複数の余弦波から、前記選択された概日振動遺伝子の発現産物量の経時的変化に類似する前記概日発現変動曲線を選択する概日発現変動曲線選択手順と、選択された前記概日発現変動曲線を特定する情報を登録する登録手順と、を行うことにより解決される。
上記課題は、請求項６に係る分子時刻表作成プログラムによれば、生物個体から採取した標準検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定するための分子時刻表を作成する情報処理装置に、所定生物種の複数個体の所定部位から採取した標準検体中の遺伝子発現産物量測定データを入力するデータ入力手順と、前記標準検体中において発現する遺伝子の中から、その遺伝子発現産物量測定データの経時的変化が、所定時間を周期とする余弦波に類似する概日振動遺伝子を選択する概日振動遺伝子選択手順と、特定時間を周期とし、位相の異なる複数の余弦波から、前記選択された概日振動遺伝子の発現産物量の経時的変化に類似する前記概日発現変動曲線を選択する概日発現変動曲線選択手順と、選択された前記概日発現変動曲線を特定する情報を登録する登録手順と、を実行させるための分子時刻表作成プログラムにより解決される。
上記課題は、請求項８に係る体内時刻推定システムによれば、生物個体から採取した標準検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定するための分子時刻表を作成すると共に、前記生物個体から採取した検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定するシステムであって、前記体内時刻情報を提供する情報センターに設置されるサーバコンピュータと、該サーバコンピュータに情報送受信可能に接続された端末コンピュータとを備え、前記サーバコンピュータは、所定生物種の複数個体の所定部位から採取した標準検体中の遺伝子発現産物量測定データを入力する標準データ入力手段と、前記標準検体中において発現する遺伝子の中から、その遺伝子発現産物量測定データの経時的変化が、所定時間を周期とする余弦波に類似する概日振動遺伝子を選択する概日振動遺伝子選択手段と、特定時間を周期とし、位相の異なる複数の余弦波から、前記選択された概日振動遺伝子の発現産物量の経時的変化に類似する前記概日発現変動曲線を選択する概日発現変動曲線選択手段と、前記概日発現変動曲線上の値が最大となる時刻を、前記概日振動遺伝子の標準分子時刻として前記体内時刻の推定のために用いられる分子時刻表に登録すると共に、前記発現産物量の前記概日振動遺伝子毎の平均値及び標準偏差を、前記概日振動遺伝子の標準発現量及び標準変動量として前記分子時刻表に登録する登録手段と、前記生物個体の前記所定部位から採取された検体に含まれる前記概日振動遺伝子の前記遺伝子発現産物量測定データを入力する測定データ入力手段と、入力された前記測定データを、前記分子時刻表で特定される前記概日発現変動曲線と照合することにより、前記生物個体が概日リズム障害であるかの判定結果及び概日リズム障害でない場合には前記生物個体の体内時刻推定結果を含む体内時刻情報を導出する体内時刻情報導出手段と、前記導出された体内時刻情報を、前記端末コンピュータに送信する体内時刻情報送信手段と、を備えることにより解決される。
このように、前記標準検体中において発現する遺伝子の中から、その遺伝子発現産物量測定データの経時的変化が、所定時間を周期とする余弦波に類似する概日振動遺伝子を選択する概日振動遺伝子選択手段と、特定時間を周期とし、位相の異なる複数の余弦波から、前記選択された概日振動遺伝子の発現産物量の経時的変化に類似する前記概日発現変動曲線を選択する概日発現変動曲線選択手段と、選択された前記概日発現変動曲線を特定する情報を登録する登録手段と、を備えているため、所定生物種の所定部位に含まれる概日振動遺伝子の概日発現変動曲線を特定する時刻表を作成することができる。つまり、その生物種の個体について、概日リズム障害であるか否かの判定、概日リズム障害でない場合の体内時刻推定の基本となる時刻表を得ることが可能になるのである。
このとき、前記登録手段は、前記概日発現変動曲線上の値が最大となる時刻を、前記概日振動遺伝子の標準分子時刻として前記体内時刻の推定のために用いられる分子時刻表に登録すると共に、前記発現産物量の前記概日振動遺伝子毎の平均値及び標準偏差を、前記概日振動遺伝子の標準発現量及び標準変動量として前記分子時刻表に登録すると好適である。
このように構成しているため、所定生物種の所定部位に含まれる概日振動遺伝子の概日発現変動曲線を特定する分子時刻表を作成することができる。つまり、その生物種の個体について、概日リズム障害であるか否かの判定、概日リズム障害でない場合の体内時刻推定に用いる分子時刻表を得ることが可能になるのである。また、概日発現変動曲線を標準分子時刻、標準発現量、標準変動量のデータのみで登録でき、概日発現変動曲線を、少ないデータで管理可能となる。
上記課題は、請求項３に係る体内時刻推定装置によれば、生物個体から採取した検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定する装置であって、所定生物種の所定部位の概日振動遺伝子の発現産物量の経時的変化を示す概日発現変動曲線を特定する分子時刻表を記憶する分子時刻表記憶手段と、前記生物個体の前記所定部位から採取された検体に含まれる前記概日振動遺伝子の前記遺伝子発現産物量測定データを入力するデータ入力手段と、入力された前記測定データを、前記分子時刻表で特定される前記概日発現変動曲線と照合することにより、前記生物個体が概日リズム障害であるかの判定結果及び概日リズム障害でない場合には前記生物個体の体内時刻推定結果を含む体内時刻情報を導出する体内時刻情報導出手段と、を備えることにより解決される。
上記課題は、請求項５に係る体内時刻推定方法によれば、生物個体から採取した検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を、情報処理装置を用いて推定する方法であって、所定生物種の所定部位の概日振動遺伝子の発現産物量の経時的変化を示す概日発現変動曲線を特定する分子時刻表を記憶する分子時刻表記憶手順と、前記生物個体の前記所定部位から採取された検体に含まれる前記概日振動遺伝子の前記遺伝子発現産物量測定データを入力するデータ入力手順と、入力された前記測定データを、前記分子時刻表で特定される前記概日発現変動曲線と照合することにより、前記生物個体が概日リズム障害であるかの判定結果及び概日リズム障害でない場合には前記生物個体の体内時刻推定結果を含む体内時刻情報を導出する体内時刻情報導出手順と、を行うことにより解決される。
上記課題は、請求項７に係る体内時刻推定プログラムによれば、生物個体から採取した検体の遺伝子発現産物量測定データに基づき前記生物個体の体内時刻を推定する情報処理装置に、所定生物種の所定部位の概日振動遺伝子の発現産物量の経時的変化を示す概日発現変動曲線を特定する分子時刻表を記憶する分子時刻表記憶手順と、前記生物個体の前記所定部位から採取された検体に含まれる前記概日振動遺伝子の前記遺伝子発現産物量測定データを入力するデータ入力手順と、入力された前記測定データを、前記分子時刻表で特定される前記概日発現変動曲線と照合することにより、前記生物個体が概日リズム障害であるかの判定結果及び概日リズム障害でない場合には前記生物個体の体内時刻推定結果を含む体内時刻情報を導出する体内時刻情報導出手順と、を実行させるための体内時刻推定プログラムにより解決される。
このように、測定データを、前記分子時刻表で特定される前記概日発現変動曲線と照合することにより、前記生物個体が概日リズム障害であるかの判定結果及び概日リズム障害でない場合には前記生物個体の体内時刻推定結果を含む体内時刻情報を導出する体内時刻情報導出手段を備えているため、１回の採取で得た検体の測定データと所定生物種の所定部位の分子時刻表により特定される前記概日発現変動曲線とを照合することにより、検体を複数回採取する必要なく、１回の採取で得た検体に基づいて簡易に概日リズム障害か否かの判定及び体内時刻の推定が可能となる。
【図面の簡単な説明】  図１は、本発明の一実施形態に係る体内時刻推定システムの全体構成を示す説明図、図２は、本発明の一実施形態に係る体内時刻推定システムを統括する体内時刻推定装置のハード構成を示す説明図、図３は、経時的遺伝子発現産物量テーブルの構成を示す説明図、図４は、分子時刻表テーブルの構成を示す説明図、図５は、体内時刻情報テーブルの構成を示す説明図、図６は、本発明の一実施形態に係る体内時刻推定システムの処理の流れを示すブロック図、図７は、本発明の一実施形態に係る概日振動遺伝子選択処理を示すフローチャート、図８は、遺伝子選択用余弦波を示す説明図、図９は、本発明の一実施形態に係る概日発現変動曲線・分子時刻表作成処理を示すフローチャート、図１０は、本発明の一実施形態に係る体内時刻の推定及び情報提供処理を示すフローチャート、図１１は、本発明の一実施形態に係る体内時刻の推定及び情報提供処理のうち初期設定の処理を示すフローチャート、図１２は、本発明の一実施形態に係る体内時刻の推定及び情報提供処理のうちサンプルの概日振動遺伝子発現産物量の正規化処理を示すフローチャート、図１３は、本発明の一実施形態に係る体内時刻の推定及び情報提供処理のうち正規化された概日振動遺伝子発現産物量と時刻ｔ時における各概日振動遺伝子産物の推定発現量との間のピアソン相関係数ｃの計算処理を示すフローチャート、図１４は、本発明の一実施形態に係る体内時刻情報報告画面を示す説明図、図１５は、検査対象者特性テーブルを示す説明図、図１６は、データ解析結果の一例を示す説明図である。
【発明を実施するための最良の形態】  以下、本発明の実施例を図面に基づいて説明する。なお、以下に説明する部材、配置等は本発明を限定するものでなく、本発明の趣旨の範囲内で種々改変することができるものである。
本発明の体内時刻推定システムとは、情報センター内のサーバコンピュータに入力・受信された人体の検体中の遺伝子発現産物量測定データを用いて、概日リズム障害か否かの判定結果と概日リズム障害でない場合の体内時刻推定結果とを含む体内時刻情報を導出すると共に、この体内時刻情報を、判定依頼者である顧客の端末コンピュータに送信するシステムである。
本実施形態では、人体の体内時刻推定を行う例について説明するが、本実施形態の体内時刻推定システムは、マウス，ラット，犬等の哺乳動物など、概日リズムを有する生物であれば、いかなる生物の体内時刻推定に用いてもよい。
また、本実施形態ではヒトの血液検体を用いて各個人の体内時刻推定を行う例について説明するが、マウスの視交叉上核，肝臓，ショウジョウバエの頭部，ヒトの皮膚，末梢血の白血球，口腔粘膜，その他各臓器など、概日リズムを有する検体であれば、いかなる検体を体内時刻推定に用いてもよい。
更に、本実施形態では、遺伝子発現産物量測定データとして、遺伝子の転写産物であるメッセンジャーＲＮＡ（以下、ｍＲＮＡとする）量の測定データを用いるが、遺伝子によりコードされる蛋白質量の測定データ、遺伝子によりコードされる蛋白質の修飾量の測定データ、遺伝子にコードされる酵素の活性量の測定データ、酵素によって代謝される化合物量の測定データ、遺伝子発現産物により制御されている体温・血圧・心拍数・血流量・呼吸数・脳波・活動量・体液ｐＨ等の物理量の測定データなど、概日リズムを有するものであれば、いかなる測定データを遺伝子発現産物量測定データとして用いても良い。
「概日振動遺伝子」は、明期と暗期が交互に繰り返す２４時間周期の明暗条件下或いは恒暗条件下で、遺伝子発現が約２４時間（２０〜２８時間）周期で振動を示すものをいう。
「標準発現量」は遺伝子の経時的発現量測定データの平均値、「標準変動量」は遺伝子の経時的発現量測定データの標準偏差、「相対発現量」は、遺伝子の遺伝子発現産物の測定値から標準発現量を減じ、標準変動量で除すことによって正規化した発現量を意味する。
「正規化」とは、各遺伝子の発現変動を同一基準で評価するための手段であって、各遺伝子の遺伝子発現産物量の測定値或いは各概日発現変動曲線上の値から標準発現量を減じ、標準変動量で除すことをいう。これによって、各概日振動遺伝子の標準変動量を１、標準発現量を０として統一的に扱うことが可能となる。
「概日発現変動曲線」とは、有限の発現データから概日振動遺伝子の経時的な発現状態を数式化したもので、経時的に測定した発現量をもとに複数の方法で作成可能である。「相対概日発現変動曲線」とは、各概日発現変動曲線上の値から標準発現量を減じ、標準変動量で除すことによって正規化した概日発現変動曲線を意味する。
「ＺＴ時刻」は、明暗条件下における点灯時をＺＴ０時とした時刻を示す。１２時間点灯した後、ＺＴ１２時に消灯する。従って、ＺＴ０時からＺＴ１２時までが明期で、その後ＺＴ１２時からＺＴ２４時までが暗期になる。
「ＣＴ時刻」は、明期・暗期の２４時間周期にリズム同調を行った後、本来ＺＴ０時に相当する時刻に点灯を行わずにＣＴ０時とする。ＣＴ０時からの恒暗条件下における２４時間周期の時刻をＣＴ時刻として表記する。
「分子時刻」とは概日振動遺伝子の概日発現変動曲線が最大値をとるＣＴ時刻或いはＺＴ時刻である。
「類似性」とは対象と対象との間の類似の度合いを表す数値である。類似性には相関係数のように値の大きい方ほど類似性が高いことを表わす「類似度」と、距離のようにその値が小さい方ほど類似性が高いことを表わす「非類似度」とがある。
非類似度を表す指標の一つである「距離」は、対象と対象の間の隔たりの程度を示す指標のことである。「ユークリッドの平方距離」は、最も良く用いられる距離のうちの一つであり、データが（ｘ１，ｙ１）、（ｘ２，ｙ２）、（ｘ３，ｙ３）…（ｘｎ，ｙｎ）で与えられた場合、変数ｘとｙの間の距離ｄｘｙは
2004012128.tif 000003
で定義される（文献３）。
「相関係数」とは類似度を表す指標の一つで相関の程度を示す指標のことである。「ピアソンの積率相関係数」は、最も良く用いられる相関係数の一つであり、
2004012128.tif 000004
で定義される（文献４）。ただし
2004012128.tif 000005
とする。
余弦波の標準偏差ｓについて説明する。余弦波上の対応する時刻ｔにおける値Ａ（ｔ）は、次のように求めることが出来る。
2004012128.tif 000006
周期時間：２４（時間）
Ｔ：余弦波の分子時刻（例えば、０／６０，１０／６０，２０／６０，３０／６０，…，１４２０／６０，１４３０／６０までの１０分刻みの１４４種類の時刻とすることができる）
ｔ：対応する時刻
余弦波上の経時的データ（サンプリング時刻における値の集合）の標準偏差ｓは、次のように求めることができる。
2004012128.tif 000007
例えば、４時間間隔で余弦波上の経時的データを取得した場合、ｔは次のような６点の時刻の組（ｎ＝６）として表記できる。（α＋０）／６０，（α＋２４０）／６０，（α＋４８０）／６０，（α＋７２０）／６０，（α＋９６０）／６０，（α＋１２００）／６０、（但し、αは、０分から２４０分までのある時刻（２４０分は含まず））。このときの標準偏差は、０．７７４５９６７となる。
本明細書中で「コンピュータ」とは、演算装置を備えた情報端末すべてを含む意味である。例えば、スーパーコンピュータ、汎用コンピュータ、オフィスコンピュータ、制御用コンピュータ、ワークステーション、パソコンのほか、携帯情報端末、演算装置を備えた携帯電話、ウェアラブルコンピュータ等をも含む。
また、本明細書中で「入力」とは、コンピュータに指示やデータを入れることをいい、キーボード、マウス等で指示やデータを入力する場合のほか、フレキシブルディスク、コンパクトディスク等の記憶媒体から指示やデータを入れる場合、通信回線を介して外部装置から指示やデータを受信する場合も含む。
（実施の形態１）
本実施形態の体内時刻推定システムは、図１に示すように、個人の体内時刻情報を顧客に提供する情報センター１と、情報センター１から個人の体内時刻情報の提供を受ける顧客としての医療機関２、スポーツセンター３、個人会員４、老人保健施設５、法人会員６、教育機関７を備えている。
情報センター１にはサーバコンピュータ１１が設置され、インターネット１３を介して顧客２〜７の各端末コンピュータ２１〜７１が接続されている。
情報センター１は、本実施形態の体内時刻推定システムを統括管理する組織であり、検査会社、製薬会社、医療機関、研究機関、分析センター、検査センター、データセンター等に含まれていてもよい。
情報センター１には、検体の遺伝子発現産物量を測定し、サーバコンピュータ１１に遺伝子発現産物量測定データを入力するＤＮＡチップ読取装置１２と、ＤＮＡチップ読取装置１２から受信した遺伝子発現産物量測定データを用いて概日リズム障害の判定処理、体内時刻推定処理を行うと共に、これらの処理で導出された体内時刻情報を端末コンピュータ２１〜７１に送信する体内時刻情報提供装置としてのサーバコンピュータ１１とが設置されている。
なお、本実施形態では、遺伝子発現産物量測定データを、ＤＮＡチップ読取装置１２から受信しているが、顧客２〜７に遺伝子発現産物量測定装置を設置し、端末コンピュータ２１〜７１から受信するようにしてもよい。
情報センター１は、顧客２〜７から検査対象者の血液サンプル１９を送付又は持込みにより受領する。公知の方法により、血液サンプル１９を、遺伝子の発現ｍＲＮＡ量を測定可能な状態に処理する。公知の方法として、例えば、血液サンプル１９から市販のＲＮＡ抽出キット、例えばＲＮＡｑｕｅｏｕｓ−Ｂｌｏｏｄ  Ｍｏｄｕｌｅ、ＲＮＡｑｕｅｏｕｓ  Ｐｈｅｎｏｌ−ｆｒｅｅ  Ｔｏｔａｌ  ＲＮＡ  ｋｉｔ（フナコシ社）を用い、添付のプロトコールに従ってｍＲＮＡを分離・精製する。なお、血液サンプル１９の処理として、グアニジン・チオシアネート塩化セシウム法、グアニジン・チオシアネート・ホット・フェノール法、グアニジン・チオシアネート−グアニジン・塩酸法等でｍＲＮＡを抽出し、ｍＲＮＡをオリゴ（ｄＴ）セルロースカラムに吸着・溶出させ、精製してもよい。
その後、処理された血液サンプル１９について、ＤＮＡチップ読取装置１２で発現ｍＲＮＡ量を測定する。測定されたｍＲＮＡ量の生データは、ＤＮＡチップ読取装置１２に蓄積される。
ＤＮＡチップ読取装置１２に蓄積された生データは、サーバコンピュータ１１に送信される。サーバコンピュータ１１では、この生データを加工して、検査対象者が概日リズム障害であるかどうかの判定を行い、概日リズム障害でない場合にはこの検査対象者の体内時刻を導出する。
導出された検査対象者の概日リズム障害の判定結果及び概日リズム障害でない場合の体内時刻の情報は、この検査対象者の体内時刻情報として、依頼した顧客の端末コンピュータ２１〜７１に送信され、報告される。
顧客２〜７は、情報センター１の顧客であって、情報センター１から本実施形態の体内時刻推定システムによる体内時刻情報の提供を受ける。
顧客２〜７にはそれぞれ、情報センター１のサーバコンピュータ１１が提供する体内時刻情報表示画面を閲覧可能な端末コンピュータ２１〜７１が設置されている。
なお、遺伝子発現産物量を測定するための遺伝子発現産物量測定装置を顧客２〜７に設置してもよい。この場合、検査対象者の血液サンプル１９の発現ｍＲＮＡ量を顧客２〜７側で測定し、その測定データを端末コンピュータ２１〜７１からサーバコンピュータ１１に送信し、その測定データに対する概日リズム障害の判定結果、体内時刻の情報を、サーバコンピュータ１１から端末コンピュータ２１〜７１で受信することになる。
医療機関２は、睡眠覚醒リズム障害、季節性うつ病、時差症候群、昼夜交代勤務労働者における睡眠障害等の患者の概日リズム障害の判定、患者の投薬時刻の決定等のために、本実施形態の体内時刻推定システムを利用する。
医療機関２にはまた、診療記録、処方箋、各検査結果等が電子データとして患者毎に登録された電子カルテサーバコンピュータ２２が設置されている。電子カルテサーバコンピュータ２２は、院内ＬＡＮで端末コンピュータ２１と接続され、サーバコンピュータ１１から端末コンピュータ２１にダウンロードした体内時刻情報を電子カルテサーバコンピュータ２２に入力可能に構成されている。
後述する検査対象者番号を、電子カルテサーバコンピュータ２２に登録された患者番号と同じ番号にしておくと共に、サーバコンピュータ１１から端末コンピュータ２１にダウンロードした検査対象者毎の体内時刻情報を自動で電子カルテサーバコンピュータ２２の不図示のカルテデータベースに取込むプログラムを格納しておけば、電子カルテに体内時刻情報を直接自動で取込むよう構成できる。
また、医療機関２の代わりに、心理カウンセラー、保健所等の保健福祉施設、診療所、児童相談所、薬局等が顧客であってもよい。
スポーツセンター３は、スポーツ選手の指導、トレーニング、育成等をする組織であって、スポーツの種類は問わない。スポーツセンター３は、国際大会参加前のスポーツ選手の練習スケジュール作成等のために、本実施形態の体内時刻推定システムを利用する。
スポーツセンター３には、端末コンピュータ３１のほかに、スポーツ選手の健康状態、練習メニュー、食事メニュー、大会参加記録等が電子データとして選手毎に登録された選手管理サーバコンピュータ３２が設置されている。選手管理サーバコンピュータ３２は、センター内ＬＡＮで端末コンピュータ３１と接続され、サーバコンピュータ１１から端末コンピュータ３１にダウンロードした体内時刻情報を選手管理サーバコンピュータ３２に入力可能に構成されている。
後述する検査対象者番号を、選手管理サーバコンピュータ３２に登録された選手番号と同じ番号にしておくと共に、サーバコンピュータ１１から端末コンピュータ３１にダウンロードした検査対象者毎の体内時刻情報を自動で選手管理サーバコンピュータ３２の不図示のデータベースに取込むプログラムを格納しておけば、データベースに体内時刻情報を直接自動で取込むよう構成できる。
また、スポーツセンター３の代わりに、スポーツジム、スポーツクラブ、クラブ活動を管理する法人やクラブ活動を管轄する教育機関等が顧客であってもよい。
個人会員４は、自分の生活改善を目的として個人的に本実施形態の体内時刻推定システムを利用する顧客である。
例えば、健康志向の高い個人、高齢者、過去に睡眠覚醒リズム障害、季節性うつ病、時差症候群等の経験のある者、昼夜交代勤務労働者等である。また、これらの者の家族や介護者が、本人から検体を採取できる場合には、検体を採取し、本人の体内時刻情報の閲覧・管理を代わりに行ってもよい。
老人保健施設５は、施設に滞在している老人の健康管理、生活管理のために本実施形態の体内時刻推定システムを利用する。
例えば、部屋割りを決定するための情報として、体内時刻情報を用いてもよい。
入所時に体内時刻情報を情報センター１から取得することにより、体内時刻の近い者を同じ部屋に滞在させることが可能である。これにより、同じ部屋の者は、体内時刻が近く、生活ペースが類似するため、起床時や就寝時の見回りなど、滞在者の管理を楽に行うことが可能になる。
また、入所時に体内時刻情報を取得しておけば、施設側も、入所後早い時期にその滞在者の生活リズムに対応できるため、新規入所者を早い時期に施設全体のペースに適合させることが可能になる。
老人保健施設５には、端末コンピュータ５１のほかに、滞在者の健康状態、食事メニュー、家族構成、リハビリメニュー等が電子データとして滞在者毎に登録された管理サーバコンピュータ５２が設置されている。管理サーバコンピュータ５２は、施設内ＬＡＮで端末コンピュータ５１と接続され、サーバコンピュータ１１から端末コンピュータ５１にダウンロードした体内時刻情報を管理サーバコンピュータ５２に入力可能に構成されている。
後述する検査対象者番号を、管理サーバコンピュータ５２に登録された滞在者通し番号と同じ番号にしておくと共に、サーバコンピュータ１１から端末コンピュータ５１にダウンロードした検査対象者毎の体内時刻情報を自動で管理サーバコンピュータ５２の不図示のデータベースに取込むプログラムを格納しておけば、データベースに体内時刻情報を直接自動で取込むよう構成できる。
また、老人保健施設５のほかに、老人病院、老人福祉施設、老人ホーム、在宅老人介護の支援業者等が顧客であってもよい。
法人会員６は、従業員の健康管理、業務効率向上等のために、本実施形態の体内時刻推定システムを利用する。
法人会員６には、端末コンピュータ６１のほかに、勤務時間、勤務時間帯、作業内容、健康状態、年齢等が電子データとして従業員毎に登録された従業員情報サーバコンピュータ６２が設置されている。従業員情報サーバコンピュータ６２は、ＬＡＮで端末コンピュータ６１と接続され、サーバコンピュータ１１から端末コンピュータ６１にダウンロードした体内時刻情報を従業員情報サーバコンピュータ６２に入力可能に構成されている。
後述する検査対象者番号を、従業員情報サーバコンピュータ６２に登録された従業員番号と同じ番号にしておくと共に、サーバコンピュータ１１から端末コンピュータ６１にダウンロードした検査対象者毎の体内時刻情報を自動で従業員情報サーバコンピュータ６２の不図示のカルテデータベースに取込むプログラムを格納しておけば、従業員情報に体内時刻情報を直接自動で取込むよう構成できる。
また、法人会員６の代わりに、官公庁及びその国内外の出先機関等が顧客であってもよい。
教育機関７は、児童・生徒・学生の健康管理、精神衛生、教育効率向上等のために、本実施形態の体内時刻推定システムを利用する。
教育機関７には、端末コンピュータ７１のほかに、成績、家庭環境、健康状態、年齢等が電子データとして毎に登録された児童・生徒・学生情報サーバコンピュータ７２が設置されている。児童・生徒・学生情報サーバコンピュータ７２は、ＬＡＮで端末コンピュータ７１と接続され、サーバコンピュータ１１から端末コンピュータ７１にダウンロードした体内時刻情報を児童・生徒・学生情報サーバコンピュータ７２に入力可能に構成されている。
後述する検査対象者番号を、児童・生徒・学生情報サーバコンピュータ７２に登録された児童・生徒・学生番号と同じ番号にしておくと共に、サーバコンピュータ１１から端末コンピュータ７１にダウンロードした検査対象者毎の体内時刻情報を自動で児童・生徒・学生情報サーバコンピュータ７２の不図示のカルテデータベースに取込むプログラムを格納しておけば、児童・生徒・学生情報に体内時刻情報を直接自動で取込むよう構成できる。
また、教育機関７の代わりに、大学、研究所等の研究機関等が顧客であってもよい。
顧客２〜７は、検査対象者の血液を採取して、情報センター１から配布された試験管に入れ、血液サンプル１９とする。この試験管には、サンプルを入れた段階で、顧客固有の顧客番号、検査対象者の番号、採取年月日・時刻を記入する。顧客２〜７は、血液サンプル１９を、送付又は持込みにより情報センター１に提出する。
血液サンプル１９提出後１〜２日で、顧客２〜７は、端末コンピュータ２１〜７１で、検査対象者が概日リズム障害であるかの判定結果と、概日リズム障害でない場合の体内時刻とを、閲覧可能である。
一方、顧客２〜７側に遺伝子発現産物量測定装置が設置されている場合、血液サンプル１９の発現ｍＲＮＡ量を顧客２〜７側で測定し、この測定データを、予め付与されている顧客番号、検査対象者の番号、採取年月日・時刻等の情報と共に端末コンピュータ２１〜７１からサーバコンピュータ１１に送信する。
測定データ送信後数時間以内で、顧客２〜７は、端末コンピュータ２１〜７１で、検査対象者が概日リズム障害であるかの判定結果、体内時刻情報とを閲覧可能となる。
次いで、各装置、コンピュータについて説明する。
ＤＮＡチップ読取装置１２は、公知のＤＮＡチップ読取装置からなり、不図示のＤＮＡチップのスキャナと、スキャナで読み取ったデータの加工を行う不図示のＣＰＵと、加工されたデータを記憶する記憶装置とを備えている。
ＤＮＡチップ読取装置１２は、検体に含まれる全ＲＮＡ中の個々のｍＲＮＡの量を測定し、測定データを不図示の記憶装置に登録する。
本実施形態では、遺伝子発現産物量測定にＤＮＡチップ読取装置１２を用いているが、定量ＰＣＲ装置，リアルタイムＰＣＲ装置，ＤＮＡマイクロアレイ装置，ＲＮａｓｅプロテクションアッセイ装置，ノザンハイブリダイゼーション装置を用いてもよい。
また、遺伝子発現産物量として蛋白質量の測定データを利用する場合には、公知の二次元泳動装置，マススペクトル装置，蛋白質チップ装置、抗体チップ装置、イムノブロット装置等を用いる。
サーバコンピュータ１１は、遺伝子発現産物量測定データを用い、検体の概日リズム障害有無の判定、検体に係る検査対象者の体内時刻推定を行う装置である。
サーバコンピュータ１１は、図２に示すように、中央演算処理装置であるＣＰＵ７０と、記憶装置としてのＲＡＭ７３、ＲＯＭ７４、ＨＤＤ７５、記憶媒体装置７６と、通信装置７７、キーボード７８、マウス７９、表示装置８０、プリンタ８１を備えている。
ＣＰＵ７０は、記憶装置、通信装置や入力装置から受け取った情報を用いて演算し、演算した結果を記憶装置、通信装置や出力装置に渡す役目をはたす。
記憶装置は各種処理を実行するためのプログラム等を記憶する。記憶装置のうち、記憶媒体装置７６は外付のＨＤＤ、ＣＤ、ＤＶＤ等の媒体装置により構成され、通信装置７７を介して受信したデータ、キーボード７８やマウス７９から入力されるデータを適宜記憶、読出し可能である。
ＲＡＭ７３には、ＣＰＵ７０が処理を実行するために必要なデータ等が記憶される。表示装置８０は、ＣＰＵ７０がＲＯＭ７４、ＨＤＤ７５等に記憶されたプログラムにより作成した画像等を表示する。プリンタ８１は各コンピュータから所定の情報を紙面上に出力する。通信装置７７は、ＤＮＡチップ読取装置１２、端末コンピュータ２１〜７１等他のコンピュータとのデータ送受信を行う。
ＨＤＤ７５には、図３の経時的遺伝子発現産物量テーブル１４と、標準検体の経時的遺伝子発現産物量測定データから選択された概日振動遺伝子候補の遺伝子発現産物量の経時的データを含む不図示の概日振動遺伝子テーブルと、図４の分子時刻表テーブル１５と、図５の体内時刻情報テーブル１６と、不図示の遺伝子選択用余弦波情報と、不図示の時刻表作成用余弦波情報とが格納されている。
図３の経時的遺伝子発現産物量テーブル１４とは、標準検体に含まれる遺伝子発現産物量の経時的データが、遺伝子毎に登録されたテーブルである。
図４の分子時刻表テーブル１５とは、サンプルの概日リズム障害判定、体内時刻推定処理に用いるテーブルであって、標準検体中から選択した各概日振動遺伝子の概日発現変動曲線を特徴付ける値、即ちＬＤ条件（１２時間明期、１２時間暗期の２４時間周期の明暗条件）における標準発現量１５３、標準変動量１５４、分子時刻１５５と、ＤＤ条件（恒暗条件）における標準発現量１５６、標準変動量１５７、分子時刻１５８が登録されている。
図５の体内時刻情報テーブル１６とは、検査対象者のサンプルの概日リズム障害判定、体内時刻推定処理の結果を登録するテーブルであって、図１０の体内時刻の推定及び情報提供処理のフローチャートを実行することにより各レコードが登録される。
体内時刻情報テーブル１６には、血液サンプル１９受領時に情報センター１で検体毎に付される受付番号１６１、顧客の会員登録時に顧客毎に付される顧客番号１６２、顧客担当者のユーザＩＤ１６３、顧客の担当者のユーザパスワード１６４、顧客側で検査対象者毎に付される検査対象者番号１６５、顧客側で血液サンプル１９を採取時にサンプル１９の試験管等に記載され、情報センター１側で入力される採取時刻１６６、サーバコンピュータ１１による概日リズム障害の有無の判定結果１６７、サーバコンピュータ１１による体内時刻の推定結果１６８、情報センター１側で必要に応じて検体毎に入力されるコメント１６９が登録される。
遺伝子選択用余弦波情報とは、標準検体に含まれる遺伝子中から概日振動遺伝子を選択するために用いられる余弦波の情報である。
遺伝子選択用余弦波とは、周期時間が２０時間から２８時間まで１時間刻みで異なると共に、６０分の１刻みで位相のずれた合計５４０種類の余弦波をいう。即ち、９種類の周期時間のそれぞれについて、位相の異なる６０種類の余弦波を作成したものである。この遺伝子選択用余弦波の一部を、図８に示す。
遺伝子選択用余弦波情報には、遺伝子選択用余弦波上の値ｚ（ｔ）を求めるために用いられるＡ（ｔ）の式
2004012128.tif 000008
周期時間：２０，２１，２２，２３，２４，２５，２６，２７，２８の９種類
Ｔ：周期時間を６０等分したときの各時刻であって且つ余弦波が最大値を示す時刻
ｔ：ある時刻
が含まれている。
時刻表作成用余弦波情報とは、標準検体に含まれる遺伝子中から選択された概日振動遺伝子の分子時刻表を作成するために用いられる余弦波の情報である。
時刻表作成用余弦波とは、２４時間周期を持ち１０分刻みで位相のずれた１４４種類の余弦波をいう。
時刻表作成用余弦波情報には、時刻表作成用余弦波上の値α（ｔ）を求めるために用いられるＡ（ｔ）の式
2004012128.tif 000009
周期時間：２４（時間）
Ｔ：余弦波の分子時刻（例えば、０／６０，１０／６０，２０／６０，３０／６０，…，１４２０／６０，１４３０／６０までの１０分刻みの１４４種類の時刻とすることができる）
ｔ：ある時刻
が含まれている。
次いで図６に基づき、本実施形態の体内時刻推定システムの処理の流れについて説明する。
情報センター１は、サーバコンピュータ１１で、インターネット１３を介して顧客向けの体内時刻情報閲覧サービスを提供する。
体内時刻推定システムの処理は、大きく分けて、図６に示す遺伝子選択・時刻表作成段階１００と、概日リズム障害判定・体内時刻推定段階２００と、顧客への報告段階３００とからなる。
遺伝子選択・時刻表作成段階１００の手順は、概日リズム障害判定・体内時刻推定の準備段階というべき手順であり、顧客への体内時刻情報提供サービス開始前に行われる手順である。
概日リズム障害判定・体内時刻推定に使用する組織・細胞等の部位、概日振動遺伝子を選択し、この概日振動遺伝子の分子時刻表を作成する。
遺伝子選択・時刻表作成段階１００では、まず、処理１０１で、標準検体の経時的遺伝子発現データを取得する。このデータ取得は、ＤＮＡチップ読取装置１２で行われる。
標準検体とは、生物の組織を複数の個体から採取したものをいう。本実施形態では、標準検体として、複数のヒトから採取した血液を用いる。
標準検体として、所定条件下の複数の人体から一定時間毎に血液を採取する。具体的には、複数の被験者を１２時間明期、１２時間暗期の２４時間周期の室内に２週間滞在させた後、２日間にわたって同様に１２時間明期、１２時間暗期の２４時間周期の明暗条件下（ＬＤ条件下）、或いは恒暗条件下（ＤＤ条件下）４時間ごとに血液を採取する。
明暗条件下における採取時刻は、点灯直前（ＺＴ０）、点灯４時間後（ＺＴ４）、点灯８時間後（ＺＴ８）、消灯直前（ＺＴ１２）、消灯４時間後（ＺＴ１６）、消灯８時間後（ＺＴ２０）、二日目の点灯直前（ＺＴ２４）、二日目の点灯４時間後（ＺＴ２８）、二日目の点灯８時間後（ＺＴ３２）、二日目の消灯直前（ＺＴ３６）、二日目の消灯４時間後（ＺＴ４０）及び二日目の消灯８時間後（ＺＴ４４）の１２点、恒暗条件下における採取時刻は、恒暗条件開始直前（ＣＴ０時）、主観的な昼の始まりから４時間後（ＣＴ４時）、主観的な昼の始まりから８時間後（ＣＴ８時）、主観的な夜の開始直前（ＣＴ１２時）、主観的な夜の始まりから４時間後（ＣＴ１６時）、主観的な夜の始まりから８時間後（ＣＴ２０時）、二日目の主観的な昼の開始直前（ＣＴ２４時）、二日目の主観的な昼の始まりから４時間後（ＣＴ２８時）、二日目の主観的な昼の始まりから８時間後（ＣＴ３２時）、二日目の主観的な夜の開始直前（ＣＴ３６時）、二日目の主観的な夜の始まりから４時間後（ＣＴ４０時）及び二日目の主観的な夜の始まりから８時間後（ＣＴ４４時）の１２点とする。
複数人数から採取した血液を公知の方法で個体毎条件毎採取時刻毎に、遺伝子発現産物量であるｍＲＮＡ量測定用の標準検体を調製する。
この個体毎条件毎採取時刻毎の標準検体中のｍＲＮＡ量を、ＤＮＡチップ読取装置１２を用いて測定し、標準検体の経時的遺伝子発現産物量測定データを得る。
この経時的遺伝子発現産物量測定データは、ＤＮＡチップ読取装置１２の不図示の記憶装置に格納されると共に、ＤＮＡチップ読取装置１２からサーバコンピュータ１１に送信される。
次いで、処理１０２で、処理１０１で受信した経時的遺伝子発現産物量測定データに基づき、サーバコンピュータ１１が、概日振動遺伝子選択を行う。
この概日振動遺伝子選択の処理を、図７の概日振動遺伝子選択フローチャートに基づき詳細に説明する。図７の処理は、サーバコンピュータ１１のＣＰＵ７０で制御される。
サーバコンピュータ１１は、ＤＮＡチップ読取装置１２から経時的遺伝子発現産物量測定データを受信すると、このデータをＨＤＤ７５に経時的遺伝子発現産物量テーブル１４として登録する。経時的遺伝子発現産物量テーブル１４を図３に示す。
経時的遺伝子発現産物量テーブル１４には、各遺伝子の各時刻における発現量であるｍＲＮＡ量が、遺伝子毎に登録されている。
ＨＤＤ７５に経時的遺伝子発現産物量テーブル１４が登録されると、図７のフローチャートがスタートする。まず、ステップＳ１０１で、経時的遺伝子発現産物量テーブル１４の遺伝子レコードを一件読込む。
次いで、ステップＳ１０２で、ＨＤＤ７５から取得した遺伝子選択用余弦波情報の式Ａ（ｔ）に、それぞれの遺伝子の標準変動量、標準発現量を適用した式ｚ（ｔ）を導き出し、ｚ（ｔ）
2004012128.tif 000010
周期時間：２０，２１，２２，２３，２４，２５，２６，２７，２８の９種類
Ａ：遺伝子Ａの標準変動量
Ｂ：上記数４及び数５により計算した余弦波の標準偏差
Ｃ：遺伝子Ａの標準発現量
Ｔ：周期時間を６０等分したときの各時刻であって且つ余弦波が最大値を示す時刻
ｔ：ある時刻
に周期時刻、Ｔをそれぞれ１件代入し、５４０種類の余弦波の中から１つの余弦波を特定する。
次いで、ステップＳ１０３で、ステップＳ１０１で読込んだ遺伝子レコードの各時刻の値と、対応する４時間刻みの時刻を式ｚ（ｔ）のｔに代入して得た値との間のピアソンの積率相関係数を算出し、算出結果をＲＡＭ７３に保存する。
即ち、遺伝子レコードの時刻４ｎ時の値ｘｎ（ＺＴ４ｎの発現量又はＣＴ４ｎの発現量（ｎは、０以上１１以下の整数））と、時刻４ｎ時の余弦波上の値ｙｎ（ＺＴ４ｎの余弦波上の値又はＣＴ４ｎの余弦波上の値（ｎは、０以上１１以下の整数））との間のピアソンの積率相関係数を算出する。
次いで、ステップＳ１０４で、まだ積率相関係数を算出していない周期時刻、Ｔがあるか判定する。
まだ積率相関係数を算出していない周期時刻、Ｔがある場合（ステップＳ１０４：ＹＥＳ）、ステップＳ１０２で、式ｚ（ｔ）に周期時刻、Ｔをそれぞれ１件代入し、５４０種類の余弦波の中から１つの余弦波を特定する。
積率相関係数を算出していない周期時刻、Ｔがない場合（ステップＳ１０４：ＮＯ）、その遺伝子レコードについては、すべての遺伝子作成用余弦波について積率相関係数の算出・保存が完了しているものとして、ステップＳ１０５で、その遺伝子レコードについて算出されたすべての積率相関係数中に、係数＞０．７５を満たすものがあるか判定する。
係数＞０．７５を満たすものがある場合（ステップＳ１０５：ＹＥＳ）、その遺伝子は、発現量が周期時刻２０時間以上２８時間以下の範囲内で振動しており、概日振動遺伝子としての特性を有しているものとして、ステップＳ１０６で、その遺伝子レコードに遺伝子候補フラグをセットする。
次いで、ステップＳ１０７で、まだ遺伝子選択用余弦波との積率相関係数算出処理を行っていない遺伝子レコードがあるか判定する。積率相関係数算出処理は、各遺伝子について、ＬＤ条件下、ＤＤ条件下における遺伝子レコード双方について行う。
まだ相関係数算出処理を行っていない遺伝子レコードがある場合（ステップＳ１０７：ＹＥＳ）、ステップＳ１０１で、遺伝子レコードを１件読込む。
相関係数算出処理を行っていない遺伝子レコードがない場合（ステップＳ１０７：ＮＯ）、すべての遺伝子について概日振動遺伝子候補の判定を完了したものとして、ステップＳ１０８で、概日振動遺伝子候補を含む不図示の概日振動遺伝子テーブルを生成する。このステップＳ１０８では、経時的遺伝子発現産物量テーブル１４を参照し、ＬＤ条件下及びＤＤ条件下の双方について遺伝子候補フラグがセットされている遺伝子レコードのみを抽出して、不図示の概日振動遺伝子テーブルを生成する。その後、処理を終了する。
以上の処理で、図６の処理１０２の概日振動遺伝子選択が完了する。
図６の処理１０２の概日振動遺伝子選択が完了すると、処理１０３で、サーバコンピュータ１１は、各遺伝子産物の概日発現変動曲線・分子時刻表作成を行う。
この処理を、図９の概日発現変動曲線・分子時刻表作成フローチャートに基づき説明する。図９の処理は、サーバコンピュータ１１のＣＰＵ７０で制御される。
図６の処理１０２の概日振動遺伝子選択が完了すると、図９のフローチャートがスタートする。
まず、ステップＳ２０１で、図７のステップＳ１０８で生成された不図示の概日振動遺伝子テーブルから、概日振動遺伝子レコードをＬＤ条件、ＤＤ条件の条件毎に１件読込む。
次いで、ステップＳ２０２で、ＨＤＤ７５から取得した時刻表作成用余弦波情報の式Ａ（ｔ）にそれぞれの遺伝子の標準変動量，標準発現量を適用した式α（ｔ）
2004012128.tif 000011
Ａ：概日振動遺伝子Ａの標準変動量
Ｂ：上記数４及び数５により計算した余弦波の標準偏差
Ｃ：遺伝子Ａの標準発現量
Ｔ：余弦波の分子時刻（例えば、０／６０，１０／６０，２０／６０，３０／６０，…，１４２０／６０，１４３０／６０までの１０分刻みの１４４種類の時刻とすることができる）
ｔ：ある時刻
を導き出し、α（ｔ）にＴを１件代入し、１４４種類の余弦波の中から１つの余弦波を特定する。
次いで、ステップＳ２０３で、ステップＳ２０１で読込んだ概日振動遺伝子レコードの各時刻の値と、対応する４時間刻みの時刻を式α（ｔ）のｔに代入して得た値との間のピアソンの積率相関係数を算出し、算出結果をＲＡＭ７３に保存する。
即ち、概日振動遺伝子レコードの時刻４ｎ時の値ｘｎ（ＺＴ４ｎの発現量又はＣＴ４ｎの発現量（ｎは、０以上１１以下の整数））と、時刻４ｎ時の余弦波上の値ｙｎ（ＺＴ４ｎの余弦波上の値又はＣＴ４ｎの余弦波上の値（ｎは、０以上１１以下の整数））との間のピアソンの積率相関係数を算出する。
次いで、ステップＳ２０４で、まだ積率相関係数を算出していないＴがあるか判定する。
まだ積率相関係数を算出していないＴがある場合（ステップＳ２０４：ＹＥＳ）、ステップＳ２０２で、式α（ｔ）にＴを１件代入し、１４４種類の余弦波の中から１つの余弦波を特定する。
積率相関係数を算出していないＴがない場合（ステップＳ２０４：ＮＯ）、その条件のその概日振動遺伝子レコードについては、すべての時刻表作成用余弦波について積率相関係数の算出・保存が完了しているものとして、ステップＳ２０５で、その条件のその概日振動遺伝子レコードについて算出されたすべての積率相関係数中で最大の値になる余弦波を抽出し、ＲＡＭ７３に保存する。この余弦波が、概日振動遺伝子の発現振動を示す概日発現変動曲線となる。
次いで、ステップＳ２０６で、抽出した余弦波である概日発現変動曲線上の値が最大となる時刻を分子時刻として、概日振動遺伝子名及び条件名、標準発現量、標準変動量と共に図４の分子時刻表テーブル１５に登録する。
この標準発現量１５３、１５６、標準変動量１５４、１５７、分子時刻１５５、１５８は、概日発現変動曲線を特徴付けるものである。
次いで、ステップＳ２０７で、まだ時刻表作成用余弦波との積率相関係数算出処理を行っていない概日振動遺伝子レコードの条件があるか判定する。
まだ時刻表作成用余弦波との積率相関係数算出処理を行っていない概日振動遺伝子レコードの条件がある場合（ステップＳ２０７：ＹＥＳ）、ステップＳ２０１で、図７のステップＳ１０８で生成された不図示の概日振動遺伝子テーブルから、概日振動遺伝子レコードをＬＤ条件、ＤＤ条件の条件毎に１件読込む。
時刻表作成用余弦波との積率相関係数算出処理を行っていない概日振動遺伝子レコードの条件がない場合（ステップＳ２０７：ＮＯ）、すべての概日振動遺伝子について概日発現変動曲線及び分子時刻表の作成が完了したものとして、処理を終了する。
以上で、図６の処理１０３の各遺伝子産物の概日発現変動曲線・分子時刻表作成が完了する。
これにより、遺伝子選択・時刻表作成段階１００が完了して、顧客への体内時刻情報提供サービスが可能な状態となる。
遺伝子選択・時刻表作成段階１００が完了すると、情報センター１は、体内時刻情報提供サービスを開始する。
顧客２〜７のユーザは、本実施形態に係る体内時刻推定システムを利用する前には、予め情報センター１に会員登録をしておく。
会員登録は、サーバコンピュータ１１から提供される会員登録画面で、会員団体名、ユーザ名、ユーザのパスワード、ユーザの電子メールアドレス、連絡先、サービス利用料金の支払方法等の情報を入力することにより行う。
本実施形態では、体内時刻情報提供サービスは、予め会員に登録した者に対するサービスとするが、一般の小売店等で、血液サンプル１９を入れるための試験管と、情報センター１のサーバコンピュータ１１に接続するための説明書をセットにしたキットを販売し、キットを購入した消費者が単発でサービスを受けるようにしてもよい。
顧客２〜７は、情報センター１から受領した試験管に患者の血液を入れ、血液サンプル１９を情報センター１に送付する。情報センター１は、ＤＮＡチップ読取装置１２で、受領した血液サンプル１９の発現ｍＲＮＡ量を公知の方法で測定する。
情報センター１は、血液サンプル１９の発現量測定が完了すると、サーバコンピュータ１１を作動し、概日リズム障害、体内時刻の推定及び情報提供処理を開始する。
概日リズム障害、体内時刻の推定及び情報提供処理において、サーバコンピュータ１１が実行する処理を図１０〜図１３のフローチャートに基づき説明する。
処理がスタートすると、まず、ステップＳ３０１で、初期設定が行われる。ステップＳ３０１の処理の詳細を、図１１のフローチャートに基づき説明する。
ステップＳ３０１の初期設定の処理では、まず、ステップＳ３５１で、サーバコンピュータ１１のユーザによりキーボード７８等で入力された閾値Ｄを入力し、ＨＤＤ７５の不図示の閾値メモリに登録する。
ここで、閾値Ｄとは、サンプル中の概日振動遺伝子発現産物量の波と概日振動遺伝子産物の推定発現量の波とのピアソン積率相関係数の最大値Ｃがどの値未満であれば概日リズム障害と判定するかを示す閾値である。
この閾値Ｄは、サンプル中の概日振動遺伝子の数によって異なる値とすると好適である。例えば、サンプル中の概日振動遺伝子の数が増加するほど閾値Ｄは小さい値に設定する。サンプル中の概日振動遺伝子の数が１００〜２００程度の場合、閾値Ｄを０．３以上、好ましくは０．５以上の値とするとよい。
本実施形態では、閾値Ｄを０．５とする。
次いで、ステップＳ３５２で、サーバコンピュータ１１のユーザによりキーボード７８等で入力された閾値Ｕを入力し、ＨＤＤ７５の不図示の閾値メモリに登録する。
ここで、閾値Ｕとは、サンプルの概日振動遺伝子について推定された体内時刻と環境時刻とがどの程度の時間ずれていた場合に概日リズム障害と判定するかを示す閾値であり、本実施形態では２時間とする。
次いで、ステップＳ３５３で、標準検体の各遺伝子産物の標準発現量・標準変動量・分子時刻を入力する。このステップでは、図６の処理１０３の各遺伝子産物の概日発現変動曲線・分子時刻表作成処理で作成された分子時刻表のうち、今回測定するサンプルの部位に対応するものを、分子時刻表テーブル１５から抽出し、ＲＯＭ７４に登録する。
次いで、ステップＳ３５４で、ステップＳ３５３で登録した今回測定するサンプルの部位に対応する分子時刻表から標準発現量・標準変動量・分子時刻の値を抽出し、標準検体の概日発現変動曲線の式
2004012128.tif 000012
Ａ：概日振動遺伝子Ａの標準変動量
Ｂ：上記数４及び数５により計算した余弦波の標準偏差
Ｃ：遺伝子Ａの標準発現量
Ｔ：余弦波の分子時刻（例えば、０／６０，１０／６０，２０／６０，３０／６０，…，１４２０／６０，１４３０／６０までの１０分刻みの１４４種類の時刻とすることができる）
ｔ：ある時刻
を導き出す。
次いで、ステップＳ３５５で、この概日発現変動曲線を正規化する。つまり、標準検体の各概日振動遺伝子の概日発現変動曲線から標準発現量を減じ、次いで標準変動量で除すことにより標準検体の各概日振動遺伝子の相対概日発現変動曲線の式
2004012128.tif 000013
Ｂ：上記数４及び数５により計算した余弦波の標準偏差
Ｔ：余弦波の分子時刻（例えば、０／６０，１０／６０，２０／６０，３０／６０，…，１４２０／６０，１４３０／６０までの１０分刻みの１４４種類の時刻とすることができる）
ｔ：ある時刻
を導き出す。この相対概日発現変動曲線をＲＯＭ７４に登録する。
以上で、図１０のステップＳ３０１の初期設定処理を終了する。
次いで、図１０のステップＳ３０２で、サーバコンピュータ１１のユーザにより入力された今回測定するサンプルのサンプリング時刻ＳをＲＡＭ７３に登録する。このサンプリング時刻Ｓは、環境時刻である。
次いで、ステップＳ３０３で、サンプルの概日振動遺伝子発現産物量測定データをＨＤＤ７５に登録する。このステップは、サーバコンピュータ１１のユーザの操作により、ＤＮＡチップ読取装置１２からサンプルの概日振動遺伝子発現産物量測定データを受信することにより行う。
次いで、ステップＳ３０４で、サンプルの概日振動遺伝子発現産物量の正規化を行い、サンプルの各概日振動遺伝子産物の相対発現量を求める。このステップＳ３０４の処理について、図１２に基づき詳細に説明する。
まず、ステップＳ３６１で、サンプルの概日振動遺伝子発現産物量から標準検体の標準発現量を減ずる。この標準発現量は、遺伝子毎に標準検体中の概日振動遺伝子産物の発現量の平均値を算出することにより求めた値である。
次いで、ステップＳ３６２で、ステップＳ３６１で得た値を、標準検体の標準変動量で除す。この標準変動量は、遺伝子毎に標準検体中の概日振動遺伝子産物の発現量の標準偏差を算出することにより求めた値である。
以上で、ステップＳ３０４のサンプルの概日振動遺伝子発現産物量の正規化処理が完了し、サンプルの各概日振動遺伝子産物の相対発現量が求められる。
次いで、ステップＳ３０５で、標準検体の相対概日発現変動曲線の式β（ｔ）のｔに０を代入し、各概日振動遺伝子の相対概日発現変動曲線上の時刻ｔ＝０における値を算出する。
次に、ステップＳ３０６で、標準検体の相対概日発現変動曲線の式β（ｔ）のｔに代入した値が２４より大きいか判定する。
標準検体の相対概日発現変動曲線の式β（ｔ）のｔに代入した値が２４より大きくない場合（ステップＳ３０６：ＮＯ）、つまり、標準検体の相対概日発現変動曲線の式β（ｔ）のｔに代入した値が０時以上２４時以下である場合、ステップＳ３０７で、正規化された概日振動遺伝子発現産物量と時刻ｔ時における各概日振動遺伝子産物の推定相対発現量との間のピアソン相関係数ｃの計算を行う。
このステップＳ３０７の処理について、図１３に基づき詳細に説明する。
まず、ステップＳ３７１で、各概日振動遺伝子の相対概日発現変動曲線上の時刻ｔ時における推定相対発現量を取得する。この処理では、分子時刻表の各概日振動遺伝子に関するレコードを１件ずつ読込み、概日振動遺伝子の相対概日発現変動曲線を求める。次いで、相対概日発現変動曲線上の時刻ｔ時における値を取得し、推定相対発現量とする。この処理をすべての概日振動遺伝子について行う。
次いで、ステップＳ３７２で、ステップＳ３０４で正規化して得たサンプリング時刻におけるサンプルの概日振動遺伝子産物の相対発現量と、各概日振動遺伝子の相対概日発現変動曲線上の時刻ｔ時における推定相対発現量との間のピアソン積率相関係数を計算し、得られた値を、時刻ｔの値と関連付けてＲＡＭ７３に登録し、図１３の処理を完了する。
次いで、図１０のステップＳ３０８に進み、標準検体の相対概日発現変動曲線の式β（ｔ）のｔに、直前のｔにｄｔを足したｔ＋ｄｔを代入し、各概日振動遺伝子の相対概日発現変動曲線上の時刻ｔ＝ｔ＋ｄｔにおける値を算出する。
なお、本実施形態では、このｄｔは、１０分、即ち１０／６０時間とする。
次いで、ステップＳ３０６に進み、標準検体の相対概日発現変動曲線の式β（ｔ）のｔに代入した値が２４より大きいか判定する。
標準検体の相対概日発現変動曲線の式β（ｔ）のｔに代入した値が２４より大きい場合（ステップＳ３０６：ＹＥＳ）、サンプルの各概日振動遺伝子の相対発現量と、０時から２４時までの間のすべての時刻についての各概日振動遺伝子の相対概日発現変動曲線上の推定相対発現量との間のピアソン積率相関係数の算出が完了しているものとして、ステップＳ３０９に進み、ピアソン積率相関係数ｃの最大値Ｃ及び最大値Ｃを与える体内時刻Ｔの探索を行う。
このステップでは、ＲＡＭ７３に登録されたピアソン積率相関係数ｃのデータ群を取得し、このデータ群中から、最大のピアソン積率相関係数値を抽出して最大値Ｃとする。次いで、ピアソン積率相関係数最大値Ｃに対応する時刻ｔを取得して、最大値Ｃを与える体内時刻Ｔとする。この最大値Ｃ及び体内時刻Ｔを、サンプルの最大値Ｃ及び体内時刻Ｔとして、ＲＡＭ７３に登録する。
次いで、ステップＳ３１０で、ＲＡＭ７３に登録した最大値Ｃが、Ｃ＜閾値Ｄを満たしているか判定する。本実施形態では、閾値Ｄが０．５として登録されているので、このステップでは、最大値Ｃが０．５未満であるか判定する。
最大値Ｃ＜閾値Ｄである場合（ステップＳ３１０：ＹＥＳ）、サンプルのサンプリング時刻における相対概日振動遺伝子発現産物量は、いずれの時刻における各概日振動遺伝子の推定相対発現量とも類似しないことから、このサンプルに係る検査対象者は、概日振動遺伝子が正常な概日振動をしていないとして、ステップＳ３１１で、このサンプルについて、概日リズム障害と判定すると共に、体内時刻情報テーブル１６の該当する受付番号の概日リズム障害の有無の判定結果１６７に、リズム障害であることを示す「１」を登録する。
最大値Ｃ＜閾値Ｄでない場合（ステップＳ３１０：ＮＯ）、即ち、最大値Ｃが閾値Ｄの値と同じかそれ以上である場合、サンプルのサンプリング時刻における相対概日振動遺伝子発現産物量と各概日振動遺伝子の推定相対発現量とが類似する時刻Ｔが存在することから、このサンプルに係る検査対象者は、概日振動遺伝子が正常な概日振動をしている可能性があるとして、ステップＳ３１２に進み、体内時刻Ｔとサンプリング時刻Ｓとの差の絶対値が
｜Ｔ−Ｓ｜＞閾値Ｕ
であるかを判定する。ここでは、閾値Ｕが２時間として登録されているので、体内時刻Ｔとサンプリング時刻Ｓの差の絶対値が２時間より大きいか判定する。
体内時刻Ｔとサンプリング時刻Ｓとの差の絶対値が閾値Ｕより大きい場合（ステップＳ３１２：ＹＥＳ）、環境時刻であるサンプリング時刻Ｓとの時間のずれが、本発明に係る体内時刻推定法の一般的な誤差を考慮しても大きすぎるとして、ステップＳ３１３で、このサンプルについて、概日リズム障害と判定すると共に、体内時刻情報テーブル１６の該当する受付番号の概日リズム障害の有無の判定結果１６７に、リズム障害であることを示す「１」を登録する。
次いで、ステップＳ３１４で、ステップＳ３０９で算出された体内時刻Ｔを体内時刻情報テーブル１６の該当する受付番号の体内時刻の推定結果１６８に登録し、処理を終了する。
体内時刻Ｔとサンプリング時刻Ｓとの差の絶対値が閾値Ｕより大きくない場合（ステップＳ３１２：ＮＯ）、即ち、体内時刻Ｔとサンプリング時刻Ｓとの差の絶対値が閾値Ｕ以下である場合、環境時刻であるサンプリング時刻Ｓとの時間のずれは正常な範囲内であり、Ｔを検査対象者の体内時刻とすることに問題がないとして、体内時刻Ｔを体内時刻情報テーブル１６の該当する受付番号の体内時刻の推定結果１６８に登録し、処理を終了する。
以上で、概日リズム障害判定処理及び体内時刻推定処理を終了する。
次に、顧客２〜７から、サーバコンピュータ１１に登録された概日リズム障害の有無及び体内時刻情報を閲覧する処理について説明する。
検査対象者の血液サンプル１９を情報センター１に送付しておいた顧客２〜７のユーザは、図６の概日リズム障害判定、体内時刻推定段階２００が完了すると、サーバコンピュータ１１から、概日リズム障害判定及び体内時刻推定処理が完了した旨の電子メールを受信する。
ユーザは、この電子メールにより概日リズム障害判定及び体内時刻推定処理が完了した旨を知ると、サーバコンピュータ１１で提供される体内時刻情報報告画面９１で、図６の体内時刻情報の報告３００を受ける。
この体内時刻情報報告画面９１で、図６の体内時刻情報の報告をする処理について説明する。
端末コンピュータ２１で、インターネットで情報センター１の体内時刻情報提供サービス画面のアドレスを入力すると、不図示の体内時刻情報提供サービスの初期画面が表示される。
この初期画面には、体内時刻や投薬時間管理等に関する一般的な情報のほか、体内時刻情報提供サービスの会員専用画面に移行するためのＩＤ，パスワード入力画面表示ボタンが表示されている。
ＩＤ，パスワード入力画面表示ボタンをクリックすると、不図示のＩＤ，パスワード入力画面が表示される。この画面でＩＤ，パスワードが入力されると、サーバコンピュータ１１は、体内時刻情報テーブル１６を、その顧客の顧客番号１６２をキーとして検索し、その顧客のレコードを抽出して、体内時刻情報報告画面９１用のデータを作成する。
その後、この体内時刻情報報告画面９１用のデータを端末コンピュータ２１側に送信する。
端末コンピュータ２１では、そのＩＤの会員専用の画面にログインし、図１４に示す体内時刻情報報告画面９１が表示される。
体内時刻情報報告画面９１には、ログインした会員がリズム障害判定及び体内時刻推定を依頼した検体の体内時刻情報リスト９１０、体内時刻情報報告画面９１を閉じるための閉じるボタン９２１、ダウンロードチェック欄９１９にチェックした検体に関する判定・推定結果をダウンロードするためのダウンロードボタン９２２、ログアウトボタン９２３が表示される。
体内時刻情報リスト９１０は、検体毎の体内時刻情報のリストであって、検体毎に、その検体の受付番号９１１、その検体の検査対象者番号９１２、その検体の採取時刻９１３、その検体を情報センター１が受領した日付を示す受領日９１４、その検体の検査結果であるリズム障害有無９１５、その検体の検査結果である体内時刻情報９１６、その検査対象者番号の体内時刻情報の履歴画面表示ボタン９１７、その検体の検査結果に関するコメント表示ボタン９１８、その検体の検査結果をダウンロードするためのダウンロードチェック欄９１９が表示されている。
ユーザは、履歴画面表示ボタン９１７をクリックすることにより、リストに表示された検体の検査対象者の過去の結果の履歴を閲覧可能である。
なお、サーバコンピュータ１１は、体内時刻推定システム提供で得た多数の検査対象者の体内時刻情報をデータ解析するプログラムを格納していてもよい。
このプログラムにより、年代毎、性別毎、人種毎、罹病疾患毎、身長・体重等の特性毎、業務形態毎等の体内時刻情報を算出し、ヒトの体内時刻と人の特性との相関を解析してもよい。解析して得たヒトの体内時刻と人の特性との相関の解析結果は、人類の健康増進や疾患への対策等のために利用可能である。
体内時刻推定システムを提供することにより得た多数の検査対象者の体内時刻情報をデータ解析する処理について説明する。
情報センター１は、血液サンプル１９受領時に、検査対象者番号１６５と共に、検査対象者の生年月日、性別、人種、罹病疾患、身長・体重等の特性データを受領しておく。
このデータは、血液サンプル１９受領時にインターネット１３を介して受信するとよいが、血液サンプル１９と共にデータが記載された用紙を受領してもよい。
情報センター１が受領した各検査対象者の特性データは、サーバコンピュータ１１の図１５の検査対象者特性テーブル１７に登録する。
体内時刻情報テーブル１６及び検査対象者特性テーブル１７に所定件数以上のレコードが蓄積されると、サーバコンピュータ１１で、データの解析を行う。
データ解析結果の一例を図１６に示す。
データ解析結果は、その後の睡眠覚醒リズム障害（睡眠相後退症候群、睡眠相前進症候群、非２４時間睡眠・覚醒症候群など）、季節性うつ病、時差症候群（ＪＥＴ−ＬＡＧ）の診断に用いてもよい。
なお、概日発現変動曲線は、有限の発現データから概日振動遺伝子の経時的な発現状態を数式化したもので経時的に測定した発現量をもとに作成したものであれば限定されない。
例えば、概日発現変動曲線は、各概日振動遺伝子の経時的な発現量測定データのフーリエ変換による周期曲線の作成法（文献１）によっても作成可能である。「フーリエ変換」とは、フランスの数学、物理学者フーリエが発表した原理（物理的過程は時間ｔの関数ｈ（ｔ）で表現することも、ある周波数ｆの関数Ｈ（ｆ）を用いて表現することもできる）に基づく変換法で２つの表現間を行き来するために使われる式である。数式で表すと次のようになる。
2004012128.tif 000014
有限個の標本値（計測値）に対するフーリエ変換は離散フーリエ変換と呼ばれる。離散フーリエ変換の対象となるＮ個の一連の測定値
2004012128.tif 000015
があるとする。標本化間隔はΔである。ここで離散フーリエ変換は次のように定義できる。
2004012128.tif 000016
このＨｎを用いて時間の関数ｈ（ｔ）を次のように推定することが可能である。
2004012128.tif 000017
また、概日発現変動曲線は、スプライン補間等の補間法を用いた補間曲線の作成法（文献１）によっても作成可能である。「スプライン補間」とは、補間曲線の生成、有限要素法、関数近似、実験データの当てはめなどに利用される補間法である（文献１、文献２）。
フーリエ変換或はスプライン補間による概日発現変動曲線も時間の関数として表記されることからある時間ｔを代入することで該概日発現変動曲線上の値を、また、その値から対応する標準発現量を減じ、標準変動量で除することによって相対発現変動曲線上の値を求めることができる。
また、体内時刻の推定は、６個以上（好ましくは３０個以上、より好ましくは５０個以上、特に好ましくは１００個以上）の任意の概日振動遺伝子群について予め作成した前記相対概日発現変動曲線上のある時刻における値と、別途測定した概日振動遺伝子群の発現産物の相対発現量との間の最も高い類似性を示す時刻を体内時刻として推定することが可能である。例えば、概日発現変動曲線を余弦波で作成し、且つ類似性をピアソンの積率相関係数として算出した場合は、最も高い類似性が０．５以上の相関係数を示したときは、２時間の誤差範囲内で体内時刻の推定が可能である。
一方、すべての時刻にわたって両者に十分な類似性が認められない場合（例えば、ピアソンの積率相関係数として０．５未満）は、概日リズム障害として検出が可能である。類似性は、対象と対象との間の類似の度合いを表す数値である限り限定されないが、ピアソンの積率相関係数が好ましい。
体内時刻の推定に用いる概日振動遺伝子は、本明細書における概日振動遺伝子としての性質を有する限り、公知の概日振動遺伝子でも良いし、ＤＮＡチップ法などで抽出した遺伝子（機能未知の遺伝子を含む）であっても良い。
また上述のようにアプリオリに類似性を定義して体内時刻を推定する方法以外にも例えば、サポートベクターマシン、遺伝的アルゴリズム、ニューラルネット等のような学習アルゴリズムなどを用いて類似性を学習させることによって体内時刻を推定することも可能である。
本明細書中で引用した文献は、文献１が、ニューメリカルレシピ・イン・シー、Ｗｉｌｌｉａｍ  Ｈ．Ｐｒｅｓｓほか、技術評論社、文献２が、Ｃによるスプライン関数、−データ解析、ＣＧ、微分方程式−監修  桜井明、著者  菅野敬祐、吉村和美、高山文雄、東京電機大学出版局、文献３が、クラスター分析とその応用、西田英朗監訳、内田老鶴圃、文献４が、基礎統計学Ｉ、統計学入門  東京大学教養学部統計学教室編  東京大学出版会である。
産業上の利用性
以上のように本発明によれば、測定データを、前記分子時刻表で特定される前記概日発現変動曲線と照合することにより、前記生物個体が概日リズム障害であるかの判定結果及び概日リズム障害でない場合には前記生物個体の体内時刻推定結果を含む体内時刻情報を導出する体内時刻情報導出手段を備えているため、１回の採取で得た検体の測定データと所定生物種の所定部位の分子時刻表により特定される前記概日発現変動曲線とを照合することにより、検体を複数回採取する必要なく、１回の採取で得た検体に基づいて簡易に概日リズム障害か否かの判定及び体内時刻の推定が可能となる。
【図１】
2004012128.tif 000018
【図２】
2004012128.tif 000019
【図３】
2004012128.tif 000020
【図４】
2004012128.tif 000021
【図５】
2004012128.tif 000022
【図６】
2004012128.tif 000023
【図７】
2004012128.tif 000024
【図８】
2004012128.tif 000025
【図９】
2004012128.tif 000026
【図１０】
2004012128.tif 000027
【図１１】
2004012128.tif 000028
【図１２】
2004012128.tif 000029
【図１３】
2004012128.tif 000030
【図１４】
2004012128.tif 000031
【図１５】
2004012128.tif 000032
【図１６】
2004012128.tif 000033
【国際調査報告】
2004012128.tif 000034
2004012128.tif 000035