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

Abstract

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

Description

本発明の実施形態は、腸内細菌叢を用いた疾病評価指標を算出する技術に関し、特に、採便を提出、かつ、属性などのアンケートを回答した被験者に対し、特定の疾病に対するリスクを示す評価レポートを提供する技術に関する。

腸内細菌叢は数百〜千種を超えるともいわれる膨大な数の細菌から構成される極めて複雑な生態系であり、腸内細菌叢を培養、観察、同定、解析する作業には莫大な労力と時間が必要である。近年、分子生物学的な手法により、これらの問題を解決することが可能になった。この手法には、細菌の遺伝子、特に１６ＳｒＲＮＡの遺伝子（１６ＳｒＤＮＡ）の配列の違いを検出する技術が用いられている。

このような技術を用いることによって、腸内細菌叢と特定の疾病との関係についての研究が進んでいる。非特許文献１では、多数の被験者に採便キットで採取した大便を提出してもらい、さらに、年齢やＢＭＩ（Ｂｏｄｙ Ｍａｓｓ Ｉｎｄｅｘ）の情報も提供してもらい、腸内細菌叢の情報を用いて疾病の罹患を評価するために、多重ロジスティック回帰分析によって特定の菌群と年齢やＢＭＩとの関係を評価する研究が行われている。
また、非特許文献２では、腸内細菌叢の情報を用いた構造方程式のモデリングが記載されている。この構造方程式モデリングは疾病ではなく、菌叢の解析手法（Ｎｅｘｔ Ｇｅｎｅｒａｔｉｏｎ Ｓｅｑｕｅｎｃｉｎｇなど）を評価することを目的に研究されている。

Ｌａｈｔｉ ｅｔ ａｌ．、Ｔｉｐｐｉｎｇ ｅｌｅｍｅｎｔｓ ｉｎ ｔｈｅ ｈｕｍａｎ ｉｎｔｅｓｔｉｎａｌ ｅｃｏｓｙｓｔｅｍ、ｎａｔｕｒｅ ＮＡＴＵＲＥ ＣＯＭＭＵＮＩＣＡＴＩＯＮＳ ５、Ａｒｔｉｃｌｅ ｎｕｍｂｅｒ：４３４４（２０１４）、Ｐｕｂｌｉｓｈｅｄ： ８ Ｊｕｌｙ ２０１４、ｈｔｔｐｓ：／／ｗｗｗ．ｎａｔｕｒｅ．ｃｏｍ／ａｒｔｉｃｌｅｓ／ｎｃｏｍｍｓ５３４４ Ｗｕ ｅｔ ａｌ．、Ｃｏｍｐａｒａｔｉｖｅ ａｎａｌｙｓｉｓ ｏｆ ｍｉｃｒｏｂｉｏｍｅ ｍｅａｓｕｒｅｍｅｎｔ ｐｌａｔｆｏｒｍｓ ｕｓｉｎｇ ｌａｔｅｎｔ ｖａｒｉａｂｌｅ ｓｔｒｕｃｔｕｒａｌ ｅｑｕａｔｉｏｎ ｍｏｄｅｌｉｎｇ、ＢＭＣ Ｂｉｏｉｎｆｏｒｍａｔｉｃｓ ２０１３ １４：７９、Ｐｕｂｌｉｓｈｅｄ ５ Ｍａｒｃｈ ２０１３、ｈｔｔｐ：／／ｗｗｗ．ｂｉｏｍｅｄｃｅｎｔｒａｌ．ｃｏｍ／１４７１−２１０５／１４／７９

このように、腸内細菌叢の情報を用いて、疾病のリスクを評価するためのモデルについては実用化に向けての研究がまだされていない。このため、せっかく採便キットで採取した大便を提出したにもかかわらず、特定の疾病に対するリスクの評価を希望する者（以下、ユーザという）がいても、ユーザ自身の腸内細菌叢の情報から、ユーザが何らかの疾病になるリスクを評価するというサービスを提供することができていない。

本発明は、このような課題に着目して鋭意研究され完成されたものであり、その目的は、腸内細菌叢を用いた疾病評価指標を算出することに関する技術を提供することにある。

上記課題を解決するために、本発明は、疾病評価指標を算出するためのモデル作成方法であって、１又は複数の疾病の入力を受け付け、複数の被験者の便を解析した結果である腸内細菌叢に関するデータを格納する第１のデータベースと、前記被験者の属性及び疾病に関するデータを格納する第２のデータベースから、健康な人及び前記疾病に罹患している人のデータを抽出し、前記属性、前記腸内細菌叢、及び、前記疾病を内部因子とする連関モデルを作成する。

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

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

