WO2020040219A1 - 疾病評価指標算出方法、装置、システム、及び、プログラム、並びに、疾病評価指標を算出するためのモデル作成方法 - Google Patents

疾病評価指標算出方法、装置、システム、及び、プログラム、並びに、疾病評価指標を算出するためのモデル作成方法 Download PDF

Info

Publication number
WO2020040219A1
WO2020040219A1 PCT/JP2019/032698 JP2019032698W WO2020040219A1 WO 2020040219 A1 WO2020040219 A1 WO 2020040219A1 JP 2019032698 W JP2019032698 W JP 2019032698W WO 2020040219 A1 WO2020040219 A1 WO 2020040219A1
Authority
WO
WIPO (PCT)
Prior art keywords
data
disease
user
attributes
intestinal flora
Prior art date
Application number
PCT/JP2019/032698
Other languages
English (en)
French (fr)
Inventor
増山博昭
和巳 蓮子
正三 本川
Original Assignee
一般社団法人 日本農業フロンティア開発機構
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
Application filed by 一般社団法人 日本農業フロンティア開発機構 filed Critical 一般社団法人 日本農業フロンティア開発機構
Priority to CN201980051821.9A priority Critical patent/CN112534513A/zh
Priority to KR1020217001054A priority patent/KR20210049081A/ko
Priority to EP19852082.7A priority patent/EP3832664A4/en
Publication of WO2020040219A1 publication Critical patent/WO2020040219A1/ja
Priority to US17/181,152 priority patent/US20210174967A1/en

Links

