第三章-BLAST原理及方法-课件

第三章BLAST:BasiclocalalignmentsearchtoolBLAST!2020/10/281Outline

SummaryofkeypointsaboutpairwisealignmentIntroductiontoBLAST:practicalguidetodatabasesearchingTheBLASTalgorithmBLASTsearchstrategies2020/10/282精品资料3Pairwisealignment:keypoints

Pairwisealignmentsallowustodescribethepercentidentitytwosequencesshare,aswellasthepercentsimilarityThescoreofapairwisealignmentincludespositivevaluesforexactmatches,andotherscoresformismatchesandgapsPAMandBLOSUMmatricesprovideasetofrulesforassigningscores.PAM10andBLOSUM80areexamplesofmatricesappropriateforthecomparisonofcloselyrelatedsequences.PAM250andBLOSUM30areexamplesofmatricesusedtoscoredistantlyrelatedproteins.Globalandlocalalignmentscanbemade.2020/10/284BLASTBLAST(BasicLocalAlignmentSearchTool)allowsrapidsequencecomparisonofaquerysequenceagainstadatabase.TheBLASTalgorithmisfast,accurate,andweb-accessible.2020/10/285WhyuseBLAST?BLAST是NCBI中用来将一个蛋白质或DNA序列和各种数据库中的其他序列进行比对的主要工具。BLAST搜索是研究一个蛋白质和基因的最基本的方法之一。2020/10/286WhyuseBLAST?BLAST具有非常广泛的应用：

确定特定的蛋白质或核酸序列有哪些已知的直系同源或旁系同源序列。除了RBP外，还有哪些其他的脂质运载蛋白是我们所知道的？当一个新的细菌的基因组被测序后，几千种蛋白质被确定，其中有多少蛋白质是同源的？从这里面测出的基因中有多少是在GenBank中找不到显著性同源物的？

确定哪些蛋白质和基因在特定的物种中出现。植物中是否也存在像RBP这样的脂质运载蛋白？鱼类中是否有反转录酶基因(HIV-1pol

基因)？2020/10/287WhyuseBLAST?BLAST具有非常广泛的应用：

确定一个DNA或蛋白质序列身份。如可能通过一个负杂交实验或芯片实验发现一个特殊的DNA序列中你所使用的实验条件下是被显著调控的，那么就可以通过将这个DNA序列在一个蛋白质数据库中进行搜索，来寻找哪些蛋白质是与你的DNA序列所编码的蛋白质相关性最高的。

发现新基因。如一个对于全基因组DNA的BLAST搜索可能会发现一个DNA所编码的蛋白质是以前所没有报道过的

确定一个特定基因或蛋白质有哪些已经发现了的变种。例如，很多病毒都具有极强的突变能力；HIV-1pol有哪些已知的变异体？2020/10/288WhyuseBLAST?BLAST具有非常广泛的应用：

研究可能存在多种剪切方式的表达序列标签。有专门用于BLAST搜索的EST数据库。实际上有许多用来进行搜索专门的数据库，如专门的包含同一个特定的物种、一种组织、一个染色体、一种DNA或一个蛋白质功能类的序列数据库。

寻找对于一个蛋白质的功能和/或结构起关键作用的氨基酸残基。一次BLAST搜索的结果可以放在一起比对，这时候，就会发现其中像半胱氨酸残基这样可能具有重要生物学功能的保守残基。2020/10/289FourcomponentstoaBLASTsearch(1)Choosethesequence(query)(2)SelecttheBLASTprogram(3)Choosethedatabasetosearch(4)ChooseoptionalparametersThenclick“BLAST”2020/10/28102020/10/2811NP_0067352020/10/28122020/10/2813Step1:Chooseyoursequence三种主要的输入方式：

