序列比对原理.ppt

三、多序列比对工具 （一）ClustalX/W ClustalX和ClustalW是两个使用最广泛的多序列比对工具，均采用渐进式多序列比对算。 ClustalX的比对步骤: 1.加载要比对的序列文件 序列格式可以支持NBRF/PIR、EMBL/Swiss-Ppot、FASTA、GDE、Clustal、MSF、RSF等。 2.多序列比对 3.比对结果输出 （二）T-Coffee工具 （三）MultAlin工具 （四）MAFFT工具 序列比对原理 第一节 序列比对相关概念 一、序列比对目的及定义 （一）序列比对目的 通过比较两条或多条序列之间是否具有足够的相似性，从而判定它们之间是否具有同源性。 （二）序列比对定义 序列比对（sequence alignment）是运用某种特定的数学模型或算法,找出两个或多个序列之间的最大匹配碱基或残基数,比对的结果反映了算法在多大程度上提供序列之间的相似性关系及它们的生物学特征。 二、序列比对类型 （一）序列比对分类 双序列比对 多序列比对 global alignment local alignment (二)编辑距离 通过编辑操作计算的两条序列的距离称为编辑距离。 （三）双序列比对 （四）全局序列比对 （五）局部序列比对 三、序列比对的相关概念 （一）同源性、同一性、相似性 相似性（similarity）是指两序列间直接的数量关系，如部分相同、相似的百分比或其他一些合适的度量。 同一性（identity）是指两序列在同一位点核苷酸或氨基酸残基完全相同的序列比例。 同源性（homology）是指从某个共同祖先经趋异进化而形成的不同序列。 （二）直系同源、旁系同源 直系同源基因（orthologous gene）是指在不同物种中有相同功能的同源基因，它是在物种形成过程中形成的。 旁系同源基因（paralogous gene）是指一个物种内的同源基因。 直系同源基因和旁系同源基因统称为同源基因（homolog）。 第二节 序列比对打分方法 一、序列比对打分 序列分析的目的是揭示核苷酸或氨基酸序列编码的高级结构或功能信息目的。 二、打分矩阵 稀疏矩阵、相似性打分矩阵 （一）DNA打分矩阵 即两个序列中相应的核苷酸相同，计1分；否则计0分。 如果考虑颠换和置换，可采用以下打分矩阵 （二） 氨基酸序列打分矩阵 1.PAM矩阵（Dayhoff突变数据矩阵） PAM250矩阵 2.BLOSUM矩阵 BLOSUM62矩阵 第三节 序列比对算法 一、dotplot算法 1.构建点阵矩阵 2.获得相似性片段 二、dynamic programming algorithm 1、计算得分矩阵 2、寻找最优的比对序列 例 s=acgctq t=catgt 算法特点： 三、BLAST算法 1、编译一个由查询序列生成的长度固定的字段编译列表； 2、在数据库中扫描获得与编译列表中的字段匹配的序列记录； 3、以编译列表中的字段对为中心向两端延伸以寻找超过阈值分数S的高分值片段对HSP。 E值计算公式： 算法特点： 第四节 序列比对工具 一、FASTA工具 二、BLAST工具 （一）基本BLAST工具 nucleotide blast Search a?nucleotide?database using a?nucleotide?queryAlgorithms:?blastn, megablast, discontiguous megablast protein blast Search?protein?database using a?protein?queryAlgorithms:?blastp, psi-blast, phi-blast, delta-blast blastx Search?protein?database using a?translated nucleotide?query tblastn Search?translated nucleotide?database using a?protein?query tblastx Search?translated nucleotide?database using a?translated nucleotide?query 比对结果 （二）高级BLAST工具 1、PSI-BLAST（position specific interated BLAST） 2、PHI-BLAST（pattern hit initiated BLAST） 3、MEGABLAST 第五节 多序列比对 一、多序列比对概述 （一）多序列比对目的 为了发现构成同一基因家族的成组序列