Clustalx 多重序列比对图解教程(By Raindy)
本帖首发于Raindy'blog,转载请保留作者信息,谢谢!欢迎有写生物学软件专长的战友,加入生信教程写作:13559330,接头暗号:你所擅长的生物学软件名称

软件简介:
  CLUSTALX-是CLUSTAL多重序列比对程序的Windows版本。Clustal X为进行多重序列和轮廓比对和分析结果提供一个整体的环境。
  序列将显示屏幕的窗口中。采用多彩的模式可以在比对中加亮保守区的特征。窗口上面的下拉菜单可让你选择传统多重比对和轮廓比对需要的所有选项。

主要功能
  你可以剪切、粘贴序列以更改比对的顺序;
  你可以选择序列子集进行比对;
  你可以选择比对的子排列(Sub-range)进行重新比对并可插入到原始比对中;
  可执行比对质量分析,低分值片段或异常残基将以高亮显示。

当前版本:1.83
PS:如果你是新手或喜欢中文界面,推荐使用本人的Clustalx 1.81
  链接地址:/index.php?Go=Show:ist&ID=7435 (请完整复制)

应用:Clustalx比对结果是构建系统发育树的前提

实例:植物呼肠孤病毒属外层衣壳蛋白P8(AA序列)为例

流程:载入序列>编辑序列>设置参数>完全比对>比对结果

  1.载入序列:运行ClustalX,主界面窗 口如下所图(图1),依次在程序上方的菜单栏选择File”-Load Sequence”载入待比对的序列,如图2所示,如果当前已载入序列,此时会提示是否替换现有序列(Replace existing sequences),根据具体情形选择操作。
                      图1

             图2
  2.编辑序列 对标尺(Ruler)上方的序列进行编辑操作,主要有Cut sequences(剪切序列)Paste sequences(粘贴)Select All sequences(选定所有序列)Clear sequence Selection(清除序列选定)Search for string(搜索字串)Remove All gaps(移除序列空位)Remove Gap-Only Columns(仅移除选定序列的空位)

             图3
  3.参数设置 可以根据分析要求设置相对的比对参数。通常情况下,我们可以使用默认参数。比对参数主要有六个,分别是Reset New Gaps before Alignment(比对前重置新的空位参数vim命令列间编辑)Reset All Gaps before Alignment(比对前重置所有空位参数),Pairwise Alignment Parameters(两两序列比对参数),Multiple Alignment Parameters(多重序列比对参数),Protein Gap Parameters(蛋白空位参数),Secondary structure Parameters(二级结构参数),如图4所示:

             图4
  修改参数只需点击相应标签,示例比对的是多序列比对,故可选择Multiple Alignment Parameters”弹出参数设置窗口,如图5所示:

             图5
  4.完全比对 返回菜单栏选择Complete Alignment”标签,此时会弹出输出文件路径的设置窗口,设置Guide Tree File(向导树或指导树文件)、Alignment File(比对文件)的保存位置(存放路径),点击Align”按钮程序自动开始序列的完全比对,比对所需时间因序列文件大小和长度、计算机性能而异, 如图6-8所示:

             图6

             图7

             图8
  当主界面的左下状态栏会提示CLUSTAL-Alignment File created []”时说明比全完毕,这时文件保存位置的目录下会生成生成两个文件,分别是*.aln*.dndaln是序列比对的文件,可以进一步用于构树系统发 育树,dnd是向导树文件(指导树),这两个文件可以用Windows系统中的记事本或第三方程序UltraEdit”等打开,如:

             图9 ALN文件

             图10 dnd文件

  5.后续分析:
  1)Clustalx比对生成的结果可读性不是太好,一般需要专业的序列着软件处理,如BoxshadeESPript,这两个工具都是在线进行,其中Boxshade图解教程详见本Blog日志。
Boxshade在线网址:/software/BOX_form.html
ESPript在线网址:espript.ibcp.fr/ESPript/i

  2)转换ALN文件,进一步构建系统发育树,转换格式依不同软件而有所不同,如在PHYLIP分析前需要将ALN格式转换为PHY格式方可...


注意事项
1)dnd是向导树文件,可以用TreeViews软件查看树图。注意:向导树不是系统发育树,两者区别敬请关注近期-系统发育分析专题。
2)多重比对文件推荐要求为规范的FASTA格式,文件扩展名不限,格式大致如下:

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。