他の本発明は、ユーザが調べて欲しいと望む１又は複数の疾病を入力する入力部と、前記ユーザの便を予め解析した結果である前記ユーザの腸内細菌叢を第１のデータベースから、及び、前記ユーザの属性を第２のデータベースから抽出する第１の抽出部と、健康な人及び前記疾病に罹患している人を含む複数の被験者の便を解析した結果である前記被験者の腸内細菌叢を前記第１のデータベースから、及び、前記被験者の属性及び疾病を前記第２のデータベースから抽出する第２の抽出部と、前記被験者の属性、前記被験者の腸内細菌叢、及び、前記疾病を内部因子とする連関モデルを作成する作成部と、前記ユーザの腸内細菌叢、前記ユーザの属性、及び、前記連関モデルを用いて、前記疾病に対するリスクを算出する算出部と、を備える疾病評価指標算出システムである。

他の本発明は、ユーザが調べて欲しいと望む１又は複数の疾病を入力する入力ステップと、前記ユーザの便を予め解析した結果である前記ユーザの腸内細菌叢を第１のデータベースから、及び、前記ユーザの属性を第２のデータベースから抽出する第１の抽出ステップと、健康な人及び前記疾病に罹患している人を含む複数の被験者の便を解析した結果である前記被験者の腸内細菌叢を前記第１のデータベースから、及び、前記被験者の属性及び疾病を前記第２のデータベースから抽出する第２の抽出ステップと、前記被験者の属性、前記被験者の腸内細菌叢、及び、前記疾病を内部因子とする連関モデルを作成する作成ステップと、前記ユーザの腸内細菌叢、前記ユーザの属性、及び、前記連関モデルを用いて、前記疾病に対するリスクを算出する算出ステップと、を備える疾病評価指標算出プログラムである。

本発明によれば、腸内細菌叢を用いた疾病評価指標を算出することに関する技術を提供することができる。

本発明の実施形態に係る疾病評価指標算出システムの全体概略図である。 本実施形態に係る腸内細菌叢解析装置２００及び腸内細菌叢ＤＢ３００の概略構成図である。 本実施形態に係る評価モデルの概念図である。 本実施形態に係るアンケート内容の説明図である。 本実施形態に係る各集団の説明図である。 実施例１に係る評価モデルの概略図である。 実施例１に係る評価モデル作成のフローチャートである。 実施例１に係る共分散構造分析の場合のモデルの概略説明図である。 実施例１に係る評価モデル作成フローチャートの一部の説明図である。 実施例１に係る評価モデルの具体例を示す図である。 実施例２に係る評価モデルの具体例を示す図である。

図面を参照しながら本発明の実施の形態を説明する。なお、各図において共通する部分には同一の符号を付し、重複した説明は省略する。

（疾病評価指標算出システム）
図１は、本発明の実施形態に係る疾病評価指標算出システムの全体概略図である。疾病評価指標算出システムは、多数の被験者及び評価希望者（ユーザ）が提出する採便キットから腸内細菌叢データベース（ＤＢ）を作成するフェーズ、ユーザが評価を希望する１又は複数の疾病（特定の疾病）の評価モデルを作成するフェーズ、及び、特定の疾病に対するリスクの評価指標を算出するフェーズを備える。

（腸内細菌叢ＤＢ作成フェーズ）
腸内細菌叢ＤＢ作成フェーズでは、数千人以上という多数の被験者が各自の大便をトイレで排泄し、それを採便キットで採取する。次に、採便キットを受け取った抽出業者は、採便キットを腸内細菌ＤＮＡ抽出装置１００に入力し、腸内細菌叢に関するＤＮＡ溶液が出力される。

ＤＮＡ溶液を受け取った解析業者は、ＤＮＡ溶液を腸内細菌叢解析装置２００に入力し、腸内細菌叢を解析する。解析業者は解析結果データを腸内細菌叢ＤＢ３００に格納する。なお、抽出業者と解析業者を分けて説明したが、同一の事業者が抽出作業と解析作業を行ってもよい。

図２は、腸内細菌叢解析装置２００及び腸内細菌叢ＤＢ３００の概略構成図である。本実施形態では、腸内細菌叢解析装置２００は、ＤＮＡ溶液から腸内細菌叢を解析する方法として２種類用意するために、ＮＧＳ解析装置２１０及びＴ−ＲＦＬＰ解析装置２２０を備えている。ただし、この２種類のうち少なくとも一方を備え、解析してもよい。

