ncbi使用技巧,spss使用指南

1、生物信息技术应用,分子序列比对分析,Sequence alignment,Contents,序列比对的应用,序列数据库,基本类型： 初级数据库 收录、存储序列的基本数据资源，如核酸（蛋白质）序列、蛋白质空间结构及基因组信息。 次级数据库 在初级库资源基础上进行整理和标注，为特定专业领域服务的派生数据库，如表达序列标记、微列阵（基因芯片）、代谢和信号途径、遗传疾病数据库、免疫数据库等等。,核酸序列数据库,以核苷酸顺序及注释信息为基本内容的数据库 世界三大核酸数据库 GenBank in USA ( ) EMBL in Europe ( htt

2、p:/www.ebi.ac.uk/embl ) DDBJ in Japan ( http:/www.ddbj.nig.ac.jp ) 1998年，GenBank, EMBL, DDBJ共同成立国际核酸序列数据库协会(INSDC, )，实现了全球范围内的核酸序列的同步更新和交换互享。,蛋白质序列数据库,以蛋白质氨基酸顺序及注释信息为基本内容的数据库 世界主要蛋白质序列数据库 (1) PIR-PSD (Protein information resource-protein sequence database) founded by NBRF of USA

3、(美国国家医学研究基金会) in 1984 ( / ) 1988年，日本国际蛋白质信息库(JIPID)和德国慕尼黑蛋白质序列信息中心(MIPS)加入PIR，合作成立了国际蛋白质信息中心（PIR-International）。PIR为较全面和权威注释的蛋白质数据库，具有非冗余、高质量注释和分类全面等特点。,蛋白质序列数据库,世界主要蛋白质序列数据库 (2) SwissProt 1986年始创于瑞士日内瓦大学 ，现由瑞士生物信息学研究所(SIB)和欧洲生物信息学研究所(EBI)共同管理和维护。 ( http:/www.expasy.ch/sprot/

4、 ) SwissProt数据库的特点：所有序列条目经过专家核实，可靠性与可信度高；注释详细，包括蛋白质的功能、序列及结构域的结构、翻译后修饰及其位点、突变体等,蛋白质序列数据库,世界主要蛋白质序列数据库 (3) TrEMBL (translation of EMBL) 建立于1996年，是从EMBL中的cDNA序列翻译得到的蛋白质数据库。相似的还有GenPept (GenBank)数据库。 TrEMBL、GenPept数据库的特点：序列条目来自核酸序列库的翻译，即时性强；但未经专家的注释、分析和核实，因而错误率和冗余度都较高。,蛋白质序列数据库,全球统一的蛋白质序列与功能数据库 UniProt

5、 (Universal Protein Resource, 通用蛋白质资源) ( / ) 2002年，PIR、SIB、EBI合并了分属其下的PIR-PSD、Swiss-Prot和TrEMBL数据库，形成了统一的蛋白质数据库UniProt 截止2008年8月，UniProt共收录蛋白质序列 6,462,751个,生物大分子结构数据库,以生物大分子各原子空间信息为基本内容的数据库 给定序列的蛋白质如何折叠为稳定、具一定生物功能的三维结构？ 信息来源：对蛋白质晶体的X射线衍射、核磁共振及冷冻电镜分析 主要数据库：美国Brookhaven国家实验室的PDB (

6、Protein Data Bank, 1971年成立) 联合MSD-EBI、PDBj，于2003年组建全球共享的worldwide PDB (wwPDB) ( / ) 截止2008年8月，共收录蛋白质结构52684个,Research Collaboratory for Structural Bioinformatics,/pdb,GenBank,National Center for Biotechnology Information (NCBI) National Library of Medicine Nation

7、al Institutes of Health /,GenBank,/,全球著名的生物信息综合数据库 GenBank（1982年成立）包含所有已知的核苷酸及蛋白质序列以及与之相关的生物学信息和参考文献，现由美国生物技术信息中心（NCBI，成立于1988年）管理维护，是世界上的权威序列数据库。 数据库序列的来源为作者直接递交或间接查寻文献所得，并与世界上其他公开发行的数据库，如EMBL，DDBJ交换每日更新的数据。,GenBank,GenBank 核酸数据库的增长（1982-2008）,截止2

8、008年2月，GenBank 共收录核酸序列近8300万条，约860亿碱基对，来自26万余种生物,GenBank,/ Home page,NCBI 数据库检索系统 Entrez,Entrez为基于万维网的NCBI数据库检索系统，通过输入关键词，运用布尔算符，可在NCBI所有数据库中进行文本搜索。,NCBI 数据库检索系统 Entrez,Entrez为基于万维网的NCBI数据库检索系统，通过输入关键词，运用布尔算符，可在NCBI所有数据库中进行文本搜索。,NCBI 数据库检索系统 Entrez,Entrez 应用举例,检索主题：小鼠(mouse

9、)转录因子(TF) -E2F mRNA的核酸序列 检索过程：NCBI主页All Databases在Search对话框输入关键词”E2F”，点”Go”输出检索结果。,Entrez 应用举例,输出相关序列3784条，有待进一步筛选：限制物种来源(小鼠 Mus musculus)、限制分子为mRNA (排除未经实验验证的预测序列)。,Entrez 应用举例,最后命中168条，从中选中NM_148952，以待详细分析。,Entrez 应用举例,NM_148952： 小鼠E2F转录因子-4。,序列文件格式：FASTA 格式,标题行：文件的第一行，以大于号“”开始，不换行。内容可自定义，包括基本信息和简

