干貨基因注釋方法

干貨基因注釋方法現(xiàn)如今隨著測序技術(shù)的發(fā)展，組裝一個完整的基因組也是越來越普遍，價格已經(jīng)不像前些年那樣昂貴（當(dāng)然超大基因組除外），待基因組組裝完成后，基因預(yù)測將會是接下來需要進(jìn)行的一個重要工作（雖然能編碼的基因相對于整個基因組只是占了很小一部分），但通過這些基因結(jié)合一些近緣物種，還是能挖掘一些比較有意思東西，如物種進(jìn)化等。那么接下來小編就給簡單介紹一些基因預(yù)測方面方法。目前，基因組預(yù)測策略大致可以分為三種：Abinitio、Homology-based和EST/Unigene。從頭預(yù)測Abinitio主要通過探索DNA序列中特異的區(qū)域，如基因的起始區(qū)域和終止區(qū)域，來進(jìn)行基因預(yù)測。目前常用的軟件有Augustus.GlimmerHMM、SNAP、GenelD、GenScan.Brak等。Augustus運用隱馬爾科夫模型，模型在DNA序列和基因結(jié)構(gòu)上定義一個概率分布，采用維特比的算法，它自身帶了一個訓(xùn)練集，如人、斑馬魚等。在進(jìn)行預(yù)測是可以選擇自帶的訓(xùn)練集，也可以用挑選轉(zhuǎn)錄組和同源預(yù)測最優(yōu)結(jié)果給它生成一個訓(xùn)練集。這里順帶說下Braker軟件，它是基于genemaker預(yù)測結(jié)果作為訓(xùn)練集，通常小編Augustus和Braker會二選一。GlimmerHMM是把一個基因看做幾種特征序列，這些特征序列包括內(nèi)含子、基因間區(qū)和四種夕卜顯子（初始、中間、最終和單一外顯子）之后進(jìn)行有序切換形成馬爾科夫鏈。示意圖如下：InitExon+ TemiExon+—*ExonSngl+ GlimmerHMM使用的模型基于以下幾個假設(shè)：-假設(shè)每個基因都開始于起始密碼子ATG-假設(shè)每個基因閱讀框內(nèi)除最后一個密碼子外沒有終止密碼子(noin-framestopcodons)。-每個外顯子與前一個外顯子在同一個閱讀框中。(翻譯閱讀時外顯子間沒有移框).它也是需要一個訓(xùn)練集，通常也是自己生成一個訓(xùn)練集的效果會略優(yōu)于已有的一些。(/software/glimmerhmm/man.shtml)SNAP通過隱馬爾科夫模型進(jìn)行預(yù)測，也是需要一個訓(xùn)練集。以上這些軟件都可以自身構(gòu)建一個訓(xùn)練集，這里小編覺得畢竟還是用自己的東西舒服，也就是自身訓(xùn)練結(jié)果要稍微優(yōu)于其他模式生物訓(xùn)練集。這里小編在做真菌時，從頭軟件一般會選取這三個，GenScan和GenelD就放棄掉了，動植物基因組通常就是多多益善嗎，能用上就都給用上。GenScan也是一款比較經(jīng)典軟件，通常在預(yù)測真核生物（人）還是有不錯的效果。GeneID可以算是元老級，第一代的基因識別軟件，這個準(zhǔn)確率不高，通常在整合是權(quán)重也不會給太高。同源預(yù)測同源預(yù)測軟件通常利用GeneWise和GeneMoMa，前者是需要同源物種的蛋白序列，后者需要同源物種基因組序列及對應(yīng)的GFF文件，目前小編已經(jīng)拋棄GeneWise，使用最多的就是GeneMoMa，但是讓小編十分頭疼的是在準(zhǔn)備GFF文件太花費精力，這個軟件真的是挑肥揀瘦，必須滿足其格式才能可以運行，目前從NCBI的Reseq和Ensemble上下載都可以，其他地方來的那就得還點時間寫個腳本改下了。轉(zhuǎn)錄組數(shù)據(jù)預(yù)測PASA軟件是基于Unigene/EST序列進(jìn)行預(yù)測軟件，這個可能就需要拿到一個混樣轉(zhuǎn)錄組數(shù)據(jù)首先進(jìn)行無參組裝，接下來根據(jù)Unigene組裝結(jié)果在進(jìn)行比對，通常用Gmap或Blat兩種方法，最好三代全長轉(zhuǎn)錄本和二代一起來進(jìn)行預(yù)測，這樣可以使得找到的結(jié)構(gòu)更為準(zhǔn)確、可靠，此外PASA還有另外的一個功能就是可以用其預(yù)測可變剪切，俗稱PASA修飾。最終結(jié)果整合這么多軟件跑出來的結(jié)果，有的可靠性高些，比如轉(zhuǎn)錄組和同源；有一些要稍微差一些如GeneID，那么就需要一個軟件將這些結(jié)果進(jìn)行一個整合，通俗些就是大家放到一起比較下，看下各個軟件預(yù)測結(jié)果分布情況，本著以少數(shù)服從多數(shù)原則（這里只是簡單比喻下莫要當(dāng)真），根據(jù)權(quán)重打分，使用EVM軟件得到一版最終結(jié)果，目前小編用到最多的就是EVM真菌、植物或動物統(tǒng)統(tǒng)可以搞定，用過一段時間Glean,感覺在整合超大基因時，容易成多個（或許是參數(shù)沒有調(diào)整合理）?？偨Y(jié)上面就是小編在進(jìn)行基因預(yù)測時的一些軟件使用心得，還有是再做一些研究比較多的物種比如水稻等，那同源權(quán)重一定要調(diào)高