ＮＧＳ解析装置２１０では、予め試薬調整分注操作が行われたＤＮＡ溶液を入力し、次世代シーケンサ（Ｎｅｘｔ Ｇｅｎｅｒａｔｉｏｎ Ｓｅｑｕｅｎｃｅｒ）で解析することによって、腸内細菌叢に関するＤＮＡ配列が出力され、腸内細菌叢（ＤＮＡ配列）ＤＢ３１０に格納される。つまり、ＮＧＳでは、配列による細菌分類を行い、菌分類毎の細菌の相対量を解析することができる。

Ｔ−ＲＦＬＰ解析装置２２０では、試薬分注と調整を施したＤＮＡ溶液を入力し、キャビラリーシーケンサがＴ−ＲＦＬＰ（Ｔｅｒｍｉｎａｌ − Ｒｅｓｔｒｉｃｔｉｏｎ Ｆｒａｇｍｅｎｔ Ｌｅｎｇｔｈ Ｐｏｌｙｍｏｒｐｈｉｓｍ）で解析することによって、腸内細菌叢に関するＤＮＡフィンガープリントが出力され、腸内細菌叢（ＤＮＡフィンガープリント）ＤＢ３２０に格納される。つまり、Ｔ−ＲＦＬＰでは、ＤＮＡ断片長による細菌分類（ＯＴＵ）を行い、ＯＴＵ毎の細菌の相対量を解析することができる。ＯＴＵは操作的分類単位（Ｏｐｅｒａｔｉｏｎａｌ Ｔａｘｏｎｏｍｉｃ Ｕｎｉｔ）といい、「Ａ１９３」などの群に分類されるが、実際の菌属の名前はわからない。なお、「Ａ１９３」の「Ａ」は制限酵素

の略称である。この制限酵素によってＤＮＡを切断している。

腸内細菌叢ＤＢ３００は、腸内細菌叢（ＤＮＡ配列）ＤＢ３１０及び腸内細菌叢（ＤＮＡフィンガープリント）ＤＢ３２０を備えている。この２種類のＤＢ（３１０及び３２０）は上述した２種類の解析装置（２１０及び２２０）の解析結果をそれぞれ格納する。なお、２種類の解析装置の一方（２１０又は２２０）が腸内細菌叢解析装置２００に含まれている場合、対応するＤＢ（３１０又は３２０）が腸内細菌叢ＤＢ３００に含まれていればよい。このようなデータベース構造により、数千人以上という多数の被験者の腸内細菌叢が異なる方法で解析された場合であっても一つのデータベースとして取り扱うことができる。そして、腸内細菌叢ＤＢ３００は各被験者のＩＤ（識別番号）情報とその被験者の腸内細菌叢とを関係付けている。

次に、採便キットを提出した多数の被験者は、アンケートを記入し、アンケート回収業者に提出する。アンケート回収業者はアンケート結果をアンケートＤＢ４００に格納する。ここで、アンケート回収の初期段階では、アンケートＤＢ４００に格納されているアンケート結果が少ないため、外部のアンケートＤＢをアンケートＤＢ４００として利用してもよい。同様に、腸内細菌叢ＤＢ３００に格納されているデータも少ないため、外部の腸内細菌叢ＤＢを腸内細菌叢ＤＢ３００として利用してもよい。

（評価モデル作成フェーズ）
評価モデル作成フェーズでは、まず、評価モデル作成業者は、評価したい属性と疾病を選定する。モデル作成装置５００は、選定された属性及び疾病の入力を受け付ける入力部と、これら属性及び疾病に関係する被験者のＩＤ情報をアンケートＤＢ４００から抽出し、該当するＩＤ情報の腸内細菌叢も腸内細菌叢ＤＢ３００から抽出する抽出部と、後述する作成方法によって、評価モデルを作成する作成部を備える。モデル作成装置５００は、作成した評価モデルを評価モデルＤＢ６００に格納する。

（疾病評価指標算出フェーズ）
疾病評価指標算出フェーズでは、評価業者は、採便キットで採取した大便を提出し、かつ、特定の疾病に対するリスクの評価を希望する者（以下、ユーザという）に対し、腸内細菌叢ＤＢ作成フェーズで説明した採便キットの提出、及び、評価モデル作成フェーズで説明したアンケートの記入を依頼する。ユーザの採便キット及びアンケートは、多数の被験者の採便キット及びアンケートと同様の処理が行われ、腸内細菌叢ＤＢ３００及びアンケートＤＢ４００に格納される。さらに、ユーザが評価を希望する疾病（特定の疾病）、及び、ユーザのＩＤ情報を疾病評価指標算出装置７００に入力する。