10、单注释； 序列行：文件第二行起至结束，中间不得有空格。 FASTA为最简单的序列（核酸或蛋白质）格式,序列文件格式： GenBank 格式,GenBank 格式注释信息丰富全面，属文本文件，包括4部分： 头部：含名称、定义、识别码、物种来源等基本信息； 引文区：含相关文献信息。 序列特征表：含序列的编码区、非编码区、功能域、修饰或突变位点、翻译序列等众多注释信息 序列区：序列本身,GenBank 格式,Locus 行信息： Locus name； Sequence length; Type of Molecule: DNA, mRNA, cDNA Taxonomy：PRI primate(灵长

11、类) 、ROD rodent(啮齿类)、 MAM - other mammalian (其它哺乳类)、VRT - other vertebrate(脊椎动物）、INV invertebrate sequences、PLN - plant, fungal, and algal; BCT - bacterial; VRL - viral, PHG - bacteriophage, SYN - synthetic; UNA - unannotated; EST - expressed sequence tags Date: 上传或最近修改时间,GenBank 格式,GenBank 格式,GenBan

12、k 格式,/,成对序列比对与BLAST工具,序列比对 概念：通过比较两个或两个以上的核酸(蛋白质)序列，显示其中相似的结构区域。 成对序列比对 (pairwise sequence alignment) 多重序列比对 (multiple sequence alignment) 功能: (1) “相似”的序列 “相似”的三维结构； (2) “相似”的序列 “相似”的功能； (3) “相似”的序列 共同的进化起源,序列比对是序列分析（结构、功能与进化关系）的基础,序列比对相关概念,相似性(Similarity)：序列间相同碱基或氨基酸残基所占比例的高低，是直观的数量关系，是序列间匹配程度的直接测度

13、。 同源性(Homology)：核酸或蛋白质序列间具有共同起源，是依据进化事实的因果判断。 直系同源(orthology)：经由物种分离事件形成的存在于不同物种中的结构及功能相似的序列 旁系同源(paralogy)：经由基因复制而演化形成在同一生物体的一系列结构及功能相似的序列 序列的相似与序列的同源之间无必然的因果关系（相似并不一定同源），但存在一定的相关性（足够的相似性往往意味着同源关系的存在）,成对序列比对,用于两个核酸（蛋白质）序列间相似性比较，是生物信息学的基本技术，是多重比对的基础。 The process of lining up two sequences to achieve

14、 maximal levels of identity (or conservation, in the case of amino acid sequences) for the purpose of assessing the degree of similarity and the possibility of homology.,DNA、蛋白质序列比对,Query: 181 catcaactacaactccaaagacacccttacacccactaggatatcaacaaacctacccac 240 | | | | | | Sbjct: 189 catcaactgcaaccccaaa

15、gccacccct-cacccactaggatatcaacaaacctacccac 247,RBP: 26 RVKENFDKARFSGTWYAMAKKDPEGLFLQDNIVA 59 + K+ + + GTW+MA + L + A glycodelin: 23 QTKQDLELPKLAGTWHSMAMA-TNNISLMATLKA 55,成对序列比对,全局比对（global alignment） 比较结果包含所比较序列全长范围内所有位点的比对，适用于整体相似水平高的同源序列，在分子系统学中常用。 局部比对（local alignment） 仅对相似水平较高的 局部片段进行比对，多用于分子结构和功

16、能域研究。,成对序列比对举例,视黄醇结合蛋白(RBP)与水解牛乳蛋白(-lactoglobulin)的比对结果,完全相同区,部分相似区,identity,Similarity,Gap,间隔,成对序列比对,完全匹配（complete match） 部分匹配（partial match） 对于蛋白质序列而言，不同但性质(size, charge, hydrophobicity, and polarity)相近的氨基酸常具有相似的功能。 碱性aa、酸性aa、中性-非极性aa、极性aa 空位（gap） 源于序列片段的插入或缺失(insertion or deletion) 序列比对不同算法对于空位的处

17、理方式有所不同，这直接影响了算法的适用性,局域比对搜索工具 BLAST,Basic Local Alignment Search Tool 基于BLAST算法（Altschul et al, 1990）的序列比对搜索工具，由NCBI研发并维护，因其快速可靠、功能全面且使用方便而被广泛运用。 现行的BLAST算法为BLAST2。,BLAST 应用示例,大肠杆菌乳糖操纵子(lactose operon)翻译的一个蛋白质序列：,MKPVTLYDVAEYAGVSYQTVSRVVNQASHVSAKTREKVEAAMAELNYIPNRVAQQLAGKQSLLIGVATSSLALHAPSQIVAAIKSRA

18、DQLGASVVVSMVERSGVEACKAAVHNLLAQRVSGLIINYPLDDQDAIAVEAACTNVPALFLDVSDQTPINSIIFSHEDGTRLGVEHLVALGHQQIALLAGPLSSVSARLRLAGWHKYLTRNQIQPIAEREGDWSAMSGFQQTMQMLNEGIVPTAMLVANDQMALGAMRAITESGLRVGADISVVGYDDTEDSSCYIPPSTTIKQDFRLLGQTSVDRLLQLSQGQAVKGNQLLPVSLVKRKTTLAPNTQTASPRALADSLMQLARQVSRLESGQ,用NCBI的BLAST工具对该序列进行数据库检索和序

19、列比对分析,/blast/Blast.cgi,BLAST 应用示例,BLAST 应用示例,NCBI翻译库(nr) Refseq专家库 SwissProt欧洲专家库 Pat专利库 Pdb三维结构库,BLAST 应用示例,BLAST参数设置（蛋白质对位）,过滤器去除待测序列中重复序列干扰,相似程度降低,期望阈值,E值= x| H0) 其服从极值分布函数 E-value P * N (设数据库中有N条库存序列，Px时由随机偶然因素产生的对位结果的次数（E=1出现一次错误对位）。 根据需要设定E-value的阈值，默认值为10，更严格的搜索需设置更小

20、的阈值，如1。,BLAST 应用示例,相似度最高的检索序列(Sbjct)与被测序列(Query)的比对：,结论：待测序列应是大肠杆菌乳糖操纵子阻遏蛋白，与库存蛋白有一个错配(LS),BLAST 应用示例,相似度第二高的检索序列与被测序列的比对：,42%完全匹配，60%相同或性质相近(+)，无间隔。两序列间可能存在很高的同源性,克氏肺炎杆菌 Klebsiella pneumoniae 乳糖操纵子阻遏蛋白,BLAST 应用示例,使用过滤器的必要性,agaggccagagagggtgtcggatcccctagacctgtagtctcggacggttgttagctgacgtgtaggtgctgggag

21、ataaacccggacgctctggaggaccagtcagtcagctcttagcccctgagccctctctcctcaattgtatctttcaacttaattttctactaaaatcttttaatactttgaaaaaaaaaaaaaaaaaaaaaaaaaa,用BlastN工具进行数据库检索和序列比对分析，若不用过滤器，则：,PolyA尾巴,可能出现大量无关的假命中,BLAST 应用示例,使用过滤器的必要性,用BlastN工具进行数据库检索和序列比对分析，若用过滤器，则：,用N替代a,过滤掉简单重复序列(low complexity),多重序列比对与Clustal工具,多重序列比

22、对 (Multiple sequence alignment) 就是把两条以上可能有系统进化关系的序列同时进行比对的方法。,如果序列总体相似程度低，则成对比对往往不能发现一些局部匹配。而通过多重比对，能大大增加发现这些相似结构域的几率，而这些结构通常具有重要的生物学意义。,多重比对工具,运用最多的多重序列比对方法 Clustal算法 By Feng “.” 号： 相似性略低的残基位点（分子大小或极性基本一致）,红色: AVFPMILW (小分子非极性氨基酸) 兰色: DE (酸性、负电氨基酸) 紫色: RHK (碱性正点氨基酸) 绿色: STYHCNGQ (带羟基、氨基的极性氨基酸),Clus

23、talX,单机版Clustal多重比对工具。 版本2.0 输入序列多 用Fasta格式,ClustalX应用示例,步骤： 载入FASTA序列 “Do Complete Alignment” 指定输出文件，包括： “指导树”文件(后缀dnd) “对位输出”文件(后缀aln),比对结果： 用彩条和符号标出相似程度高的位点 对位结果可通过调整参数进行调整，以寻找最佳对位，帮助挖掘其中的生物学意义,序列比对的应用,两两比对 在序列库中对某目标序列进行比较认定 用EST或蛋白序列与核酸序列的比对，确定编码区，用于基因预测 多重比对 通过比对寻找、分析保守区域（conserved region）： 发现新

24、的基序(motif) 用于基因家族或蛋白质家族分析 建立分子系统进化树 通过与已知基因的比对，对未知基因结构进行预测 通过与结构已知的蛋白质的比对，对未知蛋白的三维结构进行预测,多重序列比对寻找保守的序列基序,球蛋白(globin)的多重比较,多重序列比对寻找保守的序列基序,多重比对发现调控序列的保守结构,如：转录启动子的保守区域-10， -35 sequence,用BLAST 进行蛋白质家族分析和寻找保守域,用BLAST 进行蛋白质家族分析和寻找保守域,Pfam (Protein Domain Families database) 包括每个蛋白质家族的多重比对信息,用BLAST 进行蛋白质家族分析和寻找保守域,结果解读： 该蛋白序列有一个保守域，属离子转运超家族(Ion_trans_2 superfamily), 位置在60-110aa之间,ClustalX应用示例,Ning Zhe