如何利用HMMER鉴定基因家族成员

如何利用HMMER鉴定基因家族成员

通常我们用的都是通过blast比对来确定我们需要的家族成员,首先是比对序列,再次是需要目标物种的蛋白序列,来进行比对,通常比对的时候我们都需要设定e-value值。今天我们来学习一下利用HMMER来鉴定基因家族成员。

1 、利用软件HMMER进行模型搜索,
软件官网:http://hmmer.org/download.html ,主要用到hmmsearch,已经来自pfam数据库的隐马尔可夫模型。记住模型下载是在pfam数据库中,选择你自己需要的模型即可,如若没有你需要的模型怎么办呢,这时我们可以自己建立模型,请看图3后一行码。(所用到的数据都上传到网盘了,查看:http://132.232.42.33/?p=59

2、在pfam数据没有模型的时候,我们可以利用clustalX2要比对的序列进行比对,然后会陈胜*.aln文件(注:该软件的目标文件必须在桌面才可以添加)。


下载桌面版本:hmmer3.0_windows,将下载下来的版本解压到桌面 文件准备,拟南芥蛋白组数据,来自pfam数据库的模型,模型是怎样确定的,请参考:http://132.232.42.33/?p=105

将蛋白文件,模型文件放到解压的hmmer文件中

图 1

修改模型文件f—>b:

图 2
win+R-->cmd>cd desktop-->cd hmmer3.0_windows

如果你有Git的话,可以直接在该文件夹右键打开,也是一样的:

图 3

建立模型:

hmmbuild new.hmm *.aln

模型搜索并保存为wrky.txt文件:

图 4

打开wrky.txt文件结果如下:

图 5

如上图就见wrky基因鉴定出来了,如果还要进行下一步分析鉴定。包括多个转录本的情况,任然需要进行手动鉴定。以防出现假。以后得分享我基本上都会把入手的数据全部上传。

欢迎关注微信公众号:蚂蚁生信

如若有问题,或侵犯您的版权请及时联系我,谢谢:smyang2018@gmail. com

smyang2018

跳至工具栏