疾病評価指標算出装置７００は、入力情報を受付ける入力部と、入力情報に関係する情報を抽出する抽出部と、評価指標を算出する算出部と、評価レポートを出力する出力部を備える。入力部は特定の疾病及びユーザのＩＤ情報の入力を受け付ける。抽出部は外部のアンケートＤＢ４００及び腸内細菌叢ＤＢ３００からユーザのＩＤ情報に関係する情報（ユーザの要因、及び、腸内細菌叢）を抽出し、評価モデルＤＢ６００から特定の疾病に関する評価モデルを抽出する。そして、疾病評価指標算出装置７００は、抽出されたユーザの要因（属性や生活習慣など）及び腸内細菌叢を、特定の疾病の評価モデルに入力し、特定の疾病に対するリスクの評価指標を算出する。算出された評価指標に基づく、評価レポートが出力される。なお、各フェーズを行う業者を分けて説明したが、同一の事業者が全てのフェーズを行ってもよい。

図３は、評価モデルＤＢ６００に格納される評価モデルの一例を説明するための概念図である。評価モデルは、特定の疾病と一定以上の相関が認められた観測変数と、観測変数と関係する潜在変数（以下、因子という）によって構築される統計解析モデルである。なお、この評価モデルでは、要因（主な原因）も因子として取り扱っている。

観測変数は実線の長方形で表し、因子は楕円形で表している。各観測変数は、図１の腸内細菌叢ＤＢ３００又はアンケートＤＢ４００の情報の一部である。

評価モデルの主要部は、観測変数としての属性（性別・年齢）６１０−１、並びに、因子としての要因１（腸内細菌叢）６２０−１、要因２（生活・環境等）６２０−２及び因子（疾病等）６２０−３を備える。属性（性別・年齢）６１０−１はアンケートＤＢ４００の情報の一部である。矢印の向きは「影響を与える向き」を定義している。また、要因１（腸内細菌叢）６２０−１から因子（疾病等）６２０−３への太い矢印は一方向になっているが、双方向の評価モデルを構築することも可能である。

次に、因子を詳細に説明する。点線で囲まれた長方形６３０−１は要因１（腸内細菌叢）６２０−１の詳細を示す。要因１（腸内細菌叢）６２０−１は腸内細菌叢因子６２０−１’及び複数の「ＯＴＵ又は菌属（菌種）」６１０−２に分解できる。複数の「ＯＴＵ又は菌属（菌種）」６１０−２は観測変数であり、図１の腸内細菌叢ＤＢ３００の情報の一部である。

点線で囲まれた長方形６３０−２は要因２（生活・環境等）６２０−２の詳細を示す。要因２（生活・環境等）６２０−２は、環境要因６２０−２’−１、生活環境（アンケート）６１０−３及び家庭環境（アンケート）６１０−４、並びに、食生活要因６２０−２’−２、食行動（アンケート）６１０−５及び食意識（アンケート）６１０−６に分解できる。生活環境（アンケート）６１０−３及び家庭環境（アンケート）６１０−４、並びに、食行動（アンケート）６１０−５及び食意識（アンケート）６１０−６は観測変数であり、図１のアンケートＤＢ４００の情報の一部である。

点線で囲まれた長方形６３０−３は因子（疾病等）６２０−３の詳細を示す。因子（疾病等）６２０−３は、疾病・健康状態因子６２０−３’、疾病（アンケート）６１０−７及び健康状態（アンケート）６１０−８に分解できる。疾病（アンケート）６１０−７及び健康状態（アンケート）６１０−８は観測変数であり、図１のアンケートＤＢ４００の情報の一部である。

この図では、ユーザを属性（性別・年齢）６１０−１、並びに、各アンケート６１０−３乃至６１０−８で分割された集団に割り当て、その集団が保有する要因１（腸内細菌叢）６２０−１の場合、ユーザが評価を希望する疾病（特定の疾病）に対するリスクは高いか低いかを算出することができることを概念的に表している。ここで、特定の疾病に対するリスクとは、特定の疾病に対して現在罹患している、又は、将来罹患する確率（可能性も含む）をいう。

すなわち、特定の疾病と一定以上の相関が認められた観測変数として「ＯＴＵ又は菌属（菌種）」６１０−２、並びに、属性（性別・年齢）６１０−１、生活環境（アンケート）６１０−３、家庭環境（アンケート）６１０−４、食行動（アンケート）６１０−５、食意識（アンケート）６１０−６、疾病（アンケート）６１０−７、及び、健康状態（アンケート）６１０−８を選択し、これらの観測変数と関係する因子によって評価モデルが構築される。

