真核生物基因結(jié)構(gòu)的預(yù)測分析方法軟件

實習(xí)一基因組數(shù)據(jù)注釋和功能分析實習(xí)二真核生物基因結(jié)構(gòu)的預(yù)測分析實習(xí)三芯片的基本數(shù)據(jù)處理和分析實習(xí)四蛋白質(zhì)結(jié)構(gòu)與功能分析實習(xí)五蛋白質(zhì)組學(xué)數(shù)據(jù)分析實習(xí)六系統(tǒng)生物學(xué)軟件實習(xí)課程內(nèi)容基因組學(xué)轉(zhuǎn)錄物組學(xué)蛋白質(zhì)組學(xué)系統(tǒng)生物學(xué)1目前一頁\總數(shù)四十七頁\編于二十二點基因組序列cDNA序列編碼區(qū)預(yù)測CodonbiasGCContent限制性酶切位點基因結(jié)構(gòu)分析選擇性剪切轉(zhuǎn)錄調(diào)控因子序列比對功能注釋KEGGGO系統(tǒng)發(fā)育樹蛋白質(zhì)序列翻譯蛋白質(zhì)理化性質(zhì)二級結(jié)構(gòu)預(yù)測結(jié)構(gòu)域分析重要信號位點分析三級結(jié)構(gòu)預(yù)測基因組功能分析2目前二頁\總數(shù)四十七頁\編于二十二點真核生物基因的主要結(jié)構(gòu)3目前三頁\總數(shù)四十七頁\編于二十二點基因結(jié)構(gòu)分析開放讀碼框GENSCANGENOMESCANCpG島CpGPlot轉(zhuǎn)錄終止信號POLYAH啟動子/轉(zhuǎn)錄起始位點PromoterScan密碼子偏好分析CodonWmRNA剪切位點NETGENE2Spidey選擇性剪切ASTD基因結(jié)構(gòu)分析常用軟件4目前四頁\總數(shù)四十七頁\編于二十二點開放讀碼框的識別開放讀碼框（openreadingframe,ORF)

是一段起始密碼子和終止密碼子之間的堿基序列ORF是潛在的蛋白質(zhì)編碼區(qū)5目前五頁\總數(shù)四十七頁\編于二十二點基因開放閱讀框/基因結(jié)構(gòu)分析識別工具ORFFinder/gorf/gorf.htmlNCBI通用BestORF/berry.phtml?topic=bestorf&group=programs&subgroup=gfindSoftberry真核GENSCAN/GENSCAN.htmlMIT脊椎、擬南芥、玉米GeneFinder/tools/genefinder/Zhanglab人、小鼠、擬南芥、酵母FGENESH/berry.phtml?topic=fgenesh&group=programs&subgroup=gfindSoftberry真核(基因結(jié)構(gòu))GeneMark/GeneMark/eukhmm.cgiGIT原核GLIMMER/genomes/MICROBES/glimmer_3.cgi/software/glimmer

Maryland原核Fgenes/berry.phtml?topic=fgenes&group=programs&subgroup=gfindSoftberry人（基因結(jié)構(gòu)）FgeneSV/berry.phtml?topic=virus&group=programs&subgroup=gfindvSoftberry病毒Generation/generation/ORNL原核FGENESB/berry.phtml?topic=fgenesb&group=programs&subgroup=gfindbSoftberry細(xì)菌（基因結(jié)構(gòu)）GenomeScan

/genomescan.html

MIT脊椎、擬南芥、玉米GeneWise2http://www.ebi.ac.uk/Wise2/EBI人GRAIL/grailexp/ORNL人、小鼠、擬南芥、果蠅6目前六頁\總數(shù)四十七頁\編于二十二點ORF識別：GENSCAN結(jié)果返回到郵箱（可選）提交序列提交序列文件運行GENSCAN顯示氨基酸或CDS序列序列名稱（可選）是否顯示非最優(yōu)外顯子選擇物種類型7目前七頁\總數(shù)四十七頁\編于二十二點8GENSCAN輸出結(jié)果：文本8目前八頁\總數(shù)四十七頁\編于二十二點9GENSCAN輸出結(jié)果：圖形9目前九頁\總數(shù)四十七頁\編于二十二點ORF識別：