剪切然后粘贴DNA或蛋白质序列使用FASTA格式的序列简单地使用索引号码(如一个RefSeq或GenBank(GI)的序号)。SequencecanbeinputinFASTAformatorasaccessionnumber2020/10/2814ExampleoftheFASTAformatforaBLASTquery一个FASTA格式的序列以一个单行的说明开始，接下来是若干个行的序列数据。在一个BLAST搜索中输入accessionnumber通常要容易些。BLAST程序可以识别和忽略出现在你的输入序列字母中间的数字。2020/10/2815Step2:ChoosetheBLASTprogram2020/10/2816Step2:ChoosetheBLASTprogramblastn(nucleotideBLAST)blastp(proteinBLAST)tblastn(translatedBLAST)blastx(translatedBLAST)tblastx(translatedBLAST)2020/10/28172020/10/2818Step2:ChoosetheBLASTprogramblastn(nucleotideBLAST)：将一个核酸的查询序列与一个核酸序列数据库相比较。blastp(proteinBLAST)：将一个氨基酸的查询序列与一个蛋白质序列数据库相比较。这类搜索有专门与蛋白质搜索相关的可选参数，如对各种PAM和BLOSUM打分矩阵的选择。tblastn(translatedBLAST)：将一个蛋白质查询序列与一个以所有阅读框动态翻译成蛋白质的核酸序列数据库进行比较。可以用此程序来判断一个DNA数据库是否编码所感兴趣的查询蛋白。用RBP查询是否可以在某个已测序的DNA数据库中找到匹配项呢？2020/10/2819Step2:ChoosetheBLASTprogramblastx(translatedBLAST)：将一个核酸的查询序列按所有可能的阅读框翻译后的序列与一个蛋白质序列数据库进行比较。如若有一个DNA序列，想知道它编码什么蛋白质，用此程序进行搜索。它会自动将DNA翻译成6种可能的蛋白质。然后此程序就会将翻译的6个蛋白质序列逐一与蛋白质序列数据库中的各个成员进行比较。tblastx(translatedBLAST)：将一个核酸查询序列的6种框架的翻译结果与一个核酸序列数据库的6种框架翻译产物进行比较。该程序不能使用BLAST网页上提供的主要的去冗余(nr)数据库，因这一操作很消耗计算机资源。2020/10/2820ChoosetheBLASTprogramProgram

Input

Database

1blastn

DNA

DNA

1blastp

protein

protein

6blastx

DNA

protein

6tblastn

protein

DNA

36tblastx

DNA

DNA2020/10/2821DNApotentiallyencodessixproteins5’CATCAA5’ATCAAC5’TCAACT5’GTGGGT5’TGGGTA5’GGGTAG5’CATCAACTACAACTCCAAAGACACCCTTACACATCAACAAACCTACCCAC3’3’GTAGTTGATGTTGAGGTTTCTGTGGGAATGTGTAGTTGTTTGGATGGGTG5’2020/10/2822Step3:choosethedatabaseBLAST搜索可使用的数据库会列在每一个BLAST页面上，对于蛋白质数据库搜索(blastp和blastx)，两个主要的选择即nr数据库和SwissProt。nr=non-redundant(mostgeneraldatabase)dbest=databaseofexpressedsequencetagsdbsts=databaseofsequencetagsitesgss=genomicsurveysequenceshtgs=highthroughputgenomicsequence2020/10/2823Step3:choosethedatabasenr数据库是合并了若干个主要的蛋白质或DNA数据库得到的。这些数据库中经常包含有相同的序列，但nr数据库只收录其中的一个序列(即使在nr数据库中出现看上去一样的序列，实际上还是具有一些细节上的区别)。nr数据库是在要搜索现有的绝大多数序列时典型和常用的数据库。2020/10/2824去冗余GenBank编码序列PDB+SwissProt+PIR+PRF2020/10/28252020/10/2826当确定了要输入的序列和要搜索的数据库之后，还有10个其他的可选参数要确定。①LimitbyEntrezQuery：任何NCBIBLAST搜索的范围都可以用在Entrez搜索中使用的任何一种范围限定词来限定。Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2827Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2828Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2829②Maxtargetsequences：比对之后显示的最大的比对序列的数目。Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2830Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2831③期望expect：期望值E是得分大于或等于某个分值S的不同的比对的数目在随机的数据库搜索中发生的可能性。这个数值表示你仅仅因为随机性造成获得这一联配结果的可能次数。对于blastn、blastp、blastxt和blastn期望值的默认设置是10。在这个E值下，随机出现得分等于或高于比对得分S的期望数为10个(这里是假设用与实际的查询序列长度相等的随机的查询序列搜索数据库)。当将期望选项值调小时，返回的数据库搜索结果将变少，匹配被搜索到的概率也会变小。增大E值将返回更多的结果。Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2832Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2833Step4a:选择可选的搜索参数Selectoptionalsearchparameters④字段长度wordsize：BLAST程序是通过比对未知序列与数据库序列中的短序列来发现最佳匹配序列的。最初进行“扫描”(scanning)就是确定匹配片段。序列的匹配程序由短序列(定义为“word”,即字)的联配得分总和来决定。联配时，“字”的每个碱基均被计分：如果碱基对完全相同(如A与A)，得某一正值；如果碱基对不很匹配(W与A或T)，则得某一略小的正值；如果两个碱基不匹配，则得一负值。总的合计得分便决定了序列间的相似程度。对于蛋白质搜索，窗口大小可以被设定为3(默认值)或者2。当用一个查询序列来进行数据库搜索时，BLAST算法首先将查询序列分割成一系列具有特定长度(字段长度)的小的序列段(字段)。对于blastp，更大的字段长度将得到更高的搜索精度。对于任意的字段长度，每个字段的匹配结果将被延伸以得到BLAST的输出结果。实际应用中对于蛋白质搜索很少需要改变字段的长度。