ここで、属性６１０−１については、統計解析を行う母集団を男女別に分割し、性別及び年齢を観測変数として用いている。また、腸内細菌叢因子６２０−１’の観測変数として、ＯＴＵ又は菌属（菌種）名６１０−２を用いている。これは、Ｔ−ＲＦＬＰによりＤＮＡフィンガープリントがわかる場合、ＯＴＵを用い、ＮＧＳによりＤＮＡ配列が分かる場合、菌属（菌種）名を用いるためである。

図４は、本実施形態に係るアンケート内容の説明図である。アンケートの構成概念は、「属性」、「生活・家庭環境（ライフスタイル）」、「食意識」、「食行動」、「心理的負荷（ストレス）」、「投薬」、「病気・未病（健康状態）」であり、これらに該当する各質問項目は一つ又は複数である。また、「７−３身体状況」には、身体状況から判断される健康状態、例えば、アトピー性皮膚炎に関するものとして、かゆみや皮膚の赤みという項目、又は、自分で測定できる血圧などを含まれる。ここで、説明の簡略化のため、「病気・未病（健康状態）」は「疾病」と述べてもよく、「生活・家庭環境（ライフスタイル）」、「食意識」、「食行動」、「心理的負荷（ストレス）」及び「投薬」はまとめて「生活及び環境」と述べてもよい。

図５は、本実施形態に係る各集団の説明図である。全被験者で構成される母集団、疾病にかかっていない健康母集団、アトピーにかかっている集団、アトピー及び喘息にかかっている集団、喘息にかかっている集団で構成されている。

（実施例１：属性と腸内細菌叢から疾病リスクを評価する場合）
実施例１では、図３で説明した評価モデルの主要部のうち、観測変数としての属性（性別・年齢）、並びに、因子としての要因１（腸内細菌叢）及び因子（疾病等）を用いて、疾病リスクを評価する場合について説明する。

図６は、実施例１に係る評価モデルの概略図である。実施例１の評価モデルは、疾病「アトピー」又は「喘息」に対するリスクの評価指標を算出するためのモデルである。観測変数「性別」及び「年齢」は被験者が提出したアンケートから抽出し、観測変数「Ａ５８」、「Ａ１３１」などは被験者が提出した採便キットから抽出している。また、観測変数「アトピー」及び「喘息」が、ユーザが評価を希望する疾病である。この評価モデルを用いることによって、ユーザが「アトピー」及び「喘息」になるリスクの評価指標を算出することができる。

実施例１の評価モデルの内部因子は「属性」、「腸内細菌叢」及び「疾病・健康状態」である。これらの内部因子同士の関係が、連関モデルとして評価モデルの概略構成で表現されている。本願では、こうした観測変数、内部因子同士の多様な関係性を統計学的に評価する分析方法として共分散構造分析を採用する。図６に示す観測変数と内部因子との関係、及び、内部因子同士の関係は構造方程式モデル（ＳＥＭ；Ｓｔｒｕｃｔｕｒａｌ Ｅｑｕａｔｉｏｎ Ｍｏｄｅｌｉｎｇ）で表現され、評価モデルが生成される。

また、内部因子と観測変数との関係は次のようになる。内部因子「属性」の観測変数が「性別」及び「年齢」であり、内部因子「腸内細菌叢」の観測変数が「Ａ５８」、「Ａ１３１」などのＯＴＵであり、さらに、内部因子「疾病・健康状態」の観測変数が「アトピー」及び「喘息」である。なお、説明の簡略化のため、内部因子「疾病・健康状態」は内部因子「疾病」と述べてもよい。

図７は、実施例１に係る評価モデル作成のフローチャートである。ここでは、疾病の一例として「アトピー」及び「喘息」について説明するが、他の疾病、例えば、「心臓病」、「大腸がん」、「胃腸病」、「肝臓病」、「腎臓病」、「糖尿病」、「高血圧」、「脂質異常」、「低体重（ＢＭＩ）」、「肥満（ＢＭＩ）」、「骨病」、「関節痛」、「うつ」などについても同様に評価モデルを作成することができる

まず、疾病として「アトピー」及び「喘息」を選択する（Ｓ１００）。評価に使用する属性として「性別」及び「年齢」を使用することを選択する（Ｓ１１０）。

