3、.. 如何获得其他物种已发表某个基因家族的所有成员呢,最简单的就是下载该物种蛋白序列文件(可以从上述数据库中下载),然后按照文章中的ID,找到对应成员。对于没有全基因组鉴定的,可以下列数据库中找: a. NCBI: nucleotide and protein db. b. EBI: http://www.ebi.ac.uk/. c. UniProtKB:http://www.uniprot.org/uniprot/2、比对工具。一般使用blast和hmmer,具体使用命令如下:nLocal BLASTformatdb–i d
4、b.fas–p F/T;blastall–p blastp(orelse) –i known.fas–d db.fas–m 8 –b 2(or else) e 1e-5 –o alignresult.txt.-b:output two different members in subject sequences (db).nHmmer (hidden Markov Model) search. Thesame as PSI-BLAST in function. It has a higher sensitivity, but the speed isl
6、sb/bwrpsb.cgi).nEST 支持n Blast and Hmmer同时检测到4、通过上述操作获得某家族的所有成员基因家族分析套路(二)本次主要讲解在基因家族分析类文章中,进化部分分析的内容。主要是进化树的构建与分析。一、构建进化树的基本步骤1、多序列比对. Muscle program.2、Model 选择. 分别针对蛋白序列和核酸序列的模型选择程序。ProtTest program for protein and ModelTest or Jmodetlest for DNA(http://user.qzone.qq.com/580017
7、04/blog).3、算法选择。三种. NJ, ML and BI.4、软件选择。 MEGA (bootstrap least 1000 replicates), phyML and Mrbayes (http://user.qzone.qq.com/58001704/main)....5、进化树修饰. MEGA: view->options and subtree-> draw options. Also can be decorated in word (http://user.qzone.qq.com/58001704/main)二、具体步骤 2