对于核酸序列，默认的字段长度是11，BLAST的字长缺省值为11，即BLASTN将扫描数据库，直到发现那些与未知序列的11个连续碱基完全匹配的11个连续碱基长度片段为止。然后这些片段(即字)被扩展。11个碱基的字长已能有效地排除中等分叉的同源性和几乎所有随机产生的显著联配。它可以被增大(15)或减小(7)。降低字段长度将会使搜索变得更准确同时也会变得更慢。2020/10/2834Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2835⑤

矩阵matrix：对于blastp的蛋白质-蛋白质搜索有5种氨基酸替代矩阵：PAM30,PAM70,BLOSUM45,BLOSUM62(默认值)以及BLOSUM80.一些其他的BLAST服务器还提供了很多其他的替代矩阵，如PAM250。通常情况下明智的选择是在一次BLAST搜索中使用几种不同的打分矩阵。

Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2836PAM1矩阵：Dayhoff和同事利用可接受点突变的数据和每个氨基酸的发现频率产生突变概率矩阵M。矩阵元素Mij表示在一给定进化时期内氨基酸j(列)替换成氨基酸i(行)的概率。进化时期为一个PAM（PAM定义为进化趋异的单位，表示两个蛋白1%氨基酸发生变化的时间）。PAM1矩阵基于紧密相关蛋白质的比对，这些蛋白质家族内的序列一致程度至少有85%。除PAM1矩阵外的其他PAM矩阵是如何得来的？Dayhoff等用PAM1矩阵乘以自身数百次，得到其他PAM矩阵。如PAM250矩阵就是PAM1矩阵乘以自身250次产生，是BLAST搜索数据库的常用矩阵之一。Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2837Dayhoff’sPAM1mutationprobabilitymatrix原始氨基酸Eachelementofthematrixshowstheprobabilitythatanoriginalaminoacid(top)willbereplacedbyanotheraminoacid(side)替代氨基酸2020/10/2838PAM250mutationprobabilitymatrixTop:originalaminoacidSide:replacementaminoacid2020/10/2839⑤PAM0矩阵：矩阵将成为单位矩阵，因没有氨基酸发生变化。PAM∝矩阵：PAM相当大(如PAM＞2000或矩阵和自己相乘无数次)。每种氨基酸等概率出现，每行的所有值都接近于一个数值，这个数值就是氨基酸的出现频率。Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2840Dayhoff’sPAM0mutationprobabilitymatrix:therulesforextremelyslowlyevolvingproteinsTop:originalaminoacidSide:replacementaminoacid2020/10/2841Dayhoff’sPAM2000mutationprobabilitymatrix:therulesforverydistantlyrelatedproteinsPAMAAlaRArgNAsnDAspCCysQGlnEGluGGlyA8.7%8.7%8.7%8.7%8.7%8.7%8.7%8.7%R4.1%4.1%4.1%4.1%4.1%4.1%4.1%4.1%N4.0%4.0%4.0%4.0%4.0%4.0%4.0%4.0%D4.7%4.7%4.7%4.7%4.7%4.7%4.7%4.7%C3.3%3.3%3.3%3.3%3.3%3.3%3.3%3.3%Q3.8%3.8%3.8%3.8%3.8%3.8%3.8%3.8%E5.0%5.0%5.0%5.0%5.0%5.0%5.0%5.0%G8.9%8.9%8.9%8.9%8.9%8.9%8.9%8.9%Top:originalaminoacidSide:replacementaminoacid2020/10/2842Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2843Step4a:Selectoptionalsearchparameters2020/10/28442020/10/2845Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2846⑥Compositionaladjustments：这个选项是默认选择的，一般来说可改善E值的统计计算和提高灵敏度(减少返回的假阳性结果的数目)。Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/28472020/10/2848⑦选择过滤条件Filter：过滤器将锁定诸如组成低复杂(lowcompositionalcomplexity)序列区(如Alu序列)，用一系列N(NNNNNN)替代这些程序。N代表任意碱基(IUB-code)。只有未知待检序列被过滤替代，而数据库的序列将不被过滤。过滤对绝大多数序列都是有益的，“Filter”项的缺省选项为ON。例如，多A碱基的尾部和脯氨酸富积的序列，会得到人为的高联配得分而误导分析。这是因为这类序列数量极大，遍布整个基因组，直至整个数据库。Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2849⑦选择过滤条件Filter：

