平成26年度ライフサイエンスデータベース統合推進事業
トーゴーの日シンポジウム2014

ゲノム・メタゲノム情報統合による微生物DBの超高度化推進
東京工業大学地球生命研究所
黒川 顕© 2014 黒川顕（東京工業大学）licensed under CC表示2.1日本

1

微生物研究を取り巻く現状• 微生物は地球上のいたる所に存在し環境と密接に関与している• 微生物研究はバイオ分野のみならず、他の多くの分野と連携可能• 既に多様なDBが多数存在する• しかし、微生物と環境との関連性を記述しているDBは存在しない• さらに、専門知識を持っていないバイオ分野以外の人は利用困難
2integrates lots of data related to microbes.
Especially, we integrates the microbial data that can be linked to genomes.

http://microbedb.jp/

Gene
Ortholog: MBGD

Taxon

Environment
Taxonomy:
NCBI Taxonomy

Genome: GTPS/RefSeq
Annotation:TogoAnnotation

Culture Collection:
NBRC/JCM

Metadata:
INSDC SRA

Metagenome:
INSDC SRA
3

Red color indicates our collaborators.

全てのデータをゲノムを核として統合する
海水

メタゲノム
(環境)
配列相同性検索ゲノム(系統)遺伝子
クラスタリング

オーソログ
(遺伝子)

土壌

ヒト腸内

微生物における各種DBを統合化し，環境情報との連携を徹底的に記述した新規DB構築を実現• 異分野データの統合化のため、セマンティックウェ
ブの技術を徹底的に導入• データ間をリンクするためのゲノム情報、オーソロ
グ遺伝子情報、メタゲノム情報の整備• 全データのRDF化、各データID間のリンク構築• 各種オントロジーの開発、各データにマッピング• アノテーション高度化システムの開発• ユーザ認証システム基盤開発• 検索結果可視化アプリケーション「Stanza」の開発• ゲノム自動アノテーションシステムMiGAPとの連携
5

各種データのトリプル数
グラフ名refseq
mbgd
gtps
taxonomymeta16S
gazetteer
srs_metadata
srs_ortholog
gobrc
gold
srs
sopdo
meomsv
mpo
mccv
その他中間データ
合計

説明作成元
RefSeq Prokaryoteゲノムデータ
DBCLS
MBGD Orthologデータ
基生研
GTPSゲノムデータ
遺伝研
SPARQLthonで作成したNCBI Taxonomyオントロジー改良版 DBCLS,遺伝研,東工大
各SRSメタ16Sの系統組成データ
東工大
地理オントロジー
外部機関
SRSメタ16S・メタゲノムの様々なメタデータ
東工大
各SRSメタゲノムのMBGD Ortholog組成
東工大,基生研
Geneオントロジー
外部機関
JCM/NBRC菌株データ with NCBI Taxonomy ID
遺伝研,東工大,DBCLS
GOLDの個別ゲノムのMEO等へのオントロジーマッピング
東工大,DBCLS
データ
SRSメタ16S・メタゲノムのMEO等へのオントロジーマッピング
東工大
データ
Sequenceオントロジー
外部機関
感染症オントロジー + 症状オントロジー + ゲノムへのオント
東工大
ロジーマッピングデータ
微生物の生息環境オントロジー
東工大
SRSメタ16S・メタゲノムのメタデータオントロジー
東工大
微生物フェノタイプオントロジー
DBCLS
菌株オントロジー
東工大,DBCLS
いくつかのデータ集計系のSPARQLクエリは遅いため、MSSが集計結果のデータを作成

トリプル数
550,273,744
291,714,037
197,069,932
10,183,714
9,831,600
7,062,536
4,982,739
2,026,746
1,211,571
903,319
150,899
53,69143,060
8,809
4,975
1,601
734
293
440,773
1,075,964,773
6

オーソログ情報のRDF化
汎用 Ortholog Ontology (OrthO)の開発
owl:Class

RDF版MBGD に対するSPARQL queryのポータルサイト

owl:ObjectProperty

owl:DatatypeProperty開発したオントロジー• FALDO (Feature Annotation Location Description Ontology)– ゲノム中の各featureの位置情報を記述するためのオントロジー（w/BioHackathon）• INSDC Ontology– INSDCエントリのfeatureとqualifierのターム記述のためのオントロジー（w/DBCLS）• MCCV (Microbial Culture Collection Vocabulary)– 菌株データを記述するためのオントロジー• MEO (Metagenome/Microbe Environmental Ontology)– 細菌の生息環境を記述するためのオントロジー• PDO/CSSO (Pathogenic Disease Ontology with Symptom)– 細菌が引き起こす感染症の情報および感染症の症状を連結したオントロジー• GMO (Growth Media Ontology)– 細菌の培地情報を記述するためのオントロジー（w/DBCLS）