GenomeScan提交待分析序列提交同源蛋白質(zhì)序列運行GenomeScan10目前十頁\總數(shù)四十七頁\編于二十二點GenomeScan輸出結(jié)果：文本預(yù)測外顯子位置、可信度等信息同源比對信息預(yù)測結(jié)果的氨基酸序列11目前十一頁\總數(shù)四十七頁\編于二十二點GenomeScan輸出結(jié)果：圖形12目前十二頁\總數(shù)四十七頁\編于二十二點課堂練習(xí)1使用GENESCAN預(yù)測序列中可能的ORF。2使用GENOMESCAN預(yù)測序列中可能的ORF。練習(xí)用的序列文件在c:\zcni\shixi2文件下，名字為clone.fasta，使用寫字板打開查看。13目前十三頁\總數(shù)四十七頁\編于二十二點轉(zhuǎn)錄調(diào)控序列分析

CpG島、轉(zhuǎn)錄終止信號和啟動子區(qū)域的預(yù)測14目前十四頁\總數(shù)四十七頁\編于二十二點CpG島的預(yù)測CpG島常位于真核生物基因轉(zhuǎn)錄起始位點，GC含>50%，長度>200bp的一段DNA序列。15目前十五頁\總數(shù)四十七頁\編于二十二點CpGIsland分析常用軟件CpGIsland/cpgislands2/cpg.aspxWebCpGPlothttp://www.ebi.ac.uk/emboss/cpgplot/index.htmlWebCpGfinder/berry.phtml?topic=cpgfinder&group=programs&subgroup=promoterWebCpGi130/CpG130.dowebCpGproDhttp://pbil.univ-lyon1.fr/software/cpgprod_query.htmlweb16目前十六頁\總數(shù)四十七頁\編于二十二點提交序列文件提交序列參數(shù)選項CpG島的預(yù)測：CpGPlot目前十七頁\總數(shù)四十七頁\編于二十二點GENESCAN預(yù)測結(jié)果

起始為532bp

終止于51783bp18目前十八頁\總數(shù)四十七頁\編于二十二點轉(zhuǎn)錄終止信號上游作用元件：AAUAAA下游作用元件：GCrich二重對稱區(qū)、UUUUUUC-GC-GG-CG-CU-AG-CG-CC-GG-CUUUUUUUUURNA5’3’AAUAAACAAAAAAAAAAAAA成熟mRNA5’3’AAUAAACAGUmRNA前體5’3’19目前十九頁\總數(shù)四十七頁\編于二十二點轉(zhuǎn)錄終止信號預(yù)測：POLYAH

提交序列文件提交序列20目前二十頁\總數(shù)四十七頁\編于二十二點polyA位置GENESCAN預(yù)測結(jié)果

PolyA位點52490bpPOLYAH輸出結(jié)果21目前二十一頁\總數(shù)四十七頁\編于二十二點啟動子區(qū)結(jié)構(gòu)啟動子（Promoter）位于結(jié)構(gòu)基因5’端上游，能活化RNA聚合酶，使之與模板DNA結(jié)合并具有轉(zhuǎn)錄起始的特異性。轉(zhuǎn)錄起始位點(Transcriptionstartsite,TSS)PYCAPY（嘧啶）核心啟動子元件(Corepromoterelement)TATAbox，Pribnowbox

（TATAA)上游啟動子元件(Upstreampromoterelement，UPE)CAATbox，GCbox，SP1，Otc增強子(Enhancer)22目前二十二頁\總數(shù)四十七頁\編于二十二點原核和真核生物基因轉(zhuǎn)錄起始位點上游區(qū)結(jié)構(gòu)原核生物真核生物TTGACATATAATAmRNA＋1－10－35PyAPyTATAATGC區(qū)CAAT區(qū)mRNA＋1－40－25－110增強子上游啟動子元件，UPE核心啟動子元件轉(zhuǎn)錄起始位點23目前二十三頁\總數(shù)四十七頁\編于二十二點PromoterScan:80/molbio/proscan/WebPromoser/zlab/PromoSer/WebNeuralNetworkPromoterPrediction/seq_tools/promoter.htmlWebSoftberry:BPROM,TSSP,TSSG,TSSW/berry.phtml?topic=index&group=programs&subgroup=promoterWebMatInspectorhttp://www.gene-regulation.de/WebRSAThttp://rsat.ulb.ac.be/rsat/WebCister/~mfrith/cister.shtmlWeb啟動子結(jié)合位點分析常用軟件24目前二十四頁\總數(shù)四十七頁\編于二十二點啟動子預(yù)測：PromoterScan提交序列25目前二十五頁\總數(shù)四十七頁\編于二十二點PromoterScan輸出結(jié)果找到的TATAbox和轉(zhuǎn)錄起始位點預(yù)測可能的轉(zhuǎn)錄因子轉(zhuǎn)錄因子在提交序列中的位置26目前二十六頁\總數(shù)四十七頁\編于二十二點課堂練習(xí)1使用CpG