Images

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/02—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving viable microorganisms
    • C12Q1/04—Determining presence or kind of microorganism; Use of selective media for testing antibiotics or bacteriocides; Compositions containing a chemical indicator therefor
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • 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
    • C12Q1/6888—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
    • C12Q1/689—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms for bacteria
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/483—Physical analysis of biological material
    • G01N33/4833—Physical analysis of biological material of solid biological material, e.g. tissue samples, cell cultures
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/20—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for electronic clinical trials or questionnaires
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/40—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for data related to laboratory analysis, e.g. patient specimen analysis
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/60—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H15/00—ICT specially adapted for medical reports, e.g. generation or transmission thereof
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/70—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H70/00—ICT specially adapted for the handling or processing of medical references
    • G16H70/60—ICT specially adapted for the handling or processing of medical references relating to pathologies

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Chemical & Material Sciences (AREA)
  • Biomedical Technology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Pathology (AREA)
  • Databases & Information Systems (AREA)
  • Organic Chemistry (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Immunology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Biochemistry (AREA)
  • Microbiology (AREA)
  • Biotechnology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • Genetics & Genomics (AREA)
  • Food Science & Technology (AREA)
  • Optics & Photonics (AREA)
  • Hematology (AREA)
  • Urology & Nephrology (AREA)
  • Medicinal Chemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Toxicology (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Medical Treatment And Welfare Office Work (AREA)

Abstract

腸内細菌叢を用いた疾病評価指標を算出することに関する技術を提供することにある。　ユーザが調べて欲しいと望む１又は複数の疾病を入力し、ユーザの便を予め解析した結果であるユーザの腸内細菌叢、ユーザの属性、並びに、複数の被験者の属性、被験者の便を解析した結果である被験者の腸内細菌叢、及び、疾病を内部因子とする連関モデルをデータベースから抽出し、ユーザの腸内細菌叢、ユーザの属性、及び、連関モデルを用いて、前記疾病に対するリスクを算出する。

Description

疾病評価指標算出方法、装置、システム、及び、プログラム、並びに、疾病評価指標を算出するためのモデル作成方法
　本発明の実施形態は、腸内細菌叢を用いた疾病評価指標を算出する技術に関し、特に、採便を提出、かつ、属性などのアンケートを回答した被験者に対し、特定の疾病に対するリスクを示す評価レポートを提供する技術に関する。
　腸内細菌叢は数百～千種を超えるともいわれる膨大な数の細菌から構成される極めて複雑な生態系であり、腸内細菌叢を培養、観察、同定、解析する作業には莫大な労力と時間が必要である。近年、分子生物学的な手法により、これらの問題を解決することが可能になった。この手法には、細菌の遺伝子、特に１６ＳｒＲＮＡの遺伝子（１６ＳｒＤＮＡ）の配列の違いを検出する技術が用いられている。
　このような技術を用いることによって、腸内細菌叢と特定の疾病との関係についての研究が進んでいる。非特許文献１では、多数の被験者に採便キットで採取した大便を提出してもらい、さらに、年齢やＢＭＩ（Ｂｏｄｙ　Ｍａｓｓ　Ｉｎｄｅｘ）の情報も提供してもらい、腸内細菌叢の情報を用いて疾病の罹患を評価するために、多重ロジスティック回帰分析によって特定の菌群と年齢やＢＭＩとの関係を評価する研究が行われている。
　また、非特許文献２では、腸内細菌叢の情報を用いた構造方程式のモデリングが記載されている。この構造方程式モデリングは疾病ではなく、菌叢の解析手法（Ｎｅｘｔ　Ｇｅｎｅｒａｔｉｏｎ　Ｓｅｑｕｅｎｃｉｎｇなど）を評価することを目的に研究されている。
Ｌａｈｔｉ　ｅｔ　ａｌ．、Ｔｉｐｐｉｎｇ　ｅｌｅｍｅｎｔｓ　ｉｎ　ｔｈｅ　ｈｕｍａｎ　ｉｎｔｅｓｔｉｎａｌ　ｅｃｏｓｙｓｔｅｍ、ｎａｔｕｒｅ　ＮＡＴＵＲＥ　ＣＯＭＭＵＮＩＣＡＴＩＯＮＳ　５、Ａｒｔｉｃｌｅ　ｎｕｍｂｅｒ：４３４４（２０１４）、Ｐｕｂｌｉｓｈｅｄ：　８　Ｊｕｌｙ　２０１４、ｈｔｔｐｓ：／／ｗｗｗ．ｎａｔｕｒｅ．ｃｏｍ／ａｒｔｉｃｌｅｓ／ｎｃｏｍｍｓ５３４４ Ｗｕ　ｅｔ　ａｌ．、Ｃｏｍｐａｒａｔｉｖｅ　ａｎａｌｙｓｉｓ　ｏｆ　ｍｉｃｒｏｂｉｏｍｅ　ｍｅａｓｕｒｅｍｅｎｔ　ｐｌａｔｆｏｒｍｓ　ｕｓｉｎｇ　ｌａｔｅｎｔ　ｖａｒｉａｂｌｅ　ｓｔｒｕｃｔｕｒａｌ　ｅｑｕａｔｉｏｎ　ｍｏｄｅｌｉｎｇ、ＢＭＣ　Ｂｉｏｉｎｆｏｒｍａｔｉｃｓ　２０１３　１４：７９、Ｐｕｂｌｉｓｈｅｄ　５　Ｍａｒｃｈ　２０１３、ｈｔｔｐ：／／ｗｗｗ．ｂｉｏｍｅｄｃｅｎｔｒａｌ．ｃｏｍ／１４７１－２１０５／１４／７９
　このように、腸内細菌叢の情報を用いて、疾病のリスクを評価するためのモデルについては実用化に向けての研究がまだされていない。このため、せっかく採便キットで採取した大便を提出したにもかかわらず、特定の疾病に対するリスクの評価を希望する者（以下、ユーザという）がいても、ユーザ自身の腸内細菌叢の情報から、ユーザが何らかの疾病になるリスクを評価するというサービスを提供することができていない。
　本発明は、このような課題に着目して鋭意研究され完成されたものであり、その目的は、腸内細菌叢を用いた疾病評価指標を算出することに関する技術を提供することにある。
　上記課題を解決するために、第１の発明は、疾病評価指標を算出するためのモデル作成方法であって、１又は複数の疾病の入力を受け付け、複数の被験者の便を解析した結果である腸内細菌叢に関するデータを格納する第１のデータベースと、前記被験者の属性に関するデータ及び疾病に関するデータを格納する第２のデータベースから、健康な人及び前記疾病に罹患している人の前記属性に関するデータ、前記腸内細菌叢に関するデータ、及び、前記疾病に関するデータを抽出し、前記属性に関するデータ、前記腸内細菌叢に関するデータ、及び、前記疾病に関するデータを内部因子とする連関モデルを作成するコンピュータによるモデル作成方法である。
　第２の発明は、前記腸内細菌叢に関するデータを１又は複数の前記内部因子とする第１の発明に記載のコンピュータによるモデル作成方法である。
　第３の発明は、前記連関モデルは構造方程式モデルである第１の発明に記載のコンピュータによるモデル作成方法である。
　第４の発明は、前記第２のデータベースは前記被験者の生活に関するデータ及び環境に関するデータをさらに格納し、前記被験者の属性に関するデータに加え、前記被験者の生活に関するデータ及び環境に関するデータも前記内部因子に加える第１の発明に記載のコンピュータによるモデル作成方法である。
　第５の発明は、前記連関モデルを作成する際に、前記疾病に関係する前記腸内細菌叢に関するデータを選定し、前記選定された腸内細菌叢に関するデータを観測変数とし、前記連関モデルの前記内部因子となる前記腸内細菌叢に関するデータを推定する第１の発明に記載のコンピュータによるモデル作成方法である。
　第６の発明は、前記連関モデルを作成する際に、疾病リスク係数表を用いる第１の発明に記載のコンピュータによるモデル作成方法である。
　第７の発明は、前記疾病リスク係数表は、所定の属性の時に、所定の腸内細菌叢と所定の疾病・健康状態とのパスに関するものである第６の発明に記載のコンピュータによるモデル作成方法である。
　第８の発明は、疾病評価指標を算出するためのモデル作成プログラムであって、１又は複数の疾病の入力を受け付ける入力ステップと、複数の被験者の便を解析した結果である腸内細菌叢に関するデータを格納する第１のデータベースと、前記被験者の属性に関するデータ及び疾病に関するデータを格納する第２のデータベースから、健康な人及び前記疾病に罹患している人の前記属性に関するデータ、前記腸内細菌叢に関するデータ、及び、前記疾病に関するデータを抽出するステップと、前記属性に関するデータ、前記腸内細菌叢に関するデータ、及び、前記疾病に関するデータを内部因子とする連関モデルを作成する作成ステップと、をコンピュータに実行させるモデル作成プログラムである。
　第９の発明は、ユーザが調べて欲しいと望む１又は複数の疾病の入力を受け付け、前記ユーザの便を予め解析した結果である前記ユーザの腸内細菌叢に関するデータ、前記ユーザの属性に関するデータ、並びに、複数の被験者の属性に関するデータ、前記被験者の便を解析した結果である前記被験者の腸内細菌叢に関するデータ、及び、前記疾病に関するデータを内部因子とする連関モデルをデータベースから抽出し、前記ユーザの腸内細菌叢に関するデータ、前記ユーザの属性に関するデータ、及び、前記連関モデルを用いて、前記疾病に対するリスクを算出するコンピュータによる疾病評価指標算出方法である。
　第１０の発明は、前記腸内細菌叢に関するデータを１又は複数の前記内部因子とする第９の発明に記載のコンピュータによる疾病評価指標算出方法である。
　第１１の発明は、前記連関モデルは構造方程式モデルである第９の発明に記載のコンピュータによる疾病評価指標算出方法である。
　第１２の発明は、前記内部因子は、前記被験者の属性に関するデータに加え、前記被験者の生活に関するデータ及び環境に関するデータも有する第９の発明に記載のコンピュータによる疾病評価指標算出方法である。
　第１３の発明は、ユーザが調べて欲しいと望む１又は複数の疾病を入力する入力部と、前記ユーザの便を予め解析した結果である前記ユーザの腸内細菌叢に関するデータ、前記ユーザの属性に関するデータ、並びに、複数の被験者の属性に関するデータ、前記被験者の便を解析した結果である前記被験者の腸内細菌叢に関するデータ、及び、前記疾病に関するデータを内部因子とする連関モデルをデータベースから抽出する抽出部と、前記ユーザの腸内細菌叢に関するデータ、前記ユーザの属性に関するデータ、及び、前記連関モデルを用いて、前記疾病に対するリスクを算出する算出部と、を備える疾病評価指標算出装置である。
　第１４の発明は、ユーザが調べて欲しいと望む１又は複数の疾病を入力する入力部と、前記ユーザの便を予め解析した結果である前記ユーザの腸内細菌叢に関するデータを第１のデータベースから、及び、前記ユーザの属性に関するデータを第２のデータベースから抽出する第１の抽出部と、健康な人及び前記疾病に罹患している人を含む複数の被験者の便を解析した結果である前記被験者の腸内細菌叢に関するデータを前記第１のデータベースから、及び、前記被験者の属性に関するデータ及び疾病に関するデータを前記第２のデータベースから抽出する第２の抽出部と、前記被験者の属性に関するデータ、前記被験者の腸内細菌叢に関するデータ、及び、前記疾病に関するデータを内部因子とする連関モデルを作成する作成部と、前記ユーザの腸内細菌叢に関するデータ、前記ユーザの属性に関するデータ、及び、前記連関モデルを用いて、前記疾病に対するリスクを算出する算出部と、を備える疾病評価指標算出システムである。
　第１５の発明は、ユーザが調べて欲しいと望む１又は複数の疾病の入力を受け付ける入力ステップと、前記ユーザの便を予め解析した結果である前記ユーザの腸内細菌叢に関するデータを第１のデータベースから、及び、前記ユーザの属性に関するデータを第２のデータベースから抽出する第１の抽出ステップと、健康な人及び前記疾病に罹患している人を含む複数の被験者の便を解析した結果である前記被験者の腸内細菌叢に関するデータを前記第１のデータベースから、及び、前記被験者の属性に関するデータ及び疾病に関するデータを前記第２のデータベースから抽出する第２の抽出ステップと、前記被験者の属性に関するデータ、前記被験者の腸内細菌叢に関するデータ、及び、前記疾病に関するデータを内部因子とする連関モデルを作成する作成ステップと、前記ユーザの腸内細菌叢に関するデータ、前記ユーザの属性に関するデータ、及び、前記連関モデルを用いて、前記疾病に対するリスクを算出する算出ステップと、をコンピュータに実行させる疾病評価指標算出プログラムである。
　本発明によれば、腸内細菌叢を用いた疾病評価指標を算出することに関する技術を提供することができる。
本発明の実施形態に係る疾病評価指標算出システムの全体概略図である。 本実施形態に係る腸内細菌叢解析装置２００及び腸内細菌叢ＤＢ３００の概略構成図である。 本実施形態に係る評価モデルの概念図である。 本実施形態に係るアンケート内容の説明図である。 本実施形態に係る各集団の説明図である。 実施例１に係る評価モデルの概略図である。 実施例１に係る評価モデル作成のフローチャートである。 実施例１に係る共分散構造分析の場合のモデルの概略説明図である。 実施例１に係る評価モデル作成フローチャートの一部の説明図である。 実施例１に係る評価モデルの具体例を示す図である。 実施例２に係る評価モデルの具体例を示す図である。
　図面を参照しながら本発明の実施の形態を説明する。なお、各図において共通する部分には同一の符号を付し、重複した説明は省略する。
　（疾病評価指標算出システム）
　図１は、本発明の実施形態に係る疾病評価指標算出システムの全体概略図である。疾病評価指標算出システムは、多数の被験者及び評価希望者（ユーザ）が提出する採便キットから腸内細菌叢データベース（ＤＢ）を作成するフェーズ、ユーザが評価を希望する１又は複数の疾病（特定の疾病）の評価モデルを作成するフェーズ、及び、特定の疾病に対するリスクの評価指標を算出するフェーズを備える。
　（腸内細菌叢ＤＢ作成フェーズ）
　腸内細菌叢ＤＢ作成フェーズでは、数千人以上という多数の被験者が各自の大便をトイレで排泄し、それを採便キットで採取する。次に、採便キットを受け取った抽出業者は、採便キットを腸内細菌ＤＮＡ抽出装置１００に入力し、腸内細菌叢に関するＤＮＡ溶液が出力される。
　ＤＮＡ溶液を受け取った解析業者は、ＤＮＡ溶液を腸内細菌叢解析装置２００に入力し、腸内細菌叢を解析する。解析業者は解析結果データを腸内細菌叢ＤＢ３００に格納する。なお、抽出業者と解析業者を分けて説明したが、同一の事業者が抽出作業と解析作業を行ってもよい。
　図２は、腸内細菌叢解析装置２００及び腸内細菌叢ＤＢ３００の概略構成図である。本実施形態では、腸内細菌叢解析装置２００は、ＤＮＡ溶液から腸内細菌叢を解析する方法として２種類用意するために、ＮＧＳ解析装置２１０及びＴ－ＲＦＬＰ解析装置２２０を備えている。ただし、この２種類のうち少なくとも一方を備え、解析してもよい。
　ＮＧＳ解析装置２１０では、予め試薬調整分注操作が行われたＤＮＡ溶液を入力し、次世代シーケンサ（Ｎｅｘｔ　Ｇｅｎｅｒａｔｉｏｎ　Ｓｅｑｕｅｎｃｅｒ）で解析することによって、腸内細菌叢に関するＤＮＡ配列が出力され、腸内細菌叢（ＤＮＡ配列）ＤＢ３１０に格納される。つまり、ＮＧＳでは、配列による細菌分類を行い、菌分類毎の細菌の相対量を解析することができる。
　Ｔ－ＲＦＬＰ解析装置２２０では、試薬分注と調整を施したＤＮＡ溶液を入力し、キャビラリーシーケンサがＴ－ＲＦＬＰ（Ｔｅｒｍｉｎａｌ　－　Ｒｅｓｔｒｉｃｔｉｏｎ　Ｆｒａｇｍｅｎｔ　Ｌｅｎｇｔｈ　Ｐｏｌｙｍｏｒｐｈｉｓｍ）で解析することによって、腸内細菌叢に関するＤＮＡフィンガープリントが出力され、腸内細菌叢（ＤＮＡフィンガープリント）ＤＢ３２０に格納される。つまり、Ｔ－ＲＦＬＰでは、ＤＮＡ断片長による細菌分類（ＯＴＵ）を行い、ＯＴＵ毎の細菌の相対量を解析することができる。ＯＴＵは操作的分類単位（Ｏｐｅｒａｔｉｏｎａｌ　Ｔａｘｏｎｏｍｉｃ　Ｕｎｉｔ）といい、「Ａ１９３」などの群に分類されるが、実際の菌属の名前はわからない。なお、「Ａ１９３」の「Ａ」は制限酵素
Figure JPOXMLDOC01-appb-I000001
の略称である。この制限酵素によってＤＮＡを切断している。
　腸内細菌叢ＤＢ３００は、腸内細菌叢（ＤＮＡ配列）ＤＢ３１０及び腸内細菌叢（ＤＮＡフィンガープリント）ＤＢ３２０を備えている。この２種類のＤＢ（３１０及び３２０）は上述した２種類の解析装置（２１０及び２２０）の解析結果をそれぞれ格納する。なお、２種類の解析装置の一方（２１０又は２２０）が腸内細菌叢解析装置２００に含まれている場合、対応するＤＢ（３１０又は３２０）が腸内細菌叢ＤＢ３００に含まれていればよい。このようなデータベース構造により、数千人以上という多数の被験者の腸内細菌叢が異なる方法で解析された場合であっても一つのデータベースとして取り扱うことができる。そして、腸内細菌叢ＤＢ３００は各被験者のＩＤ（識別番号）情報とその被験者の腸内細菌叢とを関係付けている。
　次に、採便キットを提出した多数の被験者は、アンケートを記入し、アンケート回収業者に提出する。アンケート回収業者はアンケート結果をアンケートＤＢ４００に格納する。ここで、アンケート回収の初期段階では、アンケートＤＢ４００に格納されているアンケート結果が少ないため、外部のアンケートＤＢをアンケートＤＢ４００として利用してもよい。同様に、腸内細菌叢ＤＢ３００に格納されているデータも少ないため、外部の腸内細菌叢ＤＢを腸内細菌叢ＤＢ３００として利用してもよい。
　（評価モデル作成フェーズ）
　評価モデル作成フェーズでは、まず、評価モデル作成業者は、評価したい属性と疾病を選定する。モデル作成装置５００は、選定された属性及び疾病の入力を受け付ける入力部と、これら属性及び疾病に関係する被験者のＩＤ情報をアンケートＤＢ４００から抽出し、該当するＩＤ情報の腸内細菌叢も腸内細菌叢ＤＢ３００から抽出する抽出部と、後述する作成方法によって、評価モデルを作成する作成部を備える。モデル作成装置５００は、作成した評価モデルを評価モデルＤＢ６００に格納する。
　（疾病評価指標算出フェーズ）
　疾病評価指標算出フェーズでは、評価業者は、採便キットで採取した大便を提出し、かつ、特定の疾病に対するリスクの評価を希望する者（以下、ユーザという）に対し、腸内細菌叢ＤＢ作成フェーズで説明した採便キットの提出、及び、評価モデル作成フェーズで説明したアンケートの記入を依頼する。ユーザの採便キット及びアンケートは、多数の被験者の採便キット及びアンケートと同様の処理が行われ、腸内細菌叢ＤＢ３００及びアンケートＤＢ４００に格納される。さらに、ユーザが評価を希望する疾病（特定の疾病）、及び、ユーザのＩＤ情報を疾病評価指標算出装置７００に入力する。
　疾病評価指標算出装置７００は、入力情報を受付ける入力部と、入力情報に関係する情報を抽出する抽出部と、評価指標を算出する算出部と、評価レポートを出力する出力部を備える。入力部は特定の疾病及びユーザのＩＤ情報の入力を受け付ける。抽出部は外部のアンケートＤＢ４００及び腸内細菌叢ＤＢ３００からユーザのＩＤ情報に関係する情報（ユーザの要因、及び、腸内細菌叢）を抽出し、評価モデルＤＢ６００から特定の疾病に関する評価モデルを抽出する。そして、疾病評価指標算出装置７００は、抽出されたユーザの要因（属性や生活習慣など）及び腸内細菌叢を、特定の疾病の評価モデルに入力し、特定の疾病に対するリスクの評価指標を算出する。算出された評価指標に基づく、評価レポートが出力される。なお、各フェーズを行う業者を分けて説明したが、同一の事業者が全てのフェーズを行ってもよい。
　図３は、評価モデルＤＢ６００に格納される評価モデルの一例を説明するための概念図である。評価モデルは、特定の疾病と一定以上の相関が認められた観測変数と、観測変数と関係する潜在変数（以下、因子という）によって構築される統計解析モデルである。なお、この評価モデルでは、要因（主な原因）も因子として取り扱っている。
　観測変数は実線の長方形で表し、因子は楕円形で表している。各観測変数は、図１の腸内細菌叢ＤＢ３００又はアンケートＤＢ４００の情報の一部である。
　評価モデルの主要部は、観測変数としての属性（性別・年齢）６１０－１、並びに、因子としての要因１（腸内細菌叢）６２０－１、要因２（生活・環境等）６２０－２及び因子（疾病等）６２０－３を備える。属性（性別・年齢）６１０－１はアンケートＤＢ４００の情報の一部である。矢印の向きは「影響を与える向き」を定義している。また、要因１（腸内細菌叢）６２０－１から因子（疾病等）６２０－３への太い矢印は一方向になっているが、双方向の評価モデルを構築することも可能である。
　次に、因子を詳細に説明する。点線で囲まれた長方形６３０－１は要因１（腸内細菌叢）６２０－１の詳細を示す。要因１（腸内細菌叢）６２０－１は腸内細菌叢因子６２０－１’及び複数の「ＯＴＵ又は菌属（菌種）」６１０－２に分解できる。複数の「ＯＴＵ又は菌属（菌種）」６１０－２は観測変数であり、図１の腸内細菌叢ＤＢ３００の情報の一部である。
　点線で囲まれた長方形６３０－２は要因２（生活・環境等）６２０－２の詳細を示す。要因２（生活・環境等）６２０－２は、環境要因６２０－２’－１、生活環境（アンケート）６１０－３及び家庭環境（アンケート）６１０－４、並びに、食生活要因６２０－２’－２、食行動（アンケート）６１０－５及び食意識（アンケート）６１０－６に分解できる。生活環境（アンケート）６１０－３及び家庭環境（アンケート）６１０－４、並びに、食行動（アンケート）６１０－５及び食意識（アンケート）６１０－６は観測変数であり、図１のアンケートＤＢ４００の情報の一部である。
　点線で囲まれた長方形６３０－３は因子（疾病等）６２０－３の詳細を示す。因子（疾病等）６２０－３は、疾病・健康状態因子６２０－３’、疾病（アンケート）６１０－７及び健康状態（アンケート）６１０－８に分解できる。疾病（アンケート）６１０－７及び健康状態（アンケート）６１０－８は観測変数であり、図１のアンケートＤＢ４００の情報の一部である。
　この図では、ユーザを属性（性別・年齢）６１０－１、並びに、各アンケート６１０－３乃至６１０－８で分割された集団に割り当て、その集団が保有する要因１（腸内細菌叢）６２０－１の場合、ユーザが評価を希望する疾病（特定の疾病）に対するリスクは高いか低いかを算出することができることを概念的に表している。ここで、特定の疾病に対するリスクとは、特定の疾病に対して現在罹患している、又は、将来罹患する確率（可能性も含む）をいう。
　すなわち、特定の疾病と一定以上の相関が認められた観測変数として「ＯＴＵ又は菌属（菌種）」６１０－２、並びに、属性（性別・年齢）６１０－１、生活環境（アンケート）６１０－３、家庭環境（アンケート）６１０－４、食行動（アンケート）６１０－５、食意識（アンケート）６１０－６、疾病（アンケート）６１０－７、及び、健康状態（アンケート）６１０－８を選択し、これらの観測変数と関係する因子によって評価モデルが構築される。
　ここで、属性６１０－１については、統計解析を行う母集団を男女別に分割し、性別及び年齢を観測変数として用いている。また、腸内細菌叢因子６２０－１’の観測変数として、ＯＴＵ又は菌属（菌種）名６１０－２を用いている。これは、Ｔ－ＲＦＬＰによりＤＮＡフィンガープリントがわかる場合、ＯＴＵを用い、ＮＧＳによりＤＮＡ配列が分かる場合、菌属（菌種）名を用いるためである。
　図４は、本実施形態に係るアンケート内容の説明図である。アンケートの構成概念は、「属性」、「生活・家庭環境（ライフスタイル）」、「食意識」、「食行動」、「心理的負荷（ストレス）」、「投薬」、「病気・未病（健康状態）」であり、これらに該当する各質問項目は一つ又は複数である。また、「７－３身体状況」には、身体状況から判断される健康状態、例えば、アトピー性皮膚炎に関するものとして、かゆみや皮膚の赤みという項目、又は、自分で測定できる血圧などを含まれる。ここで、説明の簡略化のため、「病気・未病（健康状態）」は「疾病」と述べてもよく、「生活・家庭環境（ライフスタイル）」、「食意識」、「食行動」、「心理的負荷（ストレス）」及び「投薬」はまとめて「生活及び環境」と述べてもよい。
　図５は、本実施形態に係る各集団の説明図である。全被験者で構成される母集団、疾病にかかっていない健康母集団、アトピーにかかっている集団、アトピー及び喘息にかかっている集団、喘息にかかっている集団で構成されている。
　（実施例１：属性と腸内細菌叢から疾病リスクを評価する場合）
　実施例１では、図３で説明した評価モデルの主要部のうち、観測変数としての属性（性別・年齢）、並びに、因子としての要因１（腸内細菌叢）及び因子（疾病等）を用いて、疾病リスクを評価する場合について説明する。
　図６は、実施例１に係る評価モデルの概略図である。実施例１の評価モデルは、疾病「アトピー」又は「喘息」に対するリスクの評価指標を算出するためのモデルである。観測変数「性別」及び「年齢」は被験者が提出したアンケートから抽出し、観測変数「Ａ５８」、「Ａ１３１」などは被験者が提出した採便キットから抽出している。また、観測変数「アトピー」及び「喘息」が、ユーザが評価を希望する疾病である。この評価モデルを用いることによって、ユーザが「アトピー」及び「喘息」になるリスクの評価指標を算出することができる。
　実施例１の評価モデルの内部因子は「属性」、「腸内細菌叢」及び「疾病・健康状態」である。これらの内部因子同士の関係が、連関モデルとして評価モデルの概略構成で表現されている。本願では、こうした観測変数、内部因子同士の多様な関係性を統計学的に評価する分析方法として共分散構造分析を採用する。図６に示す観測変数と内部因子との関係、及び、内部因子同士の関係は構造方程式モデル（ＳＥＭ；Ｓｔｒｕｃｔｕｒａｌ　Ｅｑｕａｔｉｏｎ　Ｍｏｄｅｌｉｎｇ）で表現され、評価モデルが生成される。
　また、内部因子と観測変数との関係は次のようになる。内部因子「属性」の観測変数が「性別」及び「年齢」であり、内部因子「腸内細菌叢」の観測変数が「Ａ５８」、「Ａ１３１」などのＯＴＵであり、さらに、内部因子「疾病・健康状態」の観測変数が「アトピー」及び「喘息」である。なお、説明の簡略化のため、内部因子「疾病・健康状態」は内部因子「疾病」と述べてもよい。
　図７は、実施例１に係る評価モデル作成のフローチャートである。ここでは、疾病の一例として「アトピー」及び「喘息」について説明するが、他の疾病、例えば、「心臓病」、「大腸がん」、「胃腸病」、「肝臓病」、「腎臓病」、「糖尿病」、「高血圧」、「脂質異常」、「低体重（ＢＭＩ）」、「肥満（ＢＭＩ）」、「骨病」、「関節痛」、「うつ」などについても同様に評価モデルを作成することができる
　まず、疾病として「アトピー」及び「喘息」を選択する（Ｓ１００）。評価に使用する属性として「性別」及び「年齢」を使用することを選択する（Ｓ１１０）。
　次に、疾病「アトピー」と関係の強いＯＴＵを選定する（Ｓ１２０）。この選定は、例えばロジスティック回帰分析に基づいて行われてもよい。ここで、所定の疾病と関係の強いＯＴＵをロジスティック回帰分析により選定する際に、所定の属性の時に、所定の腸内細菌叢と所定の疾病・健康状態とのパス、そのパスの係数、オッズ比、及び、ｐ値（有意確率）で構成される疾病リスク係数表を用いる。
　次いで、選定されたＯＴＵを観測変数として、評価モデルの内部因子「腸内細菌叢」を推定する（Ｓ１３０）。この推定は、例えば因子分析に基づいて行われてもよく、あるいは主成分分析、主座標分析等、他の分析方法に基づいてもよい。また、疾病「喘息」についても同様の処理（Ｓ１２０及びＳ１３０）を行う。
　Ｓ１２０と平行して（時間の前後は問わない）、疾病「アトピー」に関係する、評価モデルの内部因子としての属性を疾病リスク係数表より選定する（Ｓ１４０）。そして、選定された属性に対し、評価モデルの内部因子「属性」を推定する（Ｓ１５０）。この推定は、例えば因子分析に基づいてもよく、あるいは主成分分析、主座標分析等、他の分析方法に基づいてもよい。また、疾病「喘息」についても同様の処理（Ｓ１４０及びＳ１５０）を行う。
　Ｓ１３０で推定された評価モデルの内部因子「腸内細菌叢」及び、Ｓ１４０で推定された内部因子「属性」、さらに内部因子「疾病・健康状態」によりＳＥＭ解析を実施する（Ｓ１６０）。
そして、図８に示す構造方程式のパラメータα、β、γを求め、評価モデルを作成する（Ｓ１７０）。
　図９は、実施例１に係る評価モデル作成フローチャートの一部（Ｓ１２０）の説明図である。特定の疾病に関して、疾病リスク係数表からＯＴＵを選定するステップについて説明する。
　まず、（１）疾病を決める。ここでは「アトピー」及び「喘息」とする。（２）係数表からＯＴＵを選定する。ここでは「Ａ５８」、「Ａ１３１」、「Ａ４７７」及び「Ａ４９８」等を記載しているが、実際には各ＯＴＵの次元は４０～１００次元程度になる。
　次に、（３）図５で説明した各集団から母集団を選定する。（４）所定の属性として「年齢」を選択する。ここでは、図５の各集団は男性のみ、女性のみ、及び、全体に予め分かれており、母集団毎の評価モデルを構築するため、「年齢」を選択する。（５）ＳＥＭ解析を実施し、（６）各パスの計数と、モデルの適合度を確認、すなわち、所定のしきい値以上であれば終了する。
　図１０は、実施例１に係るアトピー及び喘息の場合の評価モデルの具体例である。各パスの矢印の向きは影響を与える向きを表し、各パスの数値が１に近いほど正の相関度が高い、すなわち、影響を与える要素として高いことを示している。一方、中黒の三角印はマイナスの数値であり、負の相関であることを示している。
　図１０の評価モデルの具体例によれば、疾病「アトピー」及び「喘息」については、ある被験者は「年齢」が高く、「Ａ１３１」に属する菌を多く保有している場合、「年齢」と「疾病因子」とは正の相関があるため、「年齢」が高いほど「アトピー」及び「喘息」にかかりやすい、と評価できる。しかしながら、「Ａ１３１」と「腸内細菌叢Ａ」は正の相関があるため、「腸内細菌叢Ａ」が高いことを意味する。そして、「アトピー」及び「喘息」と「疾病因子」は正の相関があり、かつ、「腸内細菌叢Ａ」と「疾病因子」は負の相関があるため、「腸内細菌叢Ａ」が高ければ、「アトピー」及び「喘息」にはかかりにくい、と評価できる。
　さらに、「年齢」から「疾病因子」へのパスは、直接的なパスと、「腸内細菌叢Ａ」を経由する間接的なパスがあることに留意する。直接的なパスによれば、上述した通り「年齢」が高いほど「アトピー」及び「喘息」にかかりやすい。一方、間接的なパスによれば、年齢が高くなると、「腸内細菌叢Ａ」も少し高くなり、「アトピー」及び「喘息」にはかかりにくい。このように間接的なパスも考慮すると、「年齢」が高いほど「アトピー」及び「喘息」にかかりにくいという影響は実は小さいと評価できる。
　評価モデルの適合度指標として３種類用いており、ＧＦＩはＧｏｏｄｎｅｓｓ　ｏｆ　Ｆｉｔ　Ｉｎｄｅｘ、ＡＧＦＩはＡｄｊｕｓｔｅｄ　ＧＦＩ、ＲＭＳＥＡはＲｏｏｔ　Ｍｅａｎ　Ｓｑｕａｒｅ　Ｅｒｒｏｒ　ｏｆ　Ａｐｐｒｏｘｉｍａｔｉｏｎの略称である。ＧＦＩ及びＡＧＦＩは値が１に近くなるほどデータへの当てはまり（以下、適合度という）が良いことを意味する。一方、ＲＭＳＥＡは値が小さいほど適合度が良いことを意味する。また、評価モデルには、誤差のバラツキの指標である誤差分散を用いてもよい。
　図１０の評価モデルの適合度指標は、ＧＦＩが０．９７８８、ＡＧＦＩが０．９７３５、ＲＭＳＥＡが０．１０９６であり、この評価モデルの適合度は高いことが示されている。
　（効果）
　実施例１によれば、属性と腸内細菌叢から疾病リスクを評価する評価モデルを作成することができる。すなわち、図３の要因１（腸内細菌叢）６２０－１から因子（疾病等）６２０－３への太い矢印の関連性を評価するモデルを作成することができる。また、図６で説明したように観測変数、内部因子同士の多様な関係性をモデル化しているため、図１０で説明したように直接的なパス及び間接的なパスの両方を判断可能な評価モデルを作成することができる。
　そして、ユーザは採便キットで大便を採取、提出し、さらにアンケートに回答することによって、自分の大便が「Ａ１３１」を多く保有していれば、年齢による影響は小さく、「腸内細菌叢Ａ」が高く、「アトピー」及び「喘息」にかかりにくい、という評価を受けることができる。
　そして、疾病評価指標算出装置７００は、「アトピー」及び「喘息」になるリスクの評価を希望するユーザから、ユーザＩＤ情報を受け付け、アンケートＤＢ４００及び腸内細菌叢ＤＢ３００からユーザのＩＤ情報に関係する情報（ユーザの属性、及び、腸内細菌叢の情報）を抽出し、評価モデルＤＢ６００からユーザが評価を希望する疾病「アトピー」及び「喘息」に関する評価モデルを抽出する。そして、疾病評価指標算出装置７００は、抽出されたユーザの属性（この場合、年齢）及び腸内細菌叢の情報（この場合、「Ａ１３１」を多く保有）を、図１０の評価モデルに入力し、「アトピー」及び「喘息」に対するリスクの評価指標を算出することができる。さらに、算出された評価指標に基づく、評価レポートが出力される。
　ユーザは評価レポートを読むことによって、ユーザ自身の腸内細菌叢に、腸内細菌叢Ａ（因子）への因子負荷量が大きい「Ａ１３１」が多く含まれているため、年齢による影響は小さく（年齢が高くなっても、将来的に）「アトピー」及び「喘息」に対するリスクは低いということを理解できる。
　次に、疾病評価指標算出システム事業者にとっての効果を述べる。ユーザの腸内細菌叢に「Ａ１３１」が多く含まれていることが「アトピー」及び「喘息」に対するリスクを低減させるのであれば、「Ａ１３１」に含まれる実在の菌を摂取すれば、「アトピー」及び「喘息」に対するリスクが低減される可能性を評価レポートは示唆している。この可能性に基づき、事業者はプロバイオティクスの設計、検討、及び提案を行うことができる。
　（実施例２；属性並びに、生活及び環境と、腸内細菌叢から疾病リスクを評価する場合）
　実施例２では、図３で説明した要因２（生活・環境等）を実施例１に加え、すなわち、図３で説明した評価モデルの主要部全てを用いて、疾病リスクを評価する場合について説明する。実施例２の疾病評価指標算出システム７００も、図１の全体概略図と同じである。評価モデルを作成する際に、要因２（生活・環境等）を用いている点が、実施例１と異なる点である。
図１１は、実施例２に係る評価モデルの具体例を示す図である。実施例２特有の評価モデルの内部因子は、図３の要因２（生活・環境等）として、「食生活」を用いている。内部因子「食生活」の観測変数が「牛乳」、「パン」、「肉」等である。
　観測変数「牛乳」、「パン」、「肉」等を評価モデルに加えることは、図４で説明したアンケートの構成概念「食行動」の各質問項目をアンケート情報に追加することによって可能になる。また、各質問項目の数を増やすことによって、より詳細な評価をすることが可能になる。
　図１１の評価モデルの具体例によれば、実施例１の解釈に加え、ある被験者が「肉」を食する頻度が高い場合、「食生活」因子が高いことを意味する。これは、「腸内細菌叢Ａ」因子が低下し、「腸内細菌叢Ｂ」因子が増加し、さらに、「疾病因子」の増加を意味する。この「疾病因子」は「アトピー」及び「喘息」になるリスクを増加させていることを示唆している。また、この「食生活」因子が高いことにより、「腸内細菌叢Ｂ」因子が高くなり、その高い「腸内細菌叢Ｂ」因子が「疾病因子」を高めている可能性をさらに示唆している、と評価できる。
　図１１の評価モデルの適合度指標は、ＧＦＩが０．９２６５であり、ＡＧＦＩが０．９０４８であり、ＲＭＳＥＡが０．０７３８である。このため、この評価モデルの適合度は高いことが示されている。
　（効果）
　実施例２によれば、属性及び生活・環境等と、腸内細菌叢から疾病リスクを評価する評価モデルを作成することができる。図１１の場合、「年齢」、「腸内細菌叢Ａ」、「腸内細菌叢Ｂ」、「食生活」から「アトピー」及び「喘息」になるリスクを評価する評価モデルを作成することができる。
　そして、疾病評価指標算出装置７００は、「アトピー」及び「喘息」になるリスクの評価を希望するユーザから、ユーザＩＤ情報を受け付け、アンケートＤＢ４００及び腸内細菌叢ＤＢ３００からユーザのＩＤ情報に関係する情報（ユーザの属性及び「食生活」に関するアンケート結果、並びに、腸内細菌叢の情報）を抽出し、評価モデルＤＢ６００からユーザが評価を希望する疾病「アトピー」及び「喘息」に関する評価モデルを抽出する。そして、疾病評価指標算出装置７００は、抽出されたユーザの要因「年齢、食生活」及び腸内細菌叢の情報を、「アトピー」及び「喘息」の評価モデルに入力し、「アトピー」及び「喘息」に対するリスクの評価指標を算出することができる。さらに、算出された評価指標に基づく、評価レポートが出力される。
　ユーザは評価レポートを読むことによって、実施例１の評価に加え、「肉」を食する頻度が高いため、「腸内細菌叢Ａ」因子が低下し、「腸内細菌叢Ｂ」因子が増加する傾向にあり、「アトピー」及び「喘息」になるリスクを増加させていることを示唆している。また、「肉」を食する頻度が高いことにより、「腸内細菌叢Ｂ」因子が高くなり、「アトピー」及び「喘息」になるリスクを高めている可能性をさらに示唆している、という詳細な評価を理解できる。
　次に、疾病評価指標算出システム事業者にとっての効果を述べる。ユーザの「肉」を食する頻度が高く、「野菜」を食する頻度が低いことが「アトピー」及び「喘息」に対するリスクを増加させるのであれは、「肉」を食する頻度を抑え、また「野菜」を多く摂取すれば、「アトピー」及び「喘息」に対するリスクが低減される可能性を評価レポートは示唆している。このような可能性に基づき、事業者はプレバイオティクスの設計、検討、及び提案を行うことができる。
　以上、本発明の実施例（変形例を含む）について説明してきたが、これらのうち、２つ以上の実施例を組み合わせて実施しても構わない。あるいは、これらのうち、１つの実施例を部分的に実施しても構わない。さらには、これらのうち、２つ以上の実施例を部分的に組み合わせて実施しても構わない。
　また、本発明は、上記発明の実施例の説明に何ら限定されるものではない。特許請求の範囲の記載を逸脱せず、当業者が容易に想到できる範囲で種々の変形態様もこの発明に含まれる。
　例えば、実施例の腸内細菌叢解析装置２００としてｑＰＣＲなどの他の解析装置を用いてもよい。ｑＰＣＲ解析装置の場合、試薬分注と調整を施したＤＮＡ溶液を入力し、リアルタイムの定量ＰＣＲ（ｑｕａｎｔｉｔａｔｉｖｅ　Ｐｏｌｙｍｅｒａｓｅ　Ｃｈａｉｎ　Ｒｅａｃｔｉｏｎ）で解析することによって、腸内細菌叢に関する指標菌群が出力され、腸内細菌叢ＤＢ３００内の腸内細菌叢（指標菌群）ＤＢに格納される。
　また、実施例では、連関モデルを解析する方法として、共分散構造分析を用いたが、パス解析などの他の解析方法を用いてもよい。
１００　腸内細菌ＤＮＡ抽出装置
２００　腸内細菌叢解析装置
３００　腸内細菌叢ＤＢ
４００　アンケートＤＢ
５００　モデル作成装置
６００　評価モデルＤＢ
７００　疾病評価指標算出装置

Claims (15)

  1. 　疾病評価指標を算出するためのモデル作成方法であって、
    　１又は複数の疾病の入力を受け付け、
    　複数の被験者の便を解析した結果である腸内細菌叢に関するデータを格納する第１のデータベースと、前記被験者の属性に関するデータ及び疾病に関するデータを格納する第２のデータベースから、健康な人及び前記疾病に罹患している人の前記属性に関するデータ、前記腸内細菌叢に関するデータ、及び、前記疾病に関するデータを抽出し、
    　前記属性に関するデータ、前記腸内細菌叢に関するデータ、及び、前記疾病に関するデータを内部因子とする連関モデルを作成するコンピュータによるモデル作成方法。
  2. 　前記腸内細菌叢に関するデータを１又は複数の前記内部因子とする請求項１に記載のコンピュータによるモデル作成方法。
  3. 　前記連関モデルは構造方程式モデルである請求項１に記載のコンピュータによるモデル作成方法。
  4. 　前記第２のデータベースは前記被験者の生活に関するデータ及び環境に関するデータをさらに格納し、
    　前記被験者の属性に関するデータに加え、前記被験者の生活に関するデータ及び環境に関するデータも前記内部因子に加える請求項１に記載のコンピュータによるモデル作成方法。
  5. 　前記連関モデルを作成する際に、
    　前記疾病に関係する前記腸内細菌叢に関するデータを選定し、
    　前記選定された腸内細菌叢に関するデータを観測変数とし、前記連関モデルの前記内部因子となる前記腸内細菌叢に関するデータを推定する請求項１に記載のコンピュータによるモデル作成方法。
  6. 　前記連関モデルを作成する際に、疾病リスク係数表を用いる請求項１に記載のコンピュータによるモデル作成方法。
  7. 　前記疾病リスク係数表は、所定の属性の時に、所定の腸内細菌叢と所定の疾病・健康状態とのパスに関するものである請求項６に記載のコンピュータによるモデル作成方法。
  8. 　疾病評価指標を算出するためのモデル作成プログラムであって、
    　１又は複数の疾病の入力を受け付ける入力ステップと、
    　複数の被験者の便を解析した結果である腸内細菌叢に関するデータを格納する第１のデータベースと、前記被験者の属性に関するデータ及び疾病に関するデータを格納する第２のデータベースから、健康な人及び前記疾病に罹患している人の前記属性に関するデータ、前記腸内細菌叢に関するデータ、及び、前記疾病に関するデータを抽出するステップと、
    　前記属性に関するデータ、前記腸内細菌叢に関するデータ、及び、前記疾病に関するデータを内部因子とする連関モデルを作成する作成ステップと、
    をコンピュータに実行させるモデル作成プログラム。
  9. 　ユーザが調べて欲しいと望む１又は複数の疾病の入力を受け付け、
    　前記ユーザの便を予め解析した結果である前記ユーザの腸内細菌叢に関するデータ、前記ユーザの属性に関するデータ、並びに、複数の被験者の属性に関するデータ、前記被験者の便を解析した結果である前記被験者の腸内細菌叢に関するデータ、及び、前記疾病に関するデータを内部因子とする連関モデルをデータベースから抽出し、
    　前記ユーザの腸内細菌叢に関するデータ、前記ユーザの属性に関するデータ、及び、前記連関モデルを用いて、前記疾病に対するリスクを算出するコンピュータによる疾病評価指標算出方法。
  10. 前記腸内細菌叢に関するデータを１又は複数の前記内部因子とする請求項７に記載のコンピュータによる疾病評価指標算出方法。
  11. 　前記連関モデルは構造方程式モデルである請求項７に記載のコンピュータによる疾病評価指標算出方法。
  12. 　前記内部因子は、前記被験者の属性に関するデータに加え、前記被験者の生活に関するデータ及び環境に関するデータも有する請求項７に記載のコンピュータによる疾病評価指標算出方法。
  13. 　ユーザが調べて欲しいと望む１又は複数の疾病を入力する入力部と、
    　前記ユーザの便を予め解析した結果である前記ユーザの腸内細菌叢に関するデータ、前記ユーザの属性に関するデータ、並びに、複数の被験者の属性に関するデータ、前記被験者の便を解析した結果である前記被験者の腸内細菌叢に関するデータ、及び、前記疾病に関するデータを内部因子とする連関モデルをデータベースから抽出する抽出部と、
    　前記ユーザの腸内細菌叢に関するデータ、前記ユーザの属性に関するデータ、及び、前記連関モデルを用いて、前記疾病に対するリスクを算出する算出部と、
    を備える疾病評価指標算出装置。
  14. 　ユーザが調べて欲しいと望む１又は複数の疾病を入力する入力部と、
    　前記ユーザの便を予め解析した結果である前記ユーザの腸内細菌叢に関するデータを第１のデータベースから、及び、前記ユーザの属性に関するデータを第２のデータベースから抽出する第１の抽出部と、
    　健康な人及び前記疾病に罹患している人を含む複数の被験者の便を解析した結果である前記被験者の腸内細菌叢に関するデータを前記第１のデータベースから、及び、前記被験者の属性に関するデータ及び疾病に関するデータを前記第２のデータベースから抽出する第２の抽出部と、
    　前記被験者の属性に関するデータ、前記被験者の腸内細菌叢に関するデータ、及び、前記疾病に関するデータを内部因子とする連関モデルを作成する作成部と、
    　前記ユーザの腸内細菌叢に関するデータ、前記ユーザの属性に関するデータ、及び、前記連関モデルを用いて、前記疾病に対するリスクを算出する算出部と、
    を備える疾病評価指標算出システム。
  15. 　ユーザが調べて欲しいと望む１又は複数の疾病の入力を受け付ける入力ステップと、
    　前記ユーザの便を予め解析した結果である前記ユーザの腸内細菌叢に関するデータを第１のデータベースから、及び、前記ユーザの属性に関するデータを第２のデータベースから抽出する第１の抽出ステップと、
    　健康な人及び前記疾病に罹患している人を含む複数の被験者の便を解析した結果である前記被験者の腸内細菌叢に関するデータを前記第１のデータベースから、及び、前記被験者の属性に関するデータ及び疾病に関するデータを前記第２のデータベースから抽出する第２の抽出ステップと、
    　前記被験者の属性に関するデータ、前記被験者の腸内細菌叢に関するデータ、及び、前記疾病に関するデータを内部因子とする連関モデルを作成する作成ステップと、
    　前記ユーザの腸内細菌叢に関するデータ、前記ユーザの属性に関するデータ、及び、前記連関モデルを用いて、前記疾病に対するリスクを算出する算出ステップと、
    をコンピュータに実行させる疾病評価指標算出プログラム。
PCT/JP2019/032698 2018-08-23 2019-08-21 疾病評価指標算出方法、装置、システム、及び、プログラム、並びに、疾病評価指標を算出するためのモデル作成方法 WO2020040219A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201980051821.9A CN112534513A (zh) 2018-08-23 2019-08-21 疾病评价指标算出方法、装置、系统、及程序、以及用于算出疾病评价指标的模型制成方法
KR1020217001054A KR20210049081A (ko) 2018-08-23 2019-08-21 질병 평가 지표 산출 방법, 장치, 시스템 및 프로그램, 그리고 질병 평가 지표를 산출하기 위한 모델 작성 방법
EP19852082.7A EP3832664A4 (en) 2018-08-23 2019-08-21 METHOD, DEVICE, SYSTEM AND PROGRAM FOR CALCULATING A DISEASE ASSESSMENT INDEX, AND MODEL CREATION METHOD FOR CALCULATING THE DISEASE ASSESSMENT INDEX
US17/181,152 US20210174967A1 (en) 2018-08-23 2021-02-22 Disease evaluation index calculation method, apparatus, system, and program, and model creation method for calculating disease evaluation index

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2018156379A JP6533930B1 (ja) 2018-08-23 2018-08-23 疾病評価指標算出方法、装置、システム、及び、プログラム、並びに、疾病評価指標を算出するためのモデル作成方法。
JP2018-156379 2018-08-23

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US17/181,152 Continuation US20210174967A1 (en) 2018-08-23 2021-02-22 Disease evaluation index calculation method, apparatus, system, and program, and model creation method for calculating disease evaluation index

Publications (1)

Publication Number Publication Date
WO2020040219A1 true WO2020040219A1 (ja) 2020-02-27

Family

ID=67023697

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2019/032698 WO2020040219A1 (ja) 2018-08-23 2019-08-21 疾病評価指標算出方法、装置、システム、及び、プログラム、並びに、疾病評価指標を算出するためのモデル作成方法

Country Status (6)

Country Link
US (1) US20210174967A1 (ja)
EP (1) EP3832664A4 (ja)
JP (1) JP6533930B1 (ja)
KR (1) KR20210049081A (ja)
CN (1) CN112534513A (ja)
WO (1) WO2020040219A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3974536A4 (en) * 2020-07-30 2022-08-17 Noster Inc. BACTERIA SUPPLY SYSTEM
WO2022203350A1 (ko) * 2021-03-26 2022-09-29 주식회사 에이치이엠파마 머신러닝 모델을 이용하여 아토피 유무를 판별하는 방법 및 진단장치

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6863633B1 (ja) * 2020-03-24 2021-04-21 一般社団法人日本農業フロンティア開発機構 菌群抽出方法、装置、及び、プログラム、並びに、腸内細菌叢ｄｂ作成システム、及び、方法。
JP6986650B1 (ja) * 2020-07-16 2021-12-22 威久 山本 川崎病またはＩｇＡ血管炎での合併症発生リスク予測方法、学習済みモデル及びその生成方法
JP7193810B1 (ja) 2021-06-25 2022-12-21 京都府公立大学法人 疾患リスク評価のための腸内細菌叢のタイプ分類方法
CN113628714B (zh) * 2021-07-30 2022-04-19 美益添生物医药(武汉)有限公司 针对疾病的营养素干预方法、系统、设备及存储介质
JP7270143B1 (ja) 2022-05-30 2023-05-10 シンバイオシス・ソリューションズ株式会社 疾病評価指標算出システム、方法、及び、プログラム
CN115472231B (zh) * 2022-10-11 2024-04-09 河南省人民医院 一种基于肠道菌群的2型糖尿病检测系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015533213A (ja) * 2012-10-05 2015-11-19 ネステク ソシエテ アノニム Ｉｂｓの診断マーカーとしての微生物バイオームに対する抗体、ストレス因子及びマスト細胞マーカー
US20160326574A1 (en) * 2013-11-01 2016-11-10 Washington University Methods to establish and restore normal gut microbiota function of subject in need thereof
JP2017029133A (ja) * 2015-07-21 2017-02-09 ＴＡＫ−Ｃｉｒｃｕｌａｔｏｒ株式会社 身体状態の評価方法、情報の提示方法、および身体状態を改善又は予防する物質のスクリーニング方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070118398A1 (en) * 2005-11-18 2007-05-24 Flicker Technologies, Llc System and method for estimating life expectancy and providing customized advice for improving life expectancy
US9289418B2 (en) * 2013-03-15 2016-03-22 Cedars-Sinai Medical Center Methods of diagnosis, selection, and treatment of diseases and conditions caused by or associated with methanogens
KR101445243B1 (ko) * 2014-03-28 2014-09-29 서울대학교산학협력단 장내 세균의 군집과 기능의 변화를 이용한 대사성 및 염증성 질환의 조기진단
EP3453396A1 (en) * 2014-08-28 2019-03-13 Yale University Compositions and methods for treating an inflammatory disease or disorder
CA2962466C (en) * 2014-10-21 2023-01-10 uBiome, Inc. Method and system for microbiome-derived diagnostics and therapeutics
US20180342322A1 (en) * 2014-10-21 2018-11-29 uBiome, Inc. Method and system for characterization for appendix-related conditions associated with microorganisms
EP3539033A4 (en) * 2016-11-14 2020-06-03 Cognoa, Inc. METHOD AND DEVICE FOR EVALUATING DEVELOPMENT FAULTS AND CHECKING COVERAGE AND RELIABILITY
CN107506582A (zh) * 2017-08-16 2017-12-22 周礼红 一种基于肠道微生物的健康风险预测系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015533213A (ja) * 2012-10-05 2015-11-19 ネステク ソシエテ アノニム Ｉｂｓの診断マーカーとしての微生物バイオームに対する抗体、ストレス因子及びマスト細胞マーカー
US20160326574A1 (en) * 2013-11-01 2016-11-10 Washington University Methods to establish and restore normal gut microbiota function of subject in need thereof
JP2017029133A (ja) * 2015-07-21 2017-02-09 ＴＡＫ−Ｃｉｒｃｕｌａｔｏｒ株式会社 身体状態の評価方法、情報の提示方法、および身体状態を改善又は予防する物質のスクリーニング方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
ADACHI KOHEI: "Passages; Handbook of multivariate analysis examples", HANDBOOK OF MULTIVARIATE ANALYSIS EXAMPLES, 25 June 2002 (2002-06-25), pages 287 - 293, 349-368, XP009526434, ISBN: 4-254-12147-4 *
LAHTI ET AL.: "Tipping elements in the human intestinal ecosystem, nature", NATURE COMMUNICATIONS, vol. 5, 8 July 2014 (2014-07-08), Retrieved from the Internet <URL:https://www.nature.com/articles/ncomms5344>
See also references of EP3832664A4 *
SHIGEKI YOKOYAMA; KIMIKO MATSUOKA: "Passage; Medical data mining that anyone can understand: Utilization of big data", MEDICAL DATA MINING THAT ANYONE CAN UNDERSTAND: UTILIZATION OF BIG DATA, 14 February 2014 (2014-02-14), pages 176 - 185, XP009526433, ISBN: 978-4-9907308-0-2 *
WU ET AL.: "Comparative analysis of microbiome measurement platforms using latent variable structural equation modeling", BMC BIOINFORMATICS 2013, vol. 14, 5 March 2013 (2013-03-05), pages 79, XP021140701, Retrieved from the Internet <URL:http://www.biomedcentral.com/1471-2105/14/79> DOI: 10.1186/1471-2105-14-79

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3974536A4 (en) * 2020-07-30 2022-08-17 Noster Inc. BACTERIA SUPPLY SYSTEM
WO2022203350A1 (ko) * 2021-03-26 2022-09-29 주식회사 에이치이엠파마 머신러닝 모델을 이용하여 아토피 유무를 판별하는 방법 및 진단장치

Also Published As

Publication number Publication date
JP2020030662A (ja) 2020-02-27
EP3832664A1 (en) 2021-06-09
KR20210049081A (ko) 2021-05-04
CN112534513A (zh) 2021-03-19
US20210174967A1 (en) 2021-06-10
EP3832664A4 (en) 2021-10-13
JP6533930B1 (ja) 2019-06-26

Similar Documents

Publication Publication Date Title
WO2020040219A1 (ja) 疾病評価指標算出方法、装置、システム、及び、プログラム、並びに、疾病評価指標を算出するためのモデル作成方法
Poussin et al. Interrogating the microbiome: experimental and computational considerations in support of study reproducibility
Meyer et al. Critical assessment of metagenome interpretation: the second round of challenges
Blanco-Míguez et al. Extending and improving metagenomic taxonomic profiling with uncharacterized species using MetaPhlAn 4
Wienke Frailty models in survival analysis
JP2020030800A (ja) 疾病評価指標算出方法、装置、システム、及び、プログラム、並びに、疾病評価指標を算出するためのモデル作成方法。
KR102258899B1 (ko) 통합적 건강 정보를 이용한 식단 및 운동 추천 방법 및 서비스 시스템
US20150199477A1 (en) Method, an arrangement and a computer program product for analysing a biological or medical sample
McMahan et al. Bayesian regression for group testing data
Zhu et al. OGUs enable effective, phylogeny-aware analysis of even shallow metagenome community structures
Zhuang et al. Toward a panoramic perspective of the association between environmental factors and cardiovascular disease: An environment-wide association study from National Health and Nutrition Examination Survey 1999–2014
Zhang et al. Selection of microbial biomarkers with genetic algorithm and principal component analysis
Timonen et al. lgpr: an interpretable non-parametric method for inferring covariate effects from longitudinal data
KR102492977B1 (ko) 마이크로바이옴을 이용한 건강 정보 제공 방법 및 분석장치
Kim Bioinformatic and Statistical Analysis of Microbiome Data
KR20190046940A (ko) 개인의 생물학적 상태를 예측하기 위한 시스템, 방법 및 유전자 시그니처
Menger et al. Using cluster ensembles to identify psychiatric patient subgroups
EP1681981A2 (en) Generation of biochemical images and methods of use
US20120265446A1 (en) Biomarkers based on sets of molecular signatures
JP6863633B1 (ja) 菌群抽出方法、装置、及び、プログラム、並びに、腸内細菌叢ｄｂ作成システム、及び、方法。
Queyrel et al. Reject and cascade classifier with subgroup discovery for interpretable metagenomic signatures
Alali et al. A data mining technique for detecting malignant mesothelioma cancer using multiple regression analysis
Pérez et al. Rationale and design of the prognostic transcriptomic signature in fibrotic hypersensitivity pneumonitis (PREDICT) study
Ji et al. Optimal distance metrics for single-cell RNA-seq populations
Gwon et al. Bayesian Modeling on Microbiome Data Analysis: Application to Subgingival Microbiome Study

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19852082

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2019852082

Country of ref document: EP

Effective date: 20210301