DNA序列比对最大似然度进化模型

摘要

序列比对是生物信息学中重要的研究课题，是发现序列的功能，结构和进化信息的重要手段。现有的很多比对算法都是基于目标函数，目标函数利用替换矩阵和空位罚分对比对过程和结果进行记分，用记分值来判断序列比对结果的好坏。基于目标函数的比对算法的缺点是记分系统的些许变动就可能导致局部或全局比对剧烈变化。由此，本文提出了DNA序列最大似然度比对算法来对DNA序列进行序列比对。 本文首先介绍序列比对的基本概念，详细叙述了替换矩阵、空位罚分和目标函数以及它们对序列比对的影响，然后深入研究了双序列各种常用的比对算法：点阵分析、动态规划算法和词或k-串方法，并给出了它们的算法思想或伪代码。接着，根据基于目标函数算法的缺点，提出最大似然度比对算法，该算法分为两个部分：参数估值算法和比对算法。最大似然度比对算法首先使用进化参数估值算法对一对DNA序列相关的进化参数进行估值，然后比对算法利用估算出来的参数值对DNA序列进行比对。它是一个独立的方法，完全避开了基于替换矩阵和空位罚分的序列比对算法由于DNA序列相似度的不同而要选择与之相适应的替换矩阵的问题。最后，通过对进化模型最大似然度DNA比对算法的序列比对结果和FASTA比对程序结果进行比较，验证了进化模型最大似然度DNA比对算法的正确性和精确性。

目录

文摘

英文文摘

声明

第一章绪论

1.1引言

1.2序列比对的现状

1.3本文所做的工作

第二章双序列比对基础

2.1序列比对的基本概念

2.1.1序列比对的定义

2.1.2替换矩阵

2.1.3空位罚分

2.1.4目标函数

2.2双序列比对算法

2.2.1点阵分析

2.2.2动态规划算法

2.2.3词或k-串方法

2.3本章小结

第三章DNA序列比对的最大似然度方法

3.1 DNA序列进化模型

3.1.1替换过程

3.1.2插入/删除过程

3.2最大似然度序列比对算法

3.2.1双序列比对的似然度表达式

3.2.2进化参数估计算法

3.2.3比对算法

3.3实验数据和结果分析

3.3.1实验说明

3.3.2实验数据

3.3.3结果分析

3.4本章小结

第四章 结论与展望

附录

致谢

参考文献

研究成果