实验三基因组序列分析

实验项目三：基因序列分析

一、实验目的和要求：掌握基因可读框的识别；掌握启动子区域的预测掌握CpG岛的预测掌握转录终止信号的预测采用mRNA序列预测基因：Spidey的使用掌握各预测服务器结果的分析1本文档共46页；当前第1页；编辑于星期日\13点3分原核生物基因结构1长开放阅读框2高基因密度3简单的基因结构4基因组中GC含量变化非常大特点：2本文档共46页；当前第2页；编辑于星期日\13点3分真核生物基因结构特点：1基因结构复杂2具有复杂的基因转录调控方式3具有丰富的可变剪接4有明显的CpG岛、密码子使用具有偏好性3本文档共46页；当前第3页；编辑于星期日\13点3分基因组序列分析4本文档共46页；当前第4页；编辑于星期日\13点3分

例：WhatisGenePrediction?

GivenanuncharacterizedDNAsequence,findout:

1.Wheredoesthegenestartsandends?

2.Whichregionscodeforaprotein?

AGCATCGAAGTTGCATGACGATGCATGACCTAGCAGCATCGAAGTTGCATGACGATGCATGACCTAGCAAGTTGCATGACGATGCATGACCTAGCAGCATCGAAGTTGCATGACGATGCATGACCTAGTGCATGACGATGCATGACCTAGCAGCATCGAAGTTGCATGACGATGCATGACCTAGCAAGTTGCATGACGATTGACCTAGTGCATGACGATGCATGACCTAGCAGCATCGAAGTTGCATGACGATGCATGACCTAGCAAGAAGTTGCATGACGATGCATGACCTAGTGCATGACGATGCATGACCTAGCAGCATCGAAGTTGCATGACGATGCATGACCTAGCAAGTTGCATGACGATTGACCTAGTGCATGACGATGCATGACCTAGCAGCATCGCGATGCATGACCTAGCAAGAAGTTGCATGACGATGCATGACCTAGTGCATGACGATGCATGACCTAGCAGCATCGAAGTTGCATGACGATGCATGACCTAGCAAGTTGCATGACGATTGACCTAGTGCATGACTGACCTAGCAGCATCGAAGTTGCATGACGATGCATGACCTAGTGCATGACGATGCATGACCTAGCAGCATCGAAGTTGCATGACGATGCATGACCTAGCAAGTTGCATGACGATTGACCTAGTGCATGACGATGCATGACCTAGCAGCATCGAAGTTGCATGACGATGCATGACCTAGCAAGAAGTTGCATGACGATGCATGACCTAATGCgene1gene2gene3exonintergenicregionintron5本文档共46页；当前第5页；编辑于星期日\13点3分6本文档共46页；当前第6页；编辑于星期日\13点3分一开放读码框的识别开放读码框（openreadingframe,ORF)

是一段起始密码子和终止密码子之间的碱基序列ORF是潜在的蛋白质编码区基因预测7本文档共46页；当前第7页；编辑于星期日\13点3分基因开放阅读框/基因结构分析识别工具Getorfhttp://bioweb.pasteur.fr/seqanal/interfaces/getorf.htmlEMBOSS通用Plotorfhttp://bioweb.pasteur.fr/seqanal/interfaces/plotorf.htmlEMBOSS通用ORFFinder/gorf/gorf.htmlNCBI通用BestORF/all.htmSoftberry真核GENSCAN/GENSCAN.htmlMIT脊椎、拟南芥、玉米GeneFinder/tools/genefinder/Zhanglab人、小鼠、拟南芥、酵母FGENESH/all.htmSoftberry真核GeneMark/GeneMark/GIT原核GLIMMER/genomes/MICROBES/glimmer_3.cgi/software/glimmer

Maryland原核FgeneSB/all.htmSoftberry细菌FgeneSV/all.htmSoftberry病毒Generation/generation/ORNL原核FGENESH+/all.htmSoftberry原核GenomeScan

/genomescan.html

MIT脊椎、拟南芥、玉米GeneWise

