1.Blast
Last updated
Was this helpful?
Last updated
Was this helpful?
同源搜索是生物信息最基本的分析之一。blast在很多年来都是用来实现这一任务的一个很受欢迎的工具。
在本节中,我们将学习如何利用网页版和本地版的blast进行蛋白和核酸的同源搜索。
fasta文件是存储生物序列最常见的一种文件格式。下面就展示了一个包含一条蛋白质序列的fasta文件的内容:
一个fasta文件可以包含多条序列。在fasta文件中,以">"开头的行记录了序列id等相关的信息,这一行之后直到下一个以">"开头的行之前,记录着该序列id对应的序列数据。
一条序列可以横跨多行,但是在这种情况下每行包含序列的最大长度通常是固定的。
fasta
contains gene name and sequence
gene name start with >
blastp
query aligned to database
-
以下步骤均在 /home/test/blast/
下进行
利用 blastp
进行蛋白质比对
VIM.fasta
与 NMD.fasta
分别是金属beta酶家族的两个蛋白的序列。
利用 blastn
进行核酸序列比对
H1N1-HA.fasta
与 H7N9-HA.fasta
是两种流感病毒的序列。
我们也可以在命令行搜索ncbi的远程数据库。在网络环境不好的时候容易报错,可以跳过。
例如:在yeast基因组序列中搜索Yeast.fasta序列
-dbtype
: 待建库的类型(核酸:nucl
,蛋白:prot
)
-in
: 待建库的序列文件
-out
: 数据库前缀
你可以利用 more
, less
等命令或者利用vi
等文本编辑工具查看结果文件。
Note: Docker 中已经装好
Ubuntu 自动安装软件方法: sudo apt-get install ncbi-blast+
这里 blast 由
ncbi-blast+
提供
寻找类似如下文件:
32位计算机(老机器)安装文件:如 ncbi-blast-2.2.28+-ia32-linux.tar.gz
64位机器安装文件: 如 ncbi-blast-2.2.28+-x64-linux.tar.gz
可以通过
uname -a
查看机器类型是64还是32位
Input the following sequences:
Select "Non-redundant Protein Sequences (nr)"
Select "blastp"
对于序列MSTRSVSSSSYRRMFGGPGTASRPSSSRSYVTTSTRTYSLGSALRPSTSRSLYASSPGGVYATRSSAVRL
:
1) 请使用网页版的 blastp, 将上面的蛋白序列只与 mouse protein database 进行比对, 设置输出结果最多保留10个, E 值最大为 0.5。将操作过程和结果截图,并解释一下 E value和 P value 的实际意义。
2) 请使用 Bash 脚本编程:将上面的蛋白序列随机打乱生成10个, 然后对这10个序列两两之间进行 blast 比对,输出并解释结果。(请上传bash脚本,注意做好重要code的注释;同时上传一个结果文件用来示例程序输出的结果以及你对这些结果的解释。)
注: blast的网站会提供多个mouse的databases,可以任选1个进行比对;也可以重复几次,每次选一个不同的database看看不同的输出结果,可以在作业中比较和讨论一下输出结果不同的原因。
3)解释blast 中除了动态规划(dynamic programming)还利用了什么方法来提高速度,为什么可以提高速度。
Figure 83. Atlas of Protein Sequence and Structure, Suppl 3, 1978, M.O. Dayhoff, ed. National Biomedical Research Foundation, 1979
Jim Kent
Jim Kent 1960年生于夏威夷,在旧金山长大。
Kent的编程生涯始于23岁,1983年Kent开始在Island Graphic Inc工作。当时Kent为Amiga家用电脑编写了一个能够结合3D渐变效果与2D简单动画的软件。
1985年他干脆自己开了一家软件公司 Dancing Flame, 将之前自己写的动画程序创建成为一个汇集各种动画和绘画功能的程序Cyber Paint,为CAD-3D开发提供便利。这是是第一个允许用户跨时间制作压缩视频的软件。之后他还开发了软件Autodesk Animator用于个人电脑中,可以为各种视频游戏创作艺术作品。 2000年Kent在攻读加州大学圣克鲁兹生物学博士学位,此时人类基因组计划进行到后期,他写出了GigAssembler, 帮助人类基因组项目能够拼装和发布人类基因组序列。
在GigAssembler之后,Kent继续编写BLAT(BLAST-like alignment tool)和参与维护UCSC Genome Browser 来帮助分析重要的基因组数据。目前,Kent仍然在UCSC工作,主要是在网络工具上帮助理解人类基因组。他帮助维护和升级浏览器,并致力于比较基因组学, Parasol,UCSC kilocluster的运营和ENCODE项目。
docker images的下载链接如所示
首先,进入NCBI官网在线,显示如下图所示界面,选择一种序列比对类型,然后根据提示进行序列比对。
4)我们常见的PAM250有如下图所示的两种(一种对称、一种不对称),请阅读一下 "Symmetry of the PAM matrices" @ ,再利用Google/wikipedia等工具查阅更多资料,然后总结和解释一下这两种(对称和不对称)PAM250不一样的原因及其在应用上的不同。
PAM matrices are also used as a scoring matrix when comparing DNA sequences or protein sequences to judge the quality of the alignment. This form of scoring system is utilized by a wide range of alignment software including . — wikipedia
:
:
see Videos in the ****