8

http://microbedb.jp/9

MicrobeDB.jpの開発で実現したこと
1. 既存のゲノム中の各遺伝子の情報 (オーソログ、系統プロファイル、環境プロファイル)
2. 菌株保存機関に存在する菌株の情報 (生育培地、表現型情報、遺伝子機能組成)
3. 様々な環境中の細菌群集の情報 (系統組成、遺伝子機能組成）
4. 上記の情報をシームレスに統合
問合せ例：
高温環境に多く存在する遺伝子はどのような遺伝子か?その遺伝子は、どの系統が主に持っているのか? 10

本研究開発の目標・ねらい
MicrobeDB.jpを• より広い微生物種を対象として拡張• データ収集や更新自動化による持続可能なシステム• 最先端解析プロトコルを実装した解析結果の可視化
研究者コミュニティだけでなく不特定多数のイノベータを対象とした利用性の向上を徹底する

単なる統計量の羅列ではなく、大規模データから新規知識発見を容易に行う事が可能な、今までのDBを超えたDBシステムを構築する事を目標とする。
11

主たる共同研究者
東京工業大学黒川 顕：微生物DBにおける研究統括
山田拓司：ゲノム・メタゲノムDBの構築
森 宙史：真菌類データの整備、DB自動更新システムの開発
山本 希：解析Stanza＆オントロジー開発
鈴木真也：高度解析Stanzaの開発
国立遺伝学研究所
中村保一：藻類データの整備
菅原秀明：MiGAPとの連携強化
神沼英里：MeGAPとの連携強化
藤澤貴智：アクセスレベルの制限システムの開発
基礎生物学研究所
内山郁夫：真核生物に対するオーソログ解析手法の開発
千葉啓和：ドラフトゲノムのオーソログ解析
西出浩世：オーソログを基軸とした各種データ統合の推進

12具体的な研究開発項目（3年間）
1.
2.
3.
4.
5.

各種オミックスデータへの対応
真菌類および藻類を対象とした拡張
各種オントロジー、ボキャブラリの開発
解析プロトコルを実装した各種Stanzaの開発
データの収集およびクオリティコントロール、更新の自動化など持続可能なシステムの構築
6. データ共有、公開におけるアクセスレベルの制御シス
テムの構築
7. 構築したシステムを幅広い分野の研究者に活用してもらうためのユーザビリティの向上
13

４. 解析プロトコルを実装した
各種Stanzaの開発
多種多様な情報が混在しているゲノムやメタゲノム等の複雑なデータから知識発見をするために、比較ゲノム解析や比較メタゲノム解析など様々な解析Stanzaを開発する
各系統と温度との相関係数リスト

SRR FastaファイルをSRS単位で結合

UCLUST
Identity > 97%, cov > 90%
OTU代表配列

UCHIME Referencemodeでキメラを検出

解析Stanzaによる結果の可視化

UCHIME De novo mode
でキメラを検出

両modeでキメラとされた
OTUをキメラと判定、そのOTUを構成する全リードを除去

環境ごとの系統組成のグラフ

キメラ除去済みOTURDP Classifierを用いてbootstrap cutoff >= 50%
で各OTUの代表配列を系統アサインメント

14

Stanza Example

・ Gene Annotation

・ Ortholog list

・ Genome Information

Stanza Example

・Taxonomic compositionof 16S rRNA gene ampliconsequencing analysis

・Functional and taxonomic
composition of ametagenome sample

検索ワード: lake
lake

MEO OWL
meo:pond is_a meo:lake

菌株ーMEOマッピングRDF
Strain_A mccv:isolation_source meo:pond

Strain_Aのようなオントロジーを介した推論検索を実行
lake由来のゲノム解読済み株

lake由来のメタゲノムサンプル
で多いオーソログ

lake由来の菌株

lake由来のメタゲ
ノムサン
プル
lake由来のメタ16Sサン
プルの平均的な系統組
成lakeと他の語彙のMEO中での位置関係

17

Stanza categories in MicrobeDB.jp
Gene Definition
Gene Publication
Ortholog Definition
Gene Annotation
Ortholog Group Members
Ortholog Cluster
Genome Information
GTPS Gene/Genome Feature
RefSeq Gene/Genome Feature

GTPS Genome
GTPS Genome Definition
Other Collection Numbers
Pathogen Information
Phenotype Information
RefSeq Genome
RefSeq Genome Definition
Strain Definition
Strain Genome
Strain Reference
Taxon Definition
Taxon Hierarchy

Genes

Sample Function
Mapping to Environment (Chromosome)
Mapping to Environment (Plasmid)
Ortholog Abundance among Environments
Ortholog Abundance in Environment

Disease Definition
Environment Definition
MEO Hierarchy
Ontology View
Taxon
Environment MEO
Meta16S Sample List
Metagenome Sample List
Numeric Metadata Histogram
Sample Definition
Sample Metadata
SRS Cross Reference
Genome-Sequenced Strains
Symptom List
Sequenced Genome List
Strain List
Taxonomic Composition of Genomes
Taxonomic Composition of Meta 16S
Human Meta Body Mapping
Strain Metadata

【研究テーマ】工業排水のメタゲノム解析
1.工業廃水における細菌群集構造は？
→細菌群集構造を明らかにする

2.サンプル間でどのような変化がある？
→サンプル（時系列等）ごとに細菌群集を比較する

3.どのような細菌種、遺伝子が環境因子と強い相関を持つ？
→環境因子と強い相関を持つ遺伝子・種の探索

4.同様の環境ではどのような細菌群集が存在する？
→同様のメタデータを持つサンプルとの比較

5.同様の細菌群集はどのような環境に存在する？
→同様の細菌群集構造を示すサンプルとの比較

細菌群集構造スタンザ
サンプルにおける細菌群集構造を分類群ごとに円グラフで可視化

系統・遺伝子機能組成の比較解析スタンザ
複数の環境やサンプルの系統組成・遺伝子機能組成を比較できるスタンザ
環境、サンプルの系統的・機能的な特徴を概観できる

土壌水
糞便サンプル

21

散布図スタンザ
数値メタデータと系統や遺伝子の存在量の散布図を表示するスタンザ（１）Thermocrinis (genus) は温度と高い相関を示す
好熱性細菌（２）Ruminococcaceae (family) は20~40℃にサンプルが集中
動物共生細菌が多い

温度（℃）

温度（℃）

Thermocrinisの存在量の散布図

温度とRuminococcaceaeの存在量の散布図

22相関解析スタンザ
数値メタデータと系統・遺伝子機能の存在量の相関を概観できる
系統ID/
遺伝子機能ID
系統名/
遺伝子機能名Pearson
相関係数

サンプル数

メタデータ値
平均23

Numeric Metadata Stanza
数値メタデータ※の分布をヒストグラムで表示するスタンザ※温度、年齢、BMI etc.
入力：MSV ID (メタデータID), SRS ID (メタゲノムサンプルID))
出力：ヒストグラム
・SRS IDも入力すると、そのサンプルが位置する部分を強調
200180160例）温度メタデータのヒストグラム