过滤选项的调整对于BLAST搜索的结果将会产生显著的影响。当使用blastpnr来查询一个人类富含脯氨酸的蛋白质(NP_036522)时，会得到8个数据库中的匹配项，而其中两个具有统计显著性。而当关掉过滤选项后，就会得到包含其他一些富含脯氨酸的蛋白质在内的23个匹配。一般来讲，过滤可以帮助避免那些假的数据库匹配，但在某些情况下一些可信的匹配也可能会被过滤掉。Step4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/28502020/10/2851filtering2020/10/2852以一种富含脯氨酸的人类唾液蛋白NP_036522作为查询项使用blasstpnr进行搜索的结果。该蛋白质的大部分没有在数据库中的匹配，这是由于该蛋白的中间区域被默认的过滤功能给过滤掉了2020/10/2853关掉过滤选项后搜索的结果。这次得到了更多的与在蛋白质的全序列范围内相匹配的数据库返回结果。第一个比对显示了该蛋白极度富含脯氨酸的特点2020/10/2854NCBIblastnowoffersmaskingaslowercase/colored2020/10/2855Entrez!FilterScoringmatrixWordsizeExpectStep4a:选择可选的搜索参数Selectoptionalsearchparameters2020/10/2856BLAST:optionalparametersYoucan...•choosetheorganismtosearch•turnfilteringon/off•changethesubstitutionmatrix•changetheexpect(e)value•changethewordsize•changetheoutputformat2020/10/2857Step4b:选择可选格式参数optionalformattingparametersBLAST搜索有很多控制输出格式的参数：

Alignmentview

Descriptions

Alignments2020/10/2858Step4b:选择可选格式参数optionalformattingparametersBLAST搜索结果的顶部

顶部提供关于该搜索的详细信息:BLAST搜索的类型、关于查询内容和所搜索的数据库的描述以及一个分类连接可以将结果按照物种进行分类。2020/10/2859databaseprogramquerytaxonomy结果网页2020/10/2860taxonomy2020/10/2861Step4b:选择可选格式参数optionalformattingparametersBLAST搜索结果的中间部分显示的是数据库中序列与查询序列相匹配的项的列表。简明图形提供了用不同颜色表示的搜索结果的概况。图下面的每一个条带表示数据库中的一个与查询序列相匹配的蛋白质或核酸序列，被标以不同颜色表示亲缘关系的远近(根据比对的分)，最接近匹配用红色表示。每条线的长度对应于该序列与查询序列比对上的区域大小。所有比对用一个被称为“描述”的单行小结列表描述。按照E值增加的顺序排列。显著性最强的匹配将位于顶端。2020/10/2862结果的图示输出数据库序列的列表Highscoreslowevalues2020/10/2863Step4b:选择可选格式参数optionalformattingparametersBLAST搜索结果的靠下面部分显示的是一系列的两两序列比对可检查查询序列（输入序列）与对象序列（如和查询序列比对的的特定的数据库匹配）之间的比对情况。