Plot預(yù)測基因的CpG

island位置。2使用PolyAH預(yù)測基因可能的轉(zhuǎn)錄終止的位置。3使用PromotorScan尋找基因上游序列里可能的轉(zhuǎn)錄因子調(diào)控區(qū)域。27目前二十七頁\總數(shù)四十七頁\編于二十二點基因密碼子偏好性1.研究蛋白質(zhì)結(jié)構(gòu)功能中的作用2.在表達(dá)外源基因方面的作用3.在生物信息學(xué)研究中的作用28目前二十八頁\總數(shù)四十七頁\編于二十二點基因密碼子偏好性:CodonW粘帖目的序列密碼子表的選擇如需計算FOP/CBI選擇相應(yīng)物種如需計算CAI選擇相應(yīng)物種輸出格式(默認(rèn)不選)匯總所有基因的信息29目前二十九頁\總數(shù)四十七頁\編于二十二點參數(shù)選擇計算所有指數(shù)選擇導(dǎo)入對應(yīng)物種CAI

FOP

CBI數(shù)據(jù)計算有效密碼子數(shù)計算GC含量計算GC3s含量計算同義密碼子數(shù)量計算同義密碼子第三位堿基組成密碼子總數(shù)30目前三十頁\總數(shù)四十七頁\編于二十二點各項指數(shù)輸出結(jié)果密碼子使用頻率CodonW結(jié)果界面31目前三十一頁\總數(shù)四十七頁\編于二十二點課堂練習(xí)使用CodonW分析基因的密碼子使用偏好，了解密碼子偏好分析中各指數(shù)的含義。32目前三十二頁\總數(shù)四十七頁\編于二十二點內(nèi)含子/外顯子剪切位點識別如何分析核酸序列中的外顯子組成？通過對特征序列(GT-AG)的分析進(jìn)行直接的預(yù)測基因預(yù)測軟件（NetGene2）與相應(yīng)的基因組序列比對，分析比對片段的分布位置（Spidey）33目前三十三頁\總數(shù)四十七頁\編于二十二點34目前三十四頁\總數(shù)四十七頁\編于二十二點剪切位點識別：NetGene2http:///提交序列選擇物種35目前三十五頁\總數(shù)四十七頁\編于二十二點NetGene2輸出結(jié)果供體位點受體位點可信度

相位36目前三十六頁\總數(shù)四十七頁\編于二十二點mRNA剪切位點識別：SpideyNCBI開發(fā)的在線預(yù)測程序用于mRNA序列同基因組序列比對分析37目前三十七頁\總數(shù)四十七頁\編于二十二點Spidey同源序列的獲得:序列比對通過BLAST進(jìn)行序列比對，找到可能同源的相似性好的一系列mRNA序列。BLAST比對到的三條mRNA序列38目前三十八頁\總數(shù)四十七頁\編于二十二點輸入基因組序列或序列數(shù)據(jù)庫號輸入相似性序列判斷用于分析的序列間的差異，并調(diào)整比對參數(shù)不受默認(rèn)內(nèi)含子長度限制，默認(rèn)長度：內(nèi)部內(nèi)含子為35kb,末端內(nèi)含子為100kb比對閾值選擇物種輸出格式選擇39目前三十九頁\總數(shù)四十七頁\編于二十二點Spidey輸出結(jié)果第一條藍(lán)色序列為基因組序列，橘黃色為外顯子外顯子對應(yīng)于基因組上的起始/結(jié)束位置外顯子對應(yīng)于mRNA/cDNA上的起始/結(jié)束位置供體、受體位點外顯子長度一致性百分比錯配和gap外顯子序號序列聯(lián)配結(jié)果40目前四十頁\總數(shù)四十七頁\編于二十二點課堂練習(xí)1練習(xí)兩種預(yù)測剪切位點的軟件的使用，NetGene2和Spidey。2

Spidey的同源序列文件保存在c:\zcni\shixi2文件下，名字為Spidey.txt，使用寫字板打開查看。41目前四十一頁\總數(shù)四十七頁\編于二十二點選擇性剪切（Alternativesplicing）分析選擇性剪接是調(diào)控基因