1401201008060402000 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210 220 230 240 250 260 270 280 290 300

Temperature (℃)

VITCOMIC2スタンザ
ユーザのメタ16S・メタゲノムデータと、類似した系統組成を持ったサンプルおよび環境カテゴリ（MEO ID）を探索し明快に視覚化Human Meta Body Map Stanza
自分の興味がある系統・遺伝子はヒトの体のどこに多いのか？

MicrobeDB.jpへのゲートウェイGenome Refine拡張開発

メタゲノム解析パイプラインMeGAPデータ入力について対応するとともに、ユーザビリティの向上を目的としたウェブリソースを整備し、データの公開・再利用の促進のための強化を図る

【研究テーマ】工業排水のメタゲノム解析
1.工業廃水における細菌群集構造は？
→細菌群集構造を明らかにする

水処理関連企業A：
2.サンプル間でどのような変化がある？
→サンプル（時系列等）ごとに細菌群集を比較する
分析ソリューション事業

3.どのような細菌種、遺伝子が環境因子と強い相関を持つ？
• 分析サービスの差別化→環境因子と強い相関を持つ遺伝子・種の探索• 膨大な分析結果を視覚化

4.同様の環境ではどのような細菌群集が存在する？
→同様のメタデータを持つサンプルとの比較
H27年度から事業化

5.同様の細菌群集はどのような環境に存在する？
→同様の細菌群集構造を示すサンプルとの比較