【実施例】
【0108】
最近の証拠によると、腸内微生物群の菌群のバランスがくずれるのは、肥満および糖尿病などのヒト疾患に肝心な作用(1、2)を果たすことを表明する。疾患/健康の表現型に影響を与える腸内細菌叢の特定メンバーは、疾患診断の有力なツール(3、4)だけでなく、標的として薬物(5)、糞便微生物移植(6)および飲食(7)などの多種方法で疾患の緩和/治療に用いることができる。しかしながら、腸内細菌叢自体の複雑さおよび個体間の差異およびそれと宿主と飲食の相互作用(8)のために、最適用なヒト健康を実現して腸内微生物群に対する正確な操作は、ゲノムおよび分子レベルにより深く理解する必要がある。
【0109】
我々の前の食事介入研究の中で、全穀類、伝統的な中国薬膳およびプレバイオティクスから構成される消化しにくいが発酵可能な炭水化物の食事(WTP食事)は、腸内細菌叢を顕著改善するだけでなく、プラダー・ウィリー症候群に罹患する遺伝性肥満児童の生物学的臨床パラメータおよび炎症状況(9)を改善する。特定の善玉菌ビフィドバクテリウム・シュードカテニュレイタムは、介入した後に顕著に濃化され、それは他の潜在的な有害菌種と負の相関であり、かつ宿主臨床パラメータの改善と正の相関である(9)。
【0110】
該グループからの一つの児童は105日の介入を完了する。他の生物学的臨床パラメータが改善され、初期の体重減少が25.8kgを超え、同時に腸内微生物コロニーが顕著に転換し、例えば、大便桿菌(Faecalibacterium)、乳酸桿菌(lactobacillus)、ビフィドバクテリウム属(Bifidobacterium)が多くなる。また、介入した後のビフィドバクテリウム・シュードカテニュレイタムはビフィドバクテリウム属における最も優位な菌種であることを発見する。興味深いことは、我々は105日目に彼の糞便サンプルから5種類菌株(10)を分離し、これらの菌株は炭水化物に富む介入に対して差異応答を有する。腸内生態系における細菌ニッチアダプテーションを参与する遺伝的特性および宿主と食物との相互作用(11)を知るために、我々はビフィドバクテリウム・シュードカテニュレイタムに対して比較ゲノミクス(comparative genomics)の研究を行う。
【0111】
結果および検討
【0112】
改善された生物学的臨床パラメータおよび転換された腸内細菌叢
【0113】
介入期間に該肥満児童の生物学的臨床変数が改善される(
図1および
図7)。体重が141.1kgから114.3kgに低下し、かつ血漿グルコースおよび脂質内恒常性の両方が正常範囲に改善する。介入した後、二種全身性炎症マーカーC-反応性タンパク質(C-reactive protein、CRP)および血清アミロイドタンパク質A(SerumamyloidAprotein、SAA)が低下する。アディポネクチンは2.17μg/mlから5.39μg/mlに向上し、レプチンは63.82ng/mlから34.47ng/mlに低下し、「リスク」表現型(12)の緩和を示す。また、血液中の細菌抗原に担持された代替マーカーとしてのリポ多糖結合タンパク質(LBP)(13)が低下した。
【0114】
我々はメタゲノムシークエンスにより7つの時点(0、15、30、45、60、75および105日)に各糞便のサンプル25.2±4.8百万(平均値±s.d.)個の高品質ペアエンドリード(paired-end reads)を取得する。介入期間の腸内細菌叢の組成が転換し(
図2A)、かつ三つの介入ステージ(ステージI:0日目から60日目までの基礎介入、ステージII:60日目から75日目までの基礎介入+100g以上の3号処方、ステージIII:75日目から105日目までの100g以下の1号処方+100g以上の3号処方(詳しくは、方法参照)に応答して、異なるパターンを示した。該肥満児童は、異なるステージにおける排便頻度が類似し(1日当たり3〜4回)、下痢が発生しない。介入期間にコロニーの多様性が低下し、それはグループの変化と一致する(9)。ベースラインにある場合、ルミノコッカス属(Ruminococcus)およびブブラウティア属(Blautia)は二つの最も豊富な属であり、それぞれ26.95%および18.41%を占める。同時に、バクテロイデス属(Bacteroides)およびプレボテイラ属(Prevotella)は低い豊度を有し、該コロニーが腸型3(14)に属する可能性があることを示す。介入した後に、ルミノコッカス属およびブブラウティア属は低い豊度まで低下し、バクテロイデス属およびプレボテイラ属(Prevotella)はほとんど介入の影響を受けない。ベースラインと比べて、抗炎症菌種(15、16)および有益な(5)共生細菌の大便桿菌はステージIで増加し、15日目の豊度は41.95%に達する。前の15日の大便桿菌の急激な増加は炎症の緩和に役立つ可能性があり、該ステージで、CRPおよびSAAがそれぞれ33.37%および55.85%低下する。より多くのオリゴ糖を提供する場合、大便桿菌は、ステージIIおよびステージIIIで、低い豊度まで低下する。一方、オリゴ糖(17)を代謝するように好適に配備されたラクトバチルス属は、ベースラインとステージIの期間に低い豊度があるが、ステージIIから最も優位属の1つようになる。複数の研究は既にインスリン抵抗に対してラクトバチルス属菌株の有益な作用(18、19)を報告した。ラクトバチルス属の急激な増加は、ステージIIからのインスリン感受性の改善に役立ち、例えば、該時期の期間OGTTインスリンAUC(Area under the curve、曲線下の面積)の顕著な低下および安定したOGTTグルコースAUCに示されるように、細胞内および細胞外のルートによって複数種類の炭水化物(20)を代謝するビフィズス菌の豊度は15日目から顕著になり、かつ介入期間に亘って依然として優位である:該児童が基礎介入を受け始まる時の27.47%から75日目の65.53%まで上昇し、次に少ない複合食物繊維を含む1号処方を提供する時に、36.41%まで低下する。これらの結果は、ビフィズス菌種群が三つのステージで提供された炭水化物源の変化に応答し、かつ各ステージで、それが優位な生態的地位を占めることを示す。我々は、ビフィズス菌菌株の肥満に対する既知の有益な作用(21、22)、およびそれが従来のグループ研究の相互作用ネットワークにおける重要な作用(9)、およびそれが介入期間の持続的な高い豊度から考えると、これらが体重、BMI、胴囲および腰回りを持続的に低下させることに役立つ可能性があると推測する。したがって、我々は該属をより深く分析する。
【0115】
介入後のサンプルにおいて、合計9種類のビフィドバクテリウム属菌種を同定した(
図8)。ここで、ビフィドバクテリウム・シュードカテニュレイタムが最も優位を占め、105日目の腸内細菌叢全体における豊度は29.36%である。ビフィドバクテリウム・ロンガム、ビフィドバクテリウム・ブレーベおよびビフィドバクテリウム・アドレセンティスはそれぞれ9.94%、7.61%および3.75%を占め、他の5つの菌種の豊度は1%未満である。
【0116】
ビフィドバクテリウム・シュードカテニュレイタム菌株の介入に対する差異応答
【0117】
ビフィドバクテリウム・シュードカテニュレイタム群を詳細に研究するために、我々は、105日目に我々の対象から収集された糞便サンプルから5種類のビフィドバクテリウム・シュードカテニュレイタム菌株(以下、菌株PERFECT-2017-0001、PERFECT-2017-0003、PERFECT-2017-0004およびC95(23)と定義する)を分離して完全なシークエンス(10)を行う。メタゲノムデータを完全なゲノムデータと比べることにより、Sigma(24)を利用して該5種類の菌株の各時点の豊度の変化を同定する(
図2B)。介入前に、全ての菌株は低い豊度(最大値=0.5%で、ビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0001)を有する。前者は試験期間全体に低い豊度を維持し、後者はステージIIおよびステージIIIの期間に小さな増加を示すため、ビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0003とビフィドバクテリウム・シュードカテニュレイタムC95の豊度は食事の介入に応答しなさそうである。比べてみると、ビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0001、ビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0004およびビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0002の菌株は食物介入に応答し、それらの豊度変化はビフィドバクテリウム属に対して観察されたものと一致し、ビフィドバクテリウム・シュードカテニュレイタム菌株は主にこれらの変化をもたらすことを示す。実際に、介入後、これらの三種類の菌株は顕著に増加し、特にビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0002である。
【0118】
従来の研究(9)における発見と一致し、生物学的臨床変数の改善はビフィドバクテリウム・シュードカテニュレイタム菌株の増加に関連する(
図3)。ビフィドバクテリウム・シュードカテニュレイタムC95は人体測定マーカー(体重、BMI、胴囲および腰回りを含む)と負の相関である。また、ビフィドバクテリウム・シュードカテニュレイタムPERFECT
-2017
-0002も、炎症マーカーの改善(レプチンの減少とアディポネクチンの増加を含む)に関する。
【0119】
ビフィドバクテリウム・シュードカテニュレイタムのパンゲノム解析
【0120】
我々の5個の完全なビフィドバクテリウム・シュードカテニュレイタムゲノムおよび利用可能な共通データ(5個のビフィドバクテリウム・シュードカテニュレイタムゲノム略図とビフィドバクテリウム・シュードカテニュレイタムJCM1200
Tの完全なゲノム(25)を含む)に基づいて、パンゲノム解析に11個ビフィドバクテリウム・シュードカテニュレイタムゲノムを含む。パンゲノム曲線は、PERFECT-2017-0001の生物学的臨床変数が漸近傾向を呈し、前の6個反復における平均成長速度はゲノムあたり100個遺伝子で、その後、これよりもっと低い速度に低下することを示す(
図9A)。曲線は最終的に2482個遺伝子に達する。これは、さらに追加のゲノムを併入することはパンゲノムの大きさのわずかな増加をもたらす可能性があることを示す。コアゲノム曲線は前の6個反復におけるより明らかな漸近傾向およびより明らかな低下を示す(
図9B)。曲線は最終的に1427遺伝子に達する。パンゲノムとコアゲノムの傾向は、ビフィドバクテリウム・シュードカテニュレイタムは、閉めたパンゲノムを示し、6個ゲノムはほとんどビフィドバクテリウム・シュードカテニュレイタムの遺伝子特徴を説明できる。これらの結果に基づき、ゲノム略図でファジィをもたらすことを回避するために、我々は、後の分析において6個の完全なゲノムを使用してビフィドバクテリウム・シュードカテニュレイタムのゲノム特徴を探索する。
【0121】
ビフィドバクテリウム・シュードカテニュレイタムの一般的な特徴
【0122】
我々の5種類の菌株およびビフィドバクテリウム・シュードカテニュレイタムJCM1200
Tのゲノムは、平均2355、185bpおよび56.63G+C%を示し、これはビフィドバクテリウム属のG+C%の範囲と一致する(26)。ビフィドバクテリウム・シュードカテニュレイタムのゲノム中において5個または6個のrRNAオペロン遺伝子座を同定し、かつ各ゲノムに5SrRNArRNA遺伝子の一つの追加コピー(表1)を有する。また、16SrRNA遺伝子の異質性は、ビフィドバクテリウム・シュードカテニュレイタムJCM1200Tのゲノム以外のすべてのゲノムに存在する(表S1)。各ビフィドバクテリウム・シュードカテニュレイタムゲノム中に平均的に54個のtRNA遺伝子を含む。
【表1】
【表S1】
【0123】
各ゲノムに平均的に1871個のオープンリーディングフレームがある(Open Reading Frame、ORF)と予測し、ここで、検出された80%ORFは、NCBIデータベース上のBLASTに基づいて、コンピュータ予測(Silico prediction)によって機能的な割り当てを行って、残りの20%の予測は仮説タンパク質である(表1)。COG(オルソログクラスタグループ)(27)に応じてオルソログ遺伝子を同定することは、ビフィドバクテリウム・シュードカテニュレイタムゲノムにおける大部分の遺伝子が様々なハウスキーピング機能(house keeping function)に参与し、特に炭水化物の輸送および代謝(12.54%)およびアミノ酸輸送および代謝(10.23%)に用いるものを表明する。これらの百分率は他のビフィズス菌ゲノムの百分率と一致する(28、29)。
【0124】
ゲノム解析に基づいて、かつ実験の証拠によって支持され、ビフィズス菌の複数種類の宿主定着因子を同定した、これは胆汁抵抗性およびアドヘシン(adhesins)に参与する機能(11)を含む。胆汁酸は生理的濃度で抗微生物活性(30)を有することができるので、胆汁抵抗性は多くの腸内細菌の定着にとって重要である。全ての6種類のビフィドバクテリウム・シュードカテニュレイタム菌株において胆汁抵抗性を付与した胆汁酸塩加水分解酵素および/または胆汁酸輸送タンパク質を同定した。粘着の面で、全てのゲノムはエノラーゼ酵素およびDnaKをコードする遺伝子を有し、DnaKはビフィドバクテリウム・アニマリス乳酸菌亜種類BI07(B.animalissubsp.lactisBI07)におけるプラスミノーゲン結合の関連タンパク質(31、32)であることが証明される。また、コードで4つのビフィズス菌(B.bifidum)菌株に存在するタムスチン結合に参与するトランスアルドラーゼ(33)、および細胞外マトリックスとの接着を促進すると報告された血管性血友病因子A(von willebrand factor A)(34)の遺伝子もゲノム中に存在する。これらの機能遺伝子の存在は、他のビフィズス菌の菌種と類似し、ビフィドバクテリウム・シュードカテニュレイタムが腸内に定着するゲノム基盤を有することを示す。
【0125】
ビフィドバクテリウム・シュードカテニュレイタムのゲノムマイクロ多様性
【0126】
6個の完全なビフィドバクテリウム・シュードカテニュレイタムゲノム中の平均ヌクレオチド同一性(average nucleotide identity,ANI)は菌種の境界の閾値(35)を満たし、最小値は97.76%である。異なる生育環境から分離したビフィドバクテリウム・シュードカテニュレイタムJCM1200
Tと比較した場合(最大ANI=97.80%)、同一生育環境から分離した菌株内でやや高い類似性(ビフィドバクテリウム・シュードカテニュレイタムPERFEC
T-2017-0002、PERFECT-2017-0003、PERFECT-2017-0004およびC95、最小ANI=99.88%)が観察される。これらの結果は、スポットダイアグラムでビフィドバクテリウム・シュードカテニュレイタムJCM1200
Tと他の5種の菌株との間に小さな共線性と若干の差異(挿入欠失を含む)が存在するが、全てのこれらのゲノムにおける高い共線性がゲノムのスポットダイアグラムの対照により実証されることを示す(
図5A)。このような変異の一つの注意する価値がある実施例は細胞外多糖(EPS)をコードするeps遺伝子クラスターである。EPSは細胞に付着した粘液層を形成してもよく、環境中に放出してもよい(36)。ビフィズス菌により生成されるいくつかのEPSが潜在的に宿主に役立つ多種の有益な活性を有すると考えられ、例えば、免疫系の調節、病原体の抵抗、捕捉剤としておよび微生物コロニーの調節を含む(37)。ビフィドバクテリウム・シュードカテニュレイタムの完全なゲノムにおいて、eps遺伝子クラスターのコピーを同定し(
図4)、かつ我々の5種類のビフィドバクテリウム・シュードカテニュレイタム菌株においてこれらのepsクラスターが同じであるが、ビフィドバクテリウム・シュードカテニュレイタムJCM1200
Tにおいて発見されたものと全く異なる。
【0127】
6種類の完全なゲノムで同定された全てのORFをBLASTPと比較し、かつMCLアルゴリズムでさらにクラスタし、2115個のゲノムが存在することを示す(
図5B)。ここで、〜72%(1520)は全ての6種類のビフィドバクテリウム・シュードカテニュレイタムゲノムに共有され、それはビフィドバクテリウム・シュードカテニュレイタムのコアゲノムを代表する。試験されたビフィドバクテリウム・シュードカテニュレイタムゲノムの部分集合のみに存在する合計312個の非必須ゲノムが同定され、61.48%を超える独特なグループはビフィドバクテリウム・シュードカテニュレイタムJCM1200
Tの特有ものである。COG割り当てに基づいて類似の結果が得られ、6個の完全なビフィドバクテリウム・シュードカテニュレイタムゲノムから1101個のCOGファミリーを同定した。59個のCOGファミリーは検査されたビフィドバクテリウム・シュードカテニュレイタムゲノムの部分集合のみに存在し、かつ37個の別のCOGは独特な単一菌株であり、35個はビフィドバクテリウム・シュードカテニュレイタムJCM1200
Tのみに同定された(
図5)。なお、46個のCOGファミリーは我々の分離株に存在するが、ビフィドバクテリウム・シュードカテニュレイタムJCM1200
Tゲノム中には存在しなかった。興味深いことは、これらの差異の少なくとも一部はモバイルゲノム(mobilom)およびDNA再編成に関連し、我々の分離株はより多くのプロファージ機能、細胞プロセスおよびシグナル伝達に関連する独特なCOGを有し、ビフィドバクテリウム・シュードカテニュレイタムJCM1200
Tはより多くの独特なCOGを有し、CRISPR/CASシステムに参与するものを含む。
【0128】
特に、我々の分離株とビフィドバクテリウム・シュードカテニュレイタムJCM1200
Tは異なる炭水化物源に暴露される。前者は介入後の児童の糞便サンプルから得られ、前記児童は、食物繊維に富み、およびオリゴ糖およびマルトオリゴ糖を含む粉末の全穀類由来とTCM食用植物由来の混合材料を受けるが、後者は相対的により簡単な炭水化物を受けると推測される乳児の糞便から分離する。それに対して、炭水化物の輸送および代謝に参与する遺伝子変異を発見し、例えばビフィドバクテリウム・シュードカテニュレイタムJCM1200
TはCOG0383(α-マンノシダーゼ)、COG3594(フコース4−O−アセチラーゼまたは関連アセチルトランスフェラーゼ)、COG4209(ABC型多糖輸送システム、透過酵素成分)、COG4214(ABC型キシロース輸送システム、透過酵素成分)を欠いているが、独特にCOG1554(トレハロースおよびマルトース加水分解酵素(ホスホリラーゼなど))を有する。
【0129】
コアCOGの異なるコピー数はビフィドバクテリウム・シュードカテニュレイタムのマイクロ多様性にも役立つ。各種の完全なビフィドバクテリウム・シュードカテニュレイタムゲノムに存在する1005個コアCOGにおいて多種機能カテゴリーに参与する51個COGはコピー数の面で少なくとも2つのコピーが異なる(
図6)。これらのCOGの分布によれば、我々の5種類の菌株と比べてビフィドバクテリウム・シュードカテニュレイタムJCM1200
Tが最も異なる。
【0130】
従来の研究において、既に同じ生育環境からの異なる菌株のマイクロ多様性(38-40)を報告した。我々の5種類のビフィドバクテリウム・シュードカテニュレイタム菌株において、差異は特定のコアCOGファミリーに割り当てられた遺伝子のコピー数の変異により引き起こされ、独特なおよび/または非必須なCOGの存在/非存在の面の差異ではない。クラスタ結果(
図6)は菌株の炭水化物の介入に対する差異性応答と一致する。この場合に、ほとんど介入の影響を受けないビフィドバクテリウム・シュードカテニュレイタムC95およびビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0003は、三種類の応答菌株と明らかに分離され、介入に対して適度な応答を有するビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0001およびビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0004の両者をともにグルーピングし、介入に対する応答が最も高い菌株であるビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0002は、他の4種類の菌株とさらに分離することができるが、ビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0001およびビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0004と最も類似する。
【0131】
炭水化物の輸送および代謝について、ビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0002は、COG2814(予測されたアラビノース流出透過酵素、MFSファミリー)およびCOG3250(β−ガラクトシダーゼ/β−グルクロニダーゼ)の最高コピー数を有する。また、5種類の菌株において、106±2(平均値±S.d.)個のORFは炭水化物活性酵素(CAZy)遺伝子であり、41個のCAZyファミリーを占めることを同定した。ビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0002ゲノムは全ての同定されたCAZyファミリーを有し、かつこれらの遺伝子の最大コピー数を有する。それはさらに炭水化物エステラーゼをコードする遺伝子の最高コピー数を有し、前記酵素は複雑さを克服するように植物多糖を脱アセチル化させ、かつ植物多糖の分解中にグリコシド加水分解酵素と協力する(41)。これらの特徴はビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0002が食物介入に対して応答する主な差異であり、従って、それは全てのビフィドバクテリウム・シュードカテニュレイタムの50%以上を占め、かつ介入期間全体にビフィズス菌群全体に占める割合が大きい。
【0132】
結論
【0133】
WTP食事がどのように遺伝性肥満者の体重を減少させるかを研究している従来の研究において、ビフィドバクテリウム・シュードカテニュレイタムは最も豊富なビフィズス菌菌種であると同定される。本明細書においてビフィドバクテリウム・シュードカテニュレイタムの特定菌株が食物介入の応答に対する変化を示し、比較ゲノミクスを利用してこれらの動態背後の可能な原因を同定した。本研究の一部として分離された5種類のビフィドバクテリウム・シュードカテニュレイタム菌株は乳児の糞便から分離されたビフィドバクテリウム・シュードカテニュレイタムJCM1200
Tに対してある差異を示し、これは異なる環境パラメータがゲノムのマイクロ多様性への影響(42)を示す。我々の5種類の分離株に観察されたマイクロ多様性は、だいたいコアCOGファミリーのコピー数の変化であり、かつこれらの差異は介入期間の5種類の菌株の種群変化のため推測性説明を提供する。具体的には、ビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0002は、遺伝的により多様化と植物多糖遺伝子のより多くのコピー数を有し、食事介入に応答する最大豊度を有す。グループ全体が同遺伝子グループより広い環境条件下で生存すること(43)をサポートするため、複数種類の菌株の共存および分布は直感的である。したがって、食事に異なる応答を有する5種類の菌株の共存は、1つのメカニズムとしてヒト腸管におけるビフィドバクテリウム・シュードカテニュレイタムのような重要な善玉菌の安定および回復を保証することができる。しかし、重要なことは、本研究において同定された5種類のビフィドバクテリウム・シュードカテニュレイタム菌株も宿主生物学的臨床パラメータと異なる相関性を有し、少なくともある腸内細菌叢の変化によって有益な機能は菌株の特異性(44)であることが示される。腸内微生物群の食事操作から得られた健康への促進する全ての潜在力を実現される確保するように、本明細書において紹介されたこのようなタイプのさらなる研究は必要である。
【0134】
材料および方法
【0135】
臨床研究
【0136】
該研究は上海交通大学の生命科学と生物技術学院倫理委員会(Ethics Committee of the School of Life Sciences and Biotechnology,Shanghai Jiao Tong University)の承認(No.2012-2016)のもとで行われる。臨床試験は中国臨床試験登録センター(Chinede Clinical Trial Registry)で登録する(No.ChiCTR-ONC-12002646)。肥満児童の後見人から書面の同意書を取得する。
【0137】
中国広東省広州市婦女児童病院(GuangDong Women and children Hosiptal、Guanzhou、China)で、該肥満児童は105日間の入院食事介入を受けた。該ボランティアはいかなる運動計画を実施していない。栄養士のアドバイスに応じて、適量の野菜、果物およびナッツと組み合わせて全穀類、伝統的な中国薬膳およびプレバイオティクスに基づく食事(WTP食事)(7)(食事における三種類のプレハブ即席食品、1号処方、2号処方および3号処方は完美(中国)株式会社により製造される)を投与した。介入は三つのステージに分けられる。ステージI(0日目から60日目)で、該児童は基礎介入(9)を受ける。ステージII(60日目から75日目)で、100g以上の3号処方を消耗する。ステージIII(75日目〜105日目)の期間、彼は依然として100g以上の3号処方を消耗するが、100g以下の1号処方が提供される。
【0138】
7つの時点(0、15、30、45、60、75および105日)に生体サンプル、人体測定データおよび臨床実験室分析を取得する。生物学的臨床パラメータの測定結果は従来の研究(9)と同じである。
【0139】
メタゲノムシークエンスおよび分析
【0140】
糞便サンプルからのDNA抽出は、以上の説明した(45)ように実施されたメタゲノムシークエンスに用いる。上海Genergyバイオテクノロジー株式会社(Shanghai Genergy Bio-technology Co.,Ltd)でIllumina Hiseq2000プラットフォームを使用して7個サンプルに対してシークエンスを行った。Illuminaの説明書に従ってDNAライブラリーを構築する。サプライヤーが指示したワークフローに応じて、クラスタ生成、鋳型ハイブリダイズ、等温増幅、線形化、ブロッキングおよび変性およびシークエンスプライマーのハイブリダイズを行う。ライブラリーを構築し、次に順方向および逆方向にハイスループットシークエンスを行い、151bpを有するペアエンドリードを取得する。
【0141】
Flexbar(46)を利用してリードからのアダプタを調整する。(a)Prinseq(47)を利用して20の品質閾値で3’端からリードを調整し、一番目のヌクレオチドまで達すこと、(b)リードは75bP以下であるまたは「N」塩基を含む場合、リードペアを除去すること、および(c)リードを重複排除すること。ヒトゲノム(H.sapiens、UCSC hg19)とアライメント可能なリードを除去する(Bowtie2(48)とアライメントする)、使用-再編成-no-hd-含んでない-ダブテール(-reorder-no-hd-no-contain-dovetail)。平均的に、各サンプルは25.3±4.1百万(平均値±S.d.)のペアエンドリードを保持し、かつさらに分析するために用いられる。
【0142】
105日目にMetaphlan(49)(--bt2_psに非常に敏感な局所)を利用してビフィズス菌種の豊度を計算する。Sigmaを利用して、我々の5種類のビフィドバクテリウム・シュードカテニュレイタム菌株の豊度(24)を計算する。
【0143】
全ゲノムシークエンスおよびデータ組み合わせ
【0144】
PacBio RSIIシークエンサーを利用してビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0001、ビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0002、ビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0003、ビフィドバクテリウム・シュードカテニュレイタムPERFECT-2017-0004およびビフィドバクテリウム・シュードカテニュレイタムC95のゲノムをシークエンスし、それぞれ約245、415、285、459および198倍の被覆率(Nextomics Biosciences、Wuhan 430000、China)を有する。HGAP/Quiver(50)を用いてサブリードを新たに組み合わせ、次にMiniums2(51)およびQuiverを組み合わせる。各菌株をその染色体に対応する一つのコンティグにアセンブルに組み合わせる。
【0145】
一般的な特徴予測
【0146】
Prodigal v2_60(52)とBLASTアライメントの組み合わせにより、オープンリーディングフレーム予測(Open Reading Frame、ORF)を行う。次にNCBInrデータベースにBLASTPを利用して同定されたORFに注釈を行う。Rnammer v1.2(53)を利用してリボソームRNA遺伝子を検出し、tRNAscan SEv.14(54)を利用してトランスファーRNA遺伝子を同定した。USEARCH v8.0.1517(55)を利用して全ての16S rRNA遺伝子の間の同一性マトリックスを計算する。COGtriangles(ftp://ftp.ncbi.nih.gov/pub/wolf/COGs/COGsoft/)を利用してCOGの割り当てを完了する。
【0147】
パンゲノム計算
【0148】
PGAPv1.12(56)を利用して、ビフィドバクテリウム・シュードカテニュレイタムのパンゲノムを計算した。GF(Gene Family、遺伝子ファミリー)方法により機能遺伝子クラスタを行い、次にパンゲノムスペクトルを構築する。
【0149】
ゲノム比較
【0150】
ソフトウェアパッケージMUMmerv3.0(57)を利用してヌクレオチドレベルで全ゲノム配列アライメントを行った。各ゲノムペアに対して、プログラムJSpeciesバージョン1.2.1(35)を利用して平均的なヌクレオチド同一性(ANI)を計算した。タンパク質レベルでBLASTP(各タンパク質に対して、最大E値1e-10、最小アライメント同一性50%および最小アライメント被覆率50%)を利用してall against all 方式(all against all way)で配列を比較し、次にPGAPv1.12の中でマルコフクラスタアルゴリズム(Markov Cluster Algorithm,MCL)を利用して遺伝子ファミリーにクラスタする。
【0151】
細胞外多糖(Exoplysaccharide,EPS)遺伝子クラスターの同定
【0152】
各ゲノム(37)で仮定されたプライミング?GTF(p?gtf)遺伝子:rfbp(NP_695444)とcpsD(NP_695447)を捜索して、手動でp-gtf周辺の遺伝子を検査することにより、Bifido-epsクラスターのコンピュータ解析を行った。
【0153】
炭水化物活性酵素(CAZyS)の同定
【0154】
dbCANv3.0(58)のローカルバージョンデータベースをダウンロードする。HMMscan(59)を利用して各ゲノムにおける遺伝子とデータベースをアライメントする。dbCANによって提供されたhmmscan-parser.shでアライメントに解析を行い、かつ最適なヒットを保留する。
【0155】
データおよび材料の入手可能性
【0156】
本研究で生成された全てのゲノムおよびメタゲノムデータのデータは登録番号PRJEB18557でヨーロッパのヌクレオチドライブラリー(European Nucleotide Archive、ENA)(WWW.ebi.ac.uk/ena)のもとに寄託される。我々の分析に用いるビフィドバクテリウム・シュードカテニュレイタムの全ての他のゲノムはいずれも以下の登録番号でNCBIデータベースからダウンロードする:AP012330.1、CDPW00000000.1、ABXX02000001、JEOD01000001、JGZF01000001,ビフィドバクテリウム・シュードカテニュレイタムD2CAは、MetaHIT(http://www.sanger.ac.uk/resources/downloads/bacteria/metahit/)からダウンロードする。
【0157】
以上の明細書に言及された全ての出版物はいずれも本明細書において参考として援用される。本発明の範囲および精神から逸脱することなく、記載された本発明の方法およびシステムの様々な修正および変更は当業者にとって明らかである。以上、具体的な好適な実施形態を挙げて本発明を説明したが、保護請求された発明は適切にこのような具体例に限定されることを理解するべきである。実際に、記載された本発明の実施形態の生化学的および生物学的技術または関連分野の当業者にとって明らかな様々な修正はいずれも下記特許請求の範囲内にある。
【0158】
引用文献
1.Zhao L. 2013. The gut microbiota and obesity: from correlation to causality. Nat Rev Microbiol 11:639-647.
2.Boulange CL, Neves AL, Chilloux J, Nicholson JK, Dumas ME. 2016. Impact of the gut microbiota on inflammation, obesity, and metabolic disease. Genome Med 8:42.
3.Qin N, Yang F, Li A, Prifti E, Chen Y, Shao L, Guo J, Le Chatelier E, Yao J, Wu L, Zhou J, Ni S, Liu L, Pons N, Batto JM, Kennedy SP, Leonard P, Yuan C, Ding W, Chen Y, Hu X, Zheng B, Qian G, Xu W, Ehrlich SD, Zheng S, Li L. 2014. Alterations of the human gut microbiome in liver cirrhosis. Nature 513:59-64.
4.Zhang X, Zhang D, Jia H, Feng Q, Wang D, Liang D, Wu X, Li J, Tang L, Li Y, Lan Z, Chen B, Li Y, Zhong H, Xie H, Jie Z, Chen W, Tang S, Xu X, Wang X, Cai X, Liu S, Xia Y, Li J, Qiao X, Al-Aama JY, Chen H, Wang L, Wu QJ, Zhang F, Zheng W, Li Y, Zhang M, Luo G, Xue W, Xiao L, Li J, Chen W, Xu X, Yin Y, Yang H, Wang J, Kristiansen K, Liu L, Li T, Huang Q, Li Y, Wang J. 2015. The oral and gut microbiomes are perturbed in rheumatoid arthritis and partly normalized after treatment. Nat Med 21:895-905.
5.Xu J, Lian F, Zhao L, Zhao Y, Chen X, Zhang X, Guo Y, Zhang C, Zhou Q, Xue Z, Pang X, Zhao L, Tong X. 2015. Structural modulation of gut microbiota during alleviation of type 2 diabetes with a Chinese herbal formula. ISME J 9:552-562.
6.Kelly CR, Kahn S, Kashyap P, Laine L, Rubin D, Atreja A, Moore T, Wu G. 2015. Update on Fecal Microbiota Transplantation 2015: Indications, Methodologies, Mechanisms, and Outlook. Gastroenterology 149:223-237.
7.Xiao S, Fei N, Pang X, Shen J, Wang L, Zhang B, Zhang M, Zhang X, Zhang C, Li M, Sun L, Xue Z, Wang J, Feng J, Yan F, Zhao N, Liu J, Long W, Zhao L. 2014. A gut microbiota-targeted dietary intervention for amelioration of chronic inflammation underlying metabolic syndrome. FEMS Microbiol Ecol 87:357-367.
8.Human Microbiome Project C. 2012. Structure, function and diversity of the healthy human microbiome. Nature 486:207-214.
9.Zhang C, Yin A, Li H, Wang R, Wu G, Shen J, Zhang M, Wang L, Hou Y, Ouyang H, Zhang Y, Zheng Y, Wang J, Lv X, Wang Y, Zhang F, Zeng B, Li W, Yan F, Zhao Y, Pang X, Zhang X, Fu H, Chen F, Zhao N, Hamaker BR, Bridgewater LC, Weinkove D, Clement K, Dore J, Holmes E, Xiao H, Zhao G, Yang S, Bork P, Nicholson JK, Wei H, Tang H, Zhang X, Zhao L. 2015. Dietary Modulation of Gut Microbiota Contributes to Alleviation of Both Genetic and Simple Obesity in Children. EBioMedicine 2:966-982.
10.Wu H, Wang R, Zhao Y, Pang X, Shen J, Zhang C. 2015. Exploring Carbohydrate Utilization Capacity of Bifidobacterium pseudocate-
nulatum Isolated from a Morbidly Obese Child after Dietary Intervention. Genomics and Applied Biology 34:1384-1391.
12.Grimm V, Westermann C, Riedel CU. 2014. Bifidobacteria-host interactions--an update on colonisation factors. Biomed Res Int 2014:960826.
13.Labruna G, Pasanisi F, Nardelli C, Caso R, Vitale DF, Contaldo F, Sacchetti L. 2011. High leptin/adiponectin ratio and serum triglycerides are associated with an "at-risk" phenotype in young severely obese patients. Obesity (Silver Spring) 19:1492-1496.
14.Zweigner J, Schumann RR, Weber JR. 2006. The role of lipopolysaccharide-binding protein in modulating the innate immune response. Microbes Infect 8:946-952.
15.Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, Fernandes GR, Tap J, Bruls T, Batto JM, Bertalan M, Borruel N, Casellas F, Fernandez L, Gautier L, Hansen T, Hattori M, Hayashi T, Kleerebezem M, Kurokawa K, Leclerc M, Levenez F, Manichanh C, Nielsen HB, Nielsen T, Pons N, Poulain J, Qin J, Sicheritz-Ponten T, Tims S, Torrents D, Ugarte E, Zoetendal EG, Wang J, Guarner F, Pedersen O, de Vos WM, Brunak S, Dore J, Meta HITC, Antolin M, Artiguenave F, Blottiere HM, Almeida M, Brechot C, Cara C, Chervaux C, Cultrone A, Delorme C, Denariaz G, et al. 2011. Enterotypes of the human gut microbiome. Nature 473:174-180.
16.Fujimoto T, Imaeda H, Takahashi K, Kasumi E, Bamba S, Fujiyama Y, Andoh A. 2013. Decreased abundance of Faecalibacterium prausnitzii in the gut microbiota of Crohn's disease. J Gastroenterol Hepatol 28:613-619.
17.Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermudez-Humaran LG, Gratadoux JJ, Blugeon S, Bridonneau C, Furet JP, Corthier G, Grangette C, Vasquez N, Pochart P, Trugnan G, Thomas G, Blottiere HM, Dore J, Marteau P, Seksik P, Langella P. 2008. Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proc Natl Acad Sci U S A 105:16731-16736.
18.Ganzle MG, Follador R. 2012. Metabolism of oligosaccharides and starch in lactobacilli: a review. Front Microbiol 3:340.
19.Andreasen AS, Larsen N, Pedersen-Skovsgaard T, Berg RM, Moller K, Svendsen KD, Jakobsen M, Pedersen BK. 2010. Effects of Lactobacillus acidophilus NCFM on insulin sensitivity and the systemic inflammatory response in human subjects. Br J Nutr 104:1831-1838.
20.Hulston CJ, Churnside AA, Venables MC. 2015. Probiotic supplementation prevents high-fat, overfeeding-induced insulin resistance in human subjects. Br J Nutr 113:596-602.
21.Charbonneau MR, Blanton LV, DiGiulio DB, Relman DA, Lebrilla CB, Mills DA, Gordon JI. 2016. A microbial perspective of human developmental biology. Nature 535:48-55.
22.Wang J, Tang H, Zhang C, Zhao Y, Derrien M, Rocher E, van-Hylckama Vlieg JE, Strissel K, Zhao L, Obin M, Shen J. 2015. Modulation of gut microbiota during probiotic-mediated attenuation of metabolic syndrome in high fat diet-fed mice. ISME J 9:1-15.
23.Stenman LK, Waget A, Garret C, Klopp P, Burcelin R, Lahtinen S. 2014. Potential probiotic Bifidobacterium animalis ssp. lactis 420 prevents weight gain and glucose intolerance in diet-induced obese mice. Benef Microbes 5:437-445.
24.STALEY J, KRIEG N. 1984. Classification of prokaryotes organisms: an overview. KRIEG, NR,HOLT, JG Bergey’s manual of systemayic bacteriology 1.
25.Ahn TH, Chai J, Pan C. 2015. Sigma: strain-level inference of genomes from metagenomic analysis for biosurveillance. Bioinformatics 31:170-177.
26.Morita H, Toh H, Oshima K, Nakano A, Arakawa K, Takayama Y, Kurokawa R, Takanashi K, Honda K, Hattori M. 2015. Complete genome sequence of Bifidobacterium pseudocatenulatum JCM 1200(T) isolated from infant feces. J Biotechnol 210:68-69.
27.Bottacini F, Medini D, Pavesi A, Turroni F, Foroni E, Riley D, Giubellini V, Tettelin H, van Sinderen D, Ventura M. 2010. Comparative genomics of the genus Bifidobacterium. Microbiology 156:3243-3254.
28.Galperin MY, Makarova KS, Wolf YI, Koonin EV. 2015. Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res 43:D261-269.
29.Bottacini F, O'Connell Motherway M, Kuczynski J, O'Connell KJ, Serafini F, Duranti S, Milani C, Turroni F, Lugli GA, Zomer A, Zhurina D, Riedel C, Ventura M, van Sinderen D. 2014. Comparative genomics of the Bifidobacterium breve taxon. BMC Genomics 15:170.
30.O'Callaghan A, Bottacini F, O'Connell Motherway M, van Sinderen D. 2015. Pangenome analysis of Bifidobacterium longum and site-directed mutagenesis through by-pass of restriction-modification systems. BMC Genomics 16:832.
31.Begley M, Gahan CGM, Hill C. 2005. The interaction between bacteria and bile. Fems Microbiology Reviews 29:625-651.
32.Candela M, Biagi E, Centanni M, Turroni S, Vici M, Musiani F, Vitali B, Bergmann S, Hammerschmidt S, Brigidi P. 2009. Bifidobacterial enolase, a cell surface receptor for human plasminogen involved in the interaction with the host. Microbiology 155:3294-3303.
33.Candela M, Centanni M, Fiori J, Biagi E, Turroni S, Orrico C, Bergmann S, Hammerschmidt S, Brigidi P. 2010. DnaK from Bifidobacterium animalis subsp. lactis is a surface-exposed human plasminogen receptor upregulated in response to bile salts. Microbiology 156:1609-1618.
34.Gonzalez-Rodriguez I, Sanchez B, Ruiz L, Turroni F, Ventura M, Ruas-Madiedo P, Gueimonde M, Margolles A. 2012. Role of extracellular transaldolase from Bifidobacterium bifidum in mucin adhesion and aggregation. Appl Environ Microbiol 78:3992-3998.
35.Pareti FI, Fujimura Y, Dent JA, Holland LZ, Zimmerman TS, Ruggeri ZM. 1986. Isolation and characterization of a collagen binding domain in human von Willebrand factor. J Biol Chem 261:15310-15315.
36.Richter M, Rossello-Mora R. 2009. Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci U S A 106:19126-19131.
37.Schmid J, Sieber V, Rehm B. 2015. Bacterial exopolysaccharides: biosynthesis pathways and engineering strategies. Front Microbiol 6:496.
38.Hidalgo-Cantabrana C, Sanchez B, Milani C, Ventura M, Margolles A, Ruas-Madiedo P. 2014. Genomic overview and biological functions of exopolysaccharide biosynthesis in Bifidobacterium spp. Appl Environ Microbiol 80:9-18.
39.Zhang XL, Gao PP, Chao QF, Wang LH, Senior E, Zhao LP. 2004. Microdiversity of phenol hydroxylase genes among phenol-degrading isolates of Alcaligenes sp from an activated sludge system. Fems Microbiology Letters 237:369-375.
40.Mao YJ, Zhang XJ, Xia X, Zhong HH, Zhao LP. 2010. Versatile aromatic compound-degrading capacity and microdiversity of Thauera strains isolated from a coking wastewater treatment bioreactor. Journal of Industrial Microbiology & Biotechnology 37:927-934.
41.Patra R, Chattopadhyay S, De R, Ghosh P, Ganguly M, Chowdhury A, Ramamurthy T, Nair GB, Mukhopadhyay AK. 2012. Multiple Infection and Microdiversity among Helicobacter pylori Isolates in a Single Host in India. Plos One 7.
42.Biely P. 2012. Microbial carbohydrate esterases deacetylating plant polysaccharides. Biotechnol Adv 30:1575-1588.
43.Jezbera J, Jezberova J, Kasalicky V, Simek K, Hahn MW. 2013. Patterns of Limnohabitans microdiversity across a large set of freshwater habitats as revealed by Reverse Line Blot Hybridization. PLoS One 8:e58527.
44.Moore LR, Rocap G, Chisholm SW. 1998. Physiology and molecular phylogeny of coexisting Prochlorococcus ecotypes. Nature 393:464-467.
45.Zhang C, Zhao L. 2016. Strain-level dissection of the contribution of the gut microbiome to human metabolic disease. Genome Med 8:41.
46.Fei N, Zhao L. 2013. An opportunistic pathogen isolated from the gut of an obese human causes obesity in germfree mice. ISME J 7:880-884.
47.Dodt M, Roehr JT, Ahmed R, Dieterich C. 2012. FLEXBAR-Flexible Barcode and Adapter Processing for Next-Generation Sequencing Platforms. Biology (Basel) 1:895-905.
48.Schmieder R, Edwards R. 2011. Quality control and preprocessing of metagenomic datasets. Bioinformatics 27:863-864.
49.Langmead B, Salzberg SL. 2012. Fast gapped-read alignment with Bowtie 2. Nat Methods 9:357-359.
50.Segata N, Waldron L, Ballarini A, Narasimhan V, Jousson O, Huttenhower C. 2012. Metagenomic microbial community profiling using unique clade-specific marker genes. Nat Methods 9:811-814.
51.Chin CS, Alexander DH, Marks P, Klammer AA, Drake J, Heiner C, Clum A, Copeland A, Huddleston J, Eichler EE, Turner SW, Korlach J. 2013. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods 10:563-569.
52.Sommer DD, Delcher AL, Salzberg SL, Pop M. 2007. Minimus: a fast, lightweight genome assembler. BMC Bioinformatics 8:64.
53.Hyatt D, Chen GL, Locascio PF, Land ML, Larimer FW, Hauser LJ. 2010. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11:119.
54.Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW. 2007. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res 35:3100-3108.
55.Lowe TM, Eddy SR. 1997. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res 25:955-964.
56.Edgar RC. 2010. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26:2460-2461.
57.Zhao Y, Wu J, Yang J, Sun S, Xiao J, Yu J. 2012. PGAP: pan-genomes analysis pipeline. Bioinformatics 28:416-418.
58.Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SL. 2004. Versatile and open software for comparing large genomes. Genome Biol 5:R12.
59.Yin Y, Mao X, Yang J, Chen X, Mao F, Xu Y. 2012. dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res 40:W445-451.
60.Finn RD, Clements J, Eddy SR. 2011. HMMER web server: interactive sequence similarity searching. Nucleic Acids Res 39:W29-37.