教你使用数据库

演示文稿教你使用数据库第一页，共61页。（优选）教你使用数据库第二页，共61页。生物信息学第三讲：序列的采集、存储和查询

第三页，共61页。本章内容提要1.DNA测序2.序列数据的存储3.序列数据的文件格式4.序列数据的查询第四页，共61页。1.DNA测序1.

DNA一次测序的长度为500bp。2.基因组的测序方案：将大的染色体打断成100kbp的片断，插入到BAC(BacterialArtificialChromosome)中。再随机打断，克隆，然后再组装成长的序列(contig)。3.EST(Expressedsequencetag)测序：细胞中mRNA反转录成cDNA，方向不定，测序。4.UniGene:为每一个基因创造一个唯一的条目，收集这个基因所有的ESTs.5.GSS(基因组测序序列)：类似于ESTs，来源基因组。6.HTG(高通量基因组序列)：高通量、尚未完工的DNA序列。第五页，共61页。基因组测序：两种方案1.DNA片段在染色体上的位置、方向已知。首先染色体被打断成150kbp左右的片段，然后克隆到BACs中，再进一步打碎，克隆，测序，组装。2.“鸟枪法”，shotgun，随机将DNA片段打碎，克隆，测序，组装。DNA片段在染色体上的位置和方向未知。第六页，共61页。第七页，共61页。2.序列数据的存储1.核酸三大数据库：GenBank,EBI,DDBJ.2.Ensembl数据库：基因组注释。3.ESTs数据库；4.UniGene数据库5.Refseq数据库；6.NCBI的Gene信息数据库;7.蛋白质序列：Swissprot/TrEMBL/UniProt数据库。第八页，共61页。三大数据库之间的联系第九页，共61页。GenBank中测序最多的20个物种161.0版，2007.08第十页，共61页。dbEST:表达序列标签数据库2007.08，总序列45,660,524条，最多的20个物种如下第十一页，共61页。UniGene:AnOrganizedViewoftheTranscriptom第十二页，共61页。Refseq数据库第十三页，共61页。第十四页，共61页。NCBIGene1./sites/entrez?db=gene2.序列从Refseq数据库中得到;3.详尽的注释信息，包括基因在基因组的定位，基因名称、蛋白质名称，基因结构，等等。第十五页，共61页。UniProt1.专家审核的蛋白质序列数据与知识库；2.UniProtKnowledgebase：Release12.1，2007.083.包括：Swiss-ProtRelease54.1of21-Aug-2007:277883entries;TrEMBLRelease37.1of21-Aug-2007:4754787entries第十六页，共61页。Swiss-ProtRelease54.1第十七页，共61页。TrEMBLRelease37.1第十八页，共61页。3.序列数据的文件格式1.DNA/RNA/氨基酸代码的标识2.GenBank数据格式3.UniProt4.FASTA第十九页，共61页。GenBank数据格式Definition:标题序列长度数据类型Accessionnumber版本号GInumber第二十页，共61页。GenBank的数据类型第二十一页，共61页。UniProt数据格式Accessionnumber第二十二页，共61页。FASTA格式第二十三页，共61页。4.序列数据的查询某天，Prof.Gene发现人的Hela细胞中，有丝分裂期间有异常情况：细胞不再分裂，而是开始凋亡(表型，phenotype)，通过实验的方法(例如，酵母双杂交)，发现了与有丝分裂期间某个蛋白可能相互作用的一个基因，测序结果如下(genotype)：CCCCTGCCTGGCAGCCCTTTCTCAAGGACCACCGCATCTCTACATTCAAGAACTGGCCCTTCTTGGAGGGCTGCGCCTGCACCCCGGAGCGGATGGCCGAGGCTGGCTTCATCCACTGCCCCACTGAGAACGAGCCAGACTTGGCCCAGTGTTTCTTCTGCTTCAAGGAGCTGGAAGGCTGGGAGCCAGATGACGACCCCATAGAGGAACATAAAAAGCATTCGTCCGGTTGCGCTTTCCTTTCTGTCAAGAAGCAGTTTGAAGAATTAACCCTTGGTGAATTTTTGAAACTGGACAGAGAAAGAGCCAAGAACAAAATTGCAAAGGAAACCAACAATAAGAAGAAAGAATTTGAGGAAACTGCGGAGAAAGTGCGCCGTGCCATCGAGCAGCTGGCTGCCATGGATTGAGGCCTCTGGC第二十四页，共61页。问题：1.这是哪个基因？2.编码的蛋白质序列是怎样的？3.有没有保守的功能结构域(domain)？4.它的功能是怎样的？5.它在真核生物中保守吗？6.有没有三级结构信息？第二十五页，共61页。NCBI:BLAST/第二十六页，共61页。nucleotideblast第二十七页，共61页。Megablast:找基因序列输入序列第二十八页，共61页。提交序列第二十九页，共61页。NM_001168.2：SurvivinUniGeneGeo:基因表达信息Geneinfo：基因信息第三十页，共61页。Survivin第三十一页，共61页。Geneinfo:17号染色体第三十二页，共61页。功能注释：GeneOntology第三十三页，共61页。结论11.该基因为人的Survivin基因，染色体定位：17号染色体，73721872-73733311；基因标识符：NM_001168.2;2.初步的功能分析：细胞周期，caspase酶的抑制因子，等等。第三十四页，共61页。NM_001168.2：Survivin第三十五页，共61页。HumanSurvivin！第三十六页，共61页。获取蛋白质的序列信息第三十七页，共61页。Survivin：142aa第三十八页，共61页。结论2人的Survivin蛋白质包含142个氨基酸，序列标识符为：NP_001159.2第三十九页，共61页。获取FASTA序列第四十页，共61页。FASTA格式的序列第四十一页，共61页。PHI-BLAST:finddomain第四十二页，共61页。填入蛋白质的FASTA序列第四十三页，共61页。BIRdomain第四十四页，共61页。结论3Survivin具有保守的功能结构域BIR第四十五页，共61页。UniProt:蛋白质数据库第四十六页，共61页。在人中做BLAST搜索第四十七页，共61页。Survivin：O15392第四十八页，共61页。Survivin的蛋白质信息第四十九页，共61页。功能注释第五十页，共61页。结论4：功能分析1.在瘤形成过程中可能起一定作用；2.阻碍G2/M期的细胞编程性凋亡；3.Chromosomalpassengercomplex(CPC)的成员之一。…细胞亚定位：胞质，核。第五十一页，共61页。人的Survivin在酵母中有同源物吗？第五十二页，共61页。提交序列第五十三页，共61页。在酵母中进行序列比对第五十四页，共61页。酵母BIR1:P47134