バイオインフォマティクスの備忘録コンテンツへスキップバイオインフォマティクスの備忘録バイオインフォマティクスやPython等プログラミング関連の情報をメモしておきます。LLCからお願いいたします。賀建奎副教授がゲノム編集の技術を人の受精卵に応用し、双子の女児を誕生させたと主張して話題になっています。遺伝子編集ベイビーヒトの受精卵をゲノム編集で操作し子供を誕生させることは、倫理、技術の両面から問題が指摘されており、今までに誕生した例はありませんでした。今回の主張が事実だとすると、世界初の誕生例となります。ポスドクという経歴です。その後、南方科技大学に就職したとされています。賀建奎副教授はビジネスにも積極的で、7社の株主、6社の代表取締役を務めているようです。設置し、調査を進める予定であるとのことです。中国の科学者の間からも批判の声が上がっており、122人が共同声明を発表し、「人体実験」を行ったことに対して反対を表明しました。Editingでの発表2018年11月28日、香港大学で開かれた国際ゲノム編集会議で賀建奎副教授が発表を行いました。この発表は、今回の騒動以前から予定されていたものであり、元の発表内容を変更して行われました。発表では、マウスにおけるCCR5ノックアウトについての説明、そして今回の件のヒトにおけるCCR5ノックアウトについての説明がありました。その後、参加者から質問を受け付けました。その中で、遺伝病を持った患者たちを助けるために、今回の成果が役に立つことを強調しました。今回生まれた双子以外にも、もう一人妊娠している女性がいることを明らかにしました。バイオインフォマティクス中国で「ゲノム編集の双子」が誕生した件についてのまとめにコメントバイオインフォマティクスにおけるパイプラインの構築にはどのフレームワークを使用すれば良いかバイオインフォマティクス解析においては、多数のOSS(オープンソースソフトウェア)を組み合わせてパイプラインを構築することがよくあります。パイプラインの構築は、自分の好きな言語を選んで自己流で行うことも可能ですが、パイプライン構築用のフレームワークを使用することで、可読性が高く、簡潔なプログラムを書くことができるようになります。多くのフレームワークではDockerコンテナに対応しているため、多数のOSSの実行環境を整備するといった手間のかかる作業を短縮することができます。Dockerコンテナを使用すると、パイプラインのポータビリティが高まりますので、他の研究者へパイプラインを渡したり、クラウド上でパイプラインを実行したりすることも容易になります(再現性の向上)。パイプライン構築用のフレームワークは、awesome-pipelineにまとめられています。ありますので、そのうち主なものを以下にまとめました。snakemakePythonベースのフレームワーク。Pythonの記法が使えるので、柔軟なワークフローが設計できる。柔軟なワークフローが設計できる分、タスクごとにDockerで実行環境を分離する考え方と相性が悪く、パイプラインのポータビリティはあまり良くない印象。Language(CWL)ワークフローを記述するための統一言語。CWLを実行するためのソフトウェアはcwltool、Arvados、toil等たくさん開発されている。[]NextflowCWLと目指す方向はほぼ一緒であるが、CWLは言語仕様と実行エンジンが分離しているのに対して、Nextflowは同一である。その分、Nextflowの方がCWLと比べて柔軟にワークフローを記述できる印象。バイオインフォマティクスバイオインフォマティクスにおけるパイプラインの構築にはどのフレームワークを使用すれば良いかにコメントバイオインフォマティクスを始めるのにどの言語を勉強したら良いかバイオインフォマティクスをこれから始める人にどの言語を使用したら良いかと質問されることがよくある。Rubyのいずれかだと思う。結論からいうと、2018年現在でバイオインフォマティクスをこれから始める人にふさわしい言語はPython一択であると考えられる。Pythonはユーザー数が多いそれぞれの言語のユーザー数がどれくらいいるか正確に把握することは難しいが、参考までにGoogleトレンドの結果を示す。全世界:日本:世界的には2011年後頃からPerlを逆転してPythonがトップになっている。日本でも2014年頃にはPythonがPerlを抜いている。ユーザー数が多い言語を選ぶべきかPythonのユーザー数が多いことは理解いただけたかと思うが、なぜユーザー数が多い言語を選ぶべきかを説明する。Overflow」の検索結果を示す。圧倒的にPythonの件数が多くなっている。WEB上に情報が多ければ、問題が起きた時にググるだけで解決する可能性が高い。BioRubyを比較してみる。そのライブラリーの開発に貢献した人数である「Contributer」の数を見てみる。30人BioPythonが最も多くの人が関わっている。次に、直近1年の更新頻度を見てみる。BioPython:BioPerl:BioRuby:縦軸に注意してもらいたい。圧倒的にBioPythonの更新頻度が多くなっている。ユーザー数の多い言語のライブラリーほど、多くの人が関わり、更新頻度が高くなるので、機能が充実している&バグが少ない可能性が高い。まとめ以上のように、バイオインフォマティクスをこれから始める人にふさわしい言語はユーザー数の多いPythonであると考えられる。過去に開発されたプログラムはPerlだからPerlを使うべきという人がいるが、進歩の早いこの分野においては、10年前のプログラムは使い物にならないことが多い。Rubyを勧める人の言い分で「Rubyは日本語ドキュメントが充実している」という人がいるが、主にWEB開発に関するドキュメントであり、バイオインフォマティクスにはほとんど役に立たない。今後他の言語が主流になる可能性は否定しないが、現状としてはPythonをお勧めする。NanoPore等)を使用することが多い。第3世代シーケンサーではリードにエラーが多く含まれており、そのリードをもとにアセンブルすると、遺伝子領域に多くのギャップが含まれた結果が得られてしまう。第2世代シーケンサーの結果を使用してアセンブル結果のエラー補正を実施した場合でも、多くのギャップが残っている。バイオインフォマティクスロングリードでのアセンブルは遺伝子予測結果にエラーが多く含まれるにコメントメタゲノム解析メタゲノム解析とはメタゲノム解析とは、微生物群を直接シーケンスし、網羅的に菌種や遺伝子を特定する解析のことです。糞便や皮膚、海洋や土壌等といったサンプルを対象として研究が進められています。rRNAメタゲノム解析」の2種類があります。現在は主にショートリードシーケンサーが使用されていますが、今後はロングリードシーケンサーを使用することで、より精度の高い解析が可能になることが期待されます。ショットガンメタゲノム解析ショットガンメタゲノム解析におけるバイオインフォマティクス解析の手法は、大きく分けて2つあります。1つはリファレンス配列を使用せずアセンブリを行う手法で、「MetaSPAdes」や「MEGAHIT」といったソフトウェアがあります。1つは、リファレンス配列に対してマッピングをする手法で、「MetaPhlAn」といったソフトウェアがあります。ヒトの糞便サンプル等といった研究が比較的進んでいる分野については、マッピングの手法が使われることが多いようです。rRNAのみをシーケンス対象とすることで、解像度は落ちる分、低コストで解析を実施することが可能です。rRNAメタゲノム解析には「QIIME」といったソフトウェアがよく使われます。バイオインフォマティクスメタゲノム解析にコメントバイオインフォマティクス解析でよく使われるファイル形式バイオインフォマティクス解析では様々なファイル形式が使用されています。ここでは、主要なファイル形式を紹介します。配列をマッピングした結果を記述するための形式。Mb1/4の染色体が1本のScaffold、66%の染色体が1本もしくは2本のScaffoldにまとまった。Systemを使用して植物ゲノムをアセンブルにコメントバイオインフォマティクス解析バイオインフォマティクスとはバイオインフォマティクス(Bioinformatics)は、生命情報科学と訳されることもあり、生物学のデータを情報科学の手法によって解析する学問および技術です。生物学において大量のデータが生み出されるようになるとともに、この分野は急速に発展してきました。近年では次世代シーケンサー(NGS)を使用したシーケンスコストの低下が進み、ますます多くのデータが出力されるようになっています。これらのビッグデータを活用できる人材「バイオインフォマティシャン」の育成が急務であると言われています。Sequencingを略してNGSとも呼ばれます。従来のシーケンサー(サンガー法、第一世代シーケンサー)と比べて、1回のシーケンスで大量の配列情報が得られることが特徴です。一口に次世代シーケンサーと言っても、様々なシーケンサーが各社から販売されています。Technologies等が挙げられます。Technologies等の第三世代シーケンサーに分類できます。第二世代シーケンサーは読み取れる塩基の配列長(リード長)は数百塩基程度と短いですが、比較的高い精度で大量の配列情報を読めることが特徴です。第三世代シーケンサーはデータ量および精度の面で第二世代シーケンサーと比べて劣りますが、長いリードを取得できることが特徴です。第三世代シーケンサーIllumina第二世代シーケンサーに分類されるIlluminaシーケンサーは現在最も多くのシェアを獲得しているシーケンサーです。NovaSeq等、多くの製品が展開されており、それぞれ出力できるデータ量や適したアプリケーションが異なります。Illuminaシーケンサーの特徴は、DNA断片を両端からシーケンスすることができることです。bp」のように表現されます。Biosiences社のシーケンサーは第三世代シーケンサーに分類されます。1分子リアルタイム(SMRT)シーケンシングと呼ばれる技術により、DNA合成の様子を観測することで塩基配列を決定します。IIはすでに販売が終了しています。Technologiesのシーケンサーは2014年に販売開始された比較的新しいシーケンサーで、第三世代シーケンサーに分類されます。ナノポア(ナノスケールの穴)の中を核酸が通過する際の電流の変化を計測することで、塩基配列を同定します。MinION、GridION、PromethION等の製品が展開されており、それぞれ出力できるデータ量が異なります。MinIONは手のひらサイズのシーケンサーであり、ノートパソコンから電源を供給できるため、実験環境の整っていない場所でもシーケンスが可能です。試薬の改良等により、リード長・データ量・精度の改善が見込まれます。novoアセンブリ解析とは、ゲノム配列がまだ決定されていない生物種について新規にゲノム配列を決定する解析です。novoアセンブリ解析は、リード配列をつなぎ合わせることでより長い配列を構築する解析です。Illuminaシーケンサー等から得られるショートリードを用いたアセンブリが多く行われていました。ショートリードを用いたアセンブリのソフトウェアには「SOAPdenovo」、「SPAdes」、「Edena」、「Velvet」等があります。novoアセンブリ解析はロングリードシーケンサーを使用して行うことが多くなってきています。アセンブリには「Canu」、「SMARTdenovo」といったソフトウェアが使われます。ロングリードのみを用いたアセンブリではエラーが多くなる傾向にあるため、アセンブリの結果についてショートリードを用いてエラー補正することが多く行われます。これはポリッシングと呼ばれるステップで「Pilon」や「Nanopolish」といったソフトウェアが使われます。