3-基因家族的鉴定-基于Windows系统上的HMMER

 

 

下面讲如何在windows系统上用HMMER做基因家族的鉴定：

1， 下载

我们电脑首先要有个迅雷，然后复制这个链接，用迅雷打开(http://hmmer.janelia.org/static/binaries/hmmer3.0_windows.zip)，下载到本地，解压。

 

2， 案例分析

HMMER是根据蛋白的结构域进行搜寻的，例如我们现在要找毛果杨的bZIP转录因子家族。在以往的文献中，我们知道拟南芥的bZIP蛋白序列，我们将拟南芥的bZIP蛋白序列提交到Pfam网站（http://pfam.xfam.org/）看看它的结构域。

我们批量提交并填写邮件，其它默认参数，可能要等半天甚至更长，结果会发至邮箱。

得到结果后，发现邮件里面有结果和结果的链接，点击链接可以去到网站上，结果可以也复制粘贴到excel中处理，在结果中，发现它们几乎都含有bZIP_1和bZIP_2结构域，结构域的编号分别是PF00170，PF07716。

点进去，下载stockholm格式的文件。

并将这两个结构域的文件，和毛果杨的全部蛋白文件都拷贝一份到HMMER的文件夹里面

 

在电脑的左下角打“cmd”，进入命令提示符

因为我的HMMER放在D盘中，路径是D:\HMMER，所以先进入D盘，输入命令

d:

再进入HMMER文件夹中，输入

cd HMMER

先对PF00170（bZIP_1）处理，输入命令

hmmbuild PF00170.hmm PF00170_seed.txt

再输入命令，搜索毛果杨中含有该结构域的蛋白序列

hmmsearch PF00170.hmm Ptrichocarpa_protein.fasta > PF00170.out

2秒算完，回到HMMER的文件夹中，发现了一些中间文件和PF00170.out的结果文件，可用excel打开处理。

PF07716（bZIP_2）也是同样的操作。

 具体的结果分析，自行百度即可，网上的教程很多。

pfam网站停用之后会转到interprot网站下，具体下载流程参考这个百度教程“https://www.omicsclass.com/article/1749”

其实window系统的HMMER已经很久没有维护了，建议用Linux系统的HMMER，而且更简单更容易理解，相关教程以后再出吧。

 

 

#花吻在上~