JP5560550B2 - 生体リズム予測方法 - Google Patents

生体リズム予測方法 Download PDF

Info

Publication number
JP5560550B2
JP5560550B2 JP2008267464A JP2008267464A JP5560550B2 JP 5560550 B2 JP5560550 B2 JP 5560550B2 JP 2008267464 A JP2008267464 A JP 2008267464A JP 2008267464 A JP2008267464 A JP 2008267464A JP 5560550 B2 JP5560550 B2 JP 5560550B2
Authority
JP
Japan
Prior art keywords
expression level
time
gene
clock
circadian
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2008267464A
Other languages
English (en)
Other versions
JP2010094072A (ja
Inventor
温彦 相馬
拓郎 山本
真 明石
孝一 野出
典之 岸井
和博 中川
志功 山下
智広 早川
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sony Corp
Original Assignee
Sony Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority to JP2008267464A priority Critical patent/JP5560550B2/ja
Application filed by Sony Corp filed Critical Sony Corp
Priority to CN2009801401016A priority patent/CN102177254A/zh
Priority to EP09820656.8A priority patent/EP2336355B1/en
Priority to PCT/JP2009/067966 priority patent/WO2010044481A1/ja
Priority to US13/123,847 priority patent/US20110196619A1/en
Priority to BRPI0920426-1A priority patent/BRPI0920426A2/ja
Priority to RU2011113819/10A priority patent/RU2512069C2/ru
Publication of JP2010094072A publication Critical patent/JP2010094072A/ja
Application granted granted Critical
Publication of JP5560550B2 publication Critical patent/JP5560550B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6809—Methods for determination or identification of nucleic acids involving differential detection
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
    • G16B25/10—Gene or protein expression profiling; Expression-ratio estimation or normalisation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B10/00—Other methods or instruments for diagnosis, e.g. instruments for taking a cell sample, for biopsy, for vaccination diagnosis; Sex determination; Ovulation-period determination; Throat striking implements
    • A61B10/02—Instruments for taking cell samples or for biopsy
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/158—Expression markers
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Organic Chemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Biophysics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Biotechnology (AREA)
  • Zoology (AREA)
  • General Health & Medical Sciences (AREA)
  • Wood Science & Technology (AREA)
  • Genetics & Genomics (AREA)
  • Analytical Chemistry (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Immunology (AREA)
  • General Engineering & Computer Science (AREA)
  • Biochemistry (AREA)
  • Evolutionary Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Microbiology (AREA)
  • Public Health (AREA)
  • Artificial Intelligence (AREA)
  • Bioethics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Software Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Epidemiology (AREA)
  • Databases & Information Systems (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Description

本発明は、被検個体の生体リズムを予測する方法に関する。より詳しくは、被検個体からの生体試料の採取を24時間以内に3回のみ行うことで生体リズムを予測する方法に関する。
生物個体の様々な生体現象は、自立的に振動する「周期的なリズム」を示すことが知られている。この周期的なリズムは「生体リズム」と呼ばれている。特に、約一日を周期とする「概日リズム（サーカディアンリズム）」は、睡眠覚醒サイクルや体温、血圧、ホルモン分泌量の日内変動などの生体現象を広く支配していることが知られている。また、概日リズムは、心身の活動度や運動能力、薬剤感受性などにも関与している。
生体リズムは、「時計遺伝子（クロックジーン）」と呼ばれる遺伝子群によって制御されている。時計遺伝子は、その発現や活性、局在等を自律的に周期変動（振動）させることにより「体内時計」として機能している。
時計遺伝子の遺伝子多型や遺伝子変異は、癌や糖尿病、血管系疾患、神経変性疾患などの発症要因となることが明らかにされている。さらに、近年、双極性障害や鬱病のような精神疾患についても、時計遺伝子の遺伝子多型や変異が発症に関与していることが指摘されている。
一方、生体リズムは、体内時計による自律的な制御だけでなく、社会生活による制約も受けている。例えば、睡眠覚醒サイクルでは、日々の就寝時刻や起床時刻の変化によって、「実生活の就寝起床サイクル」と「体内時計による睡眠覚醒サイクル」との間にリズムのずれ（位相のずれ）が生じる場合がある。
このような生体リズムのずれは、いわゆる「時差ぼけ」や睡眠障害を引き起こし、さらには上述のような精神疾患の原因ともなると考えられている。これらの疾患を治療するため、変調した体内時計を光照射によってリセットする試みもなされ始めている。
また、生体リズムを利用して、薬剤治療効果の最大化を図る試みも始まっている。薬剤の標的となる分子（薬剤標的分子）の発現量や薬剤を代謝する酵素（薬物代謝酵素）の活性の概日リズムに起因して、薬剤による治療効果も日内変動することが考えられる。そこで、薬剤ごとに最適な投薬時刻を定めて、治療効果を最大化しようとする「時間医療」という考え方が提唱されてきている。
さらに、より身近には、心身の活動度や運動能力の概日リズムを利用して、学習やトレーニングにおいて自己の能力を最大限に引き出すための活動時刻や、太りにくい（又は、太りやすい）摂食時刻が検討され始めている。
以上のことから、生体リズムを正確に評価することは、時差ぼけなどの体調不良の改善、種々の疾患の予防、時間医療の実現、自己能力の発揮、ダイエットなどに非常に有益と考えられる。
特許文献１には、生物個体から採取した標準検体の遺伝子発現産物量測定データに基づき体内時刻を推定する方法などが開示されている。この体内時刻推定方法では、遺伝子発現産物量（すなわち、mRNA）の発現量に基づいて、体内時刻を推定するための分子時計表を作成するものである。
国際公開第２００４／０１２１２８号
上記特許文献１には、具体的な測定対象遺伝子は記載されていないが、従来、時計遺伝子を測定対象として、その経時的な発現量変化に基づいて生体リズムを推定することが行われている。
しかし、時計遺伝子の経時的な発現量変化を測定するためには、被検個体から継続的に生体試料を採取することが必要とされていた。すなわち、精度良く生体リズムの推定を行うためには、24時間にわたって数時間おきに被検個体から生体試料を採取し、時計遺伝子発現量の時系列データを得る必要があった。
このように一日に何回も生体試料の採取を行うことは、被検者（被検個体）にとって大きな負担となる。また、夜中に被検者を起こして生体試料の採取を行うこともあり、これによって被検者の睡眠覚醒サイクルが影響を受け、生体リズムそのものが変化してしまうおそれもあった。
生体試料の採取のための被検者の負担を軽減することは、時差ぼけなどの体調不良の改善や種々の疾患の予防等を目的とした生体リズムの評価を広く普及させるための重要な課題と考えられる。
そこで、本発明は、被検個体からの生体試料の採取回数を最小限にして、高い精度で生体リズムを予測できる方法を提供することを主な目的とする。
上記課題解決のため、本発明は、（１）発現量変化の概日周期の位相が異なる２つの時計遺伝子について、被検個体から２４時間以内に３回採取された生体試料中の発現量を測定する手順と、（２）前記（１）の手順で得られた時系列発現量データから、前記概日周期を算出する手順と、を含み、前記生体試料は、前記２つの時計遺伝子の時系列発現量データをコサインフィッティングして得たコサインカーブから算出された、１時間毎の発現量変化を示すモデルデータと、このモデルデータから抽出された任意の３点の時刻及び該時刻における発現量、並びに下記式（I）及び式（II）から算出された３点サンプリングデータと、の間で、発現量の最大値を与える時刻の時間差の平均値が０．６未満であり、かつ標準誤差が０．４未満である３点の時刻の時間間隔を採取時間間隔として、２４時間以内に３回採取されたものである、被検個体の生体リズムを予測する方法を提供する。
（式（I）中、Ea(t)、Aa、ω、Caは、一の時計遺伝子の時刻tにおける発現量、振幅、初期位相、オフセット値を示す。また、式（II）中、Eb(t)、Ab、Cbは、他の一の時計遺伝子の時刻tにおける発現量、振幅、オフセット値を示す。さらに、θは両時計遺伝子の位相差を示す。）
この生体リズム予測方法では、前記（２）の手順において、前記時系列発現量データから上記式（I）及び式（II）によって前記概日周期を算出することができる。
前記生体試料は、被検個体から２４時間以内に８時間間隔で３回採取されたものを用いることにより、高精度に生体リズムの予測を行うことができる。
この生体リズム予測方法において、前記時計遺伝子は、Per3遺伝子及びNr1d2遺伝子とすることができる。
一般に「位相」とは、周期的な変動における山や谷といった点の位置・状態をひとつの周期の中に特徴付ける量をいう。時計遺伝子の発現量変化の概日周期は、おおむね下記式（III）で示される余弦波関数として観測することができる。ここで、式中、E(t)は時刻tにおける時計遺伝子の発現量、Aは発現量の振幅、Cはオフセット値を表している。
本発明において「位相」とは、この式（III）のコサインの中を意味するものとし、具体的には式（III）の「2π(t+ω)/24」を指すものとする。また「初期位相」とは、時刻t=0における位相を決定する「ω」を指すものとする。さらに「位相差」とは、時計遺伝子間での初期位相ωの差を指す。
本発明により、被検個体からの生体試料の採取回数を最小限にして、高い精度で生体リズムを予測可能な方法が提供される。
１．生体リズム予測方法
本発明者らは、複数の時計遺伝子間での発現量変化周期のずれを利用して生体リズムの予測を行うことで、被検個体からの生体試料の採取回数を少なくできるのではないかと考え、検討を行った。そして、発現量変化の概日周期の位相が異なる2つの時計遺伝子を測定対象とすることにより、被検個体からの生体試料の採取を24時間以内に3回行うのみで、高い精度で生体リズムを予測できることを見出した。すなわち、本発明に係る生体リズム予測方法は、被検個体から24時間以内に3回採取した生体試料について、発現量変化の概日周期の位相が異なる2つの時計遺伝子の発現量を測定し、得られた時系列発現量データに基づいて、被検個体の生体リズムを予測するものである。
2つの時計遺伝子a, bが、異なる発現量変化の概日周期の位相を有する場合、その位相差を「θ」とすると、各時計遺伝子の発現量変化の概日周期はそれぞれ下記コサインカーブ式（I）及び（II）によってモデル化することができる。
ここで、式（I）中、Ea(t)、Aa、ω、Caは、時計遺伝子aの時刻tにおける発現量、振幅、初期位相、オフセット値を示す。また、式（II）中、Ea(t)、Ab、Cbは、時計遺伝子bの時刻tにおける発現量、振幅、オフセット値を示している。
時計遺伝子a, bの位相差θが既知である場合、上記コサインカーブ式（I）及び（II）は、5つの未知の定数（ω、Aa、Ab、Ca、Cb）を含む式となる。被検個体から時刻tに採取した生体試料中について、時計遺伝子a, bの発現量Ea(t), Ea(t)を測定することにより、コサインカーブ式（I）及び（II）から2つの等式を得ることができる。従って、このモデルを用いることにより、異なる3つの時刻tに生体試料の採取（サンプリング）を行って、コサインカーブ式（I）及び（II）から6つの等式を得れば、時計遺伝子a, bの発現量変化の概日周期を算出することが可能となる。
数学的には5つの未知数を明らかにするには5つの等式があれば足りる。ここでは、位相差θが既知である時計遺伝子a, bを用いて、発現量変化の概日周期を上記コサインカーブ式（I）及び（II）によりモデル化したことによって、6つの等式に基づき5つの未知数を算出する。すなわち、時計遺伝子a, bの位相差θを制約条件とすることで、未知数ω、Aa、Ab、Ca、Cbをより正確に求めることを可能にしている。これにより、一般に正確な評価を行うことが難しい生体現象である時計遺伝子発現量の概日周期を高精度に算出することが可能となる。
この生体リズム予測方法は、より具体的には、（１）被検個体から生体試料を24時間以内に3回採取する手順と、（２）発現量変化の概日周期の位相が異なる2つの時計遺伝子について、生体試料中の発現量を測定する手順と、（３）これらの手順で得られた時系列発現量データから、各時計遺伝子の発現量変化の概日周期を算出する手順と、を含む。
本発明に係る生体リズム予測方法において、対象とする被検個体には、ヒトの他、マウス・ラット・サル等の実験動物などが広く含まれる。
２．生体試料の採取
被検個体から採取する生体試料は、時計遺伝子の遺伝子発現産物（mRNA）が含まれる生体組織であれば特に限定されない。採取作業の簡便さの観点からは、毛髪や口腔粘膜、皮膚等の体表面からの採取が可能な生体組織とすることが好ましい。
毛髪のサンプリングは抜去により行うことができる。抜去された毛髪の毛根部には、毛包細胞が付着しており、この細胞中の時計遺伝子mRNAを測定できる。ここで、「毛包細胞」とは、抜去された体毛の毛根部に付着する内毛根鞘（inner root sheath）、外毛根鞘（outer root sheath）及び毛乳頭（papilla）を形成する細胞群がいうものとする。
毛髪のサンプリング部位は、例えばヒトを対象とする場合、頭髪や髭、腕や足の毛などを用いることができ、特に限定されない。測定のばらつきを抑えるため、各回の採取は近傍の部位で行なうことが望ましい。
一回にサンプリングする毛髪の本数は、ヒトの場合、頭髪では5〜10本、髭では3〜5本、腕毛では10〜20本程度である。この本数以上を使用することで、時計遺伝子の発現定量のために十分な量のmRNAを抽出することができる。
口腔粘膜のサンプリングは、例えば、ブラシやスパーテル等で口腔粘膜表面から掻き取ることによって行うことができる。これにより、掻き取られた口腔粘膜細胞中の時計遺伝子mRNAを測定できる。
口腔粘膜のサンプリング部位は、頬の裏側の粘膜が好適であり、測定のばらつきを抑えるため左右両側の粘膜から採取することが望ましい。
本発明者らは、高精度の生体リズムの予測を可能にするため、3回のサンプリング時刻の時間間隔について検討を行った。その結果、被検個体から24時間以内に8時間間隔で3回生体試料を採取することにより、高い精度でできることを見出した（実施例2参照）。
すなわち、本発明に係る生体リズム予測方法では、第1回目と第2回目のサンプリング時刻との時間間隔と、第2回目と第3回目のサンプリング時刻との時間間隔と、がともに8時間となるように3回のサンプリングを行うことで、最も高精度に生体リズムを予測することができる。
３．発現量の測定
生体試料中の時計遺伝子の発現量は、従来公知の方法によって測定することができる。例えば、市販のRNA抽出キットを使用して生体試料からRNAを抽出し、抽出したRNAを鋳型とした逆転写反応によってcDNAを合成する。そして、このcDNAを用いてDNAマイクロアレイ（DNAチップ）のような網羅的解析手法や、リアルタイムPCRのような個別的解析手法によって、発現量の定量を行う。
測定を行う時計遺伝子は、現在までに同定されている一群の時計遺伝子であってよい。代表的な時計遺伝子としては、Per3遺伝子（NCBI Accession No. NM_016831）、Per2遺伝子（NM_022817）、Bmal1遺伝子（NM_001030272）、Npas2遺伝子（NM_002518）、Nr1d2遺伝子（NM_021724）、Nr1d2遺伝子（NM_005126）、Dbp遺伝子（NM_001352）、Cry1遺伝子(NM_004075)等がある。
本発明に係る生体リズム予測方法では、これらの時計遺伝子のうち、発現量変化の日内変動の位相が異なる2つの遺伝子について発現量を測定し、時系列発現量データを得る。なお、ヒト以外を対象生物とする場合には、上記ヒト時計遺伝子の対象生物におけるホモログ（相同遺伝子）について発現量の測定を行う。
選択される2つの時計遺伝子は、発現量変化の日内変動の位相が異なれば、任意の組合せとできる。選択された任意の2つの時計遺伝子の発現量変化の日内変動の位相は、例えば実施例1において後述する方法によって決定することができる。
好適な2つの時計遺伝子の組合せには、Per3遺伝子とNr1d2遺伝子を用いることができる。Per3遺伝子とNr1d2遺伝子は、生体試料中に安定して発現し、発現量の日内変動（振幅）が大きい。そのため、これらの時系列発現量データから、Per3遺伝子とNr1d2遺伝子の発現量変化の概日周期を算出することで、生体リズムを正確に予測することができる。
４．概日周期の算出
時計遺伝子の発現量変化の概日周期の算出は、遺伝子発現量の経時的変化を示す時系列発現量データから、下記式（I）及び式（II）によって算出する。
ここで、式（I）中、Ea(t)、Aa、ω、Caは、時計遺伝子aの時刻tにおける発現量、振幅、初期位相、オフセット値を示す。また、式（II）中、Ea(t)、Ab、Cbは、時計遺伝子bの時刻tにおける発現量、振幅、オフセット値を示している。さらに、θは両時計遺伝子の位相差を示す。
3回のサンプリングで得られた2つの時計遺伝子の時系列発現量データから、上記式（I）及び（II）により6つの等式を得ることができる。この6つの等式から、共役勾配法等によって5つの未知数ω、Aa、Ab、C a、Cbを求める。これにより、被検個体の生体リズムを反映した、時計遺伝子の発現量変化の概日周期を算出することができる。
以上のように、本発明に係る生体リズム予測方法では、発現量変化の概日周期の位相が異なる2つの時計遺伝子について、3点サンプリングによって発現量の測定を行うことによって、発現量変化の概日周期を高精度に算出して、被検個体の生体リズムを予測することが可能である。
被検個体からの生体試料の採取を24時間以内に3回のみとすることで、生体試料の採取のための被検者の負担を軽減することができる。さらに、3点サンプリングを8時間間隔で行うことで、被検者が起きている時間帯のみに生体試料を採取することができるため、被検者の睡眠覚醒サイクルに影響を与えず、正確な生体リズムの予測を行うことが可能となる。
１．Per3遺伝子とNr1d2遺伝子の位相差の決定
本実施例では、発現量変化の概日周期の位相が異なる2つの時計遺伝子として、Per3遺伝子（以下、単に「Per3」という）とNr1d2遺伝子（以下、「Nr1d2」という）を選択し、両遺伝子の発現量変化の概日周期の位相差を決定した。
20〜50歳までの男女15人から頭髪を採取した。毛包細胞が付着した毛根部を素早く細胞溶解バッファー（RNeasy Microkit: QIAGEN）に浸すことにより、細胞溶解液を調製した。毛髪の採取は3〜4時間間隔で行い、各採取時刻において5〜20本を採取した。
細胞溶解バッファーに添付のプロトコールに従って、-70℃で保存しておいた細胞溶解液からトータルRNAを抽出し、逆転写反応を行なった。逆転写産物の１／２０量を用いてリアルタイムPCRを行い、Per3及びNr1d2の発現量を定量した。リアルタイムPCRは、SYBR Green（ABI）又はTaqMan MGB probe（ABI）を用いて、PRISM7300（ABI）により行った。Per3及びNr1d2の発現量は、内部標準とした18S-rRNAの発現量によって補正を行い、時系列発現量データとした。
得られた時系列発現量データに、24時間周期の下記コサインカーブ式（IV）を、非線形最小二乗法によりコサインフィッティングして、Per3とNr1d2の発現量変化の概日周期の位相差を求めた。
（式中、E(t)は時刻tにおける発現量、Aは発現量の振幅、ωは初期位相、Cはオフセット値を示す）。
結果を図1に示す。図は、コサインフィッティングを行った後のコサインカーブ式（IV）において、最大の発現量E(t)を与える時刻tをプロットしたものである。Per3の発現量E(t)が最大値となる時刻tをX軸に、Nr1d2の発現量E(t)が最大値となる時刻tをY軸にプロットしている。
図1に示されるように、15名の被検者について、Per3の発現量E(t)が最大値となる時刻tは、0時〜12時まで広く散らばっていた。同様に、Nr1d2についても最大の発現量E(t)を与える時刻tは、0時〜12時まで広く散らばっていた。
しかし、Per3の発現量E(t)が最大値となる時刻tと、Nr1d2の発現量E(t)が最大値となる時刻tとの時間差（位相差）は、各被検者において2時間に近い値を示した。図1中、Per3とNr1d2の位相差が2時間である場合のプロット位置を点線で示す。
15名の被検者について得られたPer3とNr1d2との位相差の平均値と標準偏差は、それぞれ２．３、０．８であった。この位相差（以下、「位相差θ」という）を用いて、Per3とNr1d2の発現量変化の概日周期を、それぞれ下記コサインカーブ式（V）及び式（VI）でモデル化した。
（式中、Eper3(t)は時刻tにおけるPer3発現量、Aper3は発現量の振幅、ω はPer3の初期位相、Cper3はオフセット値を示す）
（式中、Enr1d2(t)は時刻tにおけるNr1d2発現量、Anr1d2は発現量の振幅、θはPer3とNr1d2の発現量変化の概日周期の位相差、Cnr1d2はオフセット値を示す）
２．サンプリング時間間隔の検討
上記コサインカーブ式（V）及び式（VI）を用いれば、最少で3回のサンプリングでPer3又はNr1d2の発現量変化の概日周期を算出して、被検者の生体リズムを推定することが可能である。そこで、本実施例では、3回のサンプリング時刻の時間間隔について、最も高精度に生体リズムを推定することができる時間間隔を決定した。3回のサンプリングは24時間以内に行うものとし、想定される全てのサンプリング時間間隔（276通り）について検討を行った。
まず、実施例1において、Per3及びNr1d2の時系列発現量データに対し、コサインフィッティングを行って得られたコサインカーブ式（IV）から、両遺伝子の1時間毎の発現量変化を算出した。そして、算出された1時間毎の発現量変化を示すモデルデータから、任意の3点の時刻tとその時刻における発現量V(t)を3点サンプリングデータとして抽出した。
この任意の3点の時刻tとその時刻における発現量v(t)を、式（V）及び式（VI）に代入し、発現量の振幅Aper3・Anr1d2、Per3の初期位相ω、オフセット値Cper3・Cnr1d2を共役勾配法によって求めた。
共役勾配法は、以下の手順により行った。すなわち、まず、式（V）及び式（VI）に代入された発現量v(t)と、1時間毎の発現量変化を示すモデルデータの発現量E(t)と、の平方和dを、3点サンプリングデータとモデルデータとの距離として下記式（VII）により定義した。
そして、この距離dが最小になるような振幅Aper3・Anr1d2、初期位相ω、オフセット値Cper3・Cnr1dを共役勾配法により求めた。なお、これらの未知定数は、非線形最小二乗法によって求めようとした場合、距離dが最小値に収束せず、求めることができなかった。
ここで、共役勾配法を適用する際、未知定数の初期値は求まる定数に大きく影響する。従って、適切な初期値を設定しないと、距離dが局所的な極値に落ちてしまい、本来の概日周期を得ることができない。
そこで、ここでは、振幅Aper3・Anr1d2の適切な初期値として、1時間毎の発現量変化を示すモデルデータの振幅Aの平均値を設定した。また、オフセット値Cper3・Cnr1d2には、3点サンプリングデータの平均値を初期値として設定した。なお、モデルデータの振幅Aは、Per3の振幅Aper3が0.8014513、Nr1d2の振幅Anr1d2が0.6402411であった。
さらに、Per3の初期位相ωは、演算解析を容易にするため、0から23までの整数に限定した。そして、0から23までの整数をひとつずつ初期値として与えて共役勾配法を適用し、最小の距離dを与え得る初期位相ωを求めた。
以上のようにして、発現量の振幅Aper3、Per3の初期位相ω、オフセット値Cper3を求めた式（V）において発現量v(t)の最大値を与える時刻tと、モデルデータにおいて発現量v(t)の最大値を与える時刻との時刻差を求めた。そして、この時刻差が最小となるようなサンプリング時刻の時間間隔を検討した。
0時から23時までの1時間おきのすべての時刻を第1回目のサンプリング時刻として設定した。そして、この第1回目のサンプリングから、24時間以内に第2回目及び第3回目のサンプリングが終了するような3回のサンプリング時刻の時間間隔の組み合せ（276通り）について検討を行った。全ての時間間隔の組み合わせについて、上記時刻差の平均値と標準誤差を求めた。標準誤差は、そのサンプリング時間間隔で算出される概日周期が、どの程度第1回目のサンプリング時刻に依存するかを示す。また、平均値は、算出された概日周期の精度を示す。従って、標準誤差及び平均値が小さい程、算出された概日周期が正確であると評価できる。
全ての時間間隔の組み合わせについて得られた平均値と標準誤差を、「表1」〜「表5」に示す。時間間隔の組み合わせは、標準誤差及び平均値が小さかったものから順に「表1」から「表5」に挙げた。時間間隔は、第1回目と第2回目のサンプリング時刻との時間間隔と、第2回目と第3回目のサンプリング時刻との時間間隔と、を組み合わせて表示した。例えば、図中「13：06」は、第1回目と第2回目のサンプリング時刻との時間間隔が13時間、第2回目と第3回目のサンプリング時刻との時間間隔が6時間であることを示している。
また、図2には、各時間間隔の組み合わせを、X軸を標準誤差、Y軸を平均値としてプロットした図を示す。
「表1」〜「表5」及び図2に示されるように、「8：08」の時間間隔が、最も標準誤差及び平均値が小さかった。このことから、第1回目と第2回目のサンプリング時刻との時間間隔と、第2回目と第3回目のサンプリング時刻との時間間隔と、がともに8時間となるような3回サンプリングが、最も高精度に概日周期を算出できることが明らかとなった。
この他、「表1」〜「表4」に示す時間間隔の組み合わせでは、標準誤差及び平均値が十分に小さく、これらのサンプリン時間間隔を採用することで、正確な概日周期の算出が可能になると考えられた。一方、「表5」に示す時間間隔の組み合わせでは、標準誤差及び平均値が大きくなり、サンプリン時間間隔として不適当であることが明らかになった。
図3には、モデルデータから「8：08」の時間間隔（A）又は「9：15」の時間間隔（B）となる全ての3点サンプリングを行って算出したPer3発現量変化の概日周期を示す。図中、符号1はモデルデータのPer3、符号2はモデルデータのNr1ｄ2の概日周期を示す。また、符号3は3点サンプリングデータのPer3、符号4は3点サンプリングデータのNr1d2の概日周期を示す。
図3（A）に示す、第1回目と第2回目のサンプリングと、第2回目と第3回目のサンプリングをともに8時間間隔で行った場合では、モデルデータと3点サンプリングデータのPer3及びNr1d2の概日周期が良く一致し、精度良く概日周期が算出されている。
一方、図3（B）に示す、第1回目と第2回目のサンプリングを9時間間隔で、第2回目と第3回目のサンプリングを15時間間隔で行った場合には、モデルデータと3点サンプリングデータのPer3及びNr1d2の概日周期が一致せず、概日周期を算出できないことが分かる。
３．3点サンプリングによる生体リズムの推定
実施例2の結果から、3回のサンプリングを8時間間隔で行うことにより、最も高精度に概日周期を算出できることが明らかとなった。そこで、実際に8時間間隔の3点サンプリングによる生体リズムの予測を試みた。3点サンプリングは、以下の「表6」に示す3パターンを行った。
3点サンプリングa〜cについて、実施例1で説明した方法に従って、各時刻におけるPer3及びNr1d2の発現量を定量した。そして、実施例2で説明した方法と同様にして、3点の時刻tとその時刻における発現量E(t)を、式（V）及び式（VI）に代入し、発現量の振幅Aper3・Anr1d2、Per3の初期位相ω、オフセット値Cper3・Cnr1d2を共役勾配法によって求め、Per3及びNr1d2の発現量変化の概日周期を算出した。
結果を図4に示す。図中、符号a1, b1, c1はそれぞれ3点サンプリングパターンa, b, cのPer3遺伝子の概日周期、符号a2, b2, c2はそれぞれ3点サンプリングパターンa, b, cのNr1d2遺伝子の概日周期を示す。また、符号d1及びd2は、12：00〜36：00までの7回のサンプリングで得られた時系列発現量データに、上記コサインカーブ式（IV）を非線形最小二乗法によりコサインフィッティングして求めた概日周期である。符号d1はPer3、符号d2はNr1d2の概日周期を示す。
図4に示されるように、3点サンプリングa〜cにおいて、Per3及びNr1d2の発現量変化の概日周期が再現性良く算出できた。これらの3点サンプリングa〜cで得られたPer3及びNr1d2の発現量変化の概日周期と、7回のサンプリングで得られた概日周期（図中、符号d1及びd2参照）との位相の誤差は、平均で0.75時間であった。
この結果から、発現量変化の概日周期の位相が異なるPer3とNr1d2の発現量の3点サンプリングによって、発現量変化の概日周期を算出することによって、被検者の生体リズムを高精度に予測できることが示された。
本発明に係る生体リズム予測方法は、時間医療の実現や、自己能力の発揮、ダイエットに役立てることができる。また、生体リズムのずれを原因とする種々の疾患の予防や、時差ぼけなどの体調不良の改善に役立てることができる。
Per3遺伝子とNr1d2遺伝子の発現量変化の概日周期において、発現量の最大値を与える時刻をプロットした図である（実施例1）。X軸（横軸）はPer3遺伝子の発現量が最大値となる時刻を、Y軸（縦軸）はNr1d2遺伝子の発現量が最大値となる時刻を示す。図中、点線は、両時刻の時間差（位相差）が2時間となるプロット位置を示す。 3回のサンプリング時刻の時間間隔について、そのサンプリング時間間隔で算出される概日周期とモデルデータとの間における発現量の最大値を与える時刻の時刻差をプロットした図である（実施例2）。図中、X軸（横軸）は、そのサンプリング時間間隔における時刻差の標準誤差を、Y軸（縦軸）は、平均値を示す。 モデルデータから「8：08」の時間間隔（A）又は「9：15」の時間間隔（B）となる全ての3点サンプリングを行って算出したPer3発現量変化の概日周期を示す図である（実施例2）。図中、符号1はモデルデータのPer3、符号2はモデルデータのNr1ｄ2の概日周期を示す。また、符号3は3点サンプリングデータのPer3、符号4は3点サンプリングデータのNr1d2の概日周期を示す。 8時間間隔の3点サンプリングによって算出したPer3遺伝子及びNr1d2遺伝子の発現量変化の概日周期を示す図である（実施例3）。図中、符号a1, b1, c1はそれぞれ3点サンプリングパターンa, b, cのPer3遺伝子の概日周期、符号a2, b2, c2はそれぞれ3点サンプリングパターンa, b, cのNr1d2遺伝子の概日周期を示す。また、符号d1及びd2は、7点サンプリングで得られたPer3遺伝子及びNr1d2の概日周期を示す。

Claims (4)

  1. （１）発現量変化の概日周期の位相が異なる２つの時計遺伝子について、被検個体から２４時間以内に３回採取された生体試料中の発現量を測定する手順と、
    （２）前記（１）の手順で得られた時系列発現量データから、前記概日周期を算出する手順と、
    を含み、
    前記生体試料は、
    前記２つの時計遺伝子の時系列発現量データをコサインフィッティングして得たコサインカーブから算出された、１時間毎の発現量変化を示すモデルデータと、このモデルデータから抽出された任意の３点の時刻及び該時刻における発現量、並びに下記式（I）及び式（II）から算出された３点サンプリングデータと、の間で、発現量の最大値を与える時刻の時間差の平均値が０．６未満であり、かつ標準誤差が０．４未満である３点の時刻の時間間隔を採取時間間隔として、２４時間以内に３回採取されたものである、
    被検個体の生体リズムを予測する方法。
    （式（I）中、Ea(t)、Aa、ω、Caは、一の時計遺伝子の時刻tにおける発現量、振幅、初期位相、オフセット値を示す。また、式（II）中、Eb(t)、Ab、Cbは、他の一の時計遺伝子の時刻tにおける発現量、振幅、オフセット値を示す。さらに、θは両時計遺伝子の位相差を示す。）
  2. 前記（２）の手順において、前記時系列発現量データから上記式（I）及び式（II）によって前記概日周期を算出する請求項１記載の方法。
  3. 前記生体試料は、被検個体から２４時間以内に８時間間隔で３回採取されたものである請求項１又は２記載の方法。
  4. 前記時計遺伝子として、Per3遺伝子及びNr1d2遺伝子を用いる請求項１〜３の何れか１項記載の方法。
JP2008267464A 2008-10-16 2008-10-16 生体リズム予測方法 Expired - Fee Related JP5560550B2 (ja)

Priority Applications (7)

Application Number Priority Date Filing Date Title
JP2008267464A JP5560550B2 (ja) 2008-10-16 2008-10-16 生体リズム予測方法
EP09820656.8A EP2336355B1 (en) 2008-10-16 2009-10-13 Biological rhythm prediction method
PCT/JP2009/067966 WO2010044481A1 (ja) 2008-10-16 2009-10-13 生体リズム予測方法
US13/123,847 US20110196619A1 (en) 2008-10-16 2009-10-13 Method for predicting biological rhythm
CN2009801401016A CN102177254A (zh) 2008-10-16 2009-10-13 生物节律预测方法
BRPI0920426-1A BRPI0920426A2 (ja) 2008-10-16 2009-10-13 Biological cycle prediction method
RU2011113819/10A RU2512069C2 (ru) 2008-10-16 2009-10-13 Способ определения биологического ритма

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008267464A JP5560550B2 (ja) 2008-10-16 2008-10-16 生体リズム予測方法

Publications (2)

Publication Number Publication Date
JP2010094072A JP2010094072A (ja) 2010-04-30
JP5560550B2 true JP5560550B2 (ja) 2014-07-30

Family

ID=42106642

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008267464A Expired - Fee Related JP5560550B2 (ja) 2008-10-16 2008-10-16 生体リズム予測方法

Country Status (7)

Country Link
US (1) US20110196619A1 (ja)
EP (1) EP2336355B1 (ja)
JP (1) JP5560550B2 (ja)
CN (1) CN102177254A (ja)
BR (1) BRPI0920426A2 (ja)
RU (1) RU2512069C2 (ja)
WO (1) WO2010044481A1 (ja)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110145740A1 (en) * 2009-12-11 2011-06-16 Mccue Geoffrey Predictive Performance Calculator Apparatus and Method
JP5644231B2 (ja) * 2010-07-23 2014-12-24 ソニー株式会社 生体リズム推定装置、生体リズム推定方法およびプログラム
CN103436607B (zh) * 2013-08-02 2015-08-05 山东农业大学 鸡clock基因不同转录本的扩增方法及其引物
US20180060515A1 (en) * 2016-08-31 2018-03-01 Ronald Leslie Monarch AMS_BioRhythms computing program
CN112423648B (zh) * 2018-07-18 2024-03-22 苏州大学 一种筛选去同步化指标的方法
CN109600887A (zh) * 2019-01-28 2019-04-09 汇云时代科技(深圳)有限公司 一种日光模拟控制器
CN113506627B (zh) * 2020-11-25 2023-03-24 中国人民解放军陆军特种作战学院 一种夜间人体生理机能评估方法及系统
CN113176301B (zh) * 2021-04-30 2024-04-30 中国科学院地球化学研究所 一种基于电生理节律的植物健康活力的检测方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1542149A1 (en) * 2002-07-30 2005-06-15 Yamanouchi Pharmaceutical Co. Ltd. Apparatus for forming molecular timetable and appratus for estimating circadian clock
US20090202659A1 (en) * 2005-06-10 2009-08-13 Gimble Jeffrey M Modulation of Peripheral Clocks in Adipose Tissue
JP2007110912A (ja) * 2005-10-18 2007-05-10 Sony Corp ストレス評価方法
WO2008105476A1 (ja) * 2007-02-27 2008-09-04 Hiroshima University 概日リズム遺伝子および時計モチーフ
JP5320864B2 (ja) * 2007-07-06 2013-10-23 ソニー株式会社 生体リズム情報取得方法
JP4858344B2 (ja) 2007-07-25 2012-01-18 ソニー株式会社 毛を用いた生体リズム情報取得方法
JP4626700B2 (ja) * 2008-10-07 2011-02-09 ソニー株式会社 生体リズム情報取得方法

Also Published As

Publication number Publication date
EP2336355B1 (en) 2014-01-22
RU2512069C2 (ru) 2014-04-10
JP2010094072A (ja) 2010-04-30
RU2011113819A (ru) 2012-10-20
EP2336355A4 (en) 2012-06-06
US20110196619A1 (en) 2011-08-11
EP2336355A1 (en) 2011-06-22
BRPI0920426A2 (ja) 2018-10-23
CN102177254A (zh) 2011-09-07
WO2010044481A1 (ja) 2010-04-22

Similar Documents

Publication Publication Date Title
JP5560550B2 (ja) 生体リズム予測方法
Parkinson et al. Clinical features of Friedreich's ataxia: classical and atypical phenotypes
Gupta et al. FBN2 mutation associated with manifestations of Marfan syndrome and congenital contractural arachnodactyly
EP4371482A2 (en) Method for providing a signal quality degree associated with an analyte value measured in a continuous monitoring system
JP4858344B2 (ja) 毛を用いた生体リズム情報取得方法
Hossain et al. Can stature be estimated from tooth crown dimensions? A study in a sample of South-East Asians
Mesin et al. Non-invasive estimation of right atrial pressure using inferior vena cava echography
KR20200005779A (ko) 단일염기다형성을 이용한 탈모 표현형 예측 방법
CN106609300A (zh) 一种冠状动脉疾病风险评估试剂盒及风险评估方法
JP5320864B2 (ja) 生体リズム情報取得方法
EP2124749A1 (en) Noninvasive method to estimate variation of blood glucose levels using metabolic measurements
Cuspidi et al. Nocturnal non‐dipping pattern in untreated hypertensives at different cardiovascular risk according to the 2003 ESH/ESC guidelines
JP2013097693A (ja) 骨粗鬆症リスク判定システム及びプログラム
Goel et al. Sampling rate influences measures of blood pressure variability in ambulatory blood pressure monitoring
US11643688B2 (en) Method for diagnosing mood disorder by using circadian rhythm
CN105200131A (zh) 基于14个snp位点评估外周动脉疾病患病风险的试剂盒
Gazizova Comparison of ejection fraction and end volumes of the left ventricle as indicators of the severity of heart failure
Haerian et al. Analysis Periodontal Health Status in Postmenopausal Women with Osteoporosis Referring to Rheumatology Clinics in Yazd and Healthy People
Shirolkar et al. Envisioning forensic odontology under a new light: A pedodontist's perspective
McKenna et al. PES16 SCALING PROPERTIES OF THE DERMATOLOGY LIFE QUALITY INDEX (DLQI)
JP2021175394A (ja) 月経周期マーカー
Misicka et al. Variation in body-mass index but not waist-hip ratio influences risk for multiple sclerosis
Greenberg et al. The clinician's guide to the literature: Interpreting results
Cho et al. Validity of estimating EuroQol (EQ5D) from a Health Assessment Questionnaire (HAQ)
CN104830962A (zh) 用于活体外测定致肥胖风险的方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110926

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130115

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130314

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20131203

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140127

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20140513

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140526

LAPS Cancellation because of no payment of annual fees