第三章序列分析

第三章序列分析1第1页，课件共92页，创作于2023年2月什么是序列分析?拿到一个基因/蛋白质序列,我能做什么?第2页，课件共92页，创作于2023年2月在数据库中进行序列相似性搜索基因结构分析/启动子序列分析Motif的寻找与序列的模式识别(含亚细胞定位,跨膜区的预测等序列分析的内容-----为了功能的分析第3页，课件共92页，创作于2023年2月序列比对第1节第4页，课件共92页，创作于2023年2月序列的相似性相似性(similarity)

是指一种很直接的数量关系。比如说，A序列和B序列的相似性是80％，或者4/5。这是个量化的关系。第5页，课件共92页，创作于2023年2月序列的相似性相似性(similarity)通常在某些位点上有一些氨基酸被另外一些化学物理特性相近的氨基酸所代替，这种突变可称为保守突变。将保守突变的因素考虑在内，就可以对两序列的相似程度打分，所得分值即代表其相似的程度。第6页，课件共92页，创作于2023年2月一致性(identity):两个蛋白质有一定数量的氨基酸在排比的位点上是相同的，即如果38个氨基酸的蛋白质中15个位点相同，我们说它们一致性为39.4%.MSDTPSTGFSIIHPTSSEGQVPPPRHLSLTHPVVAKRISFYKSG-------------PRNGTIKIYENPARTFTRPYSAKNITIYKEND第7页，课件共92页，创作于2023年2月所以,相似性的数值一定比一致性的要()大or相等or小第8页，课件共92页，创作于2023年2月同源性(homology)：

指从一些数据中推断出的两个基因或蛋白质序列具有共同祖先的结论，属于质的判断。就是说A和B的关系上，只有是同源序列，或者非同源序列两种关系。而说A和B的同源性为80％都是不科学的。生物序列的同源性第9页，课件共92页，创作于2023年2月相似性和同源性关系一般来说，序列间的相似性越高的话，它们是同源序列的可能性就更高。注意不要等价混用这两个名词。A序列和B序列的同源性为80％，记住这种说法是错误的！！第10页，课件共92页，创作于2023年2月PorcineandhumanBmal1gene第11页，课件共92页，创作于2023年2月PorcineandhumanBmal1gene第12页，课件共92页，创作于2023年2月复习-第1节序列比对判断题：A序列和B序列的相似性是80％.A序列和B序列的一致性为39.4%.A序列和B序列的同源性是80％.相似性的数值一定比一致性的要大.两序列间的相似性越高，它们的同源性就越高.第13页，课件共92页，创作于2023年2月第2节Blast的应用第14页，课件共92页，创作于2023年2月Blast简介

BLAST是基于序列相似性的数据库搜索程序。BLAST是“局部相似性基本查询工具”(BasicLocalAlignmentSearchTool)的缩写。第15页，课件共92页，创作于2023年2月Blast程序评价序列相似性的两个数据Score：使用打分矩阵对匹配的片段进行打分求和的结果，一般来说，匹配片段越长、相似性越高则Score值越大。Evalue:氨基酸残基（或碱基）随机排列得到上述Score值的概率的大小。E值越小表示越匹配。E=0表示完全配对，不存在随机配对。第16页，课件共92页，创作于2023年2月主要的blast程序程序名查询序列数据库搜索方法Blastn核酸核酸核酸序列搜索逐一核酸数据库中的序列Blastp蛋白质蛋白质蛋白质序列搜索逐一蛋白质数据库中的序列Blastx核酸蛋白质核酸序列翻译成蛋白质序列后和蛋白质数据库中的序列逐一搜索。Tblastn蛋白质核酸蛋白质序列和核酸数据库中的核酸序列翻译后的蛋白质序列逐一比对。TBlastx核酸核酸核酸序列翻译成蛋白质序列，再和核酸数据库中的核酸序列翻译成的蛋白质序列逐一进行比对。第17页，课件共92页，创作于2023年2月/Blast.cgi第18页，课件共92页，创作于2023年2月BLAST检索中采用的数据库类别：蛋白数据库：nr:无冗余数据库，汇集了GenBank中所有的蛋白序列+PDB+Swissprot+PIR等汇集；第19页，课件共92页，创作于2023年2月BLAST检索中采用的数据库类别：核酸数据库：nr/nt:无冗余数据库，汇集了GenBank+DDBJ+EMBL中所有的核甘酸序列（不包括EST)第20页，课件共92页，创作于2023年2月BLAST的应用第21页，课件共92页，创作于2023年2月具体步骤1.登陆blast主页

/BLAST/2.根据数据类型，选择合适的程序3.填写表单信息4.提交任务5.查看和分析结果第22页，课件共92页，创作于2023年2月第23页，课件共92页，创作于2023年2月第24页，课件共92页，创作于2023年2月第25页，课件共92页，创作于2023年2月第26页，课件共92页，创作于2023年2月第27页，课件共92页，创作于2023年2月提高期望阈值（Expectthreshold）;降低延伸种子序列的长度(wordsize)主要用途：搜索短的基序，如验证PCR引物的特异性第28页，课件共92页，创作于2023年2月Blastn---1第29页，课件共92页，创作于2023年2月Blastn---1第30页，课件共92页，创作于2023年2月Blastn---1Blastn1的作用：①对于已知的基因，可以分析其相似基因；②对于未知的基因片段，可以分析其属于什么基因。第31页，课件共92页，创作于2023年2月Blastn---1第32页，课件共92页，创作于2023年2月Blastn---1第33页，课件共92页，创作于2023年2月Blastn---2WhatisEST?第34页，课件共92页，创作于2023年2月Blastn---2第35页，课件共92页，创作于2023年2月Blastn---2有什么作用?对于新基因，可以了解基因的结构特征,5-UTR和3-UTR的大小.