次に、疾病「アトピー」と関係の強いＯＴＵを選定する（Ｓ１２０）。この選定は、例えばロジスティック回帰分析に基づいて行われてもよい。ここで、所定の疾病と関係の強いＯＴＵをロジスティック回帰分析により選定する際に、所定の属性の時に、所定の腸内細菌叢と所定の疾病・健康状態とのパス、そのパスの係数、オッズ比、及び、ｐ値（有意確率）で構成される疾病リスク係数表を用いる。

次いで、選定されたＯＴＵを観測変数として、評価モデルの内部因子「腸内細菌叢」を推定する（Ｓ１３０）。この推定は、例えば因子分析に基づいて行われてもよく、あるいは主成分分析、主座標分析等、他の分析方法に基づいてもよい。また、疾病「喘息」についても同様の処理（Ｓ１２０及びＳ１３０）を行う。

Ｓ１２０と平行して（時間の前後は問わない）、疾病「アトピー」に関係する、評価モデルの内部因子としての属性を疾病リスク係数表より選定する（Ｓ１４０）。そして、選定された属性に対し、評価モデルの内部因子「属性」を推定する（Ｓ１５０）。この推定は、例えば因子分析に基づいてもよく、あるいは主成分分析、主座標分析等、他の分析方法に基づいてもよい。また、疾病「喘息」についても同様の処理（Ｓ１４０及びＳ１５０）を行う。

Ｓ１３０で推定された評価モデルの内部因子「腸内細菌叢」及び、Ｓ１４０で推定された内部因子「属性」、さらに内部因子「疾病・健康状態」によりＳＥＭ解析を実施する（Ｓ１６０）。

そして、図８に示す構造方程式のパラメータα、β、γを求め、評価モデルを作成する（Ｓ１７０）。

図９は、実施例１に係る評価モデル作成フローチャートの一部（Ｓ１２０）の説明図である。特定の疾病に関して、疾病リスク係数表からＯＴＵを選定するステップについて説明する。

まず、（１）疾病を決める。ここでは「アトピー」及び「喘息」とする。（２）係数表からＯＴＵを選定する。ここでは「Ａ５８」、「Ａ１３１」、「Ａ４７７」及び「Ａ４９８」等を記載しているが、実際には各ＯＴＵの次元は４０〜１００次元程度になる。

次に、（３）図５で説明した各集団から母集団を選定する。（４）所定の属性として「年齢」を選択する。ここでは、図５の各集団は男性のみ、女性のみ、及び、全体に予め分かれており、母集団毎の評価モデルを構築するため、「年齢」を選択する。（５）ＳＥＭ解析を実施し、（６）各パスの計数と、モデルの適合度を確認、すなわち、所定のしきい値以上であれば終了する。

図１０は、実施例１に係るアトピー及び喘息の場合の評価モデルの具体例である。各パスの矢印の向きは影響を与える向きを表し、各パスの数値が１に近いほど正の相関度が高い、すなわち、影響を与える要素として高いことを示している。一方、中黒の三角印はマイナスの数値であり、負の相関であることを示している。

図１０の評価モデルの具体例によれば、疾病「アトピー」及び「喘息」については、ある被験者は「年齢」が高く、「Ａ１３１」に属する菌を多く保有している場合、「年齢」と「疾病因子」とは正の相関があるため、「年齢」が高いほど「アトピー」及び「喘息」にかかりやすい、と評価できる。しかしながら、「Ａ１３１」と「腸内細菌叢Ａ」は正の相関があるため、「腸内細菌叢Ａ」が高いことを意味する。そして、「アトピー」及び「喘息」と「疾病因子」は正の相関があり、かつ、「腸内細菌叢Ａ」と「疾病因子」は負の相関があるため、「腸内細菌叢Ａ」が高ければ、「アトピー」及び「喘息」にはかかりにくい、と評価できる。

さらに、「年齢」から「疾病因子」へのパスは、直接的なパスと、「腸内細菌叢Ａ」を経由する間接的なパスがあることに留意する。直接的なパスによれば、上述した通り「年齢」が高いほど「アトピー」及び「喘息」にかかりやすい。一方、間接的なパスによれば、年齢が高くなると、「腸内細菌叢Ａ」も少し高くなり、「アトピー」及び「喘息」にはかかりにくい。このように間接的なパスも考慮すると、「年齢」が高いほど「アトピー」及び「喘息」にかかりにくいという影響は実は小さいと評価できる。

