(57)【要約】【課題】腸内菌叢について老齢タイプと若齢タイプのいずれであるかを判定する方法の提供。【解決手段】腸内菌叢について老齢タイプと若齢タイプのいずれであるかを判定する方法であって、糞便検体に含まれる全て菌の数に対するLachnospiraceae科に属する菌群の割合に基づき前記腸内菌叢を判定する判定工程、を有する、腸内菌叢の老若判定方法を提供する。【選択図】なし【特許請求の範囲】
【請求項１】
腸内菌叢について老齢タイプと若齢タイプのいずれであるかを判定する方法であって、
糞便検体に含まれる全て菌の数に対する（ａ）Lachnospiraceae科に属する菌群の割合に基づき前記腸内菌叢を判定する判定工程、を有する、腸内菌叢の老若判定方法。
【請求項２】
前記判定工程では、前記（ａ）Lachnospiraceae科に属する菌群の前記割合に加えて、下記（ｂ）〜（ｊ）の科に属する菌群のうち少なくともいずれか一つの菌群の、前記糞便検体に含まれる全ての菌の数に対する割合に基づき、前記腸内菌叢を判定する、
請求項１に記載の腸内菌叢の老若判定方法。
（ｂ）Bifidobacteriaceae科
（ｃ）Peptococcaceae科
（ｄ）Enterococcaceae科
（ｅ）Verrucomicrobiaceae科
（ｆ）Catabacteriaceae科
（ｇ）Clostridiaceae科
（ｈ）Caulobacteraceae科
（ｉ）Propionibacteriaceae科
（ｊ）Enterobacteriaceae科
【請求項３】
前記判定工程では、前記菌群ごとの前記割合を独立変数とし、前記糞便検体を得た個体の年齢又は月齢に関する回帰分析によって回帰式を算出して、当該回帰式に基づき前記腸内菌叢を判定する、
請求項１又は２に記載の腸内菌叢の老若判定方法。
【請求項４】
前記判定工程では、前記回帰式の解と、予め定められた基準値と、を比較して、前記解が前記基準値を超える場合には前記腸内菌叢を老齢タイプと判定し、前記解が前記基準値以下である場合には前記腸内菌叢を若齢タイプと判定する
請求項３に記載の腸内菌叢の老若判定方法。
【請求項５】
前記判定工程では、前記回帰式として下記式（１）を算出し、
下記式（１）の解ｙが０を超える場合には、前記腸内菌叢を老齢タイプと判定する、
請求項４に記載の腸内菌叢の老若判定方法。
2015188380.tif 000016
（式中、ｎは前記独立変数の数であり、ｋは前記菌群が属する科のうち第ｋ番目の科を意味し、αは前記回帰式の切片であり、βは前記回帰式の回帰係数であり、ｘは前記割合である。）
【請求項６】
前記判定工程の前に前記割合を測定する測定工程を有する、請求項１〜５のいずれか一項に記載の腸内菌叢の老若判定方法。
【請求項７】
請求項１〜６のいずれか一項に記載の腸内菌叢の老若判定方法の結果に基づいて糞便検体を複数の群に分ける選別工程と、
選別された前記複数の群のうち少なくともいずれか一つの群の糞便検体に含まれる菌群を同定する菌群同定工程と、を有する、腸内菌叢の老若判定マーカーとなる菌群のスクリーニング方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、腸内菌叢の老若判定方法及び老若判定マーカーとなる菌群のスクリーニング方法に関する。
【背景技術】
【０００２】
ヒトなど動物の腸内には、多種類の細菌が常在している。これらの常在細菌によって形成される腸内菌叢においては、乳児期にビフィズス菌が増加したり、老人でウェルシュ菌が増加し、ビフィズス菌が減少するなど、成長や加齢に伴い特定の菌群が増減することが知られている。また、宿主の健康状態や宿主が摂取する薬剤、食物等によって腸内菌叢を構成する菌群の種類やその数が変動することも知られている。
【０００３】
このような腸内菌叢については、近年、そのバランスが崩れると腸内環境を悪化させ、便秘症状や下痢症状を発症したり、宿主の免疫機能が低下することなどが明らかとなってきている。これに対し、良好な腸内菌叢の維持や悪化した腸内菌叢を改善する薬剤や飲食物も開発されてきている。例えば、特許文献１には、「ラクトパーオキシダーゼを有効成分とする腸内フローラ改善剤」が開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００３−２４６７５３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上記特許文献１に記載された腸内フローラ改善剤のような腸内菌叢に作用する成分の効能を評価するためには、腸内菌叢の状態が改善しているか否かを判定する必要がある。一方、腸内菌叢を構成する菌群は多様であり、個体差も大きいため、腸内菌叢に含まれる菌群については、未だに同定されていない種もある。このため、腸内菌叢の状態に関してより詳細に判定するために、腸内菌叢の状態と関連する菌群の同定や当該菌群の数の変動を調べることが求められている。特に、老化現象は、多様な症状を含む現象であり、個体差も大きいため、腸内菌叢の状態について「老齢の状態」にあるか又は前記老齢の状態より若い「若齢の状態」にあるかを判定するためには、適切な指標を用いることが求められている。そこで、本開示は、腸内菌叢について老齢タイプと若齢タイプのいずれであるかを判定する方法等を提供することを主な目的とする。
【課題を解決するための手段】
【０００６】
本発明者は鋭意検討した結果、糞便検体から同定された複数種類の菌群のうち、Lachnospiraceae科に属する菌群の、菌叢全体の菌数に対する割合に基づいて、糞便検体を高齢タイプと若齢タイプに、精度高く分類できることを見出し、本発明を完成させるに至った。
【０００７】
すなわち本開示は、腸内菌叢について老齢タイプと若齢タイプのいずれであるかを判定する方法であって、糞便検体に含まれる全て菌の数に対する（ａ）Lachnospiraceae科に属する菌群の割合に基づき前記腸内菌叢を判定する判定工程、を有する、腸内菌叢の老若判定方法を提供する。
前記判定工程では、前記（ａ）Lachnospiraceae科に属する菌群の前記割合に加えて、下記（ｂ）〜（ｊ）の科に属する菌群のうち少なくともいずれか一つの菌群の、前記糞便検体に含まれる全ての菌の数に対する割合に基づき、前記腸内菌叢を判定してもよい。
（ｂ）Bifidobacteriaceae科
（ｃ）Peptococcaceae科
（ｄ）Enterococcaceae科
（ｅ）Verrucomicrobiaceae科
（ｆ）Catabacteriaceae科
（ｇ）Clostridiaceae科
（ｈ）Caulobacteraceae科
（ｉ）Propionibacteriaceae科
（ｊ）Enterobacteriaceae科
また、前記判定工程では、前記菌群ごとの前記割合を独立変数とし、前記糞便検体を得た個体の年齢又は月齢に関する回帰分析によって回帰式を算出して、当該回帰式に基づき前記腸内菌叢を判定することもできる。
さらに、前記判定工程では、前記回帰式の解と、予め定められた基準値と、を比較して、前記解が前記基準値を超える場合には前記腸内菌叢を老齢タイプと判定し、前記解が前記基準値以下である場合には前記腸内菌叢を若齢タイプと判定することもできる。
本開示は、また、前記判定工程では、前記回帰式として下記式（１）を算出し、下記式（１）の解ｙが０を超える場合には、前記腸内菌叢を老齢タイプと判定する判定方法とすることもできる。
【０００８】
【数１】
2015188380.tif 000002
【０００９】
（式中、ｎは前記独立変数の数であり、ｋは前記菌群が属する科のうち第ｋ番目の科を意味し、αは前記回帰式の切片であり、βは前記回帰式の回帰係数であり、ｘは前記割合である。）
前記腸内菌叢の老若判定方法は、前記判定工程の前に前記割合を測定する測定工程を有していてもよい。
【００１０】
本開示はまた、上記腸内菌叢の老若判定方法の結果に基づいて糞便検体を複数の群に分ける選別工程と、選別された前記複数の群のうち少なくともいずれか一つの群の糞便検体に含まれる菌群を同定する菌群同定工程と、を有する、腸内菌叢の老若判定マーカーとなる菌群のスクリーニング方法を提供する。
【発明の効果】
【００１１】
本開示により、腸内菌叢について老齢タイプと若齢タイプのいずれであるかを判定する方法等が提供される。
【図面の簡単な説明】
【００１２】
【図１】本開示に係る腸内菌叢の老若判定方法の工程の一例を示すフローチャートである。
【図２】本開示に係る菌群のスクリーニング方法の工程の一例を示すフローチャートである。
【図３】選別工程の一例を示すフローチャートである。
【図４】本開示に係る物質のスクリーニング方法の工程を示すフローチャートである。
【発明を実施するための形態】
【００１３】
以下、本開示を実施するための好適な形態について説明する。なお、以下に説明する実施形態は、本開示の代表的な実施形態を示したものであり、これにより本開示の範囲が狭く解釈されることはない。
【００１４】
（１）腸内菌叢の老若判定方法
本開示に係る腸内菌叢の老若判定方法では、腸内菌叢について老齢タイプと若齢タイプのいずれであるかを判定する。図１は、本開示に係る腸内菌叢の老若判定方法の一例を示すフローチャートである。腸内菌叢の老若判定方法には、少なくとも判定工程Ｓ１２が含まれる。また、本開示に係る腸内菌叢の老若判定方法では、判定工程Ｓ１２の前に、後述する、糞便検体に含まれる全ての菌の数に対する特定の菌群の割合を測定する測定工程Ｓ１１を有していてもよい。
以下、測定工程Ｓ１１と判定工程Ｓ１２について順に説明する。なお、本開示に係る腸内菌叢の老若判定方法において、腸内細菌叢を構成する菌群の「菌」とは、細菌（ｂａｃｔｅｒｉａ）及び古細菌（ａｒｃｈａｅａ）を含み、かつ、真菌（ｆｕｎｇｉ）を含まないものである。また、これは、後述する菌群のスクリーニング方法においても同様である。
【００１５】
本開示に係る腸内菌叢の老若判定方法では、糞便検体として動物から得られたものを用いることができる。動物としては、例えば、サル、チンパンジー、マウス、ラット、ハムスター、モルモット、ウシ、ブタ、ヒツジ、ヤギ、イヌ、ネコ等が挙げられる。ヒトから得られる糞便検体が好ましく、ヒトについては、２０歳以上の成人から得られる糞便検体がより好ましい。乳幼児など、腸内菌叢に関して成人とは異なる特徴を有することが知られている年齢層を除外することによって、所定の年齢未満の若齢集団における腸内菌叢を構成する菌群の種類や菌数についてのばらつきが、より抑制される。
【００１６】
測定工程Ｓ１１は、糞便検体に含まれる菌群のうち少なくとも（ａ）Lachnospiraceae科に属する菌群について、糞便検体に含まれる全ての菌の数に対する割合を測定する工程である。また、測定対象とする菌群は、（ａ）Lachnospiraceae科に加えて、（ｂ）Bifidobacteriaceae科、（ｃ）Peptococcaceae科、（ｄ）Enterococcaceae科、（ｅ）Verrucomicrobiaceae科、（ｆ）Catabacteriaceae科、（ｇ）Clostridiaceae科、（ｈ）Caulobacteraceae科、（ｉ）Propionibacteriaceae科、（ｊ）Enterobacteriaceae科のうちいずれかの科に属する菌群も測定対象とすることができる。
【００１７】
本工程Ｓ１１では、糞便検体に含まれる全ての菌の数に対する、上記（ａ）〜（ｊ）科に属する菌群のうち少なくとも（ａ）Lachnospiraceae科に属する菌群の割合を求める。全菌数及び上記割合は、例えば、下記の手法により求めることができる。
【００１８】
＜（ａ）〜（ｊ）の科に属する菌群の測定＞
糞便検体に含まれる菌のうち、上記（ａ）〜（ｊ）の科に属する菌群に分類される菌の数は、公知手法を用いて求めることができ、その手法は特に限定されない。例えば、菌の分類については、リボソームを構成する小型のサブユニット（１６Ｓ）に存在するＲＮＡをコードする遺伝子（１６ＳｒＲＮＡ遺伝子）の塩基配列に基づき、各々の菌を、上記（ａ）〜（ｊ）の科に属する各々の菌群又は、その他の菌群に分けることができる。また、糞便検体に含まれる各々の菌群の菌数の算出についても、例えば、糞便検体中に含まれる、各々の菌に由来する１６ＳｒＲＮＡ遺伝子のＤＮＡ量やコピー数を定量することにより行うこともできる。菌数の算出について、以下に例を挙げる。
【００１９】
塩基配列に基づき菌数を算出するためには、先ず、糞便検体に含まれる菌群のＤＮＡを糞便検体から抽出する。糞便検体からのＤＮＡの抽出は、公知の手法を用いて行うことができ、その手法は特に制限されない。例えば、フェノール・クロロホルムを用いる抽出法により、糞便検体からＤＮＡを抽出することができる。本開示に係る腸内菌叢の老若判定方法においては、フェノール・クロロホルムを用いる抽出方法によりＤＮＡを抽出することが好ましい。
【００２０】
また、フィルターを備えるスピンカラムなどを用いて、糞便検体からＤＮＡを抽出してもよい。なお、本工程Ｓ１１では、ＤＮＡの抽出は必須の工程ではなく、予め糞便検体から抽出されたＤＮＡ溶液を試料とすることもできる。
【００２１】
糞便検体から得られるＤＮＡについては、菌数の測定に必要なＤＮＡ量又はコピー数を得るために、核酸増幅反応を行うことが好ましい。ＤＮＡの増幅については、ＰＣＲ（Ｐｏｌｙｍｅｒａｓｅ  Ｃｈａｉｎ  Ｒｅａｃｔｉｏｎ）等、公知の手法を用いて行うことができ、核酸増幅反応において用いるプライマーについても、公知の設計手法を用いて設計することができる。さらに、核酸増幅反応の際に、ＤＮＡの解析に必要な配列を増幅核酸鎖に組み込ませるように、プライマーを設計してもよい。
【００２２】
上記（ａ）〜（ｊ）の科に属する菌群に分類される菌の数は、例えば、上記ＤＮＡ溶液に含まれるＤＮＡを鋳型とする定量的ＰＣＲ法やＴ−ＲＦＬＰ（Ｔｅｒｍｉｎａｌ−Ｒｅｓｔｒｉｃｔｉｏｎ  Ｆｒａｇｍｅｎｔ  Ｌｅｎｇｔｈ  Ｐｏｌｙｍｏｒｐｈｉｓｍ）法によって測定することができる。定量的ＰＣＲ法では、各々の菌群の１６ＳｒＲＮＡ遺伝子のうち特異的な配列部分に結合するプライマーを用い、核酸増幅反応を行うことにより、所望の菌群の１６ＳｒＲＮＡ遺伝子のＤＮＡ量又はコピー数を測定することができる。またＴ−ＲＦＬＰ法では、末端に標識が付されたプライマーを用いて、１６ＳｒＲＮＡ遺伝子の各菌群に特異的な配列部分を含む増幅核酸鎖を生成し、得られた増幅核酸鎖を制限酵素で処理し断片化してサイズの異なるＤＮＡ断片とする。このＤＮＡ断片を電気泳動により分離し、泳動パターンにおける各々のピークの同定とピーク強度から、所望の菌群の１６ＳｒＲＮＡ遺伝子のＤＮＡ量又はコピー数を測定することができる。
【００２３】
また、上記ＤＮＡ溶液に含まれるＤＮＡについては、ＦＩＳＨ（Ｆｌｕｏｒｅｓｃｅｎｃｅ  ｉｎ  ｓｉｔｕ  ｈｙｂｒｉｄｉｚａｔｉｏｎ）法により、所望の菌群の１６ＳｒＲＮＡ遺伝子のＤＮＡ量又はコピー数を算出してもよい。この他、ＤＮＡマイクロアレイを用いて１６ＳｒＲＮＡ遺伝子のＤＮＡ量又はコピー数を算出することもできる。
【００２４】
１６ＳｒＲＮＡ遺伝子のＤＮＡ量又はコピー数は、上述した方法の他に、シーケンサーを用いて糞便検体から抽出されたＤＮＡの塩基配列を決定し、特定の配列を解読した回数から求めることもできる。また、シーケンサーを用いる方法の場合、同じ操作によって、糞便検体に含まれる上記（ａ）〜（ｊ）の科に属する菌群以外の菌の１６ＳｒＲＮＡ遺伝子についても解読することができる。このため、糞便検体に含まれる上記（ａ）〜（ｊ）の科に属する菌群以外の菌の１６ＳｒＲＮＡ遺伝子のコピー数も併せて求めることができる。
【００２５】
シーケンサーを用いる方法としては、例えば、クローンライブラリー法、ＤＧＧＥ／ＴＧＧＥ（Ｄｅｎａｔｕｒｉｎｇ／Ｔｅｍｐｅｒａｔｕｒｅ  Ｇｒａｄｉｅｎｔ  Ｇｅｌ  Ｅｌｅｃｔｒｏｐｈｏｒｅｓｉｓ）法、メタゲノム解析等がある。クローンライブラリー法では、糞便検体から抽出されたＤＮＡの１６ＳｒＲＮＡ遺伝子の、各々の菌群に特異的な配列部分を増幅して、得られた増幅核酸鎖をプラスミドに挿入して、増幅核酸鎖を含むプラスミドの塩基配列を解読する。ＤＧＧＥ／ＴＧＧＥ法では、クローンライブラリー法と同様にＤＮＡを増幅し、得られた核酸鎖を電気泳動によって核酸鎖の配列に応じて分離し、バンドを切り出して、ゲル内のＤＮＡを回収して当該ＤＮＡの塩基配列を決定する。メタゲノム解析では、糞便検体から得られたＤＮＡを断片化して、必要に応じて核酸増幅反応等を行った後、各々のＤＮＡ断片又は増幅核酸鎖の塩基配列を決定する。
【００２６】
上述した塩基配列の決定は、サンガー法の原理を利用したシーケンサーで行うこともでき、また、次世代シーケンサーで行うこともできる。特にメタゲノム解析は、解読性能がより高い次世代シーケンサーで行うことが好ましい。
【００２７】
シーケンサーを用いて解読された塩基配列から菌群を同定する場合には、解析ソフトウェアやデータベース等を利用する公知の解析方法を利用することができる。データベースとしては、例えば、Genbank、ENA、DDBJ等が挙げられる。また、１６ＳｒＲＮＡ遺伝子に特化したデータベースであるGreengenes database 等も利用することができる。
糞便検体から得られるＤＮＡの解析結果を、これらのデータベースに蓄積された菌群の塩基配列に対して相同性検索することにより、糞便検体に含まれる菌群を同定することができる。いずれかの菌群に同定された配列のうち、上記（ａ）〜（ｊ）の科に属する菌群の１６ＳｒＲＮＡ遺伝子として同定された配列のコピー数を、上記（ａ）〜（ｊ）の科に属する菌群の数とすることができる。
【００２８】
＜糞便検体に含まれる全ての菌の数の測定＞
糞便検体に含まれる全ての菌の数の測定は、公知の手法を用いて求めることができ、その手法は特に限定されない。例えば、上述したシークエンサーを用いる方法においては、１６ＳｒＲＮＡ遺伝子の配列中、糞便検体に含まれる菌で共通の塩基配列部分に結合するプライマーを用いて１６ＳｒＲＮＡ遺伝子の配列を解読することにより、上記（ａ）〜（ｊ）の科に属する菌群以外の菌の１６ＳｒＲＮＡ遺伝子のコピー数も算出することができる。このため、糞便検体に含まれる全菌数は、例えば、糞便検体で測定される全ての菌の１６ＳｒＲＮＡ遺伝子のコピー数の和としてもよい。また、顕微鏡等を用いて、糞便検体に含まれる菌を数えることによって、糞便検体に含まれる全ての菌の数を求めることもできる。本開示に係る腸内菌叢の老若判定方法においては、測定工程Ｓ１１に、１６ＳｒＲＮＡ遺伝子の塩基配列を解読する工程を有する方法が好ましい。なお、本工程Ｓ１１に用いられる糞便検体は、採取後に培養操作が行われていない検体が望ましい。培養操作において、各々の菌の増殖率が異なることによって、採取時の糞便検体における各菌の割合が変化してしまうおそれがあるためである。
【００２９】
＜割合の算出＞
糞便検体に含まれる全ての菌の数に対する上記（ａ）〜（ｊ）の科に属する菌群の各々の割合とは、糞便検体に含まれる全菌数を「１」としたときの相対値である。従って、上記割合は、上記（ａ）〜（ｊ）の科に属する菌群の各々に分類された菌の数と糞便検体に含まれる全菌数から求めることができる。また、上述したように、上記（ａ）〜（ｊ）の科に属する菌群の１６ＳｒＲＮＡ遺伝子のコピー数と、糞便検体に含まれる全て菌の１６ＳｒＲＮＡ遺伝子のコピー数から上記割合を求めることもできる。なお、上記（ａ）〜（ｊ）の科に属する菌群が検出されなかった場合には、その菌群の割合は０とする。
【００３０】
判定工程Ｓ１２は、少なくとも、糞便検体に含まれる全て菌の数に対する（ａ）Lachnospiraceae科に属する菌群の割合に基づき、腸内菌叢を老齢タイプと若齢タイプのいずれかに判定する工程である。また、本開示に係る腸内菌叢の老若判定方法において、上述した菌群の割合は、測定工程Ｓ１１によって求めることもできる。
【００３１】
「老齢タイプ」とは、任意の年齢以上のヒトを高齢集団とし、この高齢集団に特徴的な腸内菌叢のタイプを意味する。一方、「若齢タイプ」とは、任意の年齢未満のヒトを若齢集団とし、この若齢集団に特徴的な腸内菌叢のタイプを意味する。任意の年齢とは、例えば、６５歳とすることができる。
【００３２】
また、例えば、マウスの糞便検体について、本開示に係る腸内細菌叢の老若判定方法を用いる場合、任意の月齢以上のマウスを高齢集団とし、任意の月齢未満のマウスを若齢集団とすることができる。
【００３３】
また、本工程Ｓ１２では、（ａ）Lachnospiraceae科に属する菌群の上記割合に加えて、下記（ｂ）〜（ｊ）の科に属する菌群のうち少なくともいずれか一つの菌群の、糞便検体に含まれる全ての菌の数に対する割合に基づき、腸内菌叢を判定することが好ましい。
（ｂ）Bifidobacteriaceae科
（ｃ）Peptococcaceae科
（ｄ）Enterococcaceae科
（ｅ）Verrucomicrobiaceae科
（ｆ）Catabacteriaceae科
（ｇ）Clostridiaceae科
（ｈ）Caulobacteraceae科
（ｉ）Propionibacteriaceae科
（ｊ）Enterobacteriaceae科
【００３４】
本工程Ｓ１２では、実施例に示すように、（ａ）Lachnospiraceae科に属する菌群の上記割合と、（ｂ）Bifidobacteriaceae科属する菌群の上記割合と、に基づき判定することが好ましく、（ａ）Lachnospiraceae科に属する菌群の上記割合と、（ｂ）Bifidobacteriaceae科に属する菌群の上記割合と、（ｃ）Peptococcaceae科に属する菌群の上記割合と、に基づき判定することがより好ましい。
【００３５】
さらに、（ａ）Lachnospiraceae科に属する菌群の上記割合と、（ｂ）Bifidobacteriaceae科に属する菌群の上記割合と、（ｃ）Peptococcaceae科に属する菌群の上記割合と、（ｄ）Enterococcaceae科に属する菌群の上記割合と、（ｅ）Verrucomicrobiaceae科に属する菌群の上記割合と、に基づき判定することがより好ましい。そして、本工程Ｓ１２では、糞便検体に含まれる全ての菌の数に対する上記（ａ）〜（ｊ）の科に属する菌群の各々の上記割合のうち、これらの割合全てに基づき判定することがより好ましい。
【００３６】
上述した菌群の割合に基づく腸内菌叢の判定には、公知の統計的手法を利用して腸内菌叢を「老齢タイプ」と「若齢タイプ」のいずれかに判定することができる。統計的手法としては、例えば、多変量解析が挙げられる。
【００３７】
また、多変量解析は、上述した菌群の割合を独立変数として、糞便検体を得た個体の年齢又は月齢を従属変数とする回帰分析とすることが好ましい。例えば、回帰分析によって、回帰式を算出して、当該回帰式の解を予め定められた基準値と比較することにより、簡便に判定を行うことができる。従って、本工程Ｓ１２では、菌群ごとの上記割合を独立変数とし、糞便検体を得た個体の年齢又は月齢に関する回帰分析によって回帰式を算出して、当該回帰式に基づき腸内菌叢を判定することもできる。
【００３８】
回帰分析については、公知の分析法の中から適宜選択することができる。回帰分析としては、例えば、単回帰分析、重回帰分析、ロジスティック回帰分析、判別分析等が挙げられる。また、上述した回帰式を判定に用いる場合には、回帰式の解と、予め定められた基準値と、を比較して判定を行ってもよい。さらに、解が基準値を超える場合には腸内菌叢を「老齢タイプ」と判定し、解が基準値以下である場合には腸内菌叢を「若齢タイプ」と判定してもよい。
【００３９】
単回帰分析又は重回帰分析によって回帰式の算出した場合には、回帰式を用いて、上述した菌群の割合から各検体の年齢又は月齢の予測値を回帰式の解として得ることができる。そして、本工程Ｓ１２では、この予測値と予め定められた基準とを比較して、腸内菌叢を「老齢タイプ」と「若齢タイプ」のいずれかに判定してもよい。また、単回帰式又は重回帰式については、例えば、最小二乗法等を用いて算出することができる。
【００４０】
判別分析を行う場合には、量的変数である「糞便検体を得た個体の年齢又は月齢」について、これを質的変数に置き換える。質的変数への置き換えとしては、例えば、所定の年齢又は月齢以上については、「１」とし、所定の年齢又は月齢未満については「−１」とすることができる。所定の年齢又は月齢とは、例えばヒトの場合には、６５歳とすることができる。また、質的変数として用いる２値は、「１」と「−１」の他、「０」と「１」や、「１」と「１００」などとすることもできる。
【００４１】
判別分析において、腸内菌叢を「老齢タイプ」と「若齢タイプ」のいずれかに判定するための回帰式には、回帰式又は重回帰式のように、直線モデル又は平面モデルの線形関数を用いることもできる。そして、これらの回帰式によって、上述した菌群の割合に基づき各検体の年齢又は月齢の予測値を回帰式の解として得ることができる。さらに、本工程Ｓ１２では、この予測値と予め定められた基準とを比較して、腸内菌叢を「老齢タイプ」と「若齢タイプ」のいずれかに判定してもよい。また、判別分析においては、マハラノビスの距離を用いて２次曲線を求めて、これを腸内菌叢の判定に用いることもできる。
【００４２】
ロジスティック回帰分析を行う場合にも、量的変数である、糞便検体を得た個体の年齢について、これを質的変数に置き換えることができる。また、ロジスティック回帰式は、最小二乗法や最尤法によって求めることができる。上述した単回帰式又は重回帰式と同様に、ロジスティック回帰式によっても、上述した菌群の割合に基づき各検体の年齢又は月齢の予測値を得ることができる。そして、この予測値と予め定められた基準とを比較して、腸内菌叢を「老齢タイプ」と「若齢タイプ」のいずれかに判定することができる。
【００４３】
なお、上述した各種回帰式の算出については、計算機等により、上述した菌群の割合と、糞便検体を得た個体の年齢又は月齢から、自動的に回帰式を算出させることもできる。
【００４４】
下記に、本工程Ｓ１２で用いることができる回帰式の一例を示す。下記式（１）中、ｎは独立変数の数であり、ｋは菌群が属する科のうち第ｋ番目の科を意味し、αは回帰式の切片であり、βは回帰式の回帰係数であり、ｘは割合である。
【００４５】
【数２】
2015188380.tif 000003
【００４６】
本工程Ｓ１２において、（ａ）Lachnospiraceae科に属する菌群の上記割合のみに基づいて腸内菌叢を判定する場合には、上記式（１）は、単回帰式となり、ｎは「１」である。また、（ａ）Lachnospiraceae科に属する菌群の上記割合に加えて、上記（ｂ）〜（ｊ）の科に属する菌群のいずれかの上記割合をも用いる場合には、これらの科に属する菌群ごとの上記割合を、各々、独立変数とする。また、「ｋ」の最大値は、独立変数の数と同じである。
【００４７】
本工程Ｓ１２では、「糞便検体を得た個体の年齢」について、これを質的変数に置き換えた後、上記式（１）を算出してもよい。例えば、ヒトから得られた糞便検体について判定する場合、糞便検体を得た個体の年齢について、６５歳以上を「１」に、６５歳未満を「−１」に、各々置き換えて、回帰式を求め、上記式（１）の解ｙが基準値「０」を超える場合には、その腸内菌叢を、６５歳以上に特徴的な「老齢タイプ」と判定することもできる。
【００４８】
本開示に係る腸内菌叢の老若判定方法においては、上述したように、少なくとも、糞便検定に含まれる全ての菌群に対するLachnospiraceae科に属する菌群の割合に基づいて、腸内菌叢が老齢タイプと若齢タイプのいずれであるか、を判定することができる。腸内菌叢のタイプを判別することによって、例えば、腸内菌叢の状態の把握が容易となる。また、腸内菌叢のタイプを判別することは、後述する菌群のスクリーニング方法や物質のスクリーニング方法において、好適に用いられ得る。
【００４９】
（２）菌群のスクリーニング方法
本開示に係る菌群のスクリーニング方法では、腸内菌叢に含まれる菌群のうち、腸内菌叢を老齢タイプ又は若齢タイプに判定するための老若判定マーカーとなる菌群を探索する。
図２は、本開示に係る菌群のスクリーニング方法の工程を示すフローチャートである。
【００５０】
菌群のスクリーニング方法には、少なくとも選別工程Ｓ２１及び菌群同定工程Ｓ２２が含まれる。なお、老若判定マーカーとは、腸内菌叢が老齢タイプであるか否かを判定するためのマーカー、腸内菌叢が若齢タイプであるか否かを判定するためのマーカー、及び腸内菌叢が老齢タイプと若齢タイプのいずれであるかを判定するためのマーカーのうち、少なくともいずれかに該当するものである。
【００５１】
選別工程Ｓ２１は、糞便検体について、上述した腸内菌叢の老若判定方法の結果に基づき糞便検体を複数の群に分ける工程である。腸内菌叢の判定については、上述した通りである。選別工程Ｓ２１の一例を図３に示す。
【００５２】
図３に示す選別工程Ｓ２１では、上述した腸内菌叢の判定により、糞便検体を老齢タイプと若齢タイプのいずれかに分ける。また、本工程Ｓ２１では、上述した（ａ）Lachnospiraceae科に属する菌群の上記割合のみに基づいて、糞便検体を選別してもよいが、Lachnospiraceae科に属する菌群の上記割合に加えて、上記（ｂ）〜（ｊ）のうち、いずれかの科に属する菌群の上記割合にも基づいて、糞便検体を選別してもよい。
【００５３】
菌群同定工程Ｓ２２は、上述した選別工程Ｓ２１によって選別された複数の群のうち、少なくともいずれか一つの群の糞便検体に含まれる菌群を同定する工程である。即ち、老齢タイプに分類された一つの群と、若齢タイプに分類された一つの群のうち、少なくともいずれか一つの群の糞便検体に含まれる菌群を同定する工程である。対象とする群は、菌群のスクリーニングの目的に応じて決めることができる。例えば、老齢タイプと判定するためのマーカー探索が目的であれば、老齢タイプに該当する糞便検体に含まれる菌群を同定することが好ましい。また、若齢タイプと判定するためのマーカー探索が目的であれば、若齢タイプに該当する糞便検体に含まれる菌群を同定することが好ましい。
【００５４】
また、複数の群の各々について、それらの群に含まれる糞便検体を菌群の同定対象としてもよい。例えば老齢タイプに該当する糞便検体と、若齢タイプに該当する糞便検体の各々から菌群を同定することもできる。なお、本開示において「菌群を同定する」とは、少なくとも門レベルまで同定されていることを指し、さらに綱レベル、目レベル、科レベル、属レベル及び種レベルのうち、いずれかのレベルまで同定することも含まれる。
【００５５】
本工程Ｓ２２において、糞便検体に含まれる菌群の同定は、公知の手法を用いて行うことができ、その手法は限定されない。例えば、糞便検体に含まれる菌群のゲノムを解読して、塩基配列に基づいて行うこともできる。この場合、菌群の１６ＳｒＲＮＡ遺伝子の塩基配列を解析して、各々の菌群を同定することができる。以下に、本工程Ｓ２２の一例を示す。
【００５６】
先ず、糞便検体に含まれる菌群のＤＮＡを抽出する。ＤＮＡの抽出については、上述した腸内菌叢の老若判定方法における測定工程Ｓ１１と同様の方法で行うことができる。
【００５７】
次に、得られたＤＮＡの塩基配列を解読する。塩基配列の解読は、公知の手法によって行うことができ、その手法は限定されない。塩基配列の決定については、上述した腸内菌叢の老若判定方法における測定工程Ｓ１１と同様に、例えばシーケンサーを用いる方法等によって行うことができる。また、上述した腸内菌叢の老若判定方法における測定工程Ｓ１１と同様の方法により、例えば、解読されたＤＮＡの塩基配列をデータベースに照会することで、菌群を同定することができる。
【００５８】
本開示に係る菌群のスクリーニング方法では、上述した腸内菌叢の老若判定方法を利用して、糞便検体を複数の群に分けて選別することができる。例えば、予め老齢タイプに該当すると判定された糞便検体は、該当しなかった他の糞便検体に比べて、加齢によって腸内菌叢がより大きく変化している可能性の高い検体である。従って、本開示に係る菌群のスクリーニング方法では、より効率的に老若判定マーカーを探索することができる。
【００５９】
また、例えば、若齢タイプに該当した糞便検体から同定された菌群やその菌数を、老齢タイプを判定するためのマーカー検索において、リファレンスとして用いることもできる。このリファレンスにより、若齢タイプであると判定された糞便検体においても同程度の割合で存在した菌群や菌種を、老齢タイプを判定するためのマーカー候補から除外することができる。
このため、より効率的に老齢タイプと判定するためのマーカーを探索することができる。この点については、老齢タイプに該当した糞便検体から同定された菌群やその菌数を、若齢タイプを判定するためのマーカー検索において、リファレンスとして用いる場合も同様である。
【００６０】
（３）物質のスクリーニング方法
本開示に係る物質のスクリーニング方法では、食餌や物質に含まれる腸内菌叢へ作用するものを探索する。図４は、本開示に係る物質のスクリーニング方法のフローチャートである。物質のスクリーニング方法には、少なくとも測定工程Ｓ３１、相関算出工程Ｓ３２及び物質同定工程Ｓ３３が含まれる。なお、本開示に係る物質のスクリーニング方法において、腸内菌叢へ作用するものとは、例えば、腸内菌叢の老化を促進したり、あるいは抑制したりするなど、腸内菌叢のバランスを変化させるものである。
【００６１】
測定工程Ｓ３１は、上述した菌群のスクリーニング方法により、老若判定マーカーとして同定された菌群について、動物から得た糞便検体に含まれる全菌数に対する割合を測定する。菌群の割合の測定については、上述した通りである。
【００６２】
本開示に係る物質のスクリーニング方法において、食餌又は物質の動物への摂取は、経口摂取であっても、非経口摂取であってもよい。また、摂取回数は、１回でも複数回でもよい。複数回の場合、同じ食餌又は物質を複数回に渡って摂取してもよい。動物から得られる糞便検体は、食餌又は物質を摂取する前と摂取した後の各々の時点で採取されたものであってもよい。
【００６３】
相関算出工程Ｓ３２は、測定工程Ｓ３１で測定された菌群の割合について、その増減と、動物が摂取した食餌又は物質との相関を求める工程である。動物が食餌又は物質を摂取した後、上述した菌群のスクリーニングによって老若判定マーカーとされた菌群の数が糞便検体において変動すれば、その食餌又は物質は腸内菌叢へ作用するものである可能性が高い。また、動物が摂取した食餌又は物質の量に応じて、菌群の割合の変動の程度も変化すれば、その食餌又は物質は腸内菌叢へ作用するものである可能性が高い。本工程Ｓ３２では、このような、動物が摂取した食餌又は物質と、糞便検体の全菌数に対する特定の菌群の割合における相関関係を解析する。
【００６４】
本工程Ｓ３２では、例えば相関係数ｒを算出してもよい。また、例えば、カイ二乗検定やマンホイットニーＵ検定等の検定方法を用いて、食餌又は物質を摂取した群と摂取していない群との間の、老若判定マーカーとした菌群の上記割合の差異について、χ２値やＵ等を求めてもよい。相関について数値で表わすことができれば、その算出方法は特に限定されず、母集団の大きさや測定工程Ｓ３１で測定する菌群の種類の数等に合わせて、公知の統計手法の中から適宜選択できる。
【００６５】
また、相関算出工程Ｓ３２では、任意の統計手法に基づき、老若判定マーカーとして用いた菌群の上記割合の増減と、動物が摂取した食餌又は物質との相関について、計算機等により自動的に相関関係を解析してもよい。
【００６６】
物質同定工程Ｓ３３は、相関算出工程Ｓ３２において算出された値が閾値を満たす場合に、その食餌又は物質を腸内菌叢へ作用するものとして同定する工程である。閾値は、相関算出工程Ｓ３２において相関の算出に用いた方法に合わせて適宜設定できる。
例えば、相関算出工程Ｓ３２において相関係数ｒ、χ２値又はＵを算出するのであれば、本工程Ｓ３３では、これらの数値に対応するｐ値に対して閾値を設定できる。例えば、閾値をｐ＜０．０５とする場合、前記ｐ値がこれに該当すれば、閾値を満たすものとする。また、前記ｐ値がこれに該当しなければ、閾値を満たさないものとする。
【００６７】
上記閾値を満たした場合、その動物が摂取した食餌又は物質は、腸内菌叢へ作用するものとして同定される。例えば、食餌又は物質の摂取後に、老齢タイプと判定するためのマーカーとされる菌群の上記割合が老齢タイプの範囲外へ変動し、相関算出工程Ｓ３２によって算出された値が閾値を満たす場合には、その食餌又は物質は、腸内菌叢の若齢タイプへの変化を促進するものとして同定される。
また、食餌又は物質の摂取後に、若齢タイプと判定するためのマーカーとされる菌群の上記割合が若齢タイプの範囲外へ変動し、その相関が閾値を満たす場合には、その食餌又は物質は、腸内菌叢の老齢タイプへの変化を促進するものとして同定される。
【００６８】
本開示に係る物質のスクリーニング方法では、上述した菌群のスクリーニング方法によって老若判定マーカーとされる菌群の上記割合を指標とすることにより、簡便に腸内菌叢へ作用するものをスクリーニングすることができる。
【００６９】
本開示は、以下の構成を採用することもできる。
〔１〕腸内菌叢について老齢タイプと若齢タイプのいずれであるかを判定する方法であって、糞便検体に含まれる全て菌の数に対する（ａ）Lachnospiraceae科に属する菌群の割合に基づき前記腸内菌叢を判定する判定工程、を有する、腸内菌叢の老若判定方法。
〔２〕前記判定工程では、前記（ａ）Lachnospiraceae科に属する菌群の前記割合に加えて、下記（ｂ）〜（ｊ）の科に属する菌群のうち少なくともいずれか一つの菌群の、前記糞便検体に含まれる全ての菌の数に対する割合に基づき、前記腸内菌叢を判定する、上記〔１〕に記載の腸内菌叢の老若判定方法。
（ｂ）Bifidobacteriaceae科
（ｃ）Peptococcaceae科
（ｄ）Enterococcaceae科
（ｅ）Verrucomicrobiaceae科
（ｆ）Catabacteriaceae科
（ｇ）Clostridiaceae科
（ｈ）Caulobacteraceae科
（ｉ）Propionibacteriaceae科
（ｊ）Enterobacteriaceae科
〔３〕前記判定工程では、前記菌群ごとの前記割合を独立変数とし、前記糞便検体を得た個体の年齢又は月齢に関する回帰分析によって回帰式を算出して、当該回帰式に基づき前記腸内菌叢を判定する、上記〔１〕又は〔２〕に記載の腸内菌叢の老若判定方法。
〔４〕前記判定工程では、前記回帰式の解と、予め定められた基準値と、を比較して、前記解が前記基準値を超える場合には前記腸内菌叢を老齢タイプと判定し、前記解が前記基準値以下である場合には前記腸内菌叢を若齢タイプと判定する、上記〔３〕に記載の腸内菌叢の老若判定方法。
〔５〕前記判定工程では、前記回帰式として下記式（１）を算出し、
下記式（１）の解ｙが０を超える場合には、前記腸内菌叢を老齢タイプと判定する、上記〔４〕に記載の腸内菌叢の老若判定方法。
【００７０】
【数３】
2015188380.tif 000004
【００７１】
（式中、ｎは前記独立変数の数であり、ｋは前記菌群が属する科のうち第ｋ番目の科を意味し、αは前記回帰式の切片であり、βは前記回帰式の回帰係数であり、ｘは前記割合である。）
〔６〕前記判定工程の前に前記割合を測定する測定工程を有する、上記〔１〕〜〔５〕のいずれかに記載の判定方法。
【００７２】
〔７〕上記〔１〕〜〔６〕のいずれかに記載の腸内菌叢の老若判定方法の結果に基づいて糞便検体を複数の群に分ける選別工程と、選別された前記複数の群のうち、少なくともいずれか一つの群の糞便検体に含まれる菌群を同定する菌群同定工程と、を有する、腸内菌叢の老若判定マーカーとなる菌群のスクリーニング方法。
【００７３】
〔８〕上記〔７〕に記載の菌群のスクリーニング方法により同定された菌群について、動物から得た糞便検体に含まれる全て菌の数に対する割合を測定する測定工程と、前記菌群の前記割合の増減と、前記動物が摂取した食餌又は物質と、の相関を求める相関算出工程と、前記相関算出工程において算出された値が閾値を満たす場合に、前記食餌又は物質を腸内菌叢へ作用するものとして同定する物質同定工程と、を含む、物質のスクリーニング方法。
【実施例】
【００７４】
１．実験例１
本実験例では、糞便検体に含まれる菌群について、老齢及び若齢の各々に特徴的な菌群の同定を試みた。
【００７５】
＜検体＞
本実験例では、健常な１７４名（２１〜１０４歳）の日本人から提供された糞便検体を用いた。この１７４名の内訳は、２０歳代：９人、３０歳代：４５人、４０歳代：２１人、５０歳代：１３人、６０歳代：５人、７０歳代：１０人、８０歳代：４６人、９０歳代：１９人、１００歳以上：６人であり、平均年齢は６１．７歳である。
【００７６】
＜ＤＮＡの抽出＞
本実験例では、先ず、糞便検体からのＤＮＡの抽出を行った。約２０ｍｇの糞便検体に対して、４５０μｌの抽出用溶液（１００ｍＭ  Ｔｒｉｓ／ＨＣｌ、４ｍＭ  ＥＤＴＡ、ｐＨ９．０）を加えて懸濁し、さらに、１０％ＳＤＳ溶液５０μｌ、０．１ｍｍ径のガラスビーズ３００ｍｇ、５００μｌのＴＥ飽和フェノール（和光純薬）を加えて混合した。得られた懸濁液に対しては、ＦａｓｔＰｒｅｐ  ＦＰ  １００Ａ（フナコシ社製）を用いてパワーレベル５で３０秒の破砕処理を実施した後、１４，０００×ｇで５分間遠心して４００μｌの上清を得た。得られた上清に２５０μｌのフェノール・クロロホルム溶液（和光純薬）を加えて混合し、１４，０００×ｇで５分間遠心した後、２５０μｌの上清を回収した。回収した上清に、２−プロパノールを２５０μｌ加え、ＤＮＡを沈殿させ、２００μｌのＴｒｉｓ−ＥＤＴＡバッファー（ｐＨ８．０）に再度溶解させて、これをＤＮＡ溶液とした。
【００７７】
＜ＤＮＡの増幅＞
次に、上記ＤＮＡ溶液に含まれる菌群のゲノムに対して核酸増幅反応を行った。下記表１〜表３に、本実験例で用いたプライマーの配列を示す。表１に示すプライマーセットは、１６ＳｒＲＮＡ遺伝子の第３〜４可変領域を増幅させるためのプライマーである。これらのプライマーについては、便宜的に第１プライマーセットとする。
【００７８】
【表１】
2015188380.tif 000005
【００７９】
表１中、「Ｒ」は、グアニン又はアデニンを意味する。「Ｗ」は、アデニン又はチミン若しくはウラシルを意味する。「Ｈ」は、アデニン又はシトシン又はチミン若しくはウラシルを意味する。「Ｖ」は、アデニン又はグアニン又はシトシンを意味する。
【００８０】
また、下記表２及び表３に示すプライマーセットは、次世代シーケンサー（Ｍｉｓｅｑ（イルミナ社製））を用いて塩基配列を解析するために必要な配列（アダプター領域及びインデックス領域）を増幅核酸鎖に組み込むためのプライマーセットである。表２に示すプライマーはフォワードプライマーであり、表３に示すプライマーは、リバースプライマーである。これらのプライマーは、便宜的に第２プライマーセットとする。なお、表１〜表３に示すプライマーは、Ｌｉｆｅ  Ｔｅｃｈｎｏｌｏｇｉｅｓ社のオリゴプライマー作成サービスにより得た。
【００８１】
【表２】
2015188380.tif 000006
【００８２】
【表３】
2015188380.tif 000007
【００８３】
上述したＤＮＡ溶液及び第１プライマーセットを含む総液量２５μｌの反応液を、ＴａＫａＲａ  ＥｘＴａｑ  ＨＳ  ｋｉｔ（タカラバイオ社製）を用いて調製した。核酸増幅反応は、サーマルサイクラー（Ｖｅｒｉｔｉ  ２００（Ｌｉｆｅ  Ｔｅｃｈｎｏｌｏｇｉｅｓ社製））を用いて行い、９４℃−３分の後、９４℃−３０秒、５０℃−３０秒、７２℃−３０秒を１サイクルとし、これを２０回繰り返し、その後７２℃−１０分として、ＰＣＲ反応を行った。得られたＰＣＲ産物を１％アガロースゲルにて電気泳動し、バンドパターンにより、増幅核酸鎖を確認した。
【００８４】
続いて、得られたＰＣＲ産物１μｌを鋳型とし、第２プライマーセットを用いて上述した条件と同様にＰＣＲ反応を行った。但し、サイクル数は１５回とした。得られたＰＣＲ産物を１％アガロースゲルにて電気泳動し、バンドパターンを確認した後、ＱＩＡｑｕｉｃｋ  ９６  ＰＣＲ  Ｐｕｒｉｆｉｃａｔｉｏｎ  Ｋｉｔ（キアゲン社製）を用いてＰＣＲ産物を精製した。精製後のＰＣＲ産物については、Ｑｕａｎｔ−ｉＴ  ＰｉｃｏＧｒｅｅｎ  ｄｓＤＮＡ  Ａｓｓａｙ  ｋｉｔ（Ｌｉｆｅ  Ｔｅｃｈｎｏｌｏｇｉｅｓ社製）によりＤＮＡ濃度を測定した。
【００８５】
＜ＤＮＡのシークエンス＞
複数の糞便検体に由来するＰＣＲ産物を同じ濃度で混合したものをＭｉｓｅｑ  ｖ２  Ｒｅａｇｅｎｔ  ｋｉｔ（イルミナ社製）に供し、Ｍｉｓｅｑにてシークエンス解析を実施した。
【００８６】
＜菌群の同定及び菌数の算出＞
上記シークエンス解析により得られたペアエンド配列をfastq-join （version.1.1.2-301）（http://code.google.com/p/ea-utils/wiki/FastqJoin）にて、ＱＩＩＭＥソフトウェア（version 1.6.0）（http://qiime.org/）にて９７％の相同性を有する配列ごとをＯＵＴ（Operational Taxonomy Unit）とした。各ＯＵＴの代表配列を下記データベースに対してＢＬＡＳＴ検索することにより、腸内菌叢の構成を解析した。
データベース：Greengenes database 12_10
(http://greengenes.secondgenome.com/downloads/database/12_10)
【００８７】
＜結果＞
上述した腸内菌叢の構成の解析の結果、個々の糞便検体について８２９５±２３６６の配列が得られた。本実験例では、この配列数を各糞便検体に含まれる全菌数とした。また、これらの配列のうち、個々の門又は科ごとの数を各々の門又は科に属する菌群の数とした。上記相同性検索の結果、表４及び表５に示す腸内菌叢構成が得られた。表４は、門レベルでの分類を示し、表５は、科レベルでの分類を示す。
【００８８】
【表４】
2015188380.tif 000008
【００８９】
【表５】
2015188380.tif 000009
【００９０】
２．実験例２
本実験例では、上記実験例１で同定された各門又は各科に属する菌群について、糞便検体に含まれる全ての菌の数に対する各々の菌群の割合が、腸内菌叢を「老齢タイプ」と「若齢タイプ」のいずれかに判定する場合に有用であるか検討した。また、本実験例では、年齢が６５歳以上の個体から得られた糞便検体を「老齢タイプ」とし、６４歳以下の個体から得られた糞便検体を「若齢タイプ」とした。
【００９１】
＜回帰式の算出＞
本実験例では、糞便検体の上記２種類のタイプへの分類に際して、回帰式を用いた。本実験例では、表４及び表５に示す門又は科に属する菌群の上記割合を、各々、独立変数とし、糞便検体を得た１７４の個体の年齢を従属変数とし、最小二乗法により、単回帰式を算出した。即ち、本実験例における回帰分析の独立変数及び従属変数は、各々、「１」である。
【００９２】
さらに、従属変数については、回帰式を算出する前に、量的変数から質的変数に置き換えた。具体的には、１７４の糞便検体を得た各々の個体の年齢について、６５歳以上を「１」、６５歳未満を「−１」に置き換えた。
【００９３】
【表６】
2015188380.tif 000010
【００９４】
上記表６に、算出された各回帰式のうち、Lachnospiraceae科、Bifidobacteriaceae科、Clostridiaceae科、Bacteroidaceae科及びEnterococcaceae科の各々に属する菌群の上記割合に基づき算出された回帰式の切片(α）と回帰係数（β）を示す。
【００９５】
また、下記式（２）は、本実験例の回帰式を示す。なお、下記式（２）において、αは切片であり、βは回帰係数であり、ｘは各科に属する菌群の上記割合である。
【００９６】
【数４】
2015188380.tif 000011
【００９７】
本実験例では、算出された切片αと回帰係数βに基づき、上記式（２）の解ｙを求め、基準と比較した。本実験例では、基準値として「０」を設定し、基準値０を超えた場合に、その糞便検体を得た個体の腸内菌叢について、６５歳以上の「老齢タイプ」に分類した。一方、回帰式の解ｙが０以下の場合には、その糞便検体を得た個体の腸内菌叢について、６５歳未満の「若齢タイプ」に分類した。
【００９８】
＜結果＞
表７に本実験例の結果を示す。表７は、Lachnospiraceae科、Bifidobacteriaceae科、Clostridiaceae科、Bacteroidaceae科及びEnterococcaceae科の、各々に属する菌群について、回帰式の解ｙと基準値に基づく「老齢タイプ」又は「若齢タイプ」への分類と、糞便検体を得た個体の実際の年齢（実年齢）と、を比較した結果を示す。また、正答率（％）とは、検体の分類結果と検体を得た個体の実年齢とが一致した検体数の、検体数全体に占める割合（％）を示す。さらに、検出率は、１７４の糞便検体のうち、各々の科に属する菌群が検出された検体の割合を示す。なお、各々の科に属する菌群が検出されなかった糞便検体については、「判定結果と実年齢が一致しない検体数」にカウントした。
【００９９】
【表７】
2015188380.tif 000012
【０１００】
Bifidobacteriaceae科に属する菌群は、多くのヒトの腸内菌叢に存在し、また、加齢に伴い腸内で菌数が減少することが知られている菌群である。Bifidobacteriaceae科に属する菌群の上記割合に基づく糞便検体の分類において、正答率は６７．２％であった。
【０１０１】
また、Clostridiaceae科、Bacteroidaceae科及びEnterococcaceae科の、各々に属する菌群の上記割合に基づく糞便検体の分類では、正答率は、各々、６４．９％、６２．１％、５３．４％であった。一方、Lachnospiraceae科に属する菌群の上記割合に基づく糞便検体の分類では、正答率は、８３．３％であった。
【０１０２】
本実験例の結果から、Lachnospiraceae科に属する菌群の上記割合に基づき、糞便検体を分類した場合、従来、加齢に伴い菌数が変動することが知られているBifidobacteriaceae科の上記割合に基づき分類した場合よりも高い正答率で、検体を分類できることが示された。また、１７４の全ての検体で、Lachnospiraceae科に属する菌群が検出された。
【０１０３】
上述した結果から、糞便検体に含まれる全ての菌の数に対するLachnospiraceae科に属する菌群の割合は、糞便検体を「老齢タイプ」と「若齢タイプ」のいずれかに判定する判定方法において、有用な値であることが示された。従って、Lachnospiraceae科に属する菌群の上記割合に基づき糞便検体を判定することにより、精度高く糞便検体を判定できることが示された。
【０１０４】
さらに、Lachnospiraceae科に属する菌群の検出率は、１７４検体において１００％であった。このことから、Lachnospiraceae科に属する菌群の上記割合を用いる腸内菌叢の老若判定方法は、複数の糞便検体において、検出される検体が限られる菌群の上記割合を用いる方法に比べて、汎用性の高い方法であることが示された。
【０１０５】
３．実験例３
上述した実験例２において、糞便検体に含まれる全ての菌の数に対するLachnospiraceae科に属する菌群の割合に基づいて判定することにより、糞便検体を得た個体の腸内菌叢を「老齢タイプ」と「若齢タイプ」のいずれかに精度高く判定できることが明らかとなった。そこで、本実験例では、Lachnospiraceae科に属する菌群の上記割合と、他の科に属する菌群の上記割合と、を組み合わせることにより、腸内菌叢の判定の精度がより高められるか検証した。
【０１０６】
本実験例では、上記実験例２で用いた（ａ）Lachnospiraceae科に属する菌群の上記割合に加えて、（ｂ）Bifidobacteriaceae科、（ｃ）Peptococcaceae科、（ｄ）Enterococcaceae科、（ｅ）Verrucomicrobiaceae科、（ｆ）Catabacteriaceae科、（ｇ）Clostridiaceae科、（ｈ）Caulobacteraceae科、（ｉ）Propionibacteriaceae科、（ｊ）Enterobacteriaceae科、の各々に属する菌群の割合について、糞便検体に含まれる全ての細菌の数に対する各々の菌群の割合を用いた。具体的には、下記の４つの組み合わせについて、検討した。
パターンＡ：（ａ）Lachnospiraceae科、（ｂ）Bifidobacteriaceae科
パターンＢ：（ａ）Lachnospiraceae科、（ｂ）Bifidobacteriaceae科、（ｃ）Peptococcaceae科
パターンＣ：（ａ）Lachnospiraceae科、（ｂ）Bifidobacteriaceae科、（ｃ）Peptococcaceae科、（ｄ）Enterococcaceae科、（ｅ）Verrucomicrobiaceae科
パターンＤ：（ａ）Lachnospiraceae科、（ｂ）Bifidobacteriaceae科、（ｃ）Peptococcaceae科、（ｄ）Enterococcaceae科、（ｅ）Verrucomicrobiaceae科、（ｆ）Catabacteriaceae科、（ｇ）Clostridiaceae科、（ｈ）Caulobacteraceae科、（ｉ）Propionibacteriaceae科、（ｊ）Enterobacteriaceae科
【０１０７】
＜回帰式の算出＞
本実験例では、糞便検体の分類に際して、実験例２と同様に、回帰式を用いた。回帰式については、各科に属する菌群の上記割合を独立変数とし、糞便検体を得た１７４の個体の年齢を従属変数として、最小二乗法により算出した。即ち、パターンＡの独立変数は「２」であり、パターンＢの独立変数は「３」であり、パターンＣの独立変数は「５」であり、パターンＤの独立変数は「１０」である。
【０１０８】
さらに、従属変数については、実験例２と同様に、回帰式を算出する前に、１７４の糞便検体を得た各々の個体の年齢について、６５歳以上を「１」、６５歳未満を「−１」に置き換えた。各科に属する菌群の上記割合と、糞便検体を得た１７４の個体の年齢ついて、重回帰分析を行い、下記式（１）に示す回帰式を得た。
【０１０９】
【数５】
2015188380.tif 000013
【０１１０】
上記式（１）において、ｎは独立変数の数であり、ｋは菌群が属する科のうち第ｋ番目の科を意味し、αは回帰式の切片であり、βは回帰式の回帰係数であり、ｘは各科に属する菌群の上記割合である。
【０１１１】
また、下記表８に、各パターンの回帰式の切片と回帰係数（偏回帰係数）を示す。なお、菌群の括弧内のアルファベットは、上述した（ａ）〜（ｊ）の各科に対応する。
【０１１２】
【表８】
2015188380.tif 000014
【０１１３】
本実験例では、算出された切片αと回帰係数（偏回帰係数）βに基づき、回帰式の解ｙを求め、基準と比較した。基準値は、実験例２と同様に、「０」を設定し、基準値０を超えた場合に、その糞便検体を得た個体の腸内菌叢について、６５歳以上の「老齢タイプ」と判定した。
【０１１４】
＜結果＞
表９に本実験例の結果を示す。表９は、回帰式の解ｙと基準値に基づく「老齢タイプ」又は「若齢タイプ」への分類と、糞便検体を得た個体の実際の年齢（実年齢）と、を比較した結果を示す。また、正答率（％）とは、検体の分類結果と検体を得た個体の実年齢とが一致した検体数の、検体数全体に占める割合（％）を示す。
【０１１５】
【表９】
2015188380.tif 000015
【０１１６】
表９のパターンＡに示すように、Lachnospiraceae科に属する菌群の上記割合に加えて、Bifidobacteriaceae科に属する菌群の上記割合に基づき糞便検体を分類すると、その正答率は、９０％を超えた。また、パターンＢに示すように、Lachnospiraceae科に属する菌群の上記割合及びBifidobacteriaceae科に属する菌群の上記割合に加えて、Peptococcaceae科に属する菌群の上記割合に基づき糞便検体を分類すると、その正答率は、９２．５％であった。さらに、パターンＣでは、正答率が９３．７％であり、パターンＤでは、正答率が９４．８％であった。
【０１１７】
本実験例の結果から、Lachnospiraceae科に属する菌群の上記割合に加えて、Bifidobacteriaceae科、Peptococcaceae科、Enterococcaceae科、Verrucomicrobiaceae科、Catabacteriaceae科、Clostridiaceae科、Caulobacteraceae科、Propionibacteriaceae科又はEnterobacteriaceae科ののうち、少なくともいずれかの一つの科に属する菌群の上記割合に基づき糞便検体を分類した場合には、正答率をさらに高められることが示された。
【０１１８】
また、実験例２において示されたように、Bifidobacteriaceae科に属する菌群の上記割合、Clostridiaceae科に属する菌群の上記割合及びEnterococcaceae科に属する菌群の上記割合は、単独で腸内細菌叢の判定に用いた場合には、正答率が５０〜６０％程度であった。しかし、これらの菌群の上記割合については、Lachnospiraceae科に属する菌群の上記割合と組み合わせて判定に用いた場合には、腸内菌叢の老若の判定に有用であることが示された。
【図１】
2015188380.tif 000017
【図２】
2015188380.tif 000018
【図３】
2015188380.tif 000019
【図４】
2015188380.tif 000020
【配列表】
2015188380000001.app