本地化BLAST+进行目标序列比对

举个栗子🌰:

最近对拟南芥中的一个叫做CBF的基因非常感兴趣,现在我想看看这个CBF基因在一个新测序的物种Am基因组上是否存在,Am中有几个CBF基因?这样的查找可以通过本地化的BLAST+来实现。

关于BLAST:
BLAST软件能够将一个目标蛋白或核苷酸序列(称为查询query)与一个数据库进行比较,并识别某个特定阈值以上的与目标序列相似的库序列。在线的BLAST功能只需要将序列输入,然后BLAST,最后会给出相似的库序列(按照同源性高低排列)。在线BLAST明显存在两个缺陷:1.目标序列BLAST后的得到的库序列是多个物种的,如果我们只是要特定物种的话,还需要再自己找。2.如果目标序列有很多,则需要执行很多次,效率低。而本地化的BLAST+能很好的解决上面两个问题。

开始吧:

1.BLAST+的下载与安装

搜索BLAST,点击网页中的下载选项,针对不同的平台(linux,mac,windows)选择性的进行安装,如图一。
下面以mac为例进行介绍:


Screenshot 2020-04-06 at 14.45.38.png
Screenshot 2020-04-06 at 14.46.53.png
Screenshot 2020-04-06 at 14.47.15.png

安装好之后打开Terminal,输入"blastn",出现下图说明安装成功。


Screenshot 2020-04-06 at 14.51.21.png
  1. 数据准备

BLAST可以搜索核酸或者蛋白序列,下面以蛋白序列为例,之所以用蛋白序列是因为后期我想要用蛋白序列画进化树~

首先准备目标序列的fasta文件,即我感兴趣的CBF基因的fasta文件。首先在NCBI上下载到拟南芥的CBF基因的蛋白序列,拟南芥共有三个CBF基因,分别是CBF1,CBF2 和CBF3。PS:fasta文件可以用Sublime Text文件查看。

Screenshot 2020-04-06 at 15.06.31.png

接下来准备目标物种的全部蛋白序列,因为文章还没发表,所以我用的物种名称用Am代替。已发表基因组的物种蛋白序列可以在NCBI上下载,同样是fasta格式。

  1. 运行BLAST

回到Terminal,开始运行BLAST。

①设置工作路径:cd /Users/mac/Desktop/CBF-Tree/ (此处根据自己文件存放的文件夹做调整)

②自定义数据库
也称为格式化数据库。在本例中,是指将Am的蛋白序列设为目标数据库:makeblastdb -in Am_pro.fa -dbtype prot
参数说明:
-in 待格式化的序列文件
-dbtype 数据库类型,prot或nucl分别表示蛋白库或核酸库

Screenshot 2020-04-06 at 15.24.01.png

③目标序列与库序列进行比对
在本例中,是指将拟南芥的CBF的蛋白序列与Am的蛋白序列进行比对:blastp -query AtCBF.fasta -db Am_pro.fa -out AtCBF_Am.txt -evalue 1e-5 -outfmt 6
参数说明:
-query:输入文件名,代表的是要进行blastp的拟南芥CBF蛋白序列文件
-db:为前面格式化了的Am的基因组蛋白序列文件
-out:输出文件名
-evalue:为筛选标准(设置输出结果的evalue值,evalue越低,相似性越高)
-outfmt:输出文件格式,总共有12种格式,6是tabular格式对应BLAST的m8格式

Screenshot 2020-04-06 at 15.39.07.png

outfmt 6生成的表格中每一列分别是:
第一列为: Query(递交序列),
第二列为: 数据库序列(目标序列subejct),
第三列为: identity
第四列为:比对长度
第五列为:错配数
第六列为:gap数
第七列和第八列为:Query开始碱基位置和结束碱基位置
第九列和第十列为:Subject开始碱基位置和结束碱基位置
第十一列为:期望值
第十二列为:比对得分

Q&A:返回的txt文件是空文件怎么办?
blastn 最后一个参数 -outfmt 6是将结果变成一个表格输出来。如果query序列在目标数据库里没有同源基因或者没有同源序列,就没有匹配结果,这样-outfmt 6就无法用表格输出来,这样输出的txt文件会显示为空文件。解决方案是将-outfmt 6这个参数去掉,这样返回的是你要的txt文件,文件中会显示no hits。

Q&A: 用blast进行短序列搜索(20nt左右)
如果query序列是很短的序列,如为了扩增某个基因设计了一个20bp的引物,要看看这个引物对目标基因是否具有唯一性的时候,用上述代码是跑不出结果的,这时候要加上另外的参数:-task blastn-short -word_size 7 -evalue 1,具体情况请参照参考文献3。

  1. R studio中提取输出文件中同源性最高的序列

输出的结果按照同源性从高到低排列(目标序列与库序列进行同源比对时根据evalue筛选出来的可能不止一个),那么如何选出每个最匹配的呢?

打开R studio,输入以下代码,即可实现。

Screenshot 2020-04-06 at 15.51.43.png

图中是生成的结果。

Screenshot 2020-04-06 at 15.52.31.png

补充资料:
BLAST+与BLAST相比,有很多改进和提高,NCBI强烈推荐放弃BLAST,使用BLAST+, 这里说的BLAST和BLAST+,都是本地的,与之前的那个批量BLAST小程序不是一回事。BLAST下载地址:NCBI BLAST+ 。BLAST+的一般用法如下:

格式化数据库
makeblastdb -in db.fasta -dbtype prot -parse_seqids -out dbname
参数说明:
-in:待格式化的序列文件
-dbtype:数据库类型,prot或nucl
-out:数据库名

蛋白序列比对蛋白数据库(blastp)
blastp -query seq.fasta -out seq.blast -db dbname -outfmt 6 -evalue 1e-5 -num_descriptions 10 -num_threads 8
参数说明:
-query: 输入文件路径及文件名
-out:输出文件路径及文件名
-db:格式化了的数据库路径及数据库名
-outfmt:输出文件格式,总共有12种格式,6是tabular格式对应BLAST的m8格式
-evalue:设置输出结果的e-value值
-num_descriptions:tabular格式输出结果的条数
-num_threads:线程数

核酸序列比对核酸数据库(blastn)以及核酸序列比对蛋白数据库(blastx)
与上面的blastp用法类似:
blastn -query seq.fasta -out seq.blast -db dbname -outfmt 6 -evalue 1e-5 -num_descriptions 10 -num_threads 8
blastx -query seq.fasta -out seq.blast -db dbname -outfmt 6 -evalue 1e-5 -num_descriptions 10 -num_threads 8

以上的参数说明只是一些常用的参数,完整的参数说明可以用-help查询。

参考资源:

  1. 本地化的NCBI BLAST+的安装及使用
    链接:https://www.jianshu.com/p/2f551c0f9060
    来源:简书

  2. BLAST+使用方法
    链接:https://blog.csdn.net/lixiangyong123/article/details/67634935
    来源:CSDN

3.用blast进行短序列(20nt左右)的搜索
链接:https://www.jianshu.com/p/a6ee3d19b832
来源:简书

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 206,126评论 6 481
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 88,254评论 2 382
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 152,445评论 0 341
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 55,185评论 1 278
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 64,178评论 5 371
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,970评论 1 284
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,276评论 3 399
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,927评论 0 259
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 43,400评论 1 300
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,883评论 2 323
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,997评论 1 333
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,646评论 4 322
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,213评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,204评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,423评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,423评论 2 352
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,722评论 2 345

推荐阅读更多精彩内容