JP7165617B2 - 高血圧のリスクを判定する方法 - Google Patents
高血圧のリスクを判定する方法 Download PDFInfo
- Publication number
- JP7165617B2 JP7165617B2 JP2019081513A JP2019081513A JP7165617B2 JP 7165617 B2 JP7165617 B2 JP 7165617B2 JP 2019081513 A JP2019081513 A JP 2019081513A JP 2019081513 A JP2019081513 A JP 2019081513A JP 7165617 B2 JP7165617 B2 JP 7165617B2
- Authority
- JP
- Japan
- Prior art keywords
- snp
- genotype
- disease
- risk
- model
- 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.)
- Active
Links
- 206010020772 Hypertension Diseases 0.000 title claims description 9
- 238000000034 method Methods 0.000 claims description 33
- 239000002773 nucleotide Substances 0.000 claims description 22
- 125000003729 nucleotide group Chemical group 0.000 claims description 22
- 239000003795 chemical substances by application Substances 0.000 claims description 10
- FWMNVWWHGCHHJJ-SKKKGAJSSA-N 4-amino-1-[(2r)-6-amino-2-[[(2r)-2-[[(2r)-2-[[(2r)-2-amino-3-phenylpropanoyl]amino]-3-phenylpropanoyl]amino]-4-methylpentanoyl]amino]hexanoyl]piperidine-4-carboxylic acid Chemical compound C([C@H](C(=O)N[C@H](CC(C)C)C(=O)N[C@H](CCCCN)C(=O)N1CCC(N)(CC1)C(O)=O)NC(=O)[C@H](N)CC=1C=CC=CC=1)C1=CC=CC=C1 FWMNVWWHGCHHJJ-SKKKGAJSSA-N 0.000 claims description 3
- 230000002596 correlated effect Effects 0.000 claims description 3
- 238000010801 machine learning Methods 0.000 claims description 3
- 210000001124 body fluid Anatomy 0.000 claims description 2
- 239000010839 body fluid Substances 0.000 claims description 2
- 201000010099 disease Diseases 0.000 description 89
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 89
- 230000000052 comparative effect Effects 0.000 description 21
- 210000000349 chromosome Anatomy 0.000 description 19
- 239000000523 sample Substances 0.000 description 12
- 108020004414 DNA Proteins 0.000 description 10
- 238000006243 chemical reaction Methods 0.000 description 10
- 238000007477 logistic regression Methods 0.000 description 8
- 238000003752 polymerase chain reaction Methods 0.000 description 7
- 239000013611 chromosomal DNA Substances 0.000 description 6
- UYTPUPDQBNUYGX-UHFFFAOYSA-N guanine Chemical compound O=C1NC(N)=NC2=C1N=CN2 UYTPUPDQBNUYGX-UHFFFAOYSA-N 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000000018 DNA microarray Methods 0.000 description 5
- 230000003321 amplification Effects 0.000 description 5
- 239000008280 blood Substances 0.000 description 5
- 210000004369 blood Anatomy 0.000 description 5
- 230000000875 corresponding effect Effects 0.000 description 5
- 239000012634 fragment Substances 0.000 description 5
- 238000003199 nucleic acid amplification method Methods 0.000 description 5
- 102000054765 polymorphisms of proteins Human genes 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- DDRJAANPRJIHGJ-UHFFFAOYSA-N creatinine Chemical compound CN1CC(=O)NC1=N DDRJAANPRJIHGJ-UHFFFAOYSA-N 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000009396 hybridization Methods 0.000 description 4
- 230000013011 mating Effects 0.000 description 4
- 102000053602 DNA Human genes 0.000 description 3
- 238000012300 Sequence Analysis Methods 0.000 description 3
- 108020004682 Single-Stranded DNA Proteins 0.000 description 3
- 210000001766 X chromosome Anatomy 0.000 description 3
- 239000003153 chemical reaction reagent Substances 0.000 description 3
- 230000001747 exhibiting effect Effects 0.000 description 3
- 108091008146 restriction endonucleases Proteins 0.000 description 3
- 210000003296 saliva Anatomy 0.000 description 3
- 229930024421 Adenine Natural products 0.000 description 2
- GFFGJBXGBJISGV-UHFFFAOYSA-N Adenine Chemical compound NC1=NC=NC2=C1N=CN2 GFFGJBXGBJISGV-UHFFFAOYSA-N 0.000 description 2
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 2
- 102000001554 Hemoglobins Human genes 0.000 description 2
- 108010054147 Hemoglobins Proteins 0.000 description 2
- 238000000585 Mann–Whitney U test Methods 0.000 description 2
- 208000008383 Wilms tumor Diseases 0.000 description 2
- 229960000643 adenine Drugs 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000013527 convolutional neural network Methods 0.000 description 2
- 229940109239 creatinine Drugs 0.000 description 2
- OPTASPLRGRRNAP-UHFFFAOYSA-N cytosine Chemical compound NC=1C=CNC(=O)N=1 OPTASPLRGRRNAP-UHFFFAOYSA-N 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 201000008026 nephroblastoma Diseases 0.000 description 2
- 108020004707 nucleic acids Proteins 0.000 description 2
- 102000039446 nucleic acids Human genes 0.000 description 2
- 150000007523 nucleic acids Chemical group 0.000 description 2
- 238000007894 restriction fragment length polymorphism technique Methods 0.000 description 2
- RWQNBRDOKXIBIV-UHFFFAOYSA-N thymine Chemical group CC1=CNC(=O)NC1=O RWQNBRDOKXIBIV-UHFFFAOYSA-N 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- 208000007848 Alcoholism Diseases 0.000 description 1
- 206010005746 Blood pressure fluctuation Diseases 0.000 description 1
- 208000000461 Esophageal Neoplasms Diseases 0.000 description 1
- 108010049003 Fibrinogen Proteins 0.000 description 1
- 102000008946 Fibrinogen Human genes 0.000 description 1
- 208000034826 Genetic Predisposition to Disease Diseases 0.000 description 1
- 206010071602 Genetic polymorphism Diseases 0.000 description 1
- 238000007397 LAMP assay Methods 0.000 description 1
- 206010030155 Oesophageal carcinoma Diseases 0.000 description 1
- 108020005187 Oligonucleotide Probes Proteins 0.000 description 1
- 206010033645 Pancreatitis Diseases 0.000 description 1
- 201000007930 alcohol dependence Diseases 0.000 description 1
- 102000015395 alpha 1-Antitrypsin Human genes 0.000 description 1
- 108010050122 alpha 1-Antitrypsin Proteins 0.000 description 1
- 229940024142 alpha 1-antitrypsin Drugs 0.000 description 1
- 230000004872 arterial blood pressure Effects 0.000 description 1
- 230000037182 bone density Effects 0.000 description 1
- 239000000872 buffer Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 208000020832 chronic kidney disease Diseases 0.000 description 1
- 238000003776 cleavage reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 208000029078 coronary artery disease Diseases 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 229940104302 cytosine Drugs 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 201000004101 esophageal cancer Diseases 0.000 description 1
- 229940012952 fibrinogen Drugs 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012252 genetic analysis Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000024924 glomerular filtration Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 208000022013 kidney Wilms tumor Diseases 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 210000002200 mouth mucosa Anatomy 0.000 description 1
- 239000002751 oligonucleotide probe Substances 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000007637 random forest analysis Methods 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 230000007017 scission Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 210000003765 sex chromosome Anatomy 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 229940113082 thymine Drugs 0.000 description 1
- 208000001072 type 2 diabetes mellitus Diseases 0.000 description 1
- 210000002700 urine Anatomy 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Description
本発明は高血圧のリスクを判定する方法に関する。
疾患のリスクの判定に用いるために、一塩基多型（以下、「ＳＮＰ」ともいう。）と疾患との関連性の特定が進められている。ＮＣＢＩ ＳＮＰ Ｄａｔａｂａｓｅは、ヒトのＳＮＰをまとめたデータベースであり、ＳＮＰにｒｓ番号を付して管理している。本明細書におけるｒｓ番号も、このＮＣＢＩ ＳＮＰ Ｄａｔａｂａｓｅにおける登録番号を意味するものとする。
本明細書においてｒｓ番号で特定されるＳＮＰと、当該ＳＮＰに関連する疾患、病態又は状態等として非特許文献等で開示されているものとの関係は、以下のとおりである。
rs5945326：２型糖尿病に関するＳＮＰ（非特許文献 1）
rs3905706：骨密度（背骨）に関するＳＮＰ（非特許文献 2）
rs34090688：平均赤血球ヘモグロビン量に関するＳＮＰ（非特許文献 3）
rs10277115：推算糸球体濾過量に関するＳＮＰ（非特許文献 4）
rs10277115：慢性腎臓病に関するＳＮＰ（非特許文献 5）
rs10277115：血液中のクレアチニン濃度に関するＳＮＰ（非特許文献 6）
rs2840381：膵炎に関するＳＮＰ（非特許文献 7）
rs5955543：腎芽細胞腫（ウィルムス腫瘍）に関するＳＮＰ（非特許文献 8）
rs7896783：血液中のフィブリノーゲン濃度に関するＳＮＰ（非特許文献 9）
rs671：アルコール依存症に関するＳＮＰ（非特許文献 10）
rs671：冠動脈疾患に関するＳＮＰ（非特許文献 11）
rs671：飲酒量に関するＳＮＰ（非特許文献 12）
rs671：アルコール代謝に関するＳＮＰ（非特許文献 13）
rs671：食道がんに関するＳＮＰ（非特許文献 14）
rs671：平均赤血球ヘモグロビン濃度に関するＳＮＰ（非特許文献 15）
rs671：血液中のα1-アンチトリプシン濃度に関するＳＮＰ（非特許文献 16）
rs671：血液中のクレアチニン濃度に関するＳＮＰ（非特許文献 17）
rs879922：冷温刺激による血圧変動に関するＳＮＰ（非特許文献 18）
rs5945326：２型糖尿病に関するＳＮＰ（非特許文献 1）
rs3905706：骨密度（背骨）に関するＳＮＰ（非特許文献 2）
rs34090688：平均赤血球ヘモグロビン量に関するＳＮＰ（非特許文献 3）
rs10277115：推算糸球体濾過量に関するＳＮＰ（非特許文献 4）
rs10277115：慢性腎臓病に関するＳＮＰ（非特許文献 5）
rs10277115：血液中のクレアチニン濃度に関するＳＮＰ（非特許文献 6）
rs2840381：膵炎に関するＳＮＰ（非特許文献 7）
rs5955543：腎芽細胞腫（ウィルムス腫瘍）に関するＳＮＰ（非特許文献 8）
rs7896783：血液中のフィブリノーゲン濃度に関するＳＮＰ（非特許文献 9）
rs671：アルコール依存症に関するＳＮＰ（非特許文献 10）
rs671：冠動脈疾患に関するＳＮＰ（非特許文献 11）
rs671：飲酒量に関するＳＮＰ（非特許文献 12）
rs671：アルコール代謝に関するＳＮＰ（非特許文献 13）
rs671：食道がんに関するＳＮＰ（非特許文献 14）
rs671：平均赤血球ヘモグロビン濃度に関するＳＮＰ（非特許文献 15）
rs671：血液中のα1-アンチトリプシン濃度に関するＳＮＰ（非特許文献 16）
rs671：血液中のクレアチニン濃度に関するＳＮＰ（非特許文献 17）
rs879922：冷温刺激による血圧変動に関するＳＮＰ（非特許文献 18）
Hara K, Fujita H, Johnson TA, Yamauchi T, Yasuda K, Horikoshi M, et al. Genome-wide association study identifies three novel loci for type 2 diabetes. Hum. Mol. Genet. 2014;23: 239-46.
Estrada K, Styrkarsdottir U, Evangelou E, Hsu YH, Duncan EL, Ntzani EE, et al. Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture. Nat. Genet. 2012;44: 491-501.
Kamatani Y, Matsuda K, Okada Y, Kubo M, Hosono N, Daigo Y, et al. Genome-wide association study of hematological and biochemical traits in a Japanese population. Nat. Genet. 2010;42: 210-5.
Okada Y, Sim X, Go MJ, Wu JY, Gu D, Takeuchi F, et al. Meta-analysis identifies multiple loci associated with kidney function-related traits in east Asian populations. Nat. Genet. 2012;44: 904-9.
Okada Y, Sim X, Go MJ, Wu JY, Gu D, Takeuchi F, et al. Meta-analysis identifies multiple loci associated with kidney function-related traits in east Asian populations. Nat. Genet. 2012;44: 904-9.
Okada Y, Sim X, Go MJ, Wu JY, Gu D, Takeuchi F, et al. Meta-analysis identifies multiple loci associated with kidney function-related traits in east Asian populations. Nat. Genet. 2012;44: 904-9.
Ota M, Ito T, Umemura T, Katsuyama Y, Yoshizawa K, Hamano H, et al. Polymorphism in the KCNA3 gene is associated with susceptibility to autoimmune pancreatitis in the Japanese population. Dis. Markers 2011;31: 223-9.
Turnbull C, Perdeaux ER, Pernet D, Naranjo A, Renwick A, Seal S, et al. A genome-wide association study identifies susceptibility loci for Wilms tumor. Nat. Genet. 2012;44: 681-4.
Sabater-Lleal M, Huang J, Chasman D, Naitza S, Dehghan A, Johnson AD, et al. Multiethnic meta-analysis of genome-wide association studies in >100 000 subjects identifies 23 fibrinogen-associated Loci but no strong evidence of a causal association between circulating fibrinogen and cardiovascular disease. Circulation 2013;128: 1310-24.
Quillen EE, Chen XD, Almasy L, Yang F, He H, Li X, et al. ALDH2 is associated to alcohol dependence and is the major genetic determinant of "daily maximum drinks" in a GWAS study of an isolated rural Chinese sample. Am. J. Med. Genet. B Neuropsychiatr. Genet. 2014;165B: 103-10.
Takeuchi F, Yokota M, Yamamoto K, Nakashima E, Katsuya T, Asano H, et al. Genome-wide association study of coronary artery disease in the Japanese. Eur. J. Hum. Genet. 2012;20: 333-40.
Higuchi S, Muramatsu T, Matsushita S, Shirakura K. Molecular Genetics of Alcoholism. Brain Science and Mental Disorders 1997; 8:273-282. [In Japanese].
Higuchi S, Muramatsu T, Matsushita S, Shirakura K. Molecular Genetics of Alcoholism. Brain Science and Mental Disorders 1997; 8:273-282. [In Japanese].
Cui R, Kamatani Y, Takahashi A, Usami M, Hosono N, Kawaguchi T, et al. Functional variants in ADH1B and ALDH2 coupled with alcohol and smoking synergistically enhance esophageal cancer risk. Gastroenterology 2009;137: 1768-75.
Kamatani Y, Matsuda K, Okada Y, Kubo M, Hosono N, Daigo Y, et al. Genome-wide association study of hematological and biochemical traits in a Japanese population. Nat. Genet. 2010;42: 210-5.
Setoh K, Terao C, Muro S, Kawaguchi T, Tabara Y, Takahashi M, et al. Three missense variants of metabolic syndrome-related genes are associated with alpha-1 antitrypsin levels. Nat Commun 2015;6: 7754.
Okada Y, Sim X, Go MJ, Wu JY, Gu D, Takeuchi F, et al. Meta-analysis identifies multiple loci associated with kidney function-related traits in east Asian populations. Nat. Genet. 2012;44: 904-9.
Huang J, Chen S, Lu X, Zhao Q, Rao DC, Jaquish CE, et al. Polymorphisms of ACE2 are associated with blood pressure response to cold pressor test: the GenSalt study. Am. J. Hypertens. 2012;25: 937-42.
本発明は、高血圧（以下、「本疾患」ともいう。）のリスクを判定する方法を提供することを課題とする。
本発明者らは上記課題を解決するために鋭意検討した。その結果、これまで一見すると本疾患との関連性がないと認められる個々の一塩基多型を、一つのまとまったセットとしてみたときに、本疾患との関連性があることを見出した。そして、その関連性を用いることで、本疾患のリスクを判定する本発明を完成するに至った。
すなわち、本発明の方法では、本疾患との関連性を見出した、rs5945326、rs3905706、rs34090688、rs10277115、rs2840381、rs5955543、rs7896783、rs671、及びrs879922を少なくとも含む一塩基多型セット（以下、「本ＳＮＰセット」ともいう。）の遺伝子型情報に基づいて、本疾患のリスクを判定する。
本発明の方法においては、「一塩基多型セット」とは、複数の一塩基多型の一つのまとまったセットを意味し、この一つのセットにより本疾患との関連性が見出されている。
また、本発明の方法における「遺伝子型情報」とは、一塩基多型における２つのホモ接合型（ＡＡ，ＢＢ）と、ヘテロ接合型（ＡＢ）に分類して示される、一塩基多型の遺伝子型（Ｇｅｎｏｔｙｐｅ）の情報を意味し、「本ＳＮＰセットの遺伝子型情報」とは、本ＳＮＰセットにおいて特定される各一塩基多型の遺伝子型情報を一まとまりとしたセットを意味し、言い換えれば、各ｒｓ番号で示される塩基配列中の各ＳＮＰの多型となる塩基に関する情報のセットである。本ＳＮＰセットの遺伝子型情報は、図１において示すとおりである。
本発明によれば、本疾患のリスクを判定することができる。
本発明の実施の形態について説明する。以下の実施形態は、本発明を説明するための例示であり、本発明をこの実施形態にのみ限定する趣旨ではない。本発明は、その要旨を逸脱しない限り、様々な形態で実施することができる。
本実施形態において、高血圧とは、動脈内の圧力が恒常的に高くなった状態をいう。
また、本実施形態において、本疾患は、一般には、本疾患に関する医学会の公表するガイドラインに沿って診断される疾患、医療用医薬品の添付文書において、効能・効果の欄に記載される疾患、あるいは、医薬・医療業界において汎用される用語として理解される疾患の少なくともいずれかを意味するものと解することができる。
本実施形態の方法においては、一見すると本疾患との関連性がないと認められる所定数の一塩基多型セットを用いて、本疾患のリスクを判定する。
本疾患のリスクとは、本疾患の罹りやすさや罹りにくさなどの本疾患に罹る可能性をいう。「リスクを判定する」とは、例えば、現在または将来において本疾患に罹る可能性をいくつかのレベルに分けて出力することや、数値により出力することを含む。本疾患のリスクの判定には、疾患に罹りやすい傾向にあるのか、罹りにくい傾向にあるのかといった、疾患に対する遺伝的要因あるいは遺伝的感受性についての評価が含まれる。
なお、本疾患のリスクを判定するにあたっては、本疾患のリスクの判定を受ける対象者が、本疾患のリスクの判定時において、実際に本疾患に罹患しているか（発症しているか）否かは問わない。
本実施形態の方法では、本ＳＮＰセットで特定される各ＳＮＰの遺伝子型を２つのホモ接合型（ＡＡ，ＢＢ）と、ヘテロ接合型（ＡＢ）に分類した遺伝子型のセットである、本ＳＮＰセットの遺伝子型情報を用いる。そして、本ＳＮＰセットの遺伝子型情報に基づいて、対象者の本疾患のリスクを判定する。
本実施形態の方法で用いる本ＳＮＰセットは、本疾患との関連性が従来認められていなかったＳＮＰを含むセットである。すなわち、通常は、本ＳＮＰセットに含まれるＳＮＰを個別に分析したとしても、本疾患のリスクを判定することはできない。しかしながら、本実施形態の方法では、本ＳＮＰセットに含まれる各ＳＮＰの遺伝子型情報を一まとまりのセットとして分析することにより、本疾患のリスクを判定することができる。また、本ＳＮＰセットを分析した場合と、比較ＳＮＰセットを分析した場合とを比較すると、本ＳＮＰセットを分析した場合の方が統計的に有意な結果が得られている。すなわち、本実施形態の方法において、本ＳＮＰセットを分析して本疾患のリスクを判定することで、高精度が高い、あるいは予測能力が高いリスクの判定方法を提供することができる。
以下、本ＳＮＰセットに含まれる各ＳＮＰに関連して、ｒｓ番号と、各ＳＮＰが存在する染色体番号（性染色体の場合には、ＸかＹで示す）と、各ＳＮＰの染色体上の位置と、ｒｓ番号に対応する塩基配列と、を列記する。なお、各ｒｓ番号で示される塩基配列中において、ＳＮＰは［］で囲って示す。例えば、［Ａ／Ｇ］と表記した場合には、その塩基配列の位置においてＡ又はＧの一塩基多型があることを示す。また、各ＳＮＰに関する塩基配列や疾患などの情報は、例えば、ｒｓ番号に基づいてＮＣＢＩ ＳＮＰ Ｄａｔａｂａｓｅを検索することで得られる。それらの情報は当該Ｄａｔａｂａｓｅにより参照可能であり、また、本明細書で援用する。なお、以下に記す染色体上の位置は、ａｓｓｅｍｂｌｙ ｇｅｎｏｍｅのバージョンＧＲＣｈ３７に対応するものである。
ｒｓ５９４５３２６
染色体番号 Ｘ
染色体上の位置 １５２８９９９２２
塩基配列 ＡＴＴＣＣＴＴＣＣＣＴＣＣＴＧＧＧＧＴＣＴＡＣＡＣ［Ａ／Ｇ］ＣＡＴＴＧＣＴＡＣＧＧＣＣＣＣＡＴＣＣＣＡＧＡＧＣ（配列番号１）
染色体番号 Ｘ
染色体上の位置 １５２８９９９２２
塩基配列 ＡＴＴＣＣＴＴＣＣＣＴＣＣＴＧＧＧＧＴＣＴＡＣＡＣ［Ａ／Ｇ］ＣＡＴＴＧＣＴＡＣＧＧＣＣＣＣＡＴＣＣＣＡＧＡＧＣ（配列番号１）
ｒｓ３９０５７０６
染色体番号 １０
染色体上の位置 ２８４７９９４２
塩基配列 ＴＡＴＣＡＴＣＴＡＡＣＴＴＡＴＣＡＡＡＴＴＴＧＡＧ［Ａ／Ｇ］ＧＣＡＴＡＧＡＡＴＴＧＴＴＣＡＴＣＡＣＡＡＴＴＣＴ（配列番号２）
染色体番号 １０
染色体上の位置 ２８４７９９４２
塩基配列 ＴＡＴＣＡＴＣＴＡＡＣＴＴＡＴＣＡＡＡＴＴＴＧＡＧ［Ａ／Ｇ］ＧＣＡＴＡＧＡＡＴＴＧＴＴＣＡＴＣＡＣＡＡＴＴＣＴ（配列番号２）
ｒｓ３４０９０６８８
染色体番号 １
染色体上の位置 １９８５９５５０６
塩基配列 ＧＡＧＡＴＧＡＧＡＴＡＧＡＧＡＡＴＡＧＴＴＧＧＧＧ［Ａ／Ｇ］ＴＧＣＴＣＴＴＴＡＧＡＡＡＧＧＡＴＧＡＴＣＡＡＧＧ（配列番号３）
染色体番号 １
染色体上の位置 １９８５９５５０６
塩基配列 ＧＡＧＡＴＧＡＧＡＴＡＧＡＧＡＡＴＡＧＴＴＧＧＧＧ［Ａ／Ｇ］ＴＧＣＴＣＴＴＴＡＧＡＡＡＧＧＡＴＧＡＴＣＡＡＧＧ（配列番号３）
ｒｓ１０２７７１１５
染色体番号 ７
染色体上の位置 １２８５１９５
塩基配列 ＧＡＣＴＣＣＣＣＡＧＴＧＣＣＣＣＴＣＣＴＡＡＣＡＣ［Ｔ／Ａ］ＧＴＣＴＴＴＣＴＴＣＣＡＡＧＴＣＣＴＴＣＣＣＴＧＣ（配列番号４）
染色体番号 ７
染色体上の位置 １２８５１９５
塩基配列 ＧＡＣＴＣＣＣＣＡＧＴＧＣＣＣＣＴＣＣＴＡＡＣＡＣ［Ｔ／Ａ］ＧＴＣＴＴＴＣＴＴＣＣＡＡＧＴＣＣＴＴＣＣＣＴＧＣ（配列番号４）
ｒｓ２８４０３８１
染色体番号 １
染色体上の位置 １１１２１７８２３
塩基配列 ＧＧＣＣＣＧＧＧＣＣＧＧＣＡＧＣＣＴＧＡＣＴＣＣＣ［Ｔ／Ｃ］ＧＴＧＧＧＣＣＧＧＣＴＣＴＧＣＴＧＧＧＴＧＧＧＣＡ（配列番号５）
染色体番号 １
染色体上の位置 １１１２１７８２３
塩基配列 ＧＧＣＣＣＧＧＧＣＣＧＧＣＡＧＣＣＴＧＡＣＴＣＣＣ［Ｔ／Ｃ］ＧＴＧＧＧＣＣＧＧＣＴＣＴＧＣＴＧＧＧＴＧＧＧＣＡ（配列番号５）
ｒｓ５９５５５４３
染色体番号 Ｘ
染色体上の位置 １７６９８３９７
塩基配列 ＧＧＣＴＧＡＴＴＡＴＣＣＣＣＡＴＧＧＧＡＧＧＡＡＧ［Ａ／Ｇ］ＧＧＣＴＧＣＴＧＡＧＧＧＡＡＧＴＧＣＡＴＧＧＧＣＣ（配列番号６）
染色体番号 Ｘ
染色体上の位置 １７６９８３９７
塩基配列 ＧＧＣＴＧＡＴＴＡＴＣＣＣＣＡＴＧＧＧＡＧＧＡＡＧ［Ａ／Ｇ］ＧＧＣＴＧＣＴＧＡＧＧＧＡＡＧＴＧＣＡＴＧＧＧＣＣ（配列番号６）
ｒｓ７８９６７８３
染色体番号 １０
染色体上の位置 ６５１６２１５３
塩基配列 ＡＡＡＡＡＴＴＧＴＴＡＧＴＣＣＴＴＴＧＴＧＴＴＧＡ［Ａ／Ｇ］ＴＧＡＡＴＡＡＡＴＧＡＡＡＴＴＧＴＧＡＣＴＧＴＧＴ（配列番号７）
染色体番号 １０
染色体上の位置 ６５１６２１５３
塩基配列 ＡＡＡＡＡＴＴＧＴＴＡＧＴＣＣＴＴＴＧＴＧＴＴＧＡ［Ａ／Ｇ］ＴＧＡＡＴＡＡＡＴＧＡＡＡＴＴＧＴＧＡＣＴＧＴＧＴ（配列番号７）
ｒｓ６７１
染色体番号 １２
染色体上の位置 １１２２４１７６６
塩基配列 ＣＧＡＧＴＡＣＧＧＧＣＴＧＣＡＧＧＣＡＴＡＣＡＣＴ［Ａ／Ｇ］ＡＡＧＴＧＡＡＡＡＣＴＧＴＧＡＧＴＧＴＧＧＧＡＣＣ（配列番号８）
染色体番号 １２
染色体上の位置 １１２２４１７６６
塩基配列 ＣＧＡＧＴＡＣＧＧＧＣＴＧＣＡＧＧＣＡＴＡＣＡＣＴ［Ａ／Ｇ］ＡＡＧＴＧＡＡＡＡＣＴＧＴＧＡＧＴＧＴＧＧＧＡＣＣ（配列番号８）
ｒｓ８７９９２２
染色体番号 Ｘ
染色体上の位置 １５５９０８０７
塩基配列 ＴＣＣＡＡＡＡＡＣＡＴＡＴＧＴＴＣＴＴＣＡＣＣＴＡ［Ｇ／Ｃ］ＴＡＡＣＣＣＣＡＧＴＣＣＴＴＧＡＡＴＴＴＧＣＴＧＧ（配列番号９）
染色体番号 Ｘ
染色体上の位置 １５５９０８０７
塩基配列 ＴＣＣＡＡＡＡＡＣＡＴＡＴＧＴＴＣＴＴＣＡＣＣＴＡ［Ｇ／Ｃ］ＴＡＡＣＣＣＣＡＧＴＣＣＴＴＧＡＡＴＴＴＧＣＴＧＧ（配列番号９）
本実施形態の方法において、本ＳＮＰセットを構成する各ＳＮＰはｒｓ番号により特定される塩基配列を参照することによって特定可能であるが、本明細書において記載するｒｓ番号が他のｒｓ番号と併合され、新たなｒｓ番号が付与された場合には、本明細書において該当するｒｓ番号は、併合後のｒｓ番号及び併合される他のｒｓ番号をも意味する。また、本明細書において記載するｒｓ番号が複数のｒｓ番号の併合により付与された番号である場合には、本明細書において該当するｒｓ番号は、その他の元となるｒｓ番号をも意味する。
また、ＳＮＰに関する各ｒｓ番号で示される上記塩基配列は、特定の塩基配列として示しているが、人種の相違等によって、当該塩基配列において該当するＳＮＰ以外の部分における塩基配列は変更されてもよい。
本実施形態の方法は、いずれの人種の被検者に対しても用いることができるが、特に、アジア人に好適に用いることができる。アジア人の中でも日本人等の東アジア人の被検者により好適に用いることができる。また、本実施形態の方法は、いずれの性別の被検者に対しても用いてもよい。
以下、本ＳＮＰセットの遺伝子型情報を分析することにより本疾患のリスクを判定する方法の一態様について説明する。但し、判定方法は、以下に限定されない。
はじめに、対象者の試料を用いて、試料中の本ＳＮＰセットに含まれる各ＳＮＰの遺伝子型を特定する。ＳＮＰの検出に用いる試料としては、染色体ＤＮＡを含む試料であれば特に制限されない。このような試料としては、例えば、唾液、血液、尿等の体液サンプル；口腔粘膜などの細胞サンプル；毛髪等の体毛などが挙げられる。ＳＮＰの検出には、これらの試料から常法により単離した染色体ＤＮＡを直接使用してもよいし、単離した染色体ＤＮＡを増幅して、増幅後の染色体ＤＮＡを使用してもよい。
ＳＮＰの検出は、通常の遺伝子多型解析方法によって行うことができる。例えば、ＤＮＡチップ法（ＤＮＡマイクロアレイ）、サンガー法を用いた従来型のシーケンサーや次世代シーケンサー（ＮＧＳ；Ｎｅｘｔ Ｇｅｎｅｒａｔｉｏｎ Ｓｅｑｕｅｎｃｅｒ）などを用いたシーケンス解析、ＰＣＲ（Ｐｏｌｙｍｅｒａｓｅ Ｃｈａｉｎ Ｒｅａｃｔｉｏｎ）、ハイブリダイゼーション、インベーダー法などが挙げられるが、これらに限定されない。
ＤＮＡチップ法では、ＳＮＰ部位を含む多数のＤＮＡ断片（プローブ）を基板上に配置したＤＮＡチップを用い、染色体ＤＮＡをチップ上のプローブとハイブリダイズさせて、結合部位を蛍光又は電流により検出することにより、染色体ＤＮＡの配列を解析する。ＳＮＰ解析に用いられるＤＮＡチップとしては、ＳＮＰ部位を含む塩基配列を検出可能なオリゴヌクレオチドプローブが配置されたチップが挙げられる。
また、シーケンス解析は通常のサンガー法により行うことができる。例えば、多型を示す塩基の５'側の数十塩基の位置に設定したプライマーを使用してシーケンス反応を行い、その解析結果から、該当する位置がどの種類の塩基であるかを決定することができる。なお、シーケンス反応の前に、あらかじめＳＮＰ部位を含む断片をＰＣＲなどによって増幅しておくことが好ましい。効率の観点からはＮＧＳ技術を使用してもよい。
また、ＳＮＰの検出は、例えば従来のＰＣＲによる増幅の有無を調べることによって行うことができる。例えば、多型を示す塩基を含む領域に対応する配列を有し、かつ、３'末端が各多型に対応するプライマーをそれぞれ用意する。それぞれのプライマーを使用してＰＣＲを行い、増幅産物の有無によってどのタイプの多型であるかを決定することができる。また、ＬＡＭＰ法（Ｌｏｏｐ－Ｍｅｄｉａｔｅｄ Ｉｓｏｔｈｅｒｍａｌ Ａｍｐｌｉｆｉｃａｔｉｏｎ；特許第３３１３３５８号明細書）、ＮＡＳＢＡ法（Ｎｕｃｌｅｉｃ Ａｃｉｄ Ｓｅｑｕｅｎｃｅ－Ｂａｓｅｄ Ａｍｐｌｉｆｉｃａｔｉｏｎ；特許２８４３５８６号明細書）、ＩＣＡＮ法（Ｉｓｏｔｈｅｒｍａｌ ａｎｄ Ｃｈｉｍｅｒｉｃ ｐｒｉｍｅｒ－ｉｎｉｔｉａｔｅｄ Ａｍｐｌｉｆｉｃａｔｉｏｎ ｏｆ Ｎｕｃｌｅｉｃ ａｃｉｄｓ；特許第３４３３９２９号公報）などによって増幅の有無を調べることもできる。その他、単鎖増幅法やＮＧＳを用いた解析法を用いてもよい。
また、ＳＮＰ部位を含むＤＮＡ断片を増幅し、増幅産物の電気泳動における移動度の違いによってどのタイプの多型であるかを決定することもできる。このような方法としては、例えば、ＰＣＲ－ＳＳＣＰ（ｓｉｎｇｌｅ－ｓｔｒａｎｄ ｃｏｎｆｏｒｍａｔｉｏｎ ｐｏｌｙｍｏｒｐｈｉｓｍ）法（Ｇｅｎｏｍｉｃｓ． １９９２ Ｊａｎ １；１２（１）：１３９－１４６．）などが挙げられる。具体的には、まず、目的のＳＮＰを含むＤＮＡを増幅し、増幅したＤＮＡを一本鎖ＤＮＡに解離させる。次いで、解離させた一本鎖ＤＮＡを非変性ゲル上で分離し、分離した一本鎖ＤＮＡのゲル上での移動度の違いによってどのタイプの多型であるかを決定することができる。
さらに、多型を示す塩基が制限酵素認識配列に含まれる場合は、制限酵素による切断の有無によって解析することもできる（ＲＦＬＰ（Ｒｅｓｔｒｉｃｔｉｏｎ Ｆｒａｇｍｅｎｔ Ｌｅｎｇｔｈ Ｐｏｌｙｍｏｒｐｈｉｓｍ）法）。この場合、まず、ＤＮＡ試料を制限酵素により切断する。次いで、ＤＮＡ断片を分離し、検出されたＤＮＡ断片の大きさによってどのタイプの多型であるかを決定することができる。
また、ハイブリダイゼーションの有無を調べることによって多型の種類を解析することも可能である。すなわち、各塩基に対応するプローブを用意し、いずれのプローブにハイブリダイズするかを調べることによってＳＮＰがいずれの塩基であるかを調べることもできる。
このようにして、本ＳＮＰセットの各ＳＮＰに関して、対象者の遺伝子型データを決定することができる。なおここで、「対象者の遺伝子型データ」とは、対象者の有する遺伝子型情報をいう。
次いで、本ＳＮＰセットの遺伝子型情報に基づいて、本疾患のリスクを判定する。リスクの判定には、任意のモデルを用いることができる。モデルとしては、特に制限されないが、例えば、本ＳＮＰセットの遺伝子型情報を用いて、対象者の遺伝子型データから算出される特徴量を入力とし、本疾患のリスクを出力とするロジスティック回帰モデルを用いることができる。当該ロジスティック回帰モデルは、予め、本疾患に罹患したヒトの遺伝子型データと、本疾患に罹患していないヒトの遺伝子型データと、を学習データとして用いてパラメータを機械学習している。
また、疾患のリスクを判定するモデルとしては、ロジスティック回帰モデルに代えて、多層パーセプトロン、ＣＮＮ（Ｃｏｎｖｏｌｕｔｉｏｎａｌ Ｎｅｕｒａｌ Ｎｅｔｗｏｒｋ）及びＲＮＮ（Ｒｅｃｕｒｒｅｎｔ Ｎｅｕｒａｌ Ｎｅｔｗｏｒｋ）などのニューラルネットワーク、ガウシアンカーネル等の任意のカーネル関数を用いるサポートベクターマシーン、回帰木としてモデル化したランダムフォレスト、重回帰分析、隠れマルコフモデルなどを利用したモデル、統計モデルや確率モデルなど種々の他のモデルを採用することもできる。また、種々のモデルを組み合わせて総合的な判定を行うモデルを採用することもできる。
次いで、モデルを用いた本疾患のリスク判定の一例について説明する。まず、本疾患のリスクの判定をする対象者の遺伝子型データを、モデルに入力可能な特徴量に変換する。本実施形態の方法における特徴量は、例えば、本ＳＮＰセットの各ＳＮＰについて、対象者の遺伝子型データがホモ接合型（ＡＡ）、ホモ接合型（ＢＢ）、又はヘテロ接合型（ＡＢ）のいずれであるかを示すパラメータである。遺伝子型は、相同染色体のＳＮＰが共にＧ（グアニン）であることを示す“ＧＧ”や、一方がＧ（グアニン）で、他方がＡ（アデニン）であることを示す“ＡＧ”などヌクレオチドにより表記されることが一般的であるため、対象者の遺伝子型データを、本ＳＮＰセットの遺伝子型情報を用いるモデルに入力可能なパラメータに変換する。しかし、モデルが、このようなパラメータへの変換の必要がないものである場合には、上記変換は必要とされない。
対象者の遺伝子型データの特徴量への変換は、例えば、本ＳＮＰセットに含まれるＳＮＰ１つ１つに関して、対象者の遺伝子型データに値を付すことにより行うことができる。例えば、各ＳＮＰについて、対象者の遺伝子型データがホモ接合型（ＡＡ）、ホモ接合型（ＢＢ）、又はヘテロ接合型（ＡＢ）のいずれに該当するのかに応じて、そのＳＮＰに値（例えば、０又は１）を対応づける。これにより、対象者の遺伝子型データを特徴量に変換することができる。なお、以下では、各ＳＮＰに対応させる値を０又は１とした場合を例に説明するが、ＳＮＰに対応させる値は０又は１の２つの値に限られるものではない。
接合型に対応づける値はＳＮＰごとに決めることができる。例えば、あるＳＮＰは、対象者の遺伝子型データがホモ接合型（ＡＡ）である場合に値１を対応付け、ホモ接合型（ＢＢ）及びヘテロ接合型（ＡＢ）である場合に値０を対応付けるようにし、他のＳＮＰは、対象者の遺伝子型データがヘテロ接合型（ＡＢ）である場合に値１を対応付け、ホモ接合型（ＡＡ）及びホモ接合型（ＢＢ）である場合に値０を対応付けるようにしてもよい。そのほか、対象者の遺伝子型データがヘテロ接合型（ＡＢ）及びホモ接合型（ＢＢ）である場合に値１を対応付け、ホモ接合型（ＡＡ）である場合に値０を対応付けるようにしてもよい。
上記のように、対象者の遺伝子型データを特徴量に変換することができる。この特徴量への変換において対応付けに使用する値は、任意に決定することができる。例えば、上記非特許文献に基づいて、もともと各ＳＮＰが関連する疾患に関係の高い遺伝子型に対して値１を対応付けるようにし、かつ、各ＳＮＰが関連する疾患に関係の低い遺伝子型に対して値０を対応付けるようにすることができる。
このような、ＳＮＰごとの接合型とその接合型に対応付ける値の関係は、図１のような本ＳＮＰセットの遺伝子型情報をもとに、例えば、図２のような変換テーブルとして表すこともできる。図２の変換テーブルでは、網掛けをした遺伝子型と一致する場合にはそのＳＮＰに対応付ける値を１とし、一致しない場合には対応付ける値を０とする。なお、図１及び２の具体的な遺伝子型の表記において、Ａはアデニン、Ｇはグアニン、Ｃはシトシン、Ｔはチミンを示す。但し、特徴量の変換テーブルの形式は、図２に限定されるものではない。
最後に、本ＳＮＰセットの遺伝子型情報に基づいて、対象者の本疾患のリスクを判定する。より具体的には、本ＳＮＰセットの遺伝子型情報に基づく変換テーブルを用いて、対象者の遺伝子型データを、モデルに入力可能に変換した特徴量として算出し、当該特徴量を所定の判定モデルに入力し、対象者の本疾患のリスクを判定することができる。
判定モデルにおいて、特徴量には、本ＳＮＰセットのＳＮＰごとに、本疾患のリスクと正の相関があることを表す重みづけや、本疾患のリスクと負の相関があることを表す重みづけをすることができる。例えば、rs3905706に対応付けた値（特徴量）ついては、本疾患のリスクと正の相関があることを表す重みづけを行い、rs2840381、rs7896783、rs879922、rs5945326、rs671、rs5955543、rs10277115、及びrs34090688に対応付けた値（特徴量）ついては、本疾患のリスクと負の相関があることを表す重みづけを行うことができる。
例えば、特徴量に対して重みづけを行う場合には、rs3905706の遺伝子型がAGである場合に、本疾患のリスクと正の相関があることを表す重みづけを行い、rs2840381の遺伝子型がTC又はTT、rs7896783の遺伝子型がGG、rs879922の遺伝子型がCC、rs5945326の遺伝子型がAG、rs671の遺伝子型がAG、rs5955543の遺伝子型がAG、rs10277115の遺伝子型がAA、及びrs34090688の遺伝子型がGGである場合に、本疾患のリスクと負の相関があることを表す重みづけを行うことができる。また、特徴量として値０を対応付けた各ＳＮＰの遺伝子型の場合には、本疾患のリスクと相関がない或いは無視しうる程度に低いものと評価することができる。
このような、本疾患のリスクとの相関を表す重みづけは、本疾患に罹患したヒトの遺伝子型データと、本疾患に罹患していないヒトの遺伝子型データと、を学習データとして用いてパラメータを機械学習することにより特定される。この際、あるモデルにおいて、あるＳＮＰが本疾患のリスクと正の相関があることを表す重み付けがなされるとした場合、他のモデルにおいても同様にそのＳＮＰは本疾患のリスクと正の相関があることを表す重み付けがなされることが通常である。すなわち、モデルの種類等によって、あるＳＮＰにおいて本疾患のリスクとの相関関係が逆になるような事態は想定し難い。なお、重みづけの具体的な値はモデルによって異なり、特に制限されるものではない。
ここで、本ＳＮＰセットのなかで本疾患のリスクと正の相関があることを表す重みづけを行うＳＮＰのまとまりを「正相関ＳＮＰセット」といい、本疾患のリスクと負の相関があることを表す重みづけを行うＳＮＰのまとまりを「負相関ＳＮＰセット」という。本ＳＮＰセットは、正相関ＳＮＰセットと負相関ＳＮＰセットとを含むものであり、このような本ＳＮＰセットの遺伝子型情報に基づくことにより、対象者の本疾患のリスクを、リスクが上昇する要因とリスクが低下する要因の両面を総合して判定することができる。
上記のようにして得られる判定結果は、本疾患の専門医が本疾患を診断する際の補助としても用いられる。また、上記のようにして判定した本疾患のリスクと、対象者からのアンケート結果とに基づいて、本疾患のリスクの判定結果は補正されしてもよい。また、本疾患のリスクと、対象者からのアンケート結果とに基づいて、対象者に対して、生活改善に関するアドバイスを出力してもよい。
本発明は、プライマーやプローブなどの検査試薬を提供することもできる。このようなプローブとしては、上記ＳＮＰ部位を含み、ハイブリダイズの有無によってＳＮＰ部位の塩基の種類を判定できるプローブが挙げられる。また、プライマーとしては、上記ＳＮＰ部位を増幅するためのＰＣＲに用いることのできるプライマー、又は上記ＳＮＰ部位をシーケンス解析するために用いることのできるプライマーが挙げられる。本実施形態の検査試薬はこれらのプライマーやプローブに加えて、ＰＣＲ用のポリメラーゼやバッファー、ハイブリダイゼーション用試薬などを含むものであってもよい。
以下、本実施形態を実施例によりさらに具体的に説明する。但し、本実施形態はこれらの実施例に限定されない。
本ＳＮＰセットと本疾患との関連性を、以下のように検証した。
遺伝子解析サービスの利用者７万３千人以上から、利用者の同意のもと、唾液試料と、各種疾患の罹患情報を収集した。罹患情報とは、例えば、本疾患に罹患している場合に１、罹患していない場合に０となる数値である。唾液試料から、利用者ごとの遺伝子型データを特定し、利用者の遺伝子型データと各種罹患情報とを対応付けたデータベースを構築した。このデータベースの中から、本疾患に罹患している被検者１５００名と、罹患していないコントロール１５００名との症例対照セットを構築した。
次いで、被検者及びコントロールの本ＳＮＰセットの各ＳＮＰの遺伝子型を、２つのホモ接合型（ＡＡ，ＢＢ）と、ヘテロ接合型（ＡＢ）に分類した。そして、遺伝子型が図２に示す網掛けをした変換テーブルの遺伝子型と一致する場合には、ｘiの値を１とし、一致しない場合には０として、ｘ1～ｘNを以下の数式（１）で表されるロジスティック回帰モデルの説明変数とした。例えば、ｒｓ５９４５３２６の場合は、遺伝子型が“ＡＧ”であるときにｘ1の値を１とし、遺伝子型が“ＡＡ”又は“ＧＧ”であるときにはｘ1の値を０とした。なお、本実施例ではＮ＝１０である。また、以下の数式で表されるロジスティック回帰モデルの目的変数は、本疾患に罹患している確率を表す０から１の間の値ｐ（罹患情報）とした。
α＝０．１
１．ＡＵＣによるモデルの検証
本ＳＮＰセットを用いた判定方法の精度について説明する。上記データベースから、テスト用に、利用者の遺伝子型情報と罹患情報とを対応付けたデータセットを作成した。データセットにおける各利用者の本ＳＮＰセットの各ＳＮＰの遺伝子型を、ホモ接合型（ＡＡ，ＢＢ）と、ヘテロ接合型（ＡＢ）に分類し、分類した各遺伝子型が図２に示す網掛けをした遺伝子型と一致する場合には、ｘiの値を１と評価し、一致しない場合には０と評価して、ｘ1～ｘNを特徴量として算出した。
本ＳＮＰセットを用いた判定方法の精度について説明する。上記データベースから、テスト用に、利用者の遺伝子型情報と罹患情報とを対応付けたデータセットを作成した。データセットにおける各利用者の本ＳＮＰセットの各ＳＮＰの遺伝子型を、ホモ接合型（ＡＡ，ＢＢ）と、ヘテロ接合型（ＡＢ）に分類し、分類した各遺伝子型が図２に示す網掛けをした遺伝子型と一致する場合には、ｘiの値を１と評価し、一致しない場合には０と評価して、ｘ1～ｘNを特徴量として算出した。
利用者毎の本ＳＮＰセットに関する特徴量を上記ロジスティック回帰モデル（以下、「判定モデル」ともいう。）に入力し、各利用者が本疾患に罹患しているか否かを予測し、その偽陽性率と真陽性率を算出し、ＲＯＣ（Ｒｅｃｅｉｖｅｒ Ｏｐｅｒａｔｉｎｇ Ｃｈａｒａｃｔｅｒｉｓｔｉｃ）曲線とＡＵＣ（Ａｒｅａ Ｕｎｄｅｒ ｔｈｅ Ｃｕｒｖｅ）をそれぞれ求めた。より具体的には、判定モデルについて５分割交差検証を行い、５つのＲＯＣ曲線（ＲＯＣ ｆｏｌｄ １からＲＯＣ ｆｏｌｄ ５）を求めて、その平均（Ｍｅａｎ ＲＯＣ）と標準偏差（±１ｓｔｄ． ｄｅｖ．）を求めた。図３中の破線（Ｌｕｃｋ）は、本疾患に罹患しているか否かをランダムに出力する場合であり、予測能力が無いモデルのＲＯＣ曲線に対応している。
また、同様にして、本ＳＮＰセットから１つのＳＮＰを除いた各比較ＳＮＰセットに対して、上記と同様にそれぞれロジスティック回帰モデル（以下、「比較判定モデル」ともいう。）を作成した。そして、各比較ＳＮＰ関する特徴量を各比較判定モデルに入力し、各利用者が本疾患に罹患しているか否かを予測し、偽陽性率と真陽性率を算出し、ＲＯＣ曲線とＡＵＣをそれぞれ求めた。その結果を図４以降に示す。
本ＳＮＰセットを用いて本疾患を判定した場合、ＡＵＣは０．７２±０．０１であり、ランダムな出力の場合（ＡＵＣ＝０．５）と比べて有意に高く、本ＳＮＰセットを用いる判定モデルの予測能力が高いことが確認できる。
一方、各比較ＳＮＰセットを用いる比較判定モデルの場合、ＡＵＣは本ＳＮＰセットを用いる場合よりも低い。従って、各比較ＳＮＰセットを用いる比較判定モデルのＡＵＣは、ランダムな出力の場合（ＡＵＣ＝０．５）よりも高いものの、本ＳＮＰセットを用いる判定モデルのＡＵＣ（０．７２±０．０１）よりも総じて低いことが確認できる。
よって、本ＳＮＰセットに含まれるＳＮＰ全てを用いて判定することで、本ＳＮＰセットから１つのＳＮＰを除いた各比較ＳＮＰセットを用いる場合よりも、本疾患に罹患しているか否かを高精度で予測できることが分かった。
２．ウィルコクソンの順位和検定による検証
本ＳＮＰセットを用いる判定モデルが、各比較ＳＮＰセットを用いる比較判定モデルよりも有意に優れたモデルであることを確かめるために、ノンパラメトリック検定の一種であるウィルコクソンの順位和検定を行った。具体的には、本ＳＮＰセットを用いる判定モデルのＡＵＣと、各比較ＳＮＰセットを用いる比較判定モデルのＡＵＣとに差が無いという帰無仮説を設定し、有意水準を０．０１としてウィルコクソンの順位和検定を行った。
本ＳＮＰセットを用いる判定モデルが、各比較ＳＮＰセットを用いる比較判定モデルよりも有意に優れたモデルであることを確かめるために、ノンパラメトリック検定の一種であるウィルコクソンの順位和検定を行った。具体的には、本ＳＮＰセットを用いる判定モデルのＡＵＣと、各比較ＳＮＰセットを用いる比較判定モデルのＡＵＣとに差が無いという帰無仮説を設定し、有意水準を０．０１としてウィルコクソンの順位和検定を行った。
その結果、ｐ値はいずれも３．９６×１０-18であり、帰無仮説が棄却されることが確認された。すなわち、本ＳＮＰセットを用いる判定モデルのＡＵＣと、各比較ＳＮＰセットを用いる比較判定モデルのＡＵＣとは、統計的に有意な差があり、本ＳＮＰセットを用いる判定モデルは、各比較ＳＮＰセットを用いる比較判定モデルよりも優れたモデルであるといえる。
上記のとおり、本実施形態の方法は、本疾患に罹患しているか否かを予測する精度が、ランダムな予測の場合の精度よりも有意に高いという効果を有する。また、本実施形態の方法は、本ＳＮＰセットの遺伝子型情報に基づく本疾患の判定の結果と、比較ＳＮＰセットの遺伝子型情報に基づく本疾患の判定の結果との間に、有意な差があるという効果を有する。当該効果は、本ＳＮＰセットの遺伝子型情報と本疾患の間に、これまで見出されていなかった潜在的な相関性があることに基づくものであると考えられる。上記で例示したロジスティック回帰モデルやその他のモデルは、本ＳＮＰセットの遺伝子型情報を前提として、本疾患に罹患したヒトと本疾患に罹患していないヒトの遺伝子型に関するデータと罹患情報を学習データとして用いてパラメータを機械学習することなどにより得られるものである。すなわち、いずれのモデルも上記潜在的な相関性を表す一つの表現型に過ぎず、本実施形態の方法の実施において使用されるモデルの種類は特に限定されるものではない。
本発明の方法は、医療やヘルスケアに関連する分野において、本疾患のリスクを判定し、その予防および／または治療に貢献するものである。
Claims (2)
- 高血圧に罹患したヒトの遺伝子型データと、高血圧に罹患していないヒトの遺伝子型データと、を学習データとして用いて機械学習したモデルを用いて、高血圧と正の相関があるrs3905706（配列番号２）と、高血圧と負の相関があるrs2840381（配列番号５）、rs7896783、rs879922、rs5945326、rs671、rs5955543、rs10277115を少なくとも含む一塩基多型セットの遺伝子型情報であって、
rs3905706（配列番号２）の遺伝子型がAG、rs2840381（配列番号５）の遺伝子型がTC又はTT、rs7896783の遺伝子型がGG、rs879922の遺伝子型がCC、rs5945326の遺伝子型がAG、rs671の遺伝子型がAG、rs5955543の遺伝子型がAG、及びrs10277115の遺伝子型がAAであるか否かに関する遺伝子型情報に基づいて、高血圧のリスクを判定する、方法。 - リスクの判定を受ける対象者の体液サンプル、細胞サンプル又は体毛を用いる、請求項１に記載の方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019081513A JP7165617B2 (ja) | 2019-04-23 | 2019-04-23 | 高血圧のリスクを判定する方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019081513A JP7165617B2 (ja) | 2019-04-23 | 2019-04-23 | 高血圧のリスクを判定する方法 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2020174639A JP2020174639A (ja) | 2020-10-29 |
| JP2020174639A5 JP2020174639A5 (ja) | 2022-03-15 |
| JP7165617B2 true JP7165617B2 (ja) | 2022-11-04 |
Family
ID=72935557
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2019081513A Active JP7165617B2 (ja) | 2019-04-23 | 2019-04-23 | 高血圧のリスクを判定する方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7165617B2 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118813782A (zh) * | 2024-08-01 | 2024-10-22 | 珠海索因医学科技有限公司 | 一种高血压预警模型、产品、计算机可读存储介质及其应用 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105002286A (zh) | 2015-07-30 | 2015-10-28 | 中国医学科学院阜外心血管病医院 | 与高血压和/或心血管疾病发病风险相关的多个单核苷酸多态性位点及相关应用 |
| JP2018143178A (ja) | 2017-03-06 | 2018-09-20 | 国立大学法人三重大学 | 循環器疾患の遺伝的リスク検出法 |
-
2019
- 2019-04-23 JP JP2019081513A patent/JP7165617B2/ja active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105002286A (zh) | 2015-07-30 | 2015-10-28 | 中国医学科学院阜外心血管病医院 | 与高血压和/或心血管疾病发病风险相关的多个单核苷酸多态性位点及相关应用 |
| JP2018143178A (ja) | 2017-03-06 | 2018-09-20 | 国立大学法人三重大学 | 循環器疾患の遺伝的リスク検出法 |
Non-Patent Citations (2)
| Title |
|---|
| dbSNP[online], RefSNP No.rs34090688 (ss41231395)，2005年7月16日（検索日：2022年5月26日）、<https://www.ncbi.nlm.nih.gov/projects/SNP/snp_ss.cgi?subsnp_id=ss41231395> |
| Hypertension Research，2019年，Vol.42，pp.681-689，Published online: 12 December 2018 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2020174639A (ja) | 2020-10-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP7165617B2 (ja) | 高血圧のリスクを判定する方法 | |
| JP2020174538A (ja) | ２型糖尿病のリスクを判定する方法 | |
| JP7165098B2 (ja) | 動脈硬化のリスクを判定する方法 | |
| JP7137525B2 (ja) | 接触性皮膚炎のリスクを判定する方法 | |
| JP7161441B2 (ja) | 慢性副鼻腔炎のリスクを判定する方法 | |
| JP7161442B2 (ja) | リウマチのリスクを判定する方法 | |
| JP7165099B2 (ja) | 心筋梗塞及び／又は狭心症のリスクを判定する方法 | |
| JP7140707B2 (ja) | 緑内障のリスクを判定する方法 | |
| JP7097853B2 (ja) | 妊娠高血圧症候群のリスクを判定する方法 | |
| JP7107882B2 (ja) | 偏頭痛のリスクを判定する方法 | |
| JP7107884B2 (ja) | 食物アレルギーのリスクを判定する方法 | |
| JP7108574B2 (ja) | アレルギー性結膜炎のリスクを判定する方法 | |
| JP7099981B2 (ja) | 痛風のリスクを判定する方法 | |
| JP7165619B2 (ja) | 慢性気管支炎のリスクを判定する方法 | |
| JP7160751B2 (ja) | 膠原病のリスクを判定する方法 | |
| JP7106490B2 (ja) | 胆石のリスクを判定する方法 | |
| JP7161443B2 (ja) | 尿管結石及び／又は腎結石のリスクを判定する方法 | |
| JP7137524B2 (ja) | 変形性膝関節症のリスクを判定する方法 | |
| JP7108573B2 (ja) | 慢性中耳炎のリスクを判定する方法 | |
| JP7106489B2 (ja) | 脂肪肝のリスクを判定する方法 | |
| JP7097849B2 (ja) | 過敏性腸症候群のリスクを判定する方法 | |
| JP7137517B2 (ja) | 鉄欠乏性貧血のリスクを判定する方法 | |
| JP7138075B2 (ja) | 甲状腺機能低下症のリスクを判定する方法 | |
| JP7138073B2 (ja) | 注意欠陥多動性症候群のリスクを判定する方法 | |
| JP7137523B2 (ja) | じんましんのリスクを判定する方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200609 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20211227 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220304 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220602 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220729 |
|
| 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: 20221003 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20221024 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7165617 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |