Name: 实用指南系统发育的非专家 (视频) | JoVE Journal
Uploaded: 2014-02-05T17:00:00.000+00:00
Duration: 12 min
Description: 了解如何使用BLAST程序识别相似序列并进行多序列比对，从而从DNA或蛋白质序列重建可靠的系统发育树。本方法为非专业人士提供完整的生物信息学分析流程。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

感谢O’Halloran实验室成员对本文稿的评论。感谢乔治华盛顿大学生物科学系以及哥伦布艺术与科学学院为D. O’Halloran提供的经费支持。

1. 基本局部比对搜索工具（BLAST）：在线界面
<ol>
	<li>点击此链接访问美国国家生物技术信息中心（NCBI）的 BLAST1 网络服务器。 - <a href="http://blast.ncbi.nlm.nih.gov/Blast.cgi">http://blast.ncbi.nlm.nih.gov/Blast.cgi</a>&#160;（图 1）。</li>
	<li>将一条 FASTA 格式的文本序列（示例见图 2）输入到查询框中。</li>
	<li>选择适当的 BLAST 程序以及用于搜索的相关数据库或特定物种，然后点击“BLAST”。 
	注意：FASTA 格式的序列以一条描述行开始，该行以“>”符号标识；描述内容必须紧随“>”符号之后，序列（即核苷酸或氨基酸）则位于下一行。BLAST 搜索结果可显示为 HTML、纯文本、XML 或命中表格（文本或 csv 格式），默认输出格式为 HTML（图 3）。</li>
</ol>
2. 基本局部比对搜索工具（BLAST）：本地可执行程序
<ol>
	<li>从以下链接下载最新的 BLAST 命令行可执行程序： 
	<a href="ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/">ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/</a> -</li>
	<li>em&#62;对于 PC 用户：双击最新的 blast win32.exe 文件，接受许可协议并点击安装。 
	注意：默认安装目录为 C:\ncbi-blast-2.2.27+。</li>
	<li>按如下方式配置 PC 环境变量：
	<ol>
		<li>点击 PC 的“开始”按钮，然后右键单击“计算机”；</li>
		<li>点击“属性”，在弹出窗口中点击“高级”选项卡；</li>
		<li>点击“环境变量”按钮，在新弹出窗口的“用户变量”部分点击“新建”按钮；</li>
		<li>在弹出窗口中添加变量名为“Path”，变量值为“C:\ncbi-blast-2.2.27+\bin。 
		注意： bin 目录包含可执行文件（即 blastp、&#160;等）。</li>
	</ol>
	</li>
	<li>em&#62;对于 Mac 用户：打开终端应用程序（只需打开“Finder”并搜索“Terminal”，即可显示“terminal”图标）。在终端窗口中输入： 
	&#62;ftp <a href="ftp://ftp.ncbi.nih.gov">ftp.ncbi.nih.gov</a> 
	注意：也可输入上述 PC 示例中使用的 URL</li>
	<li>访问 NCBI ftp 站点时，在 Name 和 Password 处输入“anonymous”，然后输入： 
	&#62;cd blast/executables/LATEST</li>
	<li>通过输入以下命令列出可执行文件： 
	&#62;ls</li>
	<li>输入以下命令（或当前最新的版本号）获取最新版本： 
	&#62;get ncbi-blast-2.2.7-macosx.tar.gz</li>
	<li>输入“exit”退出 NCBI ftp 服务器站点。</li>
	<li>通过输入以下命令解压已下载的文件： 
	&#62;tar -xzf ncbi-blast-2.2.7-macosx.tar.gz</li>
	<li>将 blast 可执行文件的二进制文件所在路径添加到您的 PATH 中，以便 shell 在查找命令时能搜索该目录，输入： 
	&#62;PATH=$PATH:new_folder_location</li>
	<li>通过输入以下命令检查是否已将路径添加至 PATH： 
	&#62;echo $PATH</li>
	<li>点击此处下载已预格式化的 BLAST 数据库（每日更新）： 
	<a href="ftp://ftp.ncbi.nlm.nih.gov/blast/db/">ftp://ftp.ncbi.nlm.nih.gov/blast/db/</a></li>
	<li>将数据库放入“db”文件夹中。</li>
	<li>em&#62;在 PC 上：打开 MS-DOS 命令提示符（点击“开始”并在搜索栏中输入“cmd”），然后输入以下命令将目录更改为 ncbi-blast 文件夹： 
	C:\Users&#62;cd ..\ [向上移动一个文件夹] 
	C:\&#62;cd ncbi-blast-2.2.27+ 
	这将把目录更改为： 
	C:\ncbi-blast-2.2.27+&#62;</li>
	<li>使用以下“makedb”命令创建数据库： 
	&#62;makedb –in db/briggsae.fasta –dbtype prot –out db/briggsae 
	注意： 在下文示例（图 4）中，数据库命名为“briggsae”，由生物体 Caenorhabditis briggsae 的一条连锁群组成。</li>
	<li>通过将一条 FASTA 格式的蛋白质序列文本插入“db”文件夹中，创建一个名为“test”的查询蛋白序列。</li>
	<li>通过输入以下命令执行 blastp 搜索以查询数据库： 
	&#62;blastp –query db/test.txt –db db/briggsae –out text.txt</li>
	<li>em&#62;在 Mac 上：通过按照上述说明（步骤 2.4）访问 NCBI ftp 网站，下载用于本地 Blast 搜索的数据库，然后输入： 
	&#62;lcd ../databases/</li>
	<li>输入以下命令下载感兴趣的基因组或序列： 
	&#62;get NC_[Accession #].fna 
	注意： “.fna”表示 FASTA 格式的核苷酸序列，“.faa”表示 FASTA 格式的氨基酸序列。</li>
	<li>输入“quit”退出 ftp 站点。</li>
	<li>输入以下命令创建数据库： 
	&#62;makeblastdb -in db/mouse.faa -out mouse -dbtype prot</li>
	<li>将一条 FASTA 格式的查询序列插入“bin”文件夹，并使用以下命令查询数据库： 
	&#62; blastp -query “your query.fasta” -db ”your database” -out results.txt</li>
</ol>
3. 生成多序列比对
<ol>
	<li>点击以下链接以访问常用的多序列比对（MSA）程序： 
	ClustalW6 <a href="http://www.clustal.org/">http://www.clustal.org/</a> 
	Kalign7 <a href="http://msa.sbc.su.se/cgi-bin/msa.cgi">http://msa.sbc.su.se/cgi-bin/msa.cgi</a> 
	MAFFT8,9 <a href="http://mafft.cbrc.jp/alignment/software/">http://mafft.cbrc.jp/alignment/software/</a> 
	MUSCLE10 <a href="http://www.drive5.com/muscle/">http://www.drive5.com/muscle/</a> 
	T-Coffee11 <a href="http://www.tcoffee.org/Projects/tcoffee/">http://www.tcoffee.org/Projects/tcoffee/</a> 
	PROBCONS12 <a href="http://toolkit.tuebingen.mpg.de/probcons">http://toolkit.tuebingen.mpg.de/probcons</a></li>
	<li>点击此链接 - <a href="http://tcoffee.crg.cat/apps/tcoffee/do:regular">http://tcoffee.crg.cat/apps/tcoffee/do:regular</a> - 并将FASTA格式的序列数据输入查询框中 
	注： T-Coffee的示例输出见图5，相似残基以颜色编码。</li>
	<li>点击此链接下载Clustal MSA的命令行版本（ClustalW）或图形化版本（ClustalX）：<a href="http://www.clustal.org/clustal2/">http://www.clustal.org/clustal2/</a> - 然后点击相应的可执行文件（例如 win、Linux、Mac OS X）。</li>
	<li>上传FASTA格式的序列文本数据并进行比对（图6）。</li>
</ol>
4. 确定最优进化模型
<ol>
	<li>单击此处下载 ProtTest13 程序： 
	<a href="http://darwin.uvigo.es/our-software/">http://darwin.uvigo.es/our-software/</a></li>
	<li>ProtTest 下载完成后，双击 ProtTest.jar 文件</li>
	<li>启动 ProtTest 后，单击“select file”并加载序列数据（图7).</li>
	<li>然后点击“start”，程序将开始运行（图8). 
	注意： 运行结束后（图8），程序将根据标准指示最佳模型 例如 根据AIC的最佳模型：WAG+I+G</li>
</ol>
5. 基于最大似然法或贝叶斯推断推断序列系统发育树
<ol>
	<li>在此处下载 PhyML4： 
	<a href="https://code.google.com/p/phyml/">https://code.google.com/p/phyml/</a></li>
	<li>双击相应的应用程序（例如 phyml Windows、phyml Linux，&#160;等等）启动可执行文件，界面窗口将弹出（图9）。</li>
	<li>通过输入以下命令加载 PHYLIP 格式的输入序列： 
	>“文件名”.phy 
	注意： 如需在不同序列格式之间转换，请使用位于以下网址的“Readseq”网络程序：<a href="http://iubio.bio.indiana.edu/cgi-bin/readseq.cgi">http://iubio.bio.indiana.edu/cgi-bin/readseq.cgi</a>。</li>
	<li>输入“Y”以启动程序。</li>
	<li>在此处下载 MrBayes5： 
	<a href="http://mrbayes.sourceforge.net/download.php">http://mrbayes.sourceforge.net/download.php</a></li>
	<li>单击可执行文件以启动程序，并通过输入以下命令将 NEXUS 格式的序列数据读入程序： 
	&#62;execute “文件名”.nex</li>
	<li>设定进化模型。</li>
	<li>通过输入以下命令设置运行的代数： 
	&#62;mcmcp ngen = 1000000 [此命令将运行代数设为 1000000] 
	&#62;sump burnin =10000 [此命令将 burnin 设为 10000]</li>
	<li>通过输入以下命令在结果文件中保存分支长度： 
	&#62;mcmcp savebrlens = yes&#160;</li>
	<li>输入以下命令运行分析： 
	&#62;mcmc</li>
	<li>使用“sumt”命令对生成的树进行汇总。</li>
</ol>
6. 系统发育树的可视化
<ol>
	<li>在此查看树形图查看程序列表： 
	<a href="http://www.treedyn.org/overview/editors.html">http://www.treedyn.org/overview/editors.html</a></li>
	<li>在此下载 TreeView14 程序： 
	<a href="http://taxonomy.zoology.gla.ac.uk/rod/treeview.html">http://taxonomy.zoology.gla.ac.uk/rod/treeview.html</a></li>
</ol>

<ol><li>Altschul, S. F., Carroll, R. J., Lipman, D. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Weights+for+data+related+by+a+tree%5BTitle%5D)+AND+%22J.+Mol.+Biol%22%5BJournal%5D)">Weights for data related by a tree.</a> J. Mol. Biol. 207 (4), 647-653 (1989).</li>
<li>Akaike, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+new+look+at+the+statistical+model+identification%5BTitle%5D)+AND+%22IEEE+Trans.+Automat.+Contr%22%5BJournal%5D)">A new look at the statistical model identification.</a> IEEE Trans. Automat. Contr. 19 (6), 706-723 (1974).</li>
<li>Schwarz, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Estimating+the+dimension+of+a+model%5BTitle%5D)+AND+%22Ann.+Stat%22%5BJournal%5D)">Estimating the dimension of a model.</a> Ann. Stat. 6 (2), 461-464 (1978).</li>
<li>Guindon, S., Gascuel, O. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+simple%2C+fast%2C+and+accurate+algorithm+to+estimate+large+phylogenies+by+maximum+likelihood%5BTitle%5D)+AND+%22Syst.+Biol%22%5BJournal%5D)">A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood.</a> Syst. Biol. 52 (5), 696-704 (2003).</li>
<li>Huelsenbeck, J. P., Ronquist, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MRBAYES%3A+Bayesian+inference+of+phylogenetic+trees%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">MRBAYES: Bayesian inference of phylogenetic trees.</a> Bioinformatics. 17 (8), 754-755 (2001).</li>
<li>Thompson, J. D., Higgins, D. G., Gibson, T. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((CLUSTAL+W%3A+Improving+the+sensitivity+of+progressive+multiple+sequence+alignment+through+sequence+weighting%2C+position-specific+gap+penalties+and+weight+matrix+choice%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice.</a> Nucleic Acids Res. 22 (22), 4673-4680 (1994).</li>
<li>Lassmann, T., Sonnhammer, E. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Kalign--an+accurate+and+fast+multiple+sequence+alignment+algorithm%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">Kalign--an accurate and fast multiple sequence alignment algorithm.</a> BMC Bioinformatics. 6, 298 (2005).</li>
<li>Katoh, K., Kuma, K., Toh, H., Miyata, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MAFFT+version+5%3A+Improvement+in+accuracy+of+multiple+sequence+alignment%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">MAFFT version 5: Improvement in accuracy of multiple sequence alignment.</a> Nucleic Acids Res. 33 (2), 511-518 (2005).</li>
<li>Katoh, K., Misawa, K., Kuma, K., Miyata, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MAFFT%3A+A+novel+method+for+rapid+multiple+sequence+alignment+based+on+fast+fourier+transform%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">MAFFT: A novel method for rapid multiple sequence alignment based on fast fourier transform.</a> Nucleic Acids Res. 30 (14), 3059-3066 (2002).</li>
<li>Edgar, R. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MUSCLE%3A+Multiple+sequence+alignment+with+high+accuracy+and+high+throughput%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">MUSCLE: Multiple sequence alignment with high accuracy and high throughput.</a> Nucleic Acids Res. 32 (5), 1792-1797 (2004).</li>
<li>Notredame, C., Higgins, D. G., Heringa, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((T-coffee%3A+A+novel+method+for+fast+and+accurate+multiple+sequence+alignment%5BTitle%5D)+AND+%22J.+Mol.+Biol%22%5BJournal%5D)">T-coffee: A novel method for fast and accurate multiple sequence alignment.</a> J. Mol. Biol. 302 (1), 205-217 (2000).</li>
<li>Do, C. B., Mahabhashyam, M. S., Brudno, M., Batzoglou, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((ProbCons%3A+Probabilistic+consistency-based+multiple+sequence+alignment%5BTitle%5D)+AND+%22Genome+Res%22%5BJournal%5D)">ProbCons: Probabilistic consistency-based multiple sequence alignment.</a> Genome Res. 15 (2), 330-340 (2005).</li>
<li>Darriba, D., Taboada, G. L., Doallo, R., Posada, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((ProtTest+3%3A+Fast+selection+of+best-fit+models+of+protein+evolution%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">ProtTest 3: Fast selection of best-fit models of protein evolution.</a> Bioinformatics. 27 (8), 1164-1165 (2011).</li>
<li>Page, R. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((TreeView%3A+An+application+to+display+phylogenetic+trees+on+personal+computers%5BTitle%5D)+AND+%22Comput.+Appl.+Biosci%22%5BJournal%5D)">TreeView: An application to display phylogenetic trees on personal computers.</a> Comput. Appl. Biosci. 12 (4), 357-358 (1996).</li>
<li>Darriba, D., Taboada, G. L., Doallo, R., Posada, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((jModelTest+2%3A+More+models%2C+new+heuristics+and+parallel+computing%5BTitle%5D)+AND+%22Nat.+Methods%22%5BJournal%5D)">jModelTest 2: More models, new heuristics and parallel computing.</a> Nat. Methods. 9 (8), 772 (2012).</li>
<li>Chevenet, F., Brun, C., Banuls, A. L., Jacq, B., Christen, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((TreeDyn%3A+Towards+dynamic+graphics+and+annotations+for+analyses+of+trees%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">TreeDyn: Towards dynamic graphics and annotations for analyses of trees.</a> BMC Bioinformatics. 7, 439 (2006).</li>
</ol>

我们没有需要披露的内容。&#160;

本文旨在为系统发育学领域的新手研究人员或学生提供一个入门起点。近年来，基因组测序项目的成本不断降低，导致对该技术的用户需求持续增长，如今小型实验室生成大规模序列数据集已成为常态。这些数据集通常为研究人员提供一系列基因，而要初步理解其功能，则需要借助系统发育框架。此外，由于系统发育学正被越来越多的研究实验室所采用，本文也希望作为对生物学研究广泛感兴趣的学生的一种教育工具。通过向用户提供常用建树工具在“为何使用”、“如何使用”以及“在何处使用”方面的信息，我们为读者熟悉这些应用及其工作原理提供了基本框架。然而，我们建议读者尝试调整每种工具中的各项设置，以理解不同参数如何影响其序列数据，并确保在每种情况下平台与软件之间的兼容性。上述分析是在一台配备 Intel Core i7 处理器的 Dell Optiplex 990 台式机和一台配备 Intel Core 2 Duo 处理器的 MacBook 笔记本电脑上完成的，但分析速度以及具体的二进制版本（例如 32 位或 64 位）将取决于用户的平台。
编写系统发育学这类用户指南面临的一个挑战是，系统发育学以及整个生物信息学...

为了理解两个（或更多）物种是如何进化的，首先需要从每个样本中获取序列或形态学数据；这些数据代表了可用于衡量它们在进化空间中相互关系的量化信息。就像测量线性距离一样，可获得的数据越多（例如 英里、英寸、微米），测量结果就越精确。因此，研究人员推断进化距离的准确性在很大程度上取决于可用于衡量关系的信息数据量。此外，由于不同样本以不同的速率和不同的机制进化，我们用来衡量两个分类单元之间关系的方法也直接影响进化测量的准确性。&#160;因此，由于进化关系并非直接观察所得，而是从序列或形态学数据中推断出来的，推断进化关系的问题就成为一个统计学问题。系统发育学是生物学的一个分支，致力于将统计模型应用于进化模式，以最优地重建分类单元之间的进化历史。这种分类单元之间的重建称为该分类单元的系统发育关系（phylogeny）。
为了帮助弥合分子生物学家与进化生物学家之间的专业知识差距，本文介绍了一套从序列集合推断系统发育关系的逐步操作流程。首先，我们详细描述了如何通过网络界面以及使用本地可执行程序，利用基本局部比对搜索工具（Basic Local Alignment Search Tool, BLAST1）算法进行数据库检索的各个步骤；这通常是获取与未知查询序列相似的一组序列列表的第一步。尽管如此，部分研究人员也可能希望通过Phylota等网络接口（http://www.phylota.net/）收集特定类群的数据。BLAST是一种用于将原始氨基酸或核苷酸序列数据与序列数据库进行比对，以搜索与查询序列相似的“匹配”序列的算法。该程序由美国国立卫生研究院（National Institutes of Health, NIH）的Stephen Altschul et al. 设计开发1。BLAST服务器包含多个不同的程序，以下是一些最常用的BLAST程序列表：
i) 核苷酸-核苷酸BLAST（blastn）：该程序需要输入一段DNA序列，并返回用户指定的DNA数据库中最相似的DNA序列（例如，针对某一特定生物体）。
ii) 蛋白质-蛋白质BLAST（blastp）： 此处用户输入一条蛋白质序列，程序将从用户指定的蛋白质数据库中返回最相似的蛋白质序列。
iii) 位点特异性迭代 BLAST (PSI-BLAST) (blastpgp)：用户输入一条蛋白质序列，程序返回一组密切相关蛋白质，并基于该数据集生成一个保守的序列谱型。随后，仅使用这些保守的“基序”生成新的查询序列，用于搜索蛋白质数据库，从而获得更大范围的蛋白质序列；从中提取出新的一组保守“基序”，再次用于数据库搜索。此过程不断重复，直至获得更广泛的蛋白质序列并生成新的谱型。通过在每一步迭代中将相关蛋白质纳入查询序列，该程序可帮助用户识别出序列差异更大的同源序列。
iv) 核苷酸六框翻译-蛋白质（blastx）： 此处用户提供核苷酸序列输入，该序列被转换为六框概念翻译产物（即两条链）并与蛋白质序列数据库进行比对。
v) 核苷酸六框翻译-核苷酸六框翻译（tblastx）：该程序以DNA核苷酸序列为输入，将输入序列翻译成所有六个阅读框的理论翻译产物，并将其与核苷酸序列数据库的六框翻译产物进行比对。
vi) 蛋白质-核苷酸六框翻译（tblastn）：该程序使用蛋白质序列作为输入，与核苷酸序列数据库的所有六个阅读框进行比对。
接下来，我们介绍从序列数据集中生成多序列比对（Multiple Sequence Alignment, MSA）的常用程序，随后提供用于确定序列数据集最优进化模型的程序使用指南。系统发育重建是一个统计学问题，因此系统发育分析方法需要纳入统计学框架。该统计学框架即为进化模型，用于描述数据集中序列的变化过程。该进化模型包含一组关于核苷酸或氨基酸替换过程的假设，针对特定数据集的最优模型可通过统计检验来选择。通过似然比检验（LRTs）或信息准则，可比较不同模型对数据的拟合优度，从而在一组可能的模型中选出最优模型。两种常用的信息准则为赤池信息准则（Akaike information criterion, AIC）2和贝叶斯信息准则（Bayesian information criterion, BIC）3。在获得最优比对结果后，可采用多种不同方法基于比对数据构建系统发育树。推断进化关系的方法众多，大致可分为两类：基于距离的方法和基于序列的方法。基于距离的方法通过计算序列间的成对距离，再利用这些距离构建系统发育树；基于序列的方法则直接使用序列比对结果，通常依据某种最优性标准搜索树空间。本文概述两种基于序列的系统发育关系重建方法：其一是实现最大似然法框架的 PhyML4，其二是采用贝叶斯马尔可夫链蒙特卡洛（Bayesian Markov Chain Monte Carlo）推断的 MrBayes5。最大似然法与贝叶斯方法为系统发育重建提供了统计学框架。通过向读者介绍常用的建树工具及其使用方法，我们阐明了推断系统发育关系所需的关键数据信息。

<table><tbody><tr><td>BLAST 网页&amp;nbsp;</td><td></td><td></td><td>&lt;a href=&quot;http://blast.ncbi.nlm.nih.gov/Blast.cgi">http://blast.ncbi.nlm.nih.gov/Blast.cgi&lt;/a&gt;</td></tr><tr><td>BLAST 可执行程序&amp;nbsp;</td><td></td><td></td><td>&lt;a href=&quot;ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/">ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/&lt;/a&gt;</td></tr><tr><td>预格式化 BLAST 数据库</td><td></td><td></td><td>&lt;a href=&quot;ftp://ftp.ncbi.nlm.nih.gov/blast/db/">ftp://ftp.ncbi.nlm.nih.gov/blast/db/&lt;/a&gt;</td></tr><tr><td>Clustal</td><td></td><td></td><td>&lt;a href=&quot;http://www.clustal.org/">http://www.clustal.org/&lt;/a&gt;</td></tr><tr><td>Kalign</td><td></td><td></td><td>&lt;a href=&quot;http://msa.sbc.su.se/cgi-bin/msa.cgi">http://msa.sbc.su.se/cgi-bin/msa.cgi&lt;/a&gt;</td></tr><tr><td>MAFFT</td><td></td><td></td><td>&lt;a href=&quot;http://mafft.cbrc.jp/alignment/software/">http://mafft.cbrc.jp/alignment/software/&lt;/a&gt;</td></tr><tr><td>MUSCLE</td><td></td><td></td><td>&lt;a href=&quot;http://www.drive5.com/muscle/">http://www.drive5.com/muscle/&lt;/a&gt;</td></tr><tr><td>T-Coffee</td><td></td><td></td><td>&lt;a href=&quot;http://www.tcoffee.org/Projects/tcoffee/">http://www.tcoffee.org/Projects/tcoffee/&lt;/a&gt;</td></tr><tr><td>PROBCONS</td><td></td><td></td><td>&lt;a href=&quot;http://toolkit.tuebingen.mpg.de/probcons">http://toolkit.tuebingen.mpg.de/probcons&amp;nbsp;&lt;/a&gt;</td></tr><tr><td>Se-Al&amp;nbsp;</td><td></td><td></td><td>&lt;a href=&quot;http://tree.bio.ed.ac.uk/software/seal/">http://tree.bio.ed.ac.uk/software/seal/&lt;/a&gt;</td></tr><tr><td>BSEdit&amp;nbsp;</td><td></td><td></td><td>&lt;a href=&quot;http://www.bsedit.org/">http://www.bsedit.org/&lt;/a&gt;</td></tr><tr><td>JalView</td><td></td><td></td><td>&lt;a href=&quot;http://www.jalview.org/">http://www.jalview.org/&lt;/a&gt;</td></tr><tr><td>SeaView</td><td></td><td></td><td>&lt;a href=&quot;http://pbil.univ-lyon1.fr/software/seaview.html">http://pbil.univ-lyon1.fr/software/seaview.html&lt;/a&gt;</td></tr><tr><td>ProtTest&amp;nbsp;</td><td></td><td></td><td>&lt;a href=&quot;https://code.google.com/p/prottest3/">https://code.google.com/p/prottest3/&lt;/a&gt;</td></tr><tr><td>Java 运行环境&amp;nbsp;</td><td></td><td></td><td>&lt;a href=&quot;http://www.java.com/en/download/chrome.jsp">http://www.java.com/en/download/chrome.jsp&lt;/a&gt;</td></tr><tr><td>Readseq</td><td></td><td></td><td>&lt;a href=&quot;http://iubio.bio.indiana.edu/cgi-bin/readseq.cgi">http://iubio.bio.indiana.edu/cgi-bin/readseq.cgi&lt;/a&gt;</td></tr><tr><td>jModelTest</td><td></td><td></td><td>&lt;a href=&quot;https://code.google.com/p/jmodeltest2/">https://code.google.com/p/jmodeltest2/&lt;/a&gt;</td></tr><tr><td>PhyML</td><td></td><td></td><td>&lt;a href=&quot;https://code.google.com/p/phyml/">https://code.google.com/p/phyml/&lt;/a&gt;</td></tr><tr><td>MrBayes</td><td></td><td></td><td>&lt;a href=&quot;http://mrbayes.sourceforge.net/download.php">http://mrbayes.sourceforge.net/download.php&lt;/a&gt;</td></tr><tr><td>TreeView</td><td></td><td></td><td>&lt;a href=&quot;http://taxonomy.zoology.gla.ac.uk/rod/treeview.html">http://taxonomy.zoology.gla.ac.uk/rod/treeview.html&lt;/a&gt;</td></tr><tr><td>TreeDyn</td><td></td><td></td><td>&lt;a href=&quot;http://www.treedyn.org/">http://www.treedyn.org/&lt;/a&gt;</td></tr></tbody></table>

适用于非专业人士的系统发育学实用指南

A Practical Guide to Phylogenetics for Nonexperts

许多研究人员在极为多样的研究领域中都将系统发育学应用于其研究问题。然而，许多研究人员对该主题尚不熟悉，因此会面临固有的困难。本文为非专业人士汇编了一份系统发育学的实用入门指南。&#160;我们以逐步说明的方式，介绍了一套从基因序列数据集生成可靠系统发育树的分析流程。首先提供通过在线界面及本地可执行程序进行相似性搜索工具的用户指南；接着探讨用于生成多序列比对的程序，随后介绍使用软件确定最优进化模型的方案；然后概述通过最大似然法和贝叶斯准则重建系统发育关系的协议&#160;，最后描述可视化系统发育树的工具。尽管本文远非对系统发育方法的全面综述，但它为读者提供了系统发育学家常用关键软件应用的实用入门信息。本文的目标是为刚开始进行系统发育研究的研究人员提供一种实用的培训工具&#160;，同时也可作为可融入课堂或教学实验室的教育资源。

寻找与查询序列的相似性可使研究人员为新序列赋予潜在身份，并推断序列之间的关系。BLAST1的文件输入类型为FASTA格式的文本序列或GenBank登录号。FASTA格式的序列以一个由“>”符号标记的描述行开头（图2）。描述内容必须紧接在“>”符号之后，序列（即核苷酸或氨基酸）则位于下一行。在保存和编辑序列文件时，最好使用文本编辑器，例如PC上的“记事本”或Mac上的TextWrangler（<a href="http://www.barebones.com/products/textwrangler/">http://www.barebones.com/products/textwrangler/</a>）。BLAST算法执行“局部”比对，即搜索短片段序列的相似性。算法首先查找查询序列中所有可能的“片段”，并尽可能延伸这些序列，然后为每对查询序列构建比对结果。理解这些比对结果的可靠性至关重要，因此BLAST对每个匹配结果应用统计学评估，包括期望值（E值）和比特分值（bit score）。E值用于指示匹配结果的统计学显著性，E值越低，匹配结果越显著。例如，E值为0....

本文介绍了一种从核苷酸或氨基酸序列数据集中生成可靠系统发育树的逐步分析流程。本指南旨在为初次接触系统发育分析的研究人员或学生提供帮助。&#160;

许多研究人员，尽管研究领域极为广泛多样，正将其研究问题应用于系统发育学分析。然而，许多研究人员对此方法尚不熟悉……

本文的总体目标是从 DNA 或蛋白质序列中重建一个可靠的系统发育树。首先，通过在 NCBI 使用 blast 程序识别相似序列来实现这一目标。第二步是比对这些相似序列。接下来，根据比对结果确定最合适的进化模型。最后一步是基于比对后的序列推断系统发育关系。最终，这一逐步分析流程展示了用户如何从序列数据获得可靠的系统发育树。该方法可通过推断新序列的身份和功能，帮助回答多个领域的关键问题。要使用基础局部比对搜索工具（BLAST）的在线版本，请访问国家生物技术信息中心（NCBI）的BLAST网络服务器，然后点击相应的BLAST程序。在查询框中输入一个FASTA格式的文本序列，例如此处所示的序列。点击适用于搜索的相应BLAST程序，然后点击BLAST。默认输出为HTML格式，显示与输入文本序列最相似的序列。下一节介绍如何在 Windows 和 Mac 系统上使用本地 blast 可执行程序。用户可跳至名为“Blast Local executables for Macs x”的下一节。要在 Windows 计算机上运行 blast 命令行程序，请从 NCBI blast 网站下载适用于 Windows 的可执行文件。安装 Blast 程序后，按如下方式配置计算机的环境变量：单击计算机的“开始”按钮，右键单击“计算机”，然后单击“属性”。在新窗口中，选择“高级系统设置”，在弹出窗口的“高级”选项卡中，单击“环境变量”按钮。然后在用户部分的用户变量下，点击新建按钮。在弹出的新窗口中，添加变量名 path 以及此处显示的变量值。接下来，从 NCBI 网站下载每日更新的预格式化 blast 数据库，或下载特定生物体的基因组。然后，通过点击“开始”并在搜索栏中键入 CMD 来打开 MS DOS 命令提示符，接着切换到 NCBI blast 文件夹。使用此处显示的 Make blast DB 命令创建数据库。通过将一条 FASTA 格式的蛋白质序列文本插入 DB 文件夹中，创建一个名为 test 的查询蛋白序列。然后，为了鉴定与测试蛋白序列最相似的序列，可通过 blast P 查询命令检索数据库。以下部分将针对 Mac 用户重复此信息。Windows 用户可跳至第五部分，生成多序列比对。要在 Mac 上运行 blast 命令行程序，需通过远程访问 N-C-B-I-F-T-P 网站下载适用于 MAC 的可执行文件。操作方法是：打开 Finder，在终端窗口中搜索 terminal，然后在终端中输入 N-C-B-I-F-T-P 网站的 FTP 地址。用户名和密码均输入 anonymous，然后输入命令 CD blast slash executables slash latest。通过输入 LS 列出可执行文件，并通过键入以下命令下载符合您系统要求的最新版本。现在解压缩已下载的文件。接着将 blast 可执行文件的二进制程序所在路径添加到您的环境变量中，以便 shell 能够搜索该目录。在查找命令时，从 NCBI 网站下载预格式化的 blast 数据库或基因组。通过输入 CD genomes 命令进入 genomes 目录。然后按照以下步骤下载感兴趣的基因组或序列，接着输入 quit 退出 FTP 站点。接下来，通过输入 Make Blast DB 指令来创建数据库。将 FASTA 格式的查询序列插入到 bin 文件夹中，并使用 blast P 查询命令对数据库进行检索，以在常用的多重序列比对（MSA）程序中找到与测试序列数据最相似的序列，其中一种常用程序是 tea coffee。在将 FASTA 格式的序列数据输入至 tea coffee 网站的查询框后，输出结果会通过颜色编码标示出相似的残基。另一种常用的多序列比对（MSA）程序是 Clustal MSA，该程序可下载为命令行版本 CLUSTERAL W，或适用于多种操作系统的图形化版本 CLUSTERAL X。接下来，通过选择文件选项卡，以 FASTA 格式的序列文本形式将数据加载到 Clustal 程序中，然后点击加载序列按钮。现在切换到比对选项卡，点击“执行完整比对”按钮，以对序列进行比对并获得最优进化模型。下载ProTest程序。ProTest下载完成后，双击启动ProTest。启动分析后，单击比对框中的“选择文件”以加载序列数据。然后单击“开始”运行程序。运行完成后，程序将根据推断序列的标准推荐最佳模型。下载并启动 Phi ML 后，通过输入文件名和 PY 来加载输入序列，将其作为文件脂质格式序列。然后通过键入 y 启动程序。从 Mr Bays 网站下载贝叶斯推断程序后，通过单击可执行文件来启动程序。然后通过输入 execute 文件名 dot NEX 将 Nexus 格式的序列数据读入程序中。接下来，设置进化模型并选择运行的世代数。使用 mc mc 命令运行分析后，通过 sum T 命令汇总生成的树以查看系统发育树。下载树状图查看程序。最后需要说明的是，不断有新的软件发布，旨在提供更优的序列比对、相似性预测或系统发育树构建。尽管本视频概述了一些常用程序，但仍鼓励观众探索其他可选工具。BLAST 算法执行局部比对，用于搜索序列中较短的相似片段。在算法查找了查询序列的所有可能片段并进行最大程度的延伸后，便会组装出比对结果。对于每一对查询序列，e 值反映了匹配结果的统计学显著性。E 值越低，表明匹配结果越显著。例如，E 值为 0.05 的序列比对意味着仅凭随机机会出现该匹配结果的可能性为 100 次中出现 5 次。BIT 分数使用特定的打分矩阵来评估比对质量的优劣，BIT 分数越高，表示比对结果越好。多序列比对（MSA）是指对三个或更多由氨基酸、DNA 或 RNA 组成的初级序列进行的序列比对。此处显示的 MSA 结果以颜色编码标记相似的残基。以下展示的是使用 Cluster X 对六个蛋白质序列进行比对后得到的氨基酸序列比对示例。protest 程序用于确定氨基酸替换的最佳拟合模型。在分析过程中，该程序会列出正在分析的模型，并在程序完成后显示最佳拟合结果。Phi ML 可根据核苷酸或氨基酸序列比对结果，估算最大似然系统发育树。该程序集成了大量替换模型，并结合多种选项用于搜索树拓扑空间。MrBayes 利用贝叶斯CMC推断方法，在多种进化模型下重建系统发育关系。程序运行后，可以按特定时间间隔查看分析进度，如图所示。系统发育树生成后，需要对拓扑结构进行可视化。在此图中，树状图窗口显示了来自 fly.Base 的蛋白质示例树。树状图包含一个树编辑器，允许用户移动分支并重新调整树的结构。在执行此操作时，重要的是要仔细阅读每个程序的用户指南。本方案为读者提供了实际的入门起点，介绍了这些程序的工作原理。然而，我鼓励读者自行尝试并熟悉每个程序相关的各种设置。

a-practical-guide-to-phylogenetics-for-nonexperts

Research

科研

JoVE Journal

Biology

生物学

实用指南系统发育的非专家

ITS2数据库是一个用于系统发育推断的工作平台，可同时分析内转录间隔区2（ITS2）的序列及其二级结构。该平台包含数据收集与精确注释、结构预测、序列-结构多序列比对以及快速建树等功能。简而言之，该工作平台将初步的系统发育分析简化为几次简单的点击操作。...

ITS2 数据库

本实验方案的目的是为缺乏统一命名和分类标准的蛋白质群组建立一个参考体系。该参考体系将有助于对整个群组进行分析和讨论，并可与已确立的名称一并使用。...

创建并应用参考体系以促进对多样化蛋白质组的讨论与分类

本文描述了一种基于真核生物SWEET蛋白与原核生物SemiSWEET蛋白序列同源性构建系统发育树的方法。系统发育分析是揭示不同生物类群间同源蛋白或基因进化关系的有力工具。...

利用系统发育分析研究真核生物基因的起源

本实验方案的目的是利用RNA测序数据研究候选基因的进化与表达。...

利用RNA-seq研究分子进化与基因表达的生物信息学分析流程

本文所述方案旨在指导用户对从甲虫中获得的细胞色素c氧化酶亚基I（COI）基因序列进行基于计算机的处理，从而能够从DNA数据中生成称为分子操作分类单元（MOTUs）的物种聚类假说。...