第36页，课件共92页，创作于2023年2月Blastp第37页，课件共92页，创作于2023年2月Blastp第38页，课件共92页，创作于2023年2月Blastp第39页，课件共92页，创作于2023年2月Blastp第40页，课件共92页，创作于2023年2月Blastp第41页，课件共92页，创作于2023年2月Blastp第42页，课件共92页，创作于2023年2月BlastpBlastp的作用：①对于已知的蛋白，可以分析其相似蛋白；②对于未知的蛋白片段，可以分析其属于什么蛋白。第43页，课件共92页，创作于2023年2月具体步骤1.登陆blast主页

/BLAST/2.根据数据类型，选择合适的程序3.填写表单信息4.提交任务5.查看和分析结果第44页，课件共92页，创作于2023年2月第45页，课件共92页，创作于2023年2月rice第46页，课件共92页，创作于2023年2月第47页，课件共92页，创作于2023年2月第48页，课件共92页，创作于2023年2月第49页，课件共92页，创作于2023年2月Blastx第50页，课件共92页，创作于2023年2月Blastx第51页，课件共92页，创作于2023年2月Blastx第52页，课件共92页，创作于2023年2月BlastxBlastx:

分析你的基因编码什么产物?分析你的基因是否是新基因？第53页，课件共92页，创作于2023年2月Blastx第54页，课件共92页，创作于2023年2月tBlastn第55页，课件共92页，创作于2023年2月tBlastn第56页，课件共92页，创作于2023年2月tBlastntBLASTn的作用:已知一种蛋白序列，在另一物种中进行其同源蛋白基因的电子克隆(insilicocloning)；寻找一个新的蛋白质序列(如双向电泳得到的)是否已有核酸序列，是否可以克隆？第57页，课件共92页，创作于2023年2月Blastx当Blastx没有结果时，可以考虑使用。第58页，课件共92页，创作于2023年2月第一节序列比对第二节Blast应用第三节序列功能分析第59页，课件共92页，创作于2023年2月序列分析的目的是什么?

--为了功能的分析--拿到一个基因/蛋白质序列,我能做什么?第60页，课件共92页，创作于2023年2月序列组成/分子量/等电点---初级分析酶切位点分析(载体构建)基因结构分析/启动子序列分析序列功能分析的内容

第61页，课件共92页，创作于2023年2月Part1.初级序列分析序列的组成/分子量/等电点分析第62页，课件共92页，创作于2023年2月/第63页，课件共92页，创作于2023年2月第64页，课件共92页，创作于2023年2月Part2.

酶切位点分析

只要进行基因工程利用必须用到各种限制性内切酶如GGATCCBamHI第65页，课件共92页，创作于2023年2月进行酶切位点分析的时候，对于构建载体，我们需要知道的信息是你的序列中有/没有某个酶的位点？为什么？如果答案是“有”，是什么情况？“没有”又是什么情况？PlasmidvectorSacIXbaIHindIIIHindIIIXbaICloningsiteSacI第66页，课件共92页，创作于2023年2月Part3.基因结构分析/启动子序列分析第67页，课件共92页，创作于2023年2月1）基因结构分析:了解基因的内含子/外显子排列方式;例子1:水稻6PGDH基因进化分析的结果表明其可能来源于内共生（基因结构分析表明其没有内含子）.例子2：NHX基因第68页，课件共92页，创作于2023年2月1）基因结构分析:GenomicDNAcDNA第69页，课件共92页，创作于2023年2月用softberry预测基因结构/berry.phtml?topic=fgenesh&group=programs&subgroup=gfind第70页，课件共92页，创作于2023年2月什么是HMM?HiddenMarkovModels(HMMs,隐马尔可夫模型)最早是在上个世纪60年代末70年代初提出来的一种概率论模型。进入80年代以后，逐渐被利用在各个领域。主要的应用领域：语音识别系统。生物学中的DNA/protein序列的分析。机器人的控制。文本文件的信息提取。第71页，课件共92页，创作于2023年2月什么是HMM?对于给定的核苷酸序列，我们在位置p处出现的概率依赖于它后核苷酸序列。即任何一个核苷酸的出现并不是完全独立的。对于基因识别而言，HMM可以分为2种：按内容搜索：即根据核苷酸和密码子在编码区内的分布规律来确定蛋白质的编码区;按信号搜索：即通过分析编码区周围的信号(启动子\终止子\各种cis-element等来确定蛋白质的编码区。大家只要知道HMM是目前生物信息学中应用非常广泛的概率论模型，广泛应用于基因识别，其基本原理基于任何一个序列的出现都不是独立的。第72页，课件共92页，创作于2023年2月一个例子:用softberry预测基因结构第73页，课件共92页，创作于2023年2月2)启动子序列分析:

什么是启动子?启动子序列，一般在TSS之前2000bp，了解哪个位点是TSS，哪个是起始ATG?TSSATGTATApromoter转录起始位点(transcriptionstartsite，TSS)第74页，课件共92页，创作于2023年2月2)启动子序列分析:

所以，我们必须得到TSS的位置.如何通过生物信息学方法确定TSS?

首先截取包括ATG之前3000bp和基因的序列采用以下两种方法1)软件预测,如Softberry;2)搜索EST数据库;第75页，课件共92页，创作于2023年2月分析的目的:2)首先找到ATG前面约3000:如何通过生物信息学方法确定TSS?以A