評価モデルの適合度指標として３種類用いており、ＧＦＩはＧｏｏｄｎｅｓｓ ｏｆ Ｆｉｔ Ｉｎｄｅｘ、ＡＧＦＩはＡｄｊｕｓｔｅｄ ＧＦＩ、ＲＭＳＥＡはＲｏｏｔ Ｍｅａｎ Ｓｑｕａｒｅ Ｅｒｒｏｒ ｏｆ Ａｐｐｒｏｘｉｍａｔｉｏｎの略称である。ＧＦＩ及びＡＧＦＩは値が１に近くなるほどデータへの当てはまり（以下、適合度という）が良いことを意味する。一方、ＲＭＳＥＡは値が小さいほど適合度が良いことを意味する。また、評価モデルには、誤差のバラツキの指標である誤差分散を用いてもよい。

図１０の評価モデルの適合度指標は、ＧＦＩが０．９７８８、ＡＧＦＩが０．９７３５、ＲＭＳＥＡが０．１０９６であり、この評価モデルの適合度は高いことが示されている。

（効果）
実施例１によれば、属性と腸内細菌叢から疾病リスクを評価する評価モデルを作成することができる。すなわち、図３の要因１（腸内細菌叢）６２０−１から因子（疾病等）６２０−３への太い矢印の関連性を評価するモデルを作成することができる。また、図６で説明したように観測変数、内部因子同士の多様な関係性をモデル化しているため、図１０で説明したように直接的なパス及び間接的なパスの両方を判断可能な評価モデルを作成することができる。

そして、ユーザは採便キットで大便を採取、提出し、さらにアンケートに回答することによって、自分の大便が「Ａ１３１」を多く保有していれば、年齢による影響は小さく、「腸内細菌叢Ａ」が高く、「アトピー」及び「喘息」にかかりにくい、という評価を受けることができる。

そして、疾病評価指標算出装置７００は、「アトピー」及び「喘息」になるリスクの評価を希望するユーザから、ユーザＩＤ情報を受け付け、アンケートＤＢ４００及び腸内細菌叢ＤＢ３００からユーザのＩＤ情報に関係する情報（ユーザの属性、及び、腸内細菌叢の情報）を抽出し、評価モデルＤＢ６００からユーザが評価を希望する疾病「アトピー」及び「喘息」に関する評価モデルを抽出する。そして、疾病評価指標算出装置７００は、抽出されたユーザの属性（この場合、年齢）及び腸内細菌叢の情報（この場合、「Ａ１３１」を多く保有）を、図１０の評価モデルに入力し、「アトピー」及び「喘息」に対するリスクの評価指標を算出することができる。さらに、算出された評価指標に基づく、評価レポートが出力される。

ユーザは評価レポートを読むことによって、ユーザ自身の腸内細菌叢に、腸内細菌叢Ａ（因子）への因子負荷量が大きい「Ａ１３１」が多く含まれているため、年齢による影響は小さく（年齢が高くなっても、将来的に）「アトピー」及び「喘息」に対するリスクは低いということを理解できる。

次に、疾病評価指標算出システム事業者にとっての効果を述べる。ユーザの腸内細菌叢に「Ａ１３１」が多く含まれていることが「アトピー」及び「喘息」に対するリスクを低減させるのであれば、「Ａ１３１」に含まれる実在の菌を摂取すれば、「アトピー」及び「喘息」に対するリスクが低減される可能性を評価レポートは示唆している。この可能性に基づき、事業者はプロバイオティクスの設計、検討、及び提案を行うことができる。

（実施例２；属性並びに、生活及び環境と、腸内細菌叢から疾病リスクを評価する場合）
実施例２では、図３で説明した要因２（生活・環境等）を実施例１に加え、すなわち、図３で説明した評価モデルの主要部全てを用いて、疾病リスクを評価する場合について説明する。実施例２の疾病評価指標算出システム７００も、図１の全体概略図と同じである。評価モデルを作成する際に、要因２（生活・環境等）を用いている点が、実施例１と異なる点である。

図１１は、実施例２に係る評価モデルの具体例を示す図である。実施例２特有の評価モデルの内部因子は、図３の要因２（生活・環境等）として、「食生活」を用いている。内部因子「食生活」の観測変数が「牛乳」、「パン」、「肉」等である。

観測変数「牛乳」、「パン」、「肉」等を評価モデルに加えることは、図４で説明したアンケートの構成概念「食行動」の各質問項目をアンケート情報に追加することによって可能になる。また、各質問項目の数を増やすことによって、より詳細な評価をすることが可能になる。

