macでインフォマティクスmacでインフォマティクススマートフォン用の表示で見るmacでインフォマティクスNGS関連のインフォマティクス情報についてまとめています。novoアセンブラについて少し誤解している人がいたので、注意喚起を兼ねて簡単にまとめておく。relatedなゲノム配列が様々なabundanceで存在している。メタゲノムのアセンブラは、この混合の(バルクの)シーケンシングリードセットから、複数の近縁な株のコンセンサスのバックボーン配列を構築するアセンブラである。この性質により、配列多少の多様性があっても、同種のゲノムはより長い領域に渡って再構成される。配列が高度に断片化していると、ゲノムを分類したり特徴を調べたりすることができなくなるため、このようなメタゲノムアセンブラの一般的性質は、下流解析を効率的に進めるためにも欠かせない特徴と言える。この特徴により、異なる生物間で保存性の高い領域も1つの配列にまとめられる可能性がある。このことを理解していないと、ターゲットの配列が1生物からのみ由来すると判断して、解析に影響を及ぼす恐れがある。ここでは、ツールに振り回されないためにも、シミュレーションリードを使ってメタゲノムアセンブラの性質を実感できる簡単なテストを行ってみる。異なるクラスタの生物間で長い領域に渡って保存性されている配列として、遺伝子クラスターが挙げられる。ここでは、よく研究されている遺伝子クラスターであるニトロゲナーゼのnifクラスターを例に説明する。GenBnakで公開されている4生物からfetchして、blastnによるアラインメント比較した結果である。ペアワイズの配列比較と作図には、GenBankから自動でblast比較して簡単にゲノム比較の図を作れるEasyFig(紹介)を使っている。clusterをfetchした。clusterの配列をリファレンスに、ARTのHiseq2500プロファイルでシミュレートして、ペアエンドのリード(150x2)を表のcoverageの深さで発生させた。novoアセンブリを実行する。metaspades3、spadesアセンブリのGFAファイルをBandageに読み込んで視覚化した。clusterはほぼ全長が再構成できている。次がmetaspadesのアセンブリ結果。metaspadesでは再構成に成功している配列もあるものの、上の交差した配列群のように、共通した配列部分でパスを共有している領域がある。OTUなどの基準で比較を行うように、メタゲノムのアセンブラも、株レベルの細かな違いは無視して、近い距離にあるゲノムを1つのコンセンサス配列として表し、下流の様々な解析に適した配列を作成しようとする。そのため、メタゲノムのアセンブラを使うと、保存された遺伝子クラスターがパスを共有したコンセンサス配列として表現される場合があることを示した。実際のメタゲノムのシーケンシングリードをspadesアセンブラなどでアセンブリすれば分かるが、アセンブリグラフの様々な部位で分岐が生じるため、進化速度が遅い領域を除き、大半の領域では、アセンブリ後の配列は高度に断片化してしまう。これではそもそもアセンブルをする意味が薄くなってしまう。メタゲノムのアセンブラは、まじめに考えれば解釈が困難な混合物のシーケンシングデータをなんとかして処理するため、また計算負荷を下げるため、アセンブリグラフをシンプルにし長い配列を出力しようと試みる。これに対するトレードオフの1つが、よく似た配列のある種誤ったコンセンサス表現である。clusterのようなよく似た長い配列を分離するには、クラスターを構成するリードを回収して精度の高いspadesなどのアセンブラで再アセンブリするか(*1)、倍数体ゲノムからのハプロタイプ再構成のような追加のテクニックが必要になる。いうことで、メタゲノムシーケンスの目的が保存された遺伝子クラスターの分析である場合、そのままwholeアセンブリするだけでは良い結果を出せない可能性があることを示した。allのANI比較とmauveのゲノム比較を行った。(ANI比較は4つのnifクラスターしか使用していない)。ANI比較ではそれぞれのnifクラスターは70%以下の塩基配列相同性になった。mauveによるゲノム比較。シンテシーブロック表現のeasyfigより細かい部分まで視覚化される。databaseに通常のゲノムとして登録することは避けるべきです。以下を読んでみてください。*1コミュニティの複雑性が低ければうまくいくかもしれない。DoriC注目記事はてなブログをはじめよう!kazumaxneoさんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。