1、第五节第五节 基因识别基因识别主讲人：孙主讲人：孙 啸啸 制作人：刘志华制作人：刘志华东南大学 吴健雄实验室基因识别基因识别v基因识别是生物信息学领域里的一个重基因识别是生物信息学领域里的一个重要研究内容要研究内容 v基因识别问题，在近几年受到广泛的重基因识别问题，在近几年受到广泛的重视视 当人类基因组研究进入一个系统测序阶段当人类基因组研究进入一个系统测序阶段时，急需可靠自动的基因组序列翻译解释时，急需可靠自动的基因组序列翻译解释技术，以处理大量已测定的但未知功能或技术，以处理大量已测定的但未知功能或未经注释的未经注释的DNA序列序列 v原核基因识别原核基因识别重点在于识别编码区域重点在于识
2、别编码区域v非翻译区域（非翻译区域（untranslated regions,UTR）编码区域两端的编码区域两端的DNA，有一部分被转录，有一部分被转录，但是不被翻译，这一部分称为非翻译区域但是不被翻译，这一部分称为非翻译区域 v5UTR-基因上游区域的非翻译区域基因上游区域的非翻译区域 v3UTR-基因下游区域的非翻译区域基因下游区域的非翻译区域v对于任何给定的核酸序列（单链DNA或mRNA），根据密码子的起始位置，可以按照三种方式进行解释。v例如，序列ATTCGATCGCAAv这三种阅读顺序称为阅读框（reading frames）CAA A ATTCGATCGATTCGATCGCAAAT
3、TCGATCGCA（1）（3）（2）一个开放阅读框（一个开放阅读框（ORF,open reading frame）是一个没有终止编码的密码子序）是一个没有终止编码的密码子序列。列。原核基因识别任务的重点是识别开放阅读原核基因识别任务的重点是识别开放阅读框，或者说识别长的编码区域。框，或者说识别长的编码区域。基于基因密码子特性的识别方法基于基因密码子特性的识别方法v辨别编码区域与非编码区域的一种方法辨别编码区域与非编码区域的一种方法是检查终止密码子的出现频率是检查终止密码子的出现频率 终止密码子出现的期望次数为：终止密码子出现的期望次数为：每每21个（个（64/3）密码子出现一次终止密）密码子出
4、现一次终止密码子码子 基本思想：基本思想：如果能够找到一个比较长的序列，其相应如果能够找到一个比较长的序列，其相应的密码子序列不含终止密码子，则这段序的密码子序列不含终止密码子，则这段序列可能就是编码区域。列可能就是编码区域。v基本算法：基本算法：扫描给定的扫描给定的DNA序列，在三个不同的阅读序列，在三个不同的阅读框中寻找较长的框中寻找较长的ORF。遇到终止密码子以。遇到终止密码子以后，回头寻找起始密码子。后，回头寻找起始密码子。这种算法过于简单，不适合于处理短的这种算法过于简单，不适合于处理短的ORF或者交叠的或者交叠的ORF。v识别编码区域的另一种方法是分析各种识别编码区域的另一种方法是
5、分析各种密码子出现的频率密码子出现的频率 将一个随机均匀分布的将一个随机均匀分布的DNA序列翻译成氨基酸序列翻译成氨基酸序列，则在氨基酸序列中上述序列，则在氨基酸序列中上述3种氨基酸出现的种氨基酸出现的比例应该为比例应该为6:4:1例如，亮氨酸、丙氨酸、色氨酸分别有例如，亮氨酸、丙氨酸、色氨酸分别有6个、个、4个和个和1个密码子个密码子但是在真实的氨基酸序列中，上述比例并不但是在真实的氨基酸序列中，上述比例并不正确正确这说明这说明DNA的编码区域并非随机的编码区域并非随机v假设在一条假设在一条DNA序列中已经找到所有的序列中已经找到所有的ORF，那么可以利用密码子频率进一步，那么可以利用密码子
6、频率进一步区分编码区分编码ORF和非编码和非编码ORFv马尔柯夫链模型马尔柯夫链模型v利用这种方法，可以计算一个利用这种方法，可以计算一个ORF成为成为编码区域的可能性。编码区域的可能性。v一个简单的统计模型一个简单的统计模型假设相继的密码子是独立的，不存在前后依假设相继的密码子是独立的，不存在前后依赖关系。赖关系。令令fabc代表密码子代表密码子abc在编码区域出现的频率在编码区域出现的频率给定序列给定序列a1,b1,c1,a2,b2,c2,an+1,bn+1从密码子从密码子a1b1c1开始的阅读框，其开始的阅读框，其n个密码子个密码子的出现概率为的出现概率为nnncbacbacbafffp
7、.2211111v第二种和第三种阅读框第二种和第三种阅读框n个密码子出现的概个密码子出现的概率分别为率分别为1322211.2nnnacbacbacbfffp11332221.3nnnbacbacbacfffpv第第i个阅读框成为编码阅读框的概率个阅读框成为编码阅读框的概率计算：计算：v算法：算法：在序列上移动长度为在序列上移动长度为n的窗口，计算的窗口，计算Pi根据根据Pi的值识别编码的阅读框的值识别编码的阅读框321ppppPii基于编码区域碱基组成特征的识别方法基于编码区域碱基组成特征的识别方法v编码序列与非编码序列在碱基组成上编码序列与非编码序列在碱基组成上有区别有区别单个碱基的组成比
8、例单个碱基的组成比例多个碱基的组成多个碱基的组成v通过统计分析识别编码序列通过统计分析识别编码序列分析实例分析实例2、真核基因识别问题、真核基因识别问题 真核基因远比原核基因复杂：真核基因远比原核基因复杂：v一方面，真核基因的编码区域是非连续一方面，真核基因的编码区域是非连续的，编码区域被分割为若干个小片段。的，编码区域被分割为若干个小片段。v另一方面，真核基因具有更加丰富的基另一方面，真核基因具有更加丰富的基因调控信息，这些信息主要分布在基因因调控信息，这些信息主要分布在基因上游区域。上游区域。基因识别基本思路基因识别基本思路 找出基因两端的功能区域找出基因两端的功能区域:转录启动区转录启动
9、区 终止区终止区 在启动区下游位置寻找翻译起始密码子在启动区下游位置寻找翻译起始密码子 识别转录剪切位点识别转录剪切位点剪切给体位点剪切给体位点剪切接受体位点剪切接受体位点v各种不同的方法有不同的适应面，而不各种不同的方法有不同的适应面，而不同的方法有时可以结合起来以提高基因同的方法有时可以结合起来以提高基因识别的准确率。识别的准确率。v关键问题是如何提高一个识别算法的敏关键问题是如何提高一个识别算法的敏感性（感性（sensitivity，Sn）和特异性）和特异性（specificity，Sp）。）。3、基因识别的主要方法、基因识别的主要方法两大类识别方法：两大类识别方法：v从头算方法（或基于
10、统计的方法）从头算方法（或基于统计的方法）根据蛋白质编码基因的一般性质和特征进行识别，根据蛋白质编码基因的一般性质和特征进行识别，通过统计值区分外显子、内含子及基因间区域通过统计值区分外显子、内含子及基因间区域 v基于同源序列比较的方法基于同源序列比较的方法利用数据库中现有与基因有关的信息（如利用数据库中现有与基因有关的信息（如EST序序列、蛋白质序列），通过同源比较，帮助发现新列、蛋白质序列），通过同源比较，帮助发现新基因。基因。v最理想的方法是综合两大类方法的优点，最理想的方法是综合两大类方法的优点，开发混合算法。开发混合算法。v基因识别方法有基因识别方法有：（1）基于规则的系统）基于规则
11、的系统（2）语义学方法语义学方法（3）线性辨别分析（线性辨别分析（LDA）（4）决策树决策树 （5）动态规划动态规划 （6）隐马尔柯夫模型隐马尔柯夫模型 （7）剪切对比排列剪切对比排列 （spliced alignment）4、编码区域识别、编码区域识别两类方法两类方法：v基于特征信号的识别基于特征信号的识别 内部外显子内部外显子剪切位点剪切位点5端的外显子一定在核心启动子的下游端的外显子一定在核心启动子的下游3端的外显子的下游包含多聚信号和终端的外显子的下游包含多聚信号和终止编码止编码 v基于统计度量的方法基于统计度量的方法 根据密码子使用倾向根据密码子使用倾向双联密码统计度量等双联密码统计
12、度量等v在一个基因中，第在一个基因中，第i个（个（i=1，64）密码子相对使用）密码子相对使用倾向倾向RSCUi的定义如下：的定义如下：vObsi是该基因中第是该基因中第i个密码子实际出现的次数个密码子实际出现的次数Expi是对应密码子期望的出现次数是对应密码子期望的出现次数 aai是统计的第是统计的第i个密码子出现的次数个密码子出现的次数 syni是所有与第是所有与第i个密码子同义密码子出现的次数个密码子同义密码子出现的次数vRSCU大于大于1表示相应密码子出现的次数比期望次数表示相应密码子出现的次数比期望次数高，而小于高，而小于1则表示出现次数相对较少。则表示出现次数相对较少。iiiExp
13、ObsRSCU iiisynaaExp（5-66）（5-65）密码子使用倾向密码子使用倾向v设一段DNA序列为S，从S的第i位到第j位的双联密码统计度量IF6（i，j）定义为：fk是从第k位开始的双联密码的频率Fk是该双联密码随机出现的频率4,.,8,5,25,.,7,4,16,.,6,3,06)/ln()/ln()/ln(max),(jkkikijkkikijkkikiFfFfFfjiIF（5-67）双联密码统计度量双联密码统计度量通过相似搜索发现编码区域或者外显子通过相似搜索发现编码区域或者外显子 EST（Expressed Sequence Tags）cDNA 蛋白质序列蛋白质序列v目前
14、大多数预测程序都将数据库相似性目前大多数预测程序都将数据库相似性搜索的信息结合进基因预测过程搜索的信息结合进基因预测过程v同时考虑序列特征信号和统计度量同时考虑序列特征信号和统计度量GRAIL用人工神经网络识别编码区域用人工神经网络识别编码区域输入是一系列反映功能位点信号特征和序列编码统计特征的参数输入是一系列反映功能位点信号特征和序列编码统计特征的参数输出就是对一段输出就是对一段DNA序列是否是编码区域的判别结果序列是否是编码区域的判别结果神经网络具有非线性映射能力，能够发现输入和输出之间的高阶相关神经网络具有非线性映射能力，能够发现输入和输出之间的高阶相关性性5、构建基因模型v基因识别最终
15、任务是建立完整的基因结构模型v一个理想的基因识别程序应该能够发现完整的基因结构（,e1,i1,in-1,en,）ATG-外显子1内含子外显子外显子n-UAGv基因剪切位点基因剪切位点 剪切给体（剪切给体（donor）位点）位点-“gt”接受体（接受体（acceptor）位点）位点-“ag”基因的可变剪切基因的可变剪切gene A基因可变剪切示意基因可变剪切示意构建基因模型方法构建基因模型方法 v剪切位点形成外显子和内含子的边界剪切位点形成外显子和内含子的边界 搜集候选外显子搜集候选外显子 候选基因候选基因v候选基因是一条候选基因是一条非相交非相交的外显子和内含的外显子和内含子的链，表示为子的链
16、，表示为 （i0,e1,i1,en,in）其中其中ij代表内含子（代表内含子（0 j n）el代表外显子（代表外显子（1 l n）i0和和in并非真实的内含子，它们分别代表并非真实的内含子，它们分别代表基因两侧的非编码序列基因两侧的非编码序列v候选基因位于给定的候选基因位于给定的DNA序列，并满足下列一序列，并满足下列一致性条件：致性条件：（1）所有外显子加起来的长度是）所有外显子加起来的长度是3的整数倍；的整数倍；（2）在各个外显子内部（除最后一个外显子的）在各个外显子内部（除最后一个外显子的最后一个密码子），没有终止编码；最后一个密码子），没有终止编码；（3）第一个内含子）第一个内含子-外
17、显子边界（外显子边界（i0,e1）是翻）是翻译起始编码，而最后一个外显子译起始编码，而最后一个外显子-内含子边界内含子边界（en,in）是终止编码。）是终止编码。位点图位点图（分层标注剪切位点）（分层标注剪切位点）另设两个特殊的顶点，即起点（另设两个特殊的顶点，即起点（source）和终点（）和终点（sink）。从起点到终点的任何一条路径代表一个可能的基因结构。从起点到终点的任何一条路径代表一个可能的基因结构。例如例如：位点图上的路径位点图上的路径v候选基因所对应的道路图中的路径候选基因所对应的道路图中的路径v求最优路径求最优路径每一条弧附加一个权值每一条弧附加一个权值外显子、内含子度量每个节
18、点附加权值每个节点附加权值剪切位点度量 综合评价综合评价 6、用于基因识别的HMM模型v隐马尔柯夫模型隐马尔柯夫模型HMM是一条状态不可见是一条状态不可见的马尔柯夫链，其当前状态的输出是可的马尔柯夫链，其当前状态的输出是可见的。见的。每个状态按照一定的概率分布随机地从字母每个状态按照一定的概率分布随机地从字母表中取出字符并释放。表中取出字符并释放。v扩展的隐藏马尔柯夫模型（扩展的隐藏马尔柯夫模型（GHMMs）对对HMM进一步抽象，产生更一般的马尔柯进一步抽象，产生更一般的马尔柯夫模型，以分析复杂的脊椎动物基因夫模型，以分析复杂的脊椎动物基因。(1)信号传感器模型信号传感器模型v将剪切位点、起始
19、编码区域或者终止编码区域看成是DNA序列上的功能位点或者信号位点，用HMM来进行分析 内含子区域内含子区域 外显子区域外显子区域 保守位点保守位点 根据对比排列，根据对比排列，形成具有形成具有19状态状态的的HMM模型。模型。v对前一节所介绍的对前一节所介绍的HMM模型进行修改，模型进行修改，可以处理双联核苷酸的问题，即将可以处理双联核苷酸的问题，即将4种种概率分布扩展为概率分布扩展为16种。种。v假设一段序列为假设一段序列为ACTGTC，则，则 P(ACTGTC)=p1(A)p2(C A)p3(T C)p4(G T)p5(T G)p6(C T)其中其中p1是状态是状态1对于对于4种核苷酸的概
20、率，种核苷酸的概率，p2(x y)状态状态2的条件概率。的条件概率。(2)编码区模型v由于密码子的长度为3，因此密码子模型的最后一个状态应该至少为2阶。v对于2阶的状态，具有64种概率分布，可根据已知编码区域进行统计计算而得到64种分布。例如：p(ACA)=c(CAA)/c(CAA)+c(CAC)+c(CAG)+c(CAT)p(CCA)=c(CAC)/c(CAA)+c(CAC)+c(CAG)+c(CAT)p(GCA)=c(CAG)/c(CAA)+c(CAC)+c(CAG)+c(CAT)p(TCA)=c(CAT)/c(CAA)+c(CAC)+c(CAG)+c(CAT)其中，c(xyz)是密码子x
21、yz的计数。v这样的模型可以检测无结束编码的区域，因为对应于三个结束编码TAA、TAG和TGA的p(ATA)、p(GTA)和p(ATG)自动为0。(3)组合模型v将上述模型扩展，使之可以识别具有多个外显子的基因。v改进后的模型见下图、基于剪切比对的基因识别方法v基本思想是：利用数据库中的同源信息进行基因识别，包括DNA、RNA和蛋白质数据库。v其方法是：首先通过分析所有可能的剪切接受体位点和剪切给体位点，构建一组候选的外显子。然后进一步分析候选外显子，探查所有可能的外显子组合，寻找一个与已知目标蛋白质或其他表达序列最匹配的组合 v一种半自动的综合方法识别基因过程：（1）选择所有长度大于50bp
22、并介于保守的剪切接受位点和给体位点之间的ORF，作为候选的外显子；预选预选（2）对于候选的外显子计算其6目编码度量值，并从大到小将它们排列起来；减小搜索范围减小搜索范围（3）对照蛋白质序列数据库进行搜索，寻找相似体。搜索，筛选搜索，筛选、基因识别程序介绍表表5.7 基因识别程序及访问地址基因识别程序及访问地址（HP主页；主页；ESE-mail服务器；服务器；WSweb服务器；服务器；CL客户客户/服务器协议；服务器协议；EX有可执行代码；有可执行代码；SC有源代码）有源代码）表表5.8 各程序的性能比较（敏感性各程序的性能比较（敏感性(1)被预测出的真实编码核酸的被预测出的真实编码核酸的%；敏
23、感性敏感性(2)被正确识别出的编码外显子的被正确识别出的编码外显子的%；特异性特异性(1)预测出的编码核酸为真实编码核酸的预测出的编码核酸为真实编码核酸的%；特异性特异性(2)预测出外显子为真实外显子的预测出外显子为真实外显子的%）v基因识别方法存在的问题和局限性：基因识别方法存在的问题和局限性：（1）关于基因的定义不明确）关于基因的定义不明确 统一定义统一定义 （2）目前的方法仅仅识别蛋白质编码基因）目前的方法仅仅识别蛋白质编码基因 转录信号转录信号 （3）现有的许多方法仅检测单个基因）现有的许多方法仅检测单个基因部分基因、多重基因部分基因、多重基因 （4）基于同源分析的方法是保守的）基于同源分析的方法是保守的不可能发现新的基因不可能发现新的基因 （5）忽视关于基因结构的生物学知识忽视关于基因结构的生物学知识 基因表达的真实分子机制基因表达的真实分子机制 谢谢！
