(57)【特許請求の範囲】
【請求項１】
被験者から採取した唾液検体中の細菌叢の１６ＳリボソームＲＮＡ遺伝子の塩基配列を無作為に決定することによってその塩基配列のデータ群を得、得られた塩基配列のデータ群を、該データ群に含まれる各塩基配列間の、閾値を９６％以上に設定した類似度に基づいたクラスタリングに供してクラスターグループ数を求め、そのクラスターグループ数が、健常者群が呈する該クラスターグループ数の基準と比べて有意に減少しているか否かによって炎症性腸疾患と評価される唾液検体を検出することを特徴とする炎症性腸疾患の検出方法。
【請求項２】
被験者から採取した唾液検体中の細菌叢の１６ＳリボソームＲＮＡ遺伝子の塩基配列を無作為に決定することによってその塩基配列のデータ群を得、得られた塩基配列のデータ群を、該データ群に含まれる各塩基配列間の、閾値を９６％以上に設定した類似度に基づいたクラスタリングに供してクラスターグループを求め、各クラスターグループに含まれる塩基配列のデータ数と、各クラスターグループに属する代表塩基配列どうしの類似度とに基づいてＵｎｉＦｒａｃ解析を行い、健常者群との群間類似距離が有意に離れているか否かによって炎症性腸疾患と評価される唾液検体を検出することを特徴とする炎症性腸疾患の検出方法。
【請求項３】
前記群間類似距離の評価を主座標分析により行う請求項２記載の炎症性腸疾患の検出方法。
【請求項４】
被験者から採取した唾液検体中の細菌叢の１６ＳリボソームＲＮＡ遺伝子の塩基配列を無作為に決定することによってその塩基配列のデータ群を得、Firmicutes門、Bacteroidetes門、Actinobacteria門、及びProteobacteria門から選ばれた１種又は２種以上の門に属する菌種について、前記データ群に含まれる該菌種に属する塩基配列のデータ数を求め、
Firmicutes門及び／又はBacteroidetes門の場合には、そのデータ数が、健常者群が呈する該データ数の基準と比べて有意に増加しているか否かによって、
Actinobacteria門及び／又はProteobacteria門の場合には、そのデータ数が、健常者群が呈する該データ数の基準と比べて有意に減少しているか否かによって、
炎症性腸疾患と評価される唾液検体を検出することを特徴とする炎症性腸疾患の検出方法。
【請求項５】
被験者から採取した唾液検体中の細菌叢の１６ＳリボソームＲＮＡ遺伝子の塩基配列を無作為に決定することによってその塩基配列のデータ群を得、Prevotella属、Streptococcus属、Veillonella属、Atopobium属、Megasphaera属、Solobacterium属、Actinomyces属、Lachnospiraceae属、Selenomonas属、Rothia属、Haemophilus属、Neisseria属、Gemella属、Porphyromonas属、Corynebacterium属、Capnocytophaga属、Bergeyella属、Aggregatibacter属、Lautropia属、Paludibacter属、Fusobacterium属、Tannerella属、及びPropionibacterium属から選ばれた１種又は２種以上の属に属する菌種について、前記データ群に含まれる該菌種に属する塩基配列のデータ数を求め、
Prevotella属、Streptococcus属、Veillonella属、Atopobium属、Megasphaera属、Solobacterium属、Actinomyces属、Lachnospiraceae属、及び／又はSelenomonas属の場合には、そのデータ数が、健常者群が呈する該データ数の基準と比べて有意に増加しているか否かによって、
Rothia属、Haemophilus属、Neisseria属、Gemella属、Porphyromonas属、Corynebacterium属、Capnocytophaga属、Bergeyella属、Aggregatibacter属、Lautropia属、Paludibacter属、Fusobacterium属、Tannerella属、及び／又はPropionibacterium属の場合には、そのデータ数が、健常者群が呈する該データ数の基準と比べて有意に減少しているか否かによって、
炎症性腸疾患と評価される唾液検体を検出することを特徴とする炎症性腸疾患の検出方法。
【請求項６】
被験者から採取した唾液検体中の細菌叢の１６ＳリボソームＲＮＡ遺伝子の塩基配列を無作為に決定することによってその塩基配列のデータ群を得、Veillonella atypical、Granulicatella adiacens、Atopobium parvulum、Streptococcus sp. EO2001-02、Veillonella dispar、Streptococcus mitis、Prevotella DO039、Megasphaera micronuciformis、Prevotella salivae、Solobacterium moorei、Prevotella melaninogenica、Actinomyces graevenitzii、Neisseria subflava、Haemophilus nbw161b08c1、Haemophilus sp. CCUG 32367、Rothia aeria、Lautropia mirabilis、Rothia mucilaginosa、Campylobacter gracilis、Corynebacterium matruchotii、Bergeyella 602D02、Capnocytophaga gingivalis、Corynebacterium durum、Gemella haemolysans、Porphyromonas CW034、Streptococcus VG051、Fusobacterium nucleatum、Prevotella IK062、Veillonella dispar、及びHaemophilus parainfluenzaeから選ばれた１種又は２種以上の菌種について、前記データ群に含まれる該菌種に属する塩基配列のデータ数を求め、
Veillonella atypical、Granulicatella adiacens、Atopobium parvulum、Streptococcus sp. EO2001-02、Veillonella dispar、Streptococcus mitis、Prevotella DO039、Megasphaera micronuciformis、Prevotella salivae、Solobacterium moorei、Prevotella melaninogenica、及び／又はActinomyces graevenitziiの場合には、そのデータ数が、健常者群が呈する該データ数の基準と比べて有意に増加しているか否かによって、
Neisseria subflava、Haemophilus nbw161b08c1、Haemophilus sp. CCUG 32367、Rothia aeria、Lautropia mirabilis、Rothia mucilaginosa、Campylobacter gracilis、Corynebacterium matruchotii、Bergeyella 602D02、Capnocytophaga gingivalis、Corynebacterium durum、Gemella haemolysans、Porphyromonas CW034、Streptococcus VG051、Fusobacterium nucleatum、Prevotella IK062、Veillonella dispar、及び／又はHaemophilus parainfluenzaeの場合には、そのデータ数が、健常者群が呈する該データ数の基準と比べて有意に減少しているか否かによって、
炎症性腸疾患と評価される唾液検体を検出することを特徴とする炎症性腸疾患の検出方法。
【請求項７】
被験者から採取した唾液検体中の細菌叢の１６ＳリボソームＲＮＡ遺伝子の塩基配列を無作為に決定することによってその塩基配列のデータ群を得、配列番号１〜３０から選ばれた１種又は２種以上の塩基配列について、前記データ群に含まれる該塩基配列のデータ数を求め、
配列番号２８、配列番号９、配列番号１３、配列番号１８、配列番号７、配列番号１６、配列番号２１、配列番号８、配列番号６、配列番号１１、配列番号１２、及び／又は配列番号１４の場合には、そのデータ数が、健常者群が呈する該データ数の基準と比べて有意に増加しているか否かによって、
配列番号１、配列番号２、配列番号３、配列番号４、配列番号５、配列番号１０、配列番号１５、配列番号１７、配列番号１９、配列番号２０、配列番号２２、配列番号２３、配列番号２４、配列番号２５、配列番号２６、配列番号２７、配列番号２９、及び／又は配列番号３０の場合には、そのデータ数が、健常者群が呈する該データ数の基準と比べて有意に減少しているか否かによって、
炎症性腸疾患と評価される唾液検体を検出することを特徴とする炎症性腸疾患の検出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、炎症性腸疾患の検出方法及びヒト唾液細菌叢の検査方法に関し、更に詳細には、ヒトの唾液に含まれる細菌叢の菌叢構造に基づく炎症性腸疾患の検出方法及びそのためのヒト唾液細菌叢の検査方法に関する。
【背景技術】
【０００２】
消化管の粘膜に慢性の炎症または潰瘍をひきおこす原因不明の疾患が知られ、総称して炎症性腸疾患(Inflammatory Bowel Disease：IBD)と呼ばれている。炎症性腸疾患のうち、例えばクローン病は、主として若年層に発症率が高く、好発部位としては小腸であるが口腔から肛門にいたるまでの消化管のあらゆる部位に炎症や潰瘍を引き起こし、非連続性の病変を呈するのが特徴とされている。また、例えば潰瘍性大腸炎は、若年層から壮年・老年層にも広くみられ、大腸の粘膜にびらんや潰瘍ができる炎症性疾患である。
【０００３】
炎症性腸疾患では、上記消化管での炎症や潰瘍等により腹痛や下痢、血便、体重減少などの症状が生じる。薬剤による内科的治療や、症状が重い場合は外科手術により病変の除去が行なわれるが、根治を得ることは難しく、いったんは寛解に向かっても、長期にわたって再発を繰り返すケースが多い。日本では炎症性腸疾患の患者は１０万人以上とも報告されており、このうちクローン病や潰瘍性大腸炎は厚生労働省の特定難治性疾患にも指定されている。
【０００４】
炎症性腸疾患の臨床的診断では、一般的な血液・糞便検査、より精密に病変部位や病態を特定するためＸ線造影検査や内視鏡検査などが行なわれるが、予防的に疾患を早期検出するためには、健康診断や集団検診などにも応用できる、効率的なスクリーニング方法が望まれている。
【０００５】
このような問題に対して、例えば下記特許文献１〜３には、炎症性腸疾患の患者群と健常者群とで、それぞれに相応した発現量を示す物質を特定して診断マーカーに用いることが記載されている。
【０００６】
一方、炎症性腸疾患の発症や病態の変化には腸内常在菌（腸内細菌叢）が関与することが報告されている（下記非特許文献１参照）。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００４−１３５５４５号公報
【特許文献２】特表２００５−５０４３２１号公報
【特許文献３】特許第４６２５７０８号公報
【非特許文献】
【０００８】
【非特許文献１】Xavier RJ and Podolsky DK: Unravelling the pathogenesis of inflammatory bowel disease. Nature 448: 427-434 (2007)
【発明の概要】
【発明が解決しようとする課題】
【０００９】
従来、Ｘ線造影検査や内視鏡検査などの侵収性検査は患者に多くの負担と苦痛を強いるものであった。また、上記特許文献１〜３のように、特定の診断マーカーを利用した場合には、かならずしも病態をよく反映した評価であるとはいえないという問題があった。
【００１０】
そこで、本発明の目的は、簡便な検査で病態をよく反映した評価を行うことができる、炎症性腸疾患の検出方法及びヒト唾液細菌叢の検査方法を提供することにある。
【課題を解決するための手段】
【００１１】
本発明者らは、上記目的を達成するために鋭意研究し、本発明を完成するに至った。
すなわち、本発明の第１は、被験者から採取した唾液検体中の細菌叢の１６ＳリボソームＲＮＡ遺伝子の塩基配列を無作為に決定することによってその塩基配列のデータ群を得、得られた塩基配列のデータ群に基づいて炎症性腸疾患と評価される唾液検体を検出することを特徴とする炎症性腸疾患の検出方法を提供するものである。
【００１２】
上記発明においては、前記データ群に基づいて細菌叢を構成する菌種数を求め、その菌種数を指標にして炎症性腸疾患と評価される唾液検体を検出することが好ましい。
【００１３】
また、前記データ群を、該データ群に含まれる各塩基配列間の類似度に基づいたクラスタリングに供してクラスターグループを求め、そのクラスターグループ数を指標にして炎症性腸疾患と評価される唾液検体を検出することが好ましい。
【００１４】
また、前記データ群を、該データ群に含まれる各塩基配列間の類似度に基づいたクラスタリングに供してクラスターグループを求め、各クラスターグループに含まれる塩基配列のデータ数と、各クラスターグループに属する代表塩基配列どうしの類似度を指標にして、各被験者から取得された該データ群との群間類似距離を求め、その群間類似距離に基づいて炎症性腸疾患と評価される唾液検体を検出することが好ましい。
【００１５】
また、Firmicutes門、Bacteroidetes門、Actinobacteria門、及びProteobacteria門から選ばれた１種又は２種以上の門に属する菌種について、前記データ群に含まれる該菌種に属する塩基配列のデータ数を求め、そのデータ数を指標にして健常者群と比較して、炎症性腸疾患と評価される唾液検体を検出することが好ましい。
【００１６】
また、Prevotella属、Streptococcus属、Veillonella属、Atopobium属、Megasphaera属、Solobacterium属、Actinomyces属、Lachnospiraceae属、Selenomonas属、Rothia属、Haemophilus属、Neisseria属、Gemella属、Porphyromonas属、Corynebacterium属、Capnocytophaga属、Bergeyella属、Aggregatibacter属、Lautropia属、Paludibacter属、Fusobacterium属、Tannerella属、及びPropionibacterium属から選ばれた１種又は２種以上の属に属する菌種について、前記データ群に含まれる該菌種に属する塩基配列のデータ数を求め、そのデータ数を指標にして健常者群と比較して、炎症性腸疾患と評価される唾液検体を検出することが好ましい。
【００１７】
また、Veillonella atypical、Granulicatella adiacens、Atopobium parvulum、Streptococcus sp. EO2001-02、Veillonella dispar、Streptococcus mitis、Prevotella DO039、Megasphaera micronuciformis、Prevotella salivae、Solobacterium moorei、Prevotella melaninogenica、Actinomyces graevenitzii、Neisseria subflava、Haemophilus nbw161b08c1、Haemophilus sp. CCUG 32367、Rothia aeria、Lautropia mirabilis、Rothia mucilaginosa、Campylobacter gracilis、Corynebacterium matruchotii、Bergeyella 602D02、Capnocytophaga gingivalis、Corynebacterium durum、Gemella haemolysans、Porphyromonas CW034、Streptococcus VG051、Fusobacterium nucleatum、Prevotella IK062、Veillonella dispar、及びHaemophilus parainfluenzaeから選ばれた１種又は２種以上の菌種について、前記データ群に含まれる該菌種に属する塩基配列のデータ数を求め、そのデータ数を指標にして健常者群と比較して、炎症性腸疾患と評価される唾液検体を検出することが好ましい。
【００１８】
また、配列番号１〜３０から選ばれた１種又は２種以上の塩基配列について、前記データ群に含まれる該塩基配列のデータ数を求め、そのデータ数を指標にして健常者群と比較して、炎症性腸疾患と評価される唾液検体を検出することが好ましい。
【００１９】
一方、本発明の第２は、被験者から採取した唾液検体中の細菌叢の１６ＳリボソームＲＮＡ遺伝子の塩基配列を無作為に決定することによってその塩基配列のデータ群を得、得られた塩基配列のデータ群を、炎症性腸疾患と評価される唾液検体を検出するために用いることを特徴とするヒト唾液細菌叢の検査方法を提供するものである。
【発明の効果】
【００２０】
本発明の炎症性腸疾患の検出方法及びヒト唾液細菌叢の検査方法によれば、収集が容易な唾液を対象にするので、多数の検体を取り扱う健康診断や集団検診などに適応でき、内視鏡等の精密診断を受けるべき炎症性腸疾患の疑いの高い被験者の効率的なスクリーニングが可能となる。また、唾液に生息する細菌から構成される唾液細菌叢の菌叢構造に基づく評価を行なうので、炎症性腸疾患の病態をよく反映した検出又は検査を行うことができる。
【図面の簡単な説明】
【００２１】
【図１】１６ＳリボソームＲＮＡ遺伝子のV1-V2領域のPCR増幅部分を示す説明図である。
【図２】フォワードプライマー27Fmod-454A（配列番号３１）のバーコード配列の役割について、４検体を解析した場合を例として説明する説明図である。
【図３】Operational Taxonomic Unit (OTU) 解析の概略説明図である。
【図４】健常者とクローン病（ＣＤ）患者の唾液細菌叢のOTU数を比較した結果を示す図表である。
【図５】OTU 解析に供される配列データ数の影響を示す図表であり、図５Ａは試行回ごとの有意差（t-testによるp値）の結果を示す図表であり、図５Ｂはその平均の結果を示す図表である。
【図６】１０名の健常者群および１０名のクローン病（CD）患者群の各被験者間の菌叢構造類似度を２次元散布図で表した結果を示す図表である。
【図７】UniFrac解析およびそれに基づく主座標分析に供される配列データ数の影響を示す図表であり、50データを用いたときの結果を示す図表である。
【図８】UniFrac解析およびそれに基づく主座標分析に供される配列データ数の影響を示す図表であり、100データを用いたときの結果を示す図表である。
【図９】UniFrac解析およびそれに基づく主座標分析に供される配列データ数の影響を示す図表であり、2,000データを用いたときの結果を示す図表である。
【発明を実施するための形態】
【００２２】
本発明は、唾液細菌叢の菌叢構造が、炎症性腸疾患の病態をよく反映していることに基づいている。即ち、被験者から採取した唾液検体を検査し、そこに含まれる細菌叢の菌叢構造が呈する様相を調べることにより、炎症性腸疾患か否かを峻別できるという発想に基づくものである。以下に具体的に説明する。
【００２３】
本発明においては、被験者から採取した唾液検体中の細菌叢の１６ＳリボソームＲＮＡ遺伝子の塩基配列を無作為に決定することによってその塩基配列のデータ群を得ることを要する。ここで「無作為に決定する」とは、塩基配列を決定すべきＤＮＡ等混合物中から、ある特定の配列のものだけ選択したり、ある配列のものを排除したりすることなく、できるだけランダムに、その混合物中に含まれる任意の塩基配列を決定するという意味である。そのようにして塩基配列を無作為に決定すると、ＤＮＡ等混合物中での濃度が高い場合には、塩基配列決定の対象となる頻度（チャンス）が多くなり、逆に低い場合には、塩基配列決定の対象となる頻度（チャンス）が少なくなる。一方、唾液細菌叢は複数種類の細菌を含み、それぞれの菌種に特有の塩基配列を１６ＳリボソームＲＮＡ遺伝子上に有しているので、唾液細菌叢から調製したＤＮＡ等混合物の塩基配列を無作為に決定することによって得られる塩基配列のデータ群中では、任意に選択され得る塩基配列に合致する塩基配列（例えばある特定の菌種に特有の塩基配列に合致するもの）のデータ数が、唾液細菌叢を構成する細菌の存在量の多寡を反映している。即ち、その塩基配列のデータ群が、唾液細菌叢の菌叢構造を反映することになる。
【００２４】
上記のように、細菌叢の１６ＳリボソームＲＮＡ遺伝子の塩基配列を無作為に決定することによって得られる塩基配列のデータ群が、唾液細菌叢の菌叢構造を反映するようにするためには、その塩基配列決定の手法が重要であることは勿論である。具体的には、塩基配列を決定すべきＤＮＡ等混合物中から、ある特定の配列のものを指向するではなく、できるだけランダムに、その混合物中に含まれる任意の塩基配列を決定する必要がある。そのような手法は、近年、メタゲノム解析のためのシークエンス手法（塩基配列決定の手法）として知られ、後述実施例でも用いたように、そのためのシークエンス装置やシークエンスキット等も市販されているので、それらを利用し又はそれらに準じて行なうことができる。
【００２５】
上記のように、細菌叢の１６ＳリボソームＲＮＡ遺伝子の塩基配列を無作為に決定することによって得られる塩基配列のデータ群が、唾液細菌叢の菌叢構造を反映するようにするためには、被験者から採取した唾液検体から細菌のゲノムＤＮＡを効率的に抽出することも重要であるが、それは当業者に周知の手法により適宜実施し得る。具体的には、細菌叢に含まれる細菌のゲノムＤＮＡを、菌種によって抽出の効率にバラツキがないように抽出して、塩基配列決定のためのサンプルとする必要がある。そのような手法は、後述実施例でも用いたように、従来周知の手法が種々知られているので、それらを利用し又はそれらに準じて行なうことができる。
【００２６】
上記のように、細菌叢の１６ＳリボソームＲＮＡ遺伝子の塩基配列を無作為に決定することによって得られる塩基配列のデータ群が、唾液細菌叢の菌叢構造を反映するようにするためには、１６ＳリボソームＲＮＡ遺伝子上のどの領域の塩基配列を決定するかも重要であるが、それは当業者に周知の知見により適宜実施し得る。例えば、真正細菌の１６ＳリボソームＲＮＡ遺伝子上には、塩基配列が保存されずに菌種間で変化に富んだ領域Ｖ１〜Ｖ９が知られているので、それらのいずれかを含むように塩基配列を決定することが好ましい（図１）。また、塩基配列の決定のための手法として周知のＰＣＲ増幅やシークエンス等に適した、１００〜５００塩基長程度の領域の塩基配列を決定することが好ましい。更に、そのＰＣＲ増幅等のためのプライマー部位は、菌種間で比較的普遍的に保存されている領域に設定することが、プライマー設計上普遍的な配列を設計し易いので、好ましい。
【００２７】
以下、被験者から採取した唾液検体を検査する手順の一例を説明する。まず、２ｍＬ程度の唾液から細菌のゲノムＤＮＡを抽出して唾液細菌叢ＤＮＡを調製する。この唾液細菌叢ＤＮＡには細菌叢を構成する各細菌のゲノムＤＮＡが混合物として含まれている。次に、唾液細菌叢ＤＮＡの溶液中の二本鎖ＤＮＡ濃度を測定し、その測定値に基づいて４０ｎｇ程度のＤＮＡを鋳型として、この例においては、ユニバーサルプライマーセット（フォワードプライマー27Fmod-454A（配列番号３１）とリバースプライマー338R-454B（配列番号３２））を用いて、１６ＳリボソームＲＮＡ遺伝子（以下、16S遺伝子）のV1-V2領域（図１）をＰＣＲ増幅する。得られたＰＣＲ産物から過剰のプラーマーや基質のヌクレオチド等を除去して精製し、二本鎖ＤＮＡ量を正確に測定して、他の被験者からの同様に処理したＰＣＲ産物とともに、次世代型超高速シークエンス装置に供する。シークエンス装置としては、例えばロシュ社製GS FLX+ System（約１００万データ／稼働、平均７００塩基長／データ）、あるいは、小規模タイプである同GS Junior System（約１０万データ／稼働、平均４００塩基長／データ）など用いることができる。各被験者由来のサンプルについて取得された粗配列データのうち、クオリティが低いものを排除して、全粗配列データのうち６０〜７０％のデータを高精度配列データとして選択し、以下の評価に用いる塩基配列のデータ群とする。なお、通常、上記シークエンサーを用いて１０万または１００万程度の粗配列データを取得するのに要する稼働時間は、１０〜２０時間である。
【００２８】
本発明においては、上記のように得られた塩基配列のデータ群に基づいて炎症性腸疾患と評価される唾液検体を検出する。その評価方法としては、以下のような方法が挙げられるが、これらは限定的なものではなく、健常者が呈する唾液細菌叢の菌叢構造との相異をよく識別できる指標であれば、適宜そのような指標を設定して評価を行えばよい。
【００２９】
（１）唾液細菌叢を構成する細菌の菌種数を指標にする方法
後述の実施例において実証されるように、炎症性腸疾患の患者では唾液細菌叢を構成する細菌の菌種数が健常者よりも有意に少なくなっている。よって、任意の唾液検体について、上記のように得られた塩基配列のデータ群に基づいて、唾液細菌叢を構成する菌種数を求め、その値が健常者の値と比較して有意に低いと判定される場合に、唾液検体が炎症性腸疾患に由来するものと評価することができる。健常者の値と比較して有意に低いかどうかの判定は、集団検診などで多数の被験者がいる場合には、その被験者どうしの統計的差異に基づいて行ってもよく、一方、被験者が単独または少数の場合には、過去のデータから基準値（例えば所定の菌種数）を予め設定しておき、その基準値を閾値にして行なってもよい。なお、この方法の場合、安定に有意な評価を得るためには、上記データ群において１被験者あたり少なくともおよそ１０００データ以上の配列データ数を有していることが好ましい（後述の実施例参照）。
【００３０】
（２）取得したデータ群に含まれる各塩基配列間の類似度に基づいたクラスターグループ数を指標にする方法
高い配列類似度を互いに有する塩基配列データは、既知配列データベース等に照合して菌種を具体的に特定するまでもなく、同じ種の細菌に由来すると推測できる。したがって、上記のように得られた塩基配列のデータ群を、そのデータ群に含まれる各塩基配列間の類似度に基づいたクラスタリング（例えば９６％塩基配列類似度を閾値として設定して）に供することによって得られるクラスターグループは、細菌叢を構成する細菌の菌種に対応し、そのクラスターグループ数は、細菌叢を構成する細菌の菌種数と等価と考えることができる。よって、上記（１）と同様に、そのクラスターグループ数の値が健常者の値と比較して有意に低いと判定される場合に、唾液検体が炎症性腸疾患に由来するものと評価することができる。健常者の値と比較して有意に低いかどうかの判定は、集団検診などで多数の被験者がいる場合には、その被験者どうしの統計的差異に基づいて行ってもよく、一方、被験者が単独または少数の場合には、過去のデータから基準値（例えば所定のクラスターグループ数）を予め設定しておき、その基準値を閾値にして行なってもよい。なお、この方法の場合も、上記（１）と同様、安定に有意な評価を得るためには、上記データ群において１被験者あたり少なくともおよそ１０００データ以上の配列データ数を有していることが好ましい（後述の実施例参照）。
【００３１】
（３）基準となるデータ群との群間類似距離を指標にする方法
後述の実施例において実証されるように、唾液細菌叢の菌叢構造は、構成細菌の種類とその存在量の多寡によって特徴づけられ、その菌叢構造の特徴を、上記のように得られた塩基配列のデータ群の類似性として評価し得る。具体的には、上記データ群に基づいて、上述のクラスターグループ（例えば９６％塩基配列類似度を閾値として設定して）を求め、各クラスターグループに含まれる塩基配列のデータ数と、各クラスターグループに属する代表塩基配列どうしの類似度を指標にして、他の複数の被験者から取得されたデータ群との系統樹上での系統距離を算出する。この系統距離は、各被験者から取得されたデータ群との類似性をあらわす群間類似距離とも定義することができる。このような群間類似距離は、塩基配列のデータ群から構成される任意の複数群ついて、各群に属する塩基配列の配列どうしの類似度と配列数から、各群間の類似度を数値化する手法であるUniFrac解析(Lozupone C and Knight R: UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol 71: 8228-8235 (2005))などで行なうことができる。そして、この群間類似距離をもとに主座標分析で２次元散布図を得ると、例えば健常者群はすべてそのｘ軸上正の領域に布置され、一方患者群はすべてそのｘ軸上負の領域に布置されるといったように、健常者群と患者群とを明確に識別することが可能である。よって、任意の唾液検体について、基準となるデータ群との群間類似距離が、健常者群に近いか患者群に近いかで、唾液検体が炎症性腸疾患に由来するものか否かを評価することができる。健常者に近いかどうかの判定は、集団検診などで多数の被験者がいる場合には、その被験者どうしの統計的差異に基づいて行ってもよく、一方、被験者が単独または少数の場合には、過去のデータから基準値（例えば健常者群との所定の群間類似距離）を予め設定しておき、その基準値を閾値にして行なってもよい。また、既にある健常者群と患者群のデータ（過去のデータ）に被験者のデータを加えて上記のUniFrac解析とそれに基づく主座標分析を行い、得られた２次元散布図上で、その被験者が健常者群中またはそれに相対的に近いところにプロットされれば「健常」と判定し、患者群中またはそれに相対的に近いところにプロットされれば「疾患」と判定することなどで行ってもよい。更には、その判定は、過去のデータから群間類似距離に基づく主座標分析での散布図を求めておき、その散布図にあてはめて行なってもよい。
なお、この方法の場合、安定に有意な評価を得るためには、上記データ群において１被験者あたり少なくともおよそ５０データ以上の配列データ数を有していることが好ましく、およそ１００データ以上の配列データ数を有していることがより好ましく、およそ２０００データ以上配列データ数を有していることが最も好ましい（後述の実施例参照）。
【００３２】
（４）特定の門、属、または種に属する細菌の塩基配列のデータ数を指標にする方法
後述の実施例において実証されるように、炎症性腸疾患の患者の唾液細菌叢と健常者の唾液細菌叢とでは、特定の門、属、または種に属する細菌の存在量の多寡において、有意に相異している。また、上記のように得られた塩基配列のデータ群中に見出される配列データ数の大小は、唾液細菌叢を構成する細菌の存在量の多寡を反映しているものと考えられる。よって、任意の唾液検体について、上記データ群に基づいて、その特定の門、属、または種に属する細菌に帰属される塩基配列のデータ数を求め、その値が健常者の値と比較して有意に増減していると判定される場合に、唾液検体が炎症性腸疾患に由来するものと評価することができる。健常者の値と比較して有意に増減しているかどうかの判定は、集団検診などで多数の被験者がいる場合には、その被験者どうしの統計的差異に基づいて行ってもよく、一方、被験者が単独または少数の場合には、過去のデータから基準値（例えば所定の増減割合）を予め設定しておき、その基準値を閾値にして行なってもよい。その判定は門レベルで行ってもよく、属レベルで行なってもよく、種レベルで行ってもよく、それらの複数レベルを併用して評価してもよい。また、２種以上の門、属、及び／又は種についての判定を併用して評価してもよい。
【００３３】
（５）配列番号１〜３０の塩基配列のデータ数を指標にする方法
後述の実施例において実証されるように、炎症性腸疾患の患者の唾液細菌叢と健常者の唾液細菌叢とでは、配列番号１〜３０の塩基配列を有する細菌の存在量の多寡において、有意に相異している。また、上記のように得られた塩基配列のデータ群中に見出される配列データ数の大小は、唾液細菌叢を構成する細菌の存在量の多寡を反映しているものと考えられる。よって、任意の唾液検体について、上記のように得られた塩基配列のデータ群に基づいて、配列番号１〜３０の塩基配列を有する細菌に帰属される塩基配列のデータ数を求め、その値が健常者の値と比較して有意に増減していると判定される場合に、唾液検体が炎症性腸疾患に由来するものと評価することができる。健常者の値と比較して有意に増減しているかどうかの判定は、集団検診などで多数の被験者がいる場合には、その被験者どうしの統計的差異に基づいて行ってもよく、一方、被験者が単独または少数の場合には、過去のデータから基準値（例えば所定の増減割合）を予め設定しておき、その基準値を閾値にして行なってもよい。また、２種以上の該塩基配列についての判定を併用して評価してもよい。
【００３４】
本発明の適用対象となる炎症性腸疾患としては、クローン病や潰瘍性大腸炎などを例示できるが、これらと診断されるものに限られずいずれとも判別し得ない炎症性腸疾患や、その前段的な病態について、本発明の方法により好適に健常者との有意な差異を検出又は検査できる。
【００３５】
以上に説明した形態は、本発明の内容を説明するために挙げたものであり、本発明の範囲を限定するものではない。即ち、例えば本技術分野の当業者によって、本発明の範囲を逸脱しない前提の下で、種々の構成の置換や変更を加えることができることは当然であり、それもまた本発明の範囲とされる。
【実施例】
【００３６】
以下に実施例を挙げて本発明を更に説明するが、本発明は以下の実施例に限定されるものではない。
【００３７】
＜試験例１＞［塩基配列データの取得］
以下の方法により、被験者（健常者１０名およびクローン患者１０名）から唾液を採取し、その唾液検体中の細菌叢の１６ＳリボソームＲＮＡ遺伝子の塩基配列を無作為に決定して、被験者ごとに塩基配列のデータ群を得た。
【００３８】
[唾液の長期保存法]
唾液は、被験者から採取した後、次のようにして長期に保存することもできる。即ち1mlの唾液に1mlのリン酸バッファー（pH7.2、以下、PBS）1mlを加え混和後、2mlの40％グリセロール溶液を加える。この懸濁液を液体窒素で急冷後、グリセロールストックとして−80℃で保存する。本方法により、細菌叢に大きな分解等が起こらない少なくとも１年間の保存が可能である。
【００３９】
[唾液細菌叢DNAの調製]   
唾液検体中の細菌叢のDNAの調製は、以下のようにして行った。
【００４０】
(a) 唾液細菌叢DNAの調製  その１
唾液DNA採取用のOrageneキット（DNAgenotek社, Ontario, Canada）を用いて、次のようにして唾液細菌叢DNAを調製した。即ちキットに付属している1.9mlのOragene・DNA溶液（保存溶液）に約2mlの唾液（または、相当量の唾液のグリセロールストック）を加えて、5秒間ほど転倒混合する（そのまま、常温保管が可能）。その唾液サンプル溶液を加温（50℃、1時間）する。その後、500μlのサンプル溶液を新しいチューブ（1.5ml）に入れ、これに20μlのOragene DNA精製溶液（キットに付属）を加え、数秒間ボルテックスし攪拌し、氷上で10分間冷却する。その後、室温で遠心し（13,000r.p.m.×5分間）、上清を回収して新しいチューブ（1.5ml）に入れ、これに5μlのグリコーゲン溶液（キットに付属）を加え、ついで500μlのエタノール（99.5％）を加えて、室温で10分間インキュベート放置する。溶液を室温で遠心し（13,000r.p.m.×2分間）、得られたDNAペレットを真空乾燥し、ついで300μlのTE溶液（10mM Tris、1mM EDTA）（以下、TE）を加えて、唾液細菌叢DNAの溶液とする。
【００４１】
(b) 唾液細菌叢DNAの調製  その２
腸内細菌叢DNAの調製法（Morita H, Kuwahara T, Ohshima K, Sasamoto H, Itoh K, Hattori M, Hayashi T, and Takami H.: An improved isolation method for metagenomic analysis of the microbial flora of the human intestine. Microbes Environ. 22, p214-222 (2007)）に準じて、次のようにして唾液細菌叢DNAを調製した。即ち唾液のグリセロールストックを氷上にてゆっくり融解する（新鮮な唾液の場合は1mlの唾液に数mlのPBSを混合して唾液のPBS溶液を調製する。）。融解した溶液（または唾液のPBS溶液）を、孔径100μmフィルター（BD社製Falconセルストレーナー）でろ過し、不ろ過物をさらに2〜3mlのPBSで数回洗浄ろ過する。ろ液を混合し、遠心（5,000 r.p.m.×10分間）し、ペレットをPBSで1回洗浄し、TE10溶液（10 mM Tris、10 mM EDTA）（以下、TE10）でさらに2回洗浄し、遠心（5,000 r.p.m.×10分間）し、菌体ペレットを得る。菌体ペレットに3mlのTE10を加えて懸濁し、Lysozyme（Sigma社）を最終濃度が15mg/ml-cell suspensionになるように加え、軽く振とうする（37℃×1時間）。ついで、Achromopeptidase（和光純薬工業社の精製品）を最終濃度が2,000units/ml-cell suspensionになるように加え、軽く振とうする（37℃×30分間）。溶液に10% SDS (pH 7.2)を最終濃度が1%になるように加え、ついで、Proteinase K（Merck社）を最終濃度1mg/ml-lysate）になるように加え、軽く振とうする（55℃×1時間）。TE10を加えて全体量を10mlにし、10mlのフェノール/クロロホルム/イソアミルアルコール（25:24:1 [vol/vol/vol]）を加え、よく混合し、遠心する（5,000r.p.m.×10分間）。上清を回収して、再度当量のフェノール/クロロホルム/イソアミルアルコールを加え、よく混合し、遠心する（5,000 r.p.m.×10分間）。遠心後、回収した上清に1/10倍量の3M 酢酸ナトリウム（pH 5.2）を加えて混合し、さらに2倍量のエタノール（99.9％）を加え、氷上にて5分間静置する。遠心（5000 r.p.m.×10分間）してDNAペレットを得、そのDNAペレットを15mlの75％ エタノールでリンスする。このリンス処理はさらに1回行う。真空乾燥したDNAペレットに600μlのTEを加えて、このDNA溶液に6μlのRNase A （10 mg/ml、Novagen社）を加え、加温する（37℃×1時間）。ついで、300μlの1.6M NaClと300μlの26(wt/vol)% PEG#6000（ナカライテスク社）を加え、氷上で1時間静置する。溶液を4℃で遠心（12,000r.p.m.×30分間）してDNAペレットを得る。得られたDNAペレットに1mlの75% ethanolを加えてリンスし、DNAペレットを真空乾燥後、300μlのTEを加えて、RNAの混入がない唾液細菌叢DNAの溶液とする。本方法により、他の既知方法に較べて高収率で高分子量の唾液細菌叢DNAを調製することができる。
【００４２】
なお、上記で示した唾液のDNA調製法（「その１」と「その２」）の違いは、以下の通り、唾液の採取環境に応じて使い分けることができる。即ち、上記「その１」のDNA調製法を用いるのは、被験者が自らの唾液を解析者があらかじめ提供したキットに備え付けの溶液に吐き出し、その溶液を室温にて解析場所まで搬送（宅配等）する、または、病院等で医師または解析者のもとで唾液を採取し、その後実験室でDNA調製を行う、などのケースが挙げられる。この場合、唾液そのものを保存できないが、多くの被験者からの唾液採取が容易である。一方、上記「その２」のDNA調製法を用いるのは、被験者が所定のチューブに吐き出した唾液そのものを解析場所まで48時間以内に搬送（冷蔵宅配等）する、または、病院等で医師または解析者のもとで唾液を採取し、その後実験室でDNA調製を行う、などのケースが挙げられる。この場合、DNA調製までの時間的な制約があるが、解析者が唾液そのものをまとめて保存でき、また、唾液の一部をここで述べた細菌叢解析に供して、残した一部をその他の解析（たとえば、免疫系たんぱく質の検出等）にも使用できる。
【００４３】
[１６ＳリボソームＲＮＡ遺伝子のV1-V2領域のPCR増幅]
唾液細菌叢DNAの溶液中の二本鎖DNA濃度を、Qubit 2.0 Fluorometer （Invitrogen社）を用いて測定した。その測定値に基づいて40ngのDNAを鋳型として、ユニバーサルプライマーセット（フォワードプライマー27Fmod-454A（配列番号３１）とリバースプライマー338R-454B（配列番号３２））を用いて、１６ＳリボソームＲＮＡ遺伝子（以下、16S遺伝子）のV1-V2領域（図１）をPCR増幅した。PCRはタカラバイオ社製の「TaKaRa Ex Taq」（登録商標）を用いて、各プライマーを0.2μmolを含む反応液を作成し、94℃で2分間のプレヒーティングを行った後、変性、アニーリング、伸長をそれぞれ94℃×30秒間、55℃×30秒間、72℃×60秒間で行い25サイクル繰り返した。サイクル終了後、増幅DNA鎖を完全に伸長させるために72℃×14分間の処理を行った。
【００４４】
下記にはフォワードプライマー27Fmod-454Aの配列の構造を示す。このフォワードプライマーは、後述するシークエンサーでの配列決定に必要なアダプターA配列（大文字で表記）を5’末端側に含み、各検体に固有の10塩基のバーコード配列（**MID**で表記）をはさんで（このバーコード配列はサンプル間の識別に利用するもので、同時にシークエンサーに供するサンプル数に対応した任意に設計した塩基配列である。）、すべての真正細菌の16S遺伝子にアニーリングするユニバーサルプライマー配列27Fmod（小文字で表記）を3’末端側に含む。
・27Fmod-454Aの配列（配列番号３１）
5’-CCATCTCATCCCTGCGTGTCTCCGACTCAG**MID**agrgtttgatymtggctcag-3’
（ただし、mはA又はC、rはG又はA、yはT又はCを意味する。）
【００４５】
なお、図２には上記バーコード配列の役割を説明する。例えば、４人の被験者から採取した４検体を同時にシークエンサーに供する場合、バーコード配列を４通りの異なった配列にして、各バーコード配列をもつPCRプライマーで各サンプルをPCR増幅する。これらのPCR産物（4検体分）を同量混合して同時にシークエンサー（この場合、１回の稼働で10万データを生産するGS Junior Systemシークエンサーを用いたとして）に供すると、平均2.5万データ／サンプルの割合で計10万データが得られる。しかし、この段階では、各データがどのサンプルに由来しているのか分からないので、この10万データをバーコード配列を指標に各サンプルに振るい分けを行う。これにより検体ごとの配列データを得ることができる（図２）。また、例えば、２０人の被験者から採取した２０検体を同様に同時解析する場合は、２０通りの異なったバーコード配列をもった27Fmod-454Aを作り、それぞれを各検体に対してPCR増幅すればよい。これらを混合してシークエンサーに供すると（GS Junior Systemシークエンサーを用いたとして）、１回の稼働で１検体当たり5000データが得られる。更に、１稼働で100万データを得ることができるGS FLX+ Systemシークエンサーを利用すると、100検体に対応する100通りのバーコード配列を用いることで、１回の稼働で1万データ／検体の配列データを得ることができる。
【００４６】
一方、下記にはリバースプライマー338R-454Bの配列の構造を示す。このリバースプライマーは、後述するシークエンサーでの配列決定に必要なアダプターB配列（大文字で表記）を5’末端側に含み、すべての真正細菌の16S遺伝子にアニーリングするユニバーサルプライマー配列338R（小文字で表記）を3’末端側に含む。
・338R-454Bの配列（配列番号３２）
5’-CCTATCCCCTGTGTGCCTTGGCAGTCTCAGtgctgcctcccgtaggagt-3’
【００４７】
上記のユニバーサルプライマーセットを用いたPCRにより、唾液細菌叢を構成する種々の細菌種の16S遺伝子のV1-V2領域を含むDNA（約400塩基）が増幅され、それらの混合物をそのPCR産物DNAとして得ることができる。
【００４８】
[PCR産物の精製およびシークエンス用サンプルの調製]
各々の唾液細菌叢DNAから得られたPCR産物DNA（その細菌叢を構成する種々の細菌種の16S遺伝子のV1-V2領域を含むDNAの混合物）を、AMPure XP kit（BECKMAN COULTER社）にて処理して、過剰のプラーマーや基質のヌクレオチド等を除去し、精製した。精製DNAは10μlのTEで溶出・回収した。各検体から回収された精製PCR産物DNAの二本鎖DNA量をQuant-iT PicoGreen dsDNA Assay Kit（Invitrogen社）を用いてそれぞれ定量後、その定量値を得る。ついで、PCR産物DNAをその定量値に基づいて同じDNA量になるように厳密に混合する。この混合溶液中の二本鎖DNA量を、再度Quant-iT PicoGreen dsDNA Assay Kitを用いて定量して、シークエンサーに供するDNA量を調整する。これを以下のシークエンスに用いるシークエンス用サンプルとした。
【００４９】
[16S遺伝子のシークエンシングと配列データの精度評価]
上記シークエンス用サンプルを、解析したい検体数に応じて、ロシュ社製GS FLX+ SystemまたはGS Junior Systemシークエンサーに供しシークエンスを行った。シークエンスの条件・工程等はメーカー所定のプロトコールに従った。なお、このシークエンサーでは、上記で調製したPCR産物DNAの１分子を１つのビーズに固定して、ついで、水（シークエンス用鋳型DNAの増幅のためのPCRプライマー、基質ヌクレオチド、DNA合成酵素を含む）と油のエマルジョン中に独立して形成された微小水滴の１つ１つに１つ１つのビーズを捕獲して、その中でPCRを行ってシークエンス用鋳型DNAを増幅して調製するようになっている。よって、この増幅した鋳型DNAが固定された各ビーズをタイタープレート上に区画した後に、その区画位置上でシークエンス反応のシグナルを読み取りことによって、上記シークエンス用サンプル中に含まれるPCR産物DNA（その細菌叢を構成する種々の細菌種の16S遺伝子のV1-V2領域を含むDNAの混合物）の塩基配列を無作為に決定することができる。また、フォワードプライマー27Fmod-454A中の上記バーコード配列を、各被験者に由来する検体ごとに特徴的な任意の配列にしておけば、GS FLX+ Systemシークエンサーを用いて少なくとも100種類の細菌叢サンプルを、小規模タイプのGS Junior Systemを用いて少なくとも10種類の細菌叢サンプルを同時解析でき、１人の被験者由来のサンプルにつき2,000〜10,000の16S遺伝子の配列データを、およそ10〜20時間で決定することができる。
【００５０】
得られた粗配列データ（〜500塩基／データ）については、配列データに含まれるサンプル固有のバーコード配列に基づき、各配列をそれぞれの固有のサンプルに分配した。その後、以下に記した精度評価工程を通して、それらの評価条件（１）〜（３）を満たさない低精度の配列データを除去することにより、高精度データを抽出した。
【００５１】
（1）配列データの両末端配列としてユニバーサルプライマー配列(27Fmodおよび338R)との配列類似度が80％以下である配列データを除去した。この工程は相同性検索プログラムのBLASTを用いて行い、両末端にユニバーサルプライマーの配列を持たない不完全な配列データを除去した。
【００５２】
（2）シークエンサーに付属のクオリティプログラムを用いて、配列決定した塩基配列の平均クオリティ値が25以下の配列データを除去した。
【００５３】
（3）上記で選択された配列データを細菌の16S配列データの公表データベースであるRDP(http://rdp.cme.msu.edu/)とCORE(http://microbiome.osu.edu./)に登録されている16S配列、及び、NCBI(http://www.ncbi.nlm.nih.gov/)とHMP(http://commonfund.nih.gov/hmp/)に登録されている細菌ゲノム配列中の16S配列とのアラインメント塩基長が90%以上を有する配列データを選択した。90%以下のアラインメント塩基長を有する配列データにはきわめてシークエンスエラーの多いデータやキメラデータが含まれる確率が高いため、これらをこの工程で除去した。
【００５４】
以上の精度評価により、各被験者由来のサンプルにおいて、〜10,000の全粗配列データのうち、60〜70%のデータが高精度配列データとして選択された。
【００５５】
＜試験例２＞［Operational Taxonomic Unit (OTU) 解析］
試験例１のようにして取得した高精度配列データからランダムに抽出した3,000データ（各被験者当り）を、クラスタリング（類似度96%の閾値）によるOperational Taxonomic Unit解析（以下、OTU解析）に供した。図３にはOTU解析の概略説明図を示す。図３に示すように、OTU解析においては、配列データの類似度を基準にして各配列データをグループ化する操作を行う。ここでは96%以上の配列類似度を互いに有する配列データのクラスターグループ（以下、OTU）を検出している。なお、配列データのクラスタリングはフリーウェアUSEARCH（http://drive5.com/usearch/usearch3.0.html）などを用いて行うことができる。各OTUは同じ種の細菌に由来すると推測できる。よって、クラスタリングによって得られるOTUの総数（OTU数）は、検出可能な範囲において、その細菌叢を構成する細菌種の数と等価と考えることができる。また、各OTU中に含まれる配列データ数からは、配列データ数全体中の各OTUの割合、つまり菌種組成比を求めることができる。さらに、各OTUの代表配列データについて上記した16S及び細菌ゲノムのデータベースへの相同性検索を行うことにより、もっとも高い配列類似度を有する既知菌種へ帰属、つまり、OTUの菌種を特定できる。
【００５６】
図４には、１０名の健常者群および１０名のクローン病（CD）患者群について、各被験者に由来する塩基配列のデータ群のOTU解析によって得られたOTU数を、健常者群とCD患者群とで比較した結果を示す。図４に示されるように、健常者群の平均OTU数（１０名の平均）とCD患者群の平均OTU数（１０名の平均）はそれぞれ132.3(標準偏差29.97)と103.2(標準偏差25.25)であり、CD患者のOTU数は健常者のそれよりも統計学的に有意に低いことが示された（t-testによるp値＜0.05）。
【００５７】
以上の結果から、CD患者では唾液細菌叢を構成する菌種数（OTU数）が健常者よりも有意に少なくなっていることが明らかとなった。よってその菌種数（OTU数）を指標にして、任意の唾液検体についてクローン病（炎症性腸疾患）と評価される唾液検体を検出できると考えられた。
【００５８】
＜試験例３＞［OTU解析に供される配列データ数の影響］
試験例２では高精度配列データからランダムに抽出した3,000データ（各被験者当り）を用いたが、用いる配列データ数を300, 500, 700, 1,000, 1,200, 1,500と代えて同様に解析を行い、健常者群の平均OTU数とCD患者群の平均OTU数との有意差（t-testによるp値）を求めた。またそれを３回試行した。図５Ａには試行回ごとのp値を示し、図５Ｂは各データ数における平均p値を示す。
【００５９】
図５Ａ，Ｂに示されるように、OTU解析において安定した結果（p値＜0.05）を得るためには、1,000データ程度以上の配列データを用いてOTU解析を行うことが必要であった。
【００６０】
＜試験例４＞［UniFrac解析およびそれに基づく主座標分析］
１０名の健常者群および１０名のクローン病（CD）患者群について、試験例１のようにして取得した高精度配列データからランダムに抽出した3,000データ（各被験者当り）を、UniFrac解析に供して各被験者間の類似度を求め、その類似度に基づく主座標分析を行った。ここでUniFrac解析は、塩基配列のデータ群から構成される任意の複数群ついて、各群に属する塩基配列の配列どうしの類似度と配列数から、各群間の類似度を数値化する手法である(Lozupone C and Knight R: UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol 71: 8228-8235 (2005))。また、主座標分析は、対象についての任意の基準の類似度を元にして、その対象をｎ次元座標上に布置する手法である。なお、UniFrac解析を行うにはフリーウェア（http://bmf.colorado.edu/unifrac/）などが利用可能である。また、主座標分析についても市販のプログラムなどが利用可能である。
【００６１】
この試験例では、各被験者に由来する塩基配列のデータ群のOTU解析を試験例２のようにして行い、得られたOTUについて、各OTUに属する代表塩基配列どうしの類似度と各OTUに含まれる塩基配列のデータ数とに基づいてUniFrac解析を行い、１０名の健常者群および１０名のクローン病（CD）患者群について、各被験者に由来する塩基配列のデータ群の類似度を、系統樹上での系統距離(UniFrac Distance)（以下、群間類似距離）として算出した。そして、その群間類似距離(UniFrac Distance)に基づいた主座標分析を行い、２次元座標上への主座標1および主座標2の値をもとに、各被験者間の菌叢構造の類似度を２次元散布図で表した。その結果を図６に示す。
【００６２】
図６に示されるように、健常者群が布置される座標領域とCD患者群が布置される座標領域とでは、それぞれ明確に異なる座標領域が形成された。このことから、CD患者の唾液細菌叢と健常者の唾液細菌叢とでは、少なくともそれらを構成する細菌の種類とその存在量の多寡の菌叢構造において、有意に相異していることが明らかとなった。また、各被験者間の菌叢構造の類似度を、上記群間類似距離またはそれに基づく主座標分析により評価でき、健常者群とCD患者群を明確に識別できることが明らかとなった。
【００６３】
以上の結果から、健常者群に対する上記群間類似距離を指標にして、任意の唾液検体からクローン病（炎症性腸疾患）と評価される唾液検体を検出できると考えられた。
【００６４】
＜試験例５＞［UniFrac解析および主座標分析に供される配列データ数の影響］
試験例４では高精度配列データからランダムに抽出した3,000データ（各被験者当り）を用いたが、用いる配列データ数を50, 100, 2,000と代えて同様にUniFrac解析を行い、それに基づいて各群間の菌叢構造類似度を２次元散布図で表した。その結果を図７〜９に示す。
【００６５】
図７〜９に示されるように、データ数が多いほどより明確に健常者群とCD患者群とを識別できたが、被験者当たりのデータ数が50のときでも、主座標分析で形成された２次元散布図上、健常者群はすべてそのｘ軸上正の領域に布置され、一方CD患者群はすべてそのｘ軸上負の領域に布置されて、健常者群とCD患者群とを明確に識別することが可能であった。
【００６６】
＜試験例６＞［健常者とCD患者の間で有意に増減するOTUの探索と菌種の特定］
１０名の健常者群および１０名のクローン病（CD）患者群について、試験例１のようにして取得した高精度配列データからランダムに抽出した3,000データ（各被験者当り）を２０名分まとめて、その60,000データをOTU解析に供した。クラスタリングのための閾値としては、上記した個別細菌叢のOTU解析と同じ96%の配列類似度を閾値として設定した。得られたOTUのうち、健常者群とCD患者群との群間での配列データ数の増減においてt-testによるp値が0.05以下（門レベルと属レベル）及び0.01以下（種レベル）を示すOTUを検出した。また、これらのOTUの門、属、種レベルでの菌種の帰属を、16S遺伝子配列のデータベースであるRDPとCORE、及びゲノム配列のデータベースであるNCBIとHMP（以下まとめてNCBI_genomeと表記）への相同性検索により行った。
【００６７】
その結果、下記の解析結果が得られた。
【００６８】
（１）門レベルでの有意に増減する菌種
OTU解析から、20名の唾液細菌叢に計12門の細菌種が検出された。その12門のうち、Firmicutes、Bacteroidetes、Actinobacteria、Proteobacteriaの4門に属する菌種の配列データ数が、健常者群とCD患者群の群間で有意に増減していた（t-testによるp値が0.05以下）。そして、Firmicutes門とBacteroidetes門に属する菌種の配列データ数において、CD患者群では健常者群に比べて有意な増加を示し、一方で、Actinobacteria門とProteobacteria門に属する菌種の配列データ数において、CD患者群では健常者群に比べて有意な減少を示した。
【００６９】
その結果を表１にまとめて示す。
【表１】
0005863035.tif 000002
【００７０】
（２）属レベルでの有意に増減する菌種
OTU解析から、20名の唾液細菌叢に計100属の細菌種が検出された。その100属のうち、Prevotella、Streptococcus、Veillonella、Atopobium、Megasphaera、Solobacterium、Actinomyces、Lachnospiraceae、Selenomonas、Rothia、Haemophilus、Neisseria、Gemella、Porphyromonas、Corynebacterium、Capnocytophaga、Bergeyella、Aggregatibacter、Lautropia、Paludibacter、Fusobacterium、Tannerella、Propionibacteriumの23属に属する菌種の配列データ数が、健常者群とCD患者群の群間で有意に増減していた（t-testによるp値が0.05以下）。この23属のうち、9属（Prevotella属、Streptococcus属、Veillonella属、Atopobium属、Megasphaera属、Solobacterium属、Actinomyces属、Uncultured Lachnospiraceae属、Selenomonas属）に属する菌種の配列データ数において、CD患者群では健常者群に比べて有意な増加を示し、一方で、14属（Rothia属、Haemophilus属、Neisseria属、Gemella属、Porphyromonas属、Corynebacterium属、Capnocytophaga属、Bergeyella属、Aggregatibacter属、Lautropia属、Paludibacter属、Fusobacterium属、Tannerella属、Propionibacterium属に属する菌種の配列データ数において、CD患者群では健常者群に比べて有意な減少を示した。
【００７１】
その結果を表２にまとめて示す。
【表２】
0005863035.tif 000003
【００７２】
（３）種レベルでの有意に増減する菌種
OTU解析から、20名の唾液細菌叢に計547種の細菌種が検出された。その547種のうち、Veillonella atypical、Granulicatella adiacens、Atopobium parvulum、Streptococcus sp. EO2001-02、Veillonella dispar、Streptococcus mitis、Prevotella DO039、Megasphaera micronuciformis、Prevotella salivae、Solobacterium moorei、Prevotella melaninogenica、Actinomyces graevenitzii、Neisseria subflava、Haemophilus nbw161b08c1、Haemophilus sp. CCUG 32367、Rothia aeria、Lautropia mirabilis、Rothia mucilaginosa、Campylobacter gracilis、Corynebacterium matruchotii、Bergeyella 602D02、Capnocytophaga gingivalis、Corynebacterium durum、Gemella haemolysans、Porphyromonas CW034、Streptococcus VG051、Fusobacterium nucleatum、Prevotella IK062、Veillonella dispar、Haemophilus parainfluenzaeの30菌種の配列データ数が、健常者群とCD患者群の群間で有意に増減していた（t-testによるp値が0.01以下）。この30菌種のうち、12菌種（Veillonella atypical、Granulicatella adiacens、Atopobium parvulum、Streptococcus sp. EO2001-02、Veillonella dispar、Streptococcus mitis、Prevotella DO039、Megasphaera micronuciformis、Prevotella salivae、Solobacterium moorei、Prevotella melaninogenica、Actinomyces graevenitzii）の配列データ数において、CD患者群では健常者群に比べて有意な増加を示し、一方で、18菌種（Neisseria subflava、Haemophilus nbw161b08c1、Haemophilus sp. CCUG 32367、Rothia aeria、Lautropia mirabilis、Rothia mucilaginosa、Campylobacter gracilis、Corynebacterium matruchotii、Bergeyella 602D02、Capnocytophaga gingivalis、Corynebacterium durum、Gemella haemolysans、Porphyromonas CW034、Streptococcus VG051、Fusobacterium nucleatum、Prevotella IK062、Veillonella dispar、Haemophilus parainfluenzae）の配列データ数において、CD患者群では健常者群に比べて有意な減少を示した。
【００７３】
その結果を表３にまとめて示す。
【表３】
0005863035.tif 000004
【００７４】
なお、上記解析における配列データ数は、被験者から採取した唾液検体中の細菌叢の16S遺伝子の塩基配列を無作為に決定することによって取得された塩基配列のデータ群中に見出される該当配列の総数（表１〜３中では各群１０名の平均値として表わされている。）であり、その大小は、唾液細菌叢を構成する細菌の存在量の多寡を反映しているものと考えられた。
【００７５】
また、上記解析においては、各OTUに属する代表塩基配列と既知菌種の16S配列との相同性検索により、各OTUがいずれの門、属、菌種に属するかを特定した。表３には、上記（３）の解析で検出された30 OTUの代表塩基配列（配列番号１〜３０）と、それぞれが帰属する既知菌種の16S配列との間の配列類似度とを合わせて示す。表３に示すように、配列類似度の範囲は96.8〜100%となり、そのうち27 OTUでは既知菌種と97%以上の類似度を示し、系統的に同一菌種であると結論づけられた。残りの3 OTUでは配列類似度は97%以下であり、これらのOTUは、それぞれ既知菌種に最も近縁の菌種に帰属するものと考えられた。
【００７６】
以上の結果から、CD患者の唾液細菌叢と健常者の唾液細菌叢とでは、特定の門、属、または種に属する細菌の存在量の多寡において、有意に相異していることが明らかとなった。よって、被験者から採取した唾液検体中の細菌叢の16S遺伝子の塩基配列を無作為に決定することによって取得された塩基配列のデータ群中に見出される、その特定の門、属、または種に属する細菌の塩基配列のデータ数、もしくは上記配列番号１〜３０の塩基配列のデータ数を指標にして、任意の唾液検体についてクローン病（炎症性腸疾患）と評価される唾液検体を検出できると考えられた。
【００７７】
「配列表フリーテキスト」
配列番号３１：１６ＳリボソームＲＮＡ遺伝子のV1-V2領域のPCR増幅のためのフォワードプライマー27Fmod-454A
配列番号３２：１６ＳリボソームＲＮＡ遺伝子のV1-V2領域のPCR増幅のためのリバースプライマー338R-454B
【図１】
0005863035.tif 000005
【図２】
0005863035.tif 000006
【図３】
0005863035.tif 000007
【図４】
0005863035.tif 000008
【図５】
0005863035.tif 000009
【図６】
0005863035.tif 000010
【図７】
0005863035.tif 000011
【図８】
0005863035.tif 000012
【図９】
0005863035.tif 000013
【配列表】
0005863035000001.app