一种用于生成河岸带与水生甲虫分子操作分类单元（MOTUs）的整合分析流程

梯子行走任务是一种新的测试方法，用于评估精细行走能力，并测量前肢和后肢的放置、迈步以及肢体间的协调性。...

梯子横杆行走任务：评分系统及其实际应用

可诱导的遗传谱系示踪（GIFM）可实现对细胞的精细时空控制下的标记与追踪 体内...

一种用于遗传诱导命运图谱绘制的实用方法：标记与追踪细胞的可视化指南 体内

我们正在描述一种从用于血浆/血清学检测的全血样本中分离基因组DNA的新方法。在采集血浆后，通常会丢弃剩余的浓缩血液。我们的新方法相较于现有技术具有显著改进，可在单次采血后同时获得DNA和血浆，无需额外采集血液样本。...

一种用于从血液采集套件中提取基因组DNA并同时保存血浆蛋白的实用且新颖的方法

单个胚胎细胞的命运可能受到继承分子和/或邻近细胞信号的影响。利用卵裂期非洲爪蟾胚胎的命运图谱，可以鉴定出单个卵裂球并进行分离培养，以评估继承分子与细胞间相互作用各自的贡献。...

用于检测胚胎发育过程中细胞命运决定的卵裂球外植实验

用于蛋白质印迹中实用、便捷且可靠的总蛋白上样对照的 V3 免染工作流程

了解如何使用BLAST程序识别相似序列并进行多序列比对，从而从DNA或蛋白质序列重建可靠的系统发育树。本方法为非专业人士提供完整的生物信息学分析流程。

适用于非专业人士的系统发育学实用指南

探索更多视频