図１１の評価モデルの具体例によれば、実施例１の解釈に加え、ある被験者が「肉」を食する頻度が高い場合、「食生活」因子が高いことを意味する。これは、「腸内細菌叢Ａ」因子が低下し、「腸内細菌叢Ｂ」因子が増加し、さらに、「疾病因子」の増加を意味する。この「疾病因子」は「アトピー」及び「喘息」になるリスクを増加させていることを示唆している。また、この「食生活」因子が高いことにより、「腸内細菌叢Ｂ」因子が高くなり、その高い「腸内細菌叢Ｂ」因子が「疾病因子」を高めている可能性をさらに示唆している、と評価できる。

図１１の評価モデルの適合度指標は、ＧＦＩが０．９２６５であり、ＡＧＦＩが０．９０４８であり、ＲＭＳＥＡが０．０７３８である。このため、この評価モデルの適合度は高いことが示されている。

（効果）
実施例２によれば、属性及び生活・環境等と、腸内細菌叢から疾病リスクを評価する評価モデルを作成することができる。図１１の場合、「年齢」、「腸内細菌叢Ａ」、「腸内細菌叢Ｂ」、「食生活」から「アトピー」及び「喘息」になるリスクを評価する評価モデルを作成することができる。

そして、疾病評価指標算出装置７００は、「アトピー」及び「喘息」になるリスクの評価を希望するユーザから、ユーザＩＤ情報を受け付け、アンケートＤＢ４００及び腸内細菌叢ＤＢ３００からユーザのＩＤ情報に関係する情報（ユーザの属性及び「食生活」に関するアンケート結果、並びに、腸内細菌叢の情報）を抽出し、評価モデルＤＢ６００からユーザが評価を希望する疾病「アトピー」及び「喘息」に関する評価モデルを抽出する。そして、疾病評価指標算出装置７００は、抽出されたユーザの要因「年齢、食生活」及び腸内細菌叢の情報を、「アトピー」及び「喘息」の評価モデルに入力し、「アトピー」及び「喘息」に対するリスクの評価指標を算出することができる。さらに、算出された評価指標に基づく、評価レポートが出力される。

ユーザは評価レポートを読むことによって、実施例１の評価に加え、「肉」を食する頻度が高いため、「腸内細菌叢Ａ」因子が低下し、「腸内細菌叢Ｂ」因子が増加する傾向にあり、「アトピー」及び「喘息」になるリスクを増加させていることを示唆している。また、「肉」を食する頻度が高いことにより、「腸内細菌叢Ｂ」因子が高くなり、「アトピー」及び「喘息」になるリスクを高めている可能性をさらに示唆している、という詳細な評価を理解できる。

次に、疾病評価指標算出システム事業者にとっての効果を述べる。ユーザの「肉」を食する頻度が高く、「野菜」を食する頻度が低いことが「アトピー」及び「喘息」に対するリスクを増加させるのであれは、「肉」を食する頻度を抑え、また「野菜」を多く摂取すれば、「アトピー」及び「喘息」に対するリスクが低減される可能性を評価レポートは示唆している。このような可能性に基づき、事業者はプレバイオティクスの設計、検討、及び提案を行うことができる。

以上、本発明の実施例（変形例を含む）について説明してきたが、これらのうち、２つ以上の実施例を組み合わせて実施しても構わない。あるいは、これらのうち、１つの実施例を部分的に実施しても構わない。さらには、これらのうち、２つ以上の実施例を部分的に組み合わせて実施しても構わない。

また、本発明は、上記発明の実施例の説明に何ら限定されるものではない。特許請求の範囲の記載を逸脱せず、当業者が容易に想到できる範囲で種々の変形態様もこの発明に含まれる。

例えば、実施例の腸内細菌叢解析装置２００としてｑＰＣＲなどの他の解析装置を用いてもよい。ｑＰＣＲ解析装置の場合、試薬分注と調整を施したＤＮＡ溶液を入力し、リアルタイムの定量ＰＣＲ（ｑｕａｎｔｉｔａｔｉｖｅ Ｐｏｌｙｍｅｒａｓｅ Ｃｈａｉｎ Ｒｅａｃｔｉｏｎ）で解析することによって、腸内細菌叢に関する指標菌群が出力され、腸内細菌叢ＤＢ３００内の腸内細菌叢（指標菌群）ＤＢに格納される。

また、実施例では、連関モデルを解析する方法として、共分散構造分析を用いたが、パス解析などの他の解析方法を用いてもよい。

１００ 腸内細菌ＤＮＡ抽出装置
２００ 腸内細菌叢解析装置
３００ 腸内細菌叢ＤＢ
４００ アンケートＤＢ
５００ モデル作成装置
６００ 評価モデルＤＢ
７００ 疾病評価指標算出装置

Claims (13)

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

Images