ログイン 概要 よくある質問

Hisat2のindexを作る作業が動きません。 このようなエラーがでます。

$ hisat2-build /Mus_musculus/UCSC/mm10/Sequence/WholeGenomeFasta/genome.fa

genome_index

dyld: Symbol not found:
__ZNSt7__cxx1112basic_stringIcSt11char_traitsIcESaIcEED1Ev

Referenced from: /usr/local/Cellar/hisat2/2.1.0/bin/hisat2-build-s Expected in: /usr/lib/libstdc++.6.0.9.dylib in /usr/local/Cellar/hisat2/2.1.0/bin/hisat2-build-s Abort trap: 6

いろいろ調べて試しましたが、意味がわかりませんでした。問題解決のヒントをご教示いただければ幸いです。 Brew update doctor upgradeは問題なかったです。

configです。

$ brew config
HOMEBREW_VERSION: 1.7.0
ORIGIN: https://github.com/Homebrew/brew
HEAD: 7e8fb9a0f8ab5e841763ec6c7fefa72a8462b594
Last commit: 5 days ago
Core tap ORIGIN: https://github.com/Homebrew/homebrew-core
Core tap HEAD: fa5253faa1d416851bbd3d957f67f5ccbcb15441
Core tap last commit: 6 hours ago
HOMEBREW_PREFIX: /usr/local
CPU: quad-core 64-bit skylake
Homebrew Ruby: 2.3.7 => /usr/local/Homebrew/Library/Homebrew/vendor/portable-ruby/2.3.7/bin/ruby
Clang: 9.1 build 902
Git: 2.18.0 => /usr/local/bin/git
Curl: 7.54.0 => /usr/bin/curl
Java: 10.0.2, 1.8.0_121
macOS: 10.13.6-x86_64
CLT: 9.4.1.0.1.1528165917
Xcode: 9.4.1
XQuartz: 2.7.11 => /opt/X11

どうぞよろしくお願いします。

質問日 Jul 21 '18 at 15:14

mikan03's gravatar image

mikan03
3135


環境がないので検証できないのですが、ひとまず公式のOS X用バイナリを使って動くか検証してみるのはどうでしょう? https://ccb.jhu.edu/software/hisat2/index.shtml ftp://ftp.ccb.jhu.edu/pub/infphilo/hisat2/downloads/hisat2-2.1.0-OSX_x86_64.zip

マウスのwhole genomeであれば公式でビルド済みのインデックスが配布されているのでそれを使うという手もあるかとおもいます。 ftp://ftp.ccb.jhu.edu/pub/infphilo/hisat2/data

回答日 Jul 21 '18 at 17:45

inutano's gravatar image

inutano
163

回答ありがとうございました。 結論からいうとたぶん解決しました。(今、ライブラリをビルドしています。) 一人で黙々やっていると思考が回らなくなるので、コメントに本当に感謝しています。

まずは結論から。 /usr/local/lib/libstdc++.6.0.9.dylib を削除しました。 (詳細は上げておきます。) これで動きました。

ありがとうございました。

(Jul 23 '18 at 14:03) mikan03 mikan03's gravatar image

マウスのindexもダウンロードしてみたのですが、最初はたくさんht1.ht2とかファイルがあったので、どれを指定したら良いかわからずに、後回しにしていました。 indexを自分で作ってみて、やっとファイルがたくさん意味がある意味がわかりました。indexは複数にわかれてできるんですね。出来上がったindex fileをフォルダ(mkdir index_mm10)にいれて、index_mm10/genome_mm10でindex指定して動きました。助かりました。

(Jul 23 '18 at 14:19) mikan03 mikan03's gravatar image

まずは結論から。 /usr/local/lib/libstdc++.6.0.9.dylib を削除しました。 ーーーーーーーーーー 以下、試したことを記述しておきます。(すみません。ログが流れてしまいましたので、部分的な情報だけです)

まず、brew uninstall hisat2をしてトラブル回避。

OSX用のバイナリからのビルドも試みましたが、makeがうまくいきませんでした。

unzip hisat2-2.1.0-OSX_x86_64.zip
cd hisat2-2.1.0
ls
/.configure
make

makeでconfigureが見つからないといわれましたので、 lsでみてconfig.というファイルがいろいろあったので、片っ端から/.configu.しましたが、特に動きなし。

続いて、他のパッケージングソフトでhisat2が入らないかを試す。 portはできなかったので、 minicondaをインストール。 sudo port install miniconda

minicondaでHIisat2をインストール which hisatでpathは通っているが、anaconda経由になっていた。 そこでもう一回hisat2でindex構築実行。

同じエラー発生。 hiast2 -hでヘルプはみられる。ということはシステム的な問題?

もういちどログ確認。 内容はあまり理解できませんが Expected in: /usr/lib/libstdc++.6.0.9.dylib どうやらこいつが悪さをしている?

下記参照サイトを元にggplotからlibstdc++.6.0.9.dylibを輸入してみるが.libstdc++.6.dylibしかない。試しに置き換えてみる。 githubからあたらしくlibstdc++.6.0.9.dylibを取ってきていれてみた。エラーログの内容が変わったが、相変わらず動かない。

libstdc++.6.0.9.dylibが消えるかを試したが、消えないことが判明。lsで表示されるのにfile not exist. rm optionいろいろ使ってみるが、どれも効果なし。finderからl/usr/local/lib//libstdc++.6.0.9.dylib(不可視ファイル)探してゴミ箱にぽいで消えた。

それからindex構築動きました。 mappingも動いて無事にでました。

以下参考にしたサイト。ありがとうございます。

https://bioconda.github.io/recipes/hisat2/README.html

http://blog.toshihiro-kaneko.tonkotsu.jp/?cid=3

https://github.com/phracker/MacOSX-SDKs/blob/master/MacOSX10.7.sdk/usr/lib/libstdc++.6.0.9.dylib

回答日 Jul 23 '18 at 14:10

mikan03's gravatar image

mikan03
3135

edited Jul 23 '18 at 14:14

動いたとのことでよかったです。(ライブラリを直接いじるとシステムに影響がないか心配ですが…)ちなみに、バイナリとはソースコードをコンパイルしたものなので、configure/make は必要ありません。unzip したあと、中に入っている hisat2 を直接実行できます。

回答日 Jul 25 '18 at 13:28

inutano's gravatar image

inutano
163

ありがとうございました。直接実行できるのも知りませんでした。記憶が曖昧なのですが、たしかそのまま実行したことがあって、そのときも同じエラーが出たので、configure/makeのせいなのか?と思い込んでいました。またわからないことがでるとかもしれませんが、そのときはどうぞよろしくお願いします。

(Jul 25 '18 at 13:44) mikan03 mikan03's gravatar image
あなたの回答
プレビューをトグルする

この質問をフォローする

By Email:

Once you sign in you will be able to subscribe for any updates here

By RSS:

回答

回答とコメント

タグ:

×1
×1
×1

質問日: Jul 21 '18 at 15:14

閲覧数: 5,548 回

最終更新日: Jul 25 '18 at 13:44

関係した質問

powered by OSQA