http://www.ebi.ac.uk/Wise2/EBI人、蠕虫GRAIL/grailexp/ORNL人、小鼠、拟南芥、果蝇8本文档共46页；当前第8页；编辑于星期日\13点3分1.ORFFinder的使用及结果分析9本文档共46页；当前第9页；编辑于星期日\13点3分1.ORFFinder的使用及结果分析10本文档共46页；当前第10页；编辑于星期日\13点3分1.ORFFinder的使用及结果分析11本文档共46页；当前第11页；编辑于星期日\13点3分1.ORFFinder的使用及结果分析12本文档共46页；当前第12页；编辑于星期日\13点3分1.ORFFinder的使用及结果分析13本文档共46页；当前第13页；编辑于星期日\13点3分1.ORFFinder的使用及结果分析Blast比对结果搜索到多个显著相似的序列，故所预测的ORF的可信度较高。如果要获取该ORF所编码的蛋白质序列，可以点击“Accept”按钮后，在“1GenBank”的下拉框中选择“3Fasta”，并点击“view”，即可获取该ORF所编码的蛋白质序列。14本文档共46页；当前第14页；编辑于星期日\13点3分1.ORFFinder的使用及结果分析15本文档共46页；当前第15页；编辑于星期日\13点3分1.ORFFinder的使用及结果分析16本文档共46页；当前第16页；编辑于星期日\13点3分1.ORFFinder的使用及结果分析17本文档共46页；当前第17页；编辑于星期日\13点3分1.ORFFinder的使用及结果分析18本文档共46页；当前第18页；编辑于星期日\13点3分提交序列提交序列文件运行GENSCAN选择物种显示氨基酸或CDS序列序列名称（可选）是否显示非最优外显子2.Genscan的使用及结果分析19本文档共46页；当前第19页；编辑于星期日\13点3分基因、外显子及类型正链、负链预测单元起始、终止及长度相位编码区打分值可信概率、得分值2.Genscan的结果分析20本文档共46页；当前第20页；编辑于星期日\13点3分3.FGENESH的使用及结果分析输入序列的Fasta文件21本文档共46页；当前第21页；编辑于星期日\13点3分3.FGENESH的使用及结果分析起始外显子中间及末端外显子PolyA位点起始碱基终止碱基打分长度22本文档共46页；当前第22页；编辑于星期日\13点3分3.FGENESH的使用及结果分析23本文档共46页；当前第23页；编辑于星期日\13点3分3.FGENESH的使用及结果分析24本文档共46页；当前第24页；编辑于星期日\13点3分二.原核和真核生物基因转录起始位点上游区结构

原核生物真核生物TTGACATATAATAmRNA＋1－10－35PyAPyTATAATGC区CAAT区mRNA＋1－40－25－110增强子上游启动子元件，UPE核心启动子元件转录起始位点25本文档共46页；当前第25页；编辑于星期日\13点3分原核生物真核生物26本文档共46页；当前第26页；编辑于星期日\13点3分二.启动子预测输入序列的Fasta文件27本文档共46页；当前第27页；编辑于星期日\13点3分启动子预测结果从预测结果可知，预测的启动子区在32564至32783之间，启动子阈值系统默认为53.00，预测的启动子分值为84.69，高于阈值，分值越高，说明预测的准确性大。与该启动子可能结合的转录因子如下所示28本文档共46页；当前第28页；编辑于星期日\13点3分三CpG岛预测CpG岛CpG岛又称为HTF岛，是DNA上的一个区域，此区域富含GC，二者以磷酸酯键相连。位于真核生物基因转录起始位点上游，GC含>50%，长度>200bpCpG岛常出现在管家基因或频繁表达的基因的启动子附近，在这些部位，CpG岛具有阻止序列甲基化的作用，因此，搜索CpG岛可以为基因及其启动子的预测提供线索。CpGIsland分析CpGIsland/cpgislands2/cpg.aspxWebCpGfinder/berry.phtml?topic=cpgfinder&group=programs&subgroup=promoterWebCpGPlot/CpGReport/Isochorehttp://www.ebi.ac.uk/emboss/cpgplot/index.htmlWeb29本文档共46页；当前第29页；编辑于星期日\13点3分输入序列的Fasta文件30本文档共46页；当前第30页；编辑于星期日\13点3分从该序列的预测结果来看，找到两个CpG岛，分别位于501-727，长度为227个碱基，54380-54691，长度为31231本文档共46页；当前第31页；编辑于星期日\13点3分四转录终止信号加polyA信号：AAUAAA转录终止信号：GCrich二重对称区、UUUUUUC-GC-GG-CG-CU-AG-CG-CC-GG-CUUUUUUUUURNA5’3’AAUAAACAAAAAAAAAAAAA成熟mRNA5’3’AAUAAACAGUmRNA前体5’3’32本文档共46页；当前第32页；编辑于星期日\13点3分33本文档共46页；当前第33页；编辑于星期日\13点3分转录终止信号预测Hcpolyar.it/~webgene/wwwHC_polya.htmlWebPOLYAH/berry.phtml?topic=polyah&group=programs&subgroup=promoterWebpolyadq/tools/polyadq/polyadq_form.htmlWeb34本文档共46页；当前第34页；编辑于星期日\13点3分POLYAH的使用及结果分析输入序列的Fasta文件35本文档共46页；当前第35页；编辑于星期日\13点3分POLYAH的使用及结果分析预测的POLYA位点，LDF为权重36本文档共46页；当前第36页；编辑于星期日\13点3分内含子/外显子剪切位点识别对基因组序列的读码框区域进行预测内含子5’端供体位点(donorsplicesite):GT内含子3’端受体位点(acceptorsplicesite):AG预测工具：GENSCAN，GENEMARKNetGene2,SpliceView37本文档共46页；当前第37页；编辑于星期日\13点3分38本文档共46页；当前第38页；编辑于星期日\13点3分mRNA剪切位点识别：spideyNCBI开发的在线预测程序用于mRNA序列同基因组序列比对分析39本文档共46页；当前第39页；编辑于星期日\13点3分40本文档共46页；当前第40页；编辑于星期日\13点3分序列在线提交形式：界面中有两个窗口：上方窗口用于输入基因组序列（直接粘贴序列或用GenbankID/AC号）下方窗口用于输入cDNA/mRNA序列（直接粘贴序列或用GenbankID/AC号）可同时输入多条cDNA/mRNA序列与同一条基因组序列进行分析Spidey序列提交页面输入基因组序列或序列数据库号AC002390.141