如何查找基因的序列(带图表)
- 格式:pdf
- 大小:2.24 MB
- 文档页数:34
gi|125806|sp|P01658|KV3F_MOUSE IG KAPPA CHAIN V-III REGION ... 33
0.30
gi|125808|sp|P01659|KV3G_MOUSE IG KAPPA CHAIN V-III REGION ... 33
0.30
gi|1172451|sp|Q05793|PGBM_MOUSE Basement membrane-specific ... 33
Query: 181 Sbjct: 175
EYEDENLYEGLNLDDCSMYEDISRGLQGTYQDVGSLNIGDVQLEKP 226 +YEDENLYEGLNLDDCSMYEDISRGLQGTYQDVG+L+IGD QLEKP DYEDENLYEGLNLDDCSMYEDISRGLQGTYQDVGNLHIGDAQLEKP 220
如果想找的基因是第一个序列即isoform a, 就可以点击NM_001025366.1, 得到如下界面:
If the gene sequence is known, and you want to find the corresponding GenBank ID, use /BLAST/ Click on the type of nucleotide search you want and enter the sequence. Results will be displayed with the GenBank ID included. To view information about a gene sequence with a given GenBank ID, go to /entrez/query.fcgi? db=Nucleotide
35
0.095
gi|1730075|sp|P01625|KV4A_HUMAN IG KAPPA CHAIN V-IV REGION LEN 34 0.26
gi|1718188|sp|P52583|VGR2_COTJA Vascular endothelial growth... 33
0.28
gi|125833|sp|P06313|KV4B_HUMAN IG KAPPA CHAIN V-IV REGION J... 33 0.30
举例说明,例如:在2003年JBC的文章(Conditional Knock-out of Integrin-
linked Kinase Demonstrates an Essential Role in Protein Kinase B/Akt Activation)中出现了“calreticulin (GenBank accession number gi 16151096)”, 那么把“16151096”输入GO前面的文本框中,点“GO”,就可以找到该基因了(当 然包括基因序列等相关信息)。
0.052
gi|114764|sp|P15530|C79B_MOUSE B-cell antigen receptor comp... 36
0.064
gi|1718161|sp|P53767|VGR1_RAT Vascular endothelial growth f...
35
0.080
gi|125735|sp|P01681|KV01_RAT Ig kappa chain V region S211
1
2 3
比对结果会在10秒左右后出现, 注意结果很多,一直下拉右侧工具条,会看到很多不同意义的比对结果。
• Search BLAST (/BLAST/) for P11912
Database: All non-redundant(非冗余) sequences 6,507,231 sequences; 2,219,987,828 total letters
如果你只想获得序列(例如去设计PCR引物的时候),那就可以选择FASTA,这样 就得到了FASTA格式的序列文件,没有其他数字和格式的干扰。
这就是FASTA格式的序列:
练习 1: 请在Entrez的蛋白质数据库中查找,gi (NCBI GenBank
ID)为13375618所登录的蛋白质的名称,并且回答该蛋白含
0.41
直接往下拉,会出现下面的画面—直接的每一组比对结果
>gi|126779|sp|P11911|C79A_MOUSE B-cell antigen receptor complex associated protein alpha-chain precursor (IG-alpha) (MB-1 membrane glycoprotein) (Surface-IGM-associated protein) (Membrane-bound immunoglobulin associated protein) (CD79A) Length = 220
e-133
gi|728993|sp|P40293|C79A_BOVIN B-cell antigen receptor comp...
312
3e-85
gi|126779|sp|P11911|C79A_MOUSE B-cell antigen receptor comp... 278
5e-75
gi|728994|sp|P40259|C79B_HUMAN B-cell antigen receptor comp... 55
Score = 278 bits (711), Expect = 5e-75 Identities = 150/226 (66%), Positives = 165/226 (73%), Gaps = 6/226 (2%)
Query: 1 Sbjct: 1
MPGGPGVLQALPATIFLLFLLSAVYLGPGCQALWMHKVPASLMVSLGEDAHFQCPHNSSN 60
1e-07
gi|125781|sp|P01618|KV1_CANFA IG KAPPA CHAIN V REGION GOM 38
0.019
gi|125361|sp|P17948|VGR1_HUMAN Vascular endothelial growth ... 37
0.042
gi|549319|sp|P35969|VGR1_MOUSE Vascular endothelial growth ... 36
Too much!!!
点击箭头所指的Limits,
出现了如下界面
Limits的意思其实就是高级检索,你可以在这里对检索词进行很多限制,这样能 大大精简查询结果。
在Limits这个界面,先选择查询的限定范围。人的VEGF基因,那就开始选择: 先选Gene name(基因名称);然后再选择Limit by Taxonomy(生物分类限定) 中的Homo sapiens(人类),然后再点击“GO”
直接往下拉,会出现下面的画面
• Distribution of Hits:
ห้องสมุดไป่ตู้
直接往下拉,会出现下面的画面
• Score E Sequences producing significant alignments:
Score E-Value
gi|547896|sp|P11912|C79A_HUMAN B-cell antigen receptor comp... 473
How to use the BLAST?
Example: NCBI Blast
Step 1: Run the Blast
NCBI BLAST /BLAST/
Step 2:Choose the program to use
Step 3: Input the data, 可以直接输入登录号,或者是gi number
可以通过查找原始数据库获得登录号。
• Search SWISSPROT for Immunoglobulin:
SWISS_PROT:C79A_HUMAN P11912
如果不知道登录号或者gi number,则可以直接粘贴已知序列
Step 4: Choose search set and program selection, and put the BLAST to start the search
MPGG
+ LL LS LGPGCQAL + P SL V+LGE+A C N+
MPGG----LEALRALPLLLFLSYACLGPGCQALRVEGGPPSLTVNLGEEARLTC-ENNGR 55
Query: 61 Sbjct: 56
NANVTWWRVLHGNYTWPPEFLGPGEDPNGTLIIQNVNKSHGGIYVCRVQEGNESYQQSCG 120 N N+TWW L N TWPP LGPG+ G L VNK+ G C+V E N ++SCG NPNITWWFSLQSNITWPPVPLGPGQGTTGQLFFPEVNKNTGACTGCQVIE-NNILKRSCG 114
下面就是点击VEGFA后出现的界面:
再往下看,可以看到Genomic regions, transcripts, and products,这里显示 了该基因在基因组中的位置,以及转录本的生成情况
就看见了目的基因的mRNA的链接(如NM_001025366.1)和蛋白质的链接 (如NP_001020537.2)
如何查找基因序列? ——在Genbank中寻找目的基因
1. 根据文献
搞reasearch肯定要读文献的,如果你曾经在文献中看到过你感兴趣的基因,而 且文中还提到了该基因在Genbank中的ID号,那就好办了,直接打开 ,在Search后的下拉框中选择Nucleotide,把 Genbank ID号输入GO前面的文本框中,点“GO”,就可以找到他了。