4种衡量的分数：比特分数、期望分数、一致性百分比、正性(相似性百分比)2020/10/28642020/10/2865Step4b:选择可选格式参数optionalformattingparameters可以不用整体地进行BLAST搜索而仅通过改变格式选项来提供一些不同的输出结果一些选项可以把比对序列显示成多序列比对的形式，这对于确定一个蛋白质或DNA家族中的保守的或趋异的氨基酸残基非常有用。

2020/10/28662020/10/2867BLASTformatoptions2020/10/2868BLASTformatoptions2020/10/28692020/10/28702020/10/2871thresholdscore=11EVDparametersBLOSUMmatrixEffectivesearchspace=mn=lengthofqueryxdblength10.0istheEvaluegappenaltiescut-offparameters2020/10/2872如何评价搜索结果的显著性如何处理过多的结果如何处理过少的结果多结构域蛋白HIV-1pol的BLAST检索BLAST检索脂质运载蛋白

lipocalins

:改变打分矩阵的作用BLAST搜索的一些策略2020/10/2873BLAST搜索策略总图.2020/10/2874如何评价搜索结果的显著性同源性是基于由搜索结果的统计评估支持的序列相似性推出的。指认基因或蛋白质之间的相似性的问题不能单靠序列分析解决，还需要使用生物学上的标准来支持同源性的推导。可以用对于蛋白质的结构和功能的评估来补充BLAST的结果。真正有亲缘关系的蛋白质的序列可以有很大的差别，即使是那些有相近三维结构的蛋白质。因此，认为数据库搜索和蛋白质比对可能会得出一定数量的假阴性匹配。如脂质运载蛋白家族很多成员如RBP4和气味分子结合蛋白OBP只有非常有限的序列一致性，但它们的三维结构关系非常近，而且它们作为疏水性配体的载体的功能也被认为是相同的。BLAST搜索的一些策略2020/10/2875SometimesarealmatchhasanEvalue>12020/10/2876SometimesasimilarEvalueoccursforashortexactmatchandlonglessexactmatch2020/10/2877决定两个蛋白质或DNA序列是否同源，可以问下面的问题：期望值究竟是不是显著？两个蛋白质是不是具有近似的大小？实际上同源蛋白质不要求具有相同的大小，两个蛋白质可能只共享有限的一个相同的结构域。但是，产生一个对两个蛋白质同源的可能性的生物学上的直觉是非常重要的。比如一个1000个氨基酸具有跨膜结构域的蛋白质不太可能与RBP蛋白同源，因绝大多数脂质运载蛋白长度都近似200aa(20-25kD)。如何评价搜索结果的显著性2020/10/2878决定两个蛋白质或DNA序列是否同源，可以问下面的问题：这两个蛋白是否有共同的基序或信号序列？是。孕激素相关子宫内膜蛋白和RBP4蛋白都含有一个GXW载脂蛋白标签。这两个蛋白质是不是一个合理的多序列比对的一部分？这两个蛋白质是否共有一个相似的生物学功能？就像所有的载脂蛋白一样，这两个蛋白都是小的、亲水性的含量丰富的分泌分子。这两个蛋白质是否具有相似的三维结构？虽然载脂蛋白序列多种多样，它们却共享一个显著的非常保守的结构。这个结构(一个杯状的盂)可使它们将疏水性配体运输通过一个水环境区间。载脂蛋白如何评价搜索结果的显著性2020/10/2879决定两个蛋白质或DNA序列是否同源，可以问下面的问题：如果BLAST搜索得到一个对另一个蛋白质的边缘匹配，以这个具有较远亲缘关系的蛋白质作为查询项再进行一次新的搜索。如用孕激素相关子宫内膜蛋白进行blastpnr搜索可得到一些在RBP4蛋白的搜索中也检测到的蛋白。这个发现增加了认为这两个蛋白实际上是一个蛋白超家族中的同源成员的信心。载脂蛋白如何评价搜索结果的显著性2020/10/2880AssessingwhetherproteinsarehomologousRBP4andPAEP:Lowbitscore,Evalue0.49,24%identity(“twilightzone”).Buttheyareindeedhomologous.TryaBLASTsearchwithPAEPasaquery,andfindmanyotherlipocalins.2020/10/2881Theuniverseoflipocalins(eachdotisaprotein)retinol-bindingproteinodorant-bindingproteinapolipoproteinD2020/10/2882BLASTsearchwithPAEPasaqueryfindsmany