新城疫病毒生物信息分析系統(tǒng)的構建及其全基因組的比較研究.pdf

1、隨著后基因組時代的到來，生物信息學的主要任務已經(jīng)從生物數(shù)據(jù)的積累發(fā)展到對生物數(shù)據(jù)的整合與處理階段。由于生物信息學分析資源多樣、生物數(shù)據(jù)海量、分析任務復雜，構建針對不同需求的生物信息分析系統(tǒng)顯得非常重要。本研究根據(jù)禽類重要傳染病病原——新城疫病毒(Newcastle disease virus，NDV)基因組分析的需求，開發(fā)針對NDV的生物信息分析系統(tǒng)，并利用生物信息學分析方法，開展NDV全基因組的比較研究，在全基因組水平上闡明NDV遺傳

2、進化規(guī)律，為防控新城疫(Newcastle disease，ND)的流行提供理論依據(jù)。
　　1.基于Z曲線理論的核酸序列可視化分析系統(tǒng)的構建。Z曲線是我國科學家發(fā)展起來的一種用幾何學方法分析比較基因組核酸序列的有效工具，它將基因組字母序列記錄的信息轉換成三維坐標的空間曲線，使研究者可以直觀地對研究對象進行分析比較。本研究根據(jù)Z曲線理論，在對系統(tǒng)需求進行分析的基礎上，采用計算機圖形學中的包圍盒抽稀算法思想，使用Visual C++編

3、程，開發(fā)了具有DNA序列的座標轉換、Z曲線顯示和比較、基本輸入輸出等功能的病毒基因組生物信息分析系統(tǒng)。通過以不同基因型NDV基因組核酸序列的顯示為例，證明該系統(tǒng)能夠被用于對病毒基因組序列進行顯示和比較。
　　2.基于Web的新城疫病毒基因分型系統(tǒng)的構建。人們根據(jù)新城疫病毒的分子生物學特性，建立了基于新城疫病毒F基因限制性酶切位點和F蛋白特定氨基酸位點的基因分型方法，由于分型過程中序列處理過程復雜、結果判定容易產(chǎn)生偏差，本研究根據(jù)N

4、DV基因型分析的規(guī)律，以瀏覽器為用戶界面，采用JSP技術，后臺數(shù)據(jù)庫采用SQL Server，以開源軟件eclipse3.2為開發(fā)平臺，構建了基于Web的新城疫病毒基因自動分型系統(tǒng)，通過實例檢驗，系統(tǒng)實現(xiàn)了快速準確地判別NDV毒株基因型的功能。
　　3.基于Web的新城疫病毒生物信息分析系統(tǒng)的構建。針對新城疫病毒基因組生物信息分析的需要，以IBM服務器作為硬件平臺，以Microsoft Windows Service2003作為操

5、作系統(tǒng)、Microsoft SQL Service2005為數(shù)據(jù)庫開發(fā)環(huán)境，采用VB+COM組件編程開發(fā)了后臺NDV生物信息數(shù)據(jù)自動獲取管理系統(tǒng)，采用ASP、JSP+COM組件編程開發(fā)了前臺基于WEB的生物信息分析系統(tǒng)。設計與實現(xiàn)過程包括系統(tǒng)多層結構模型的設計、后臺數(shù)據(jù)更新系統(tǒng)的設計與實現(xiàn)、前臺數(shù)據(jù)分析系統(tǒng)的設計與實現(xiàn)。系統(tǒng)實現(xiàn)了對NDV核酸和蛋白序列數(shù)據(jù)的提交和檢索、核酸序列的分析、蛋白質序列的分析以及蛋白質結構的預測等功能。

6、　　4.新城疫病毒分離株全基因組序列測定與分子生物學特性鑒定。根據(jù)GenBank發(fā)布的NDV全基因組序列，設計了10對引物，對由本室分離鑒定的3個鴨源弱毒株和國內(nèi)相關研究單位提供的4個鴿源、2個雞源毒株進行了全基因組序列測定，結果表明，3個鴨源毒株全基因組序列長度為15186nt，基因組AT含量為53％，GC含量為47％，4個鴿源和2個雞源株全基因組的長度為15192nt，基因組AT含量為54％，GC含量為46％。所測毒株各基岡的起始位

7、置與均已經(jīng)發(fā)表的毒株相一致。各毒株F蛋白裂解位點區(qū)域的序列符合其毒力強弱的特征，除鴨源毒株的HN蛋白長度為616和577個氨基酸外，其它均為571個氨基酸。與GenBank進行的同源比對以及與各基因型代表株的序列比對結果表明，所測毒株與已經(jīng)公布的毒株之間同源性最高的為95％左右，最低的為84％左右。基于F基因片段構建的系統(tǒng)進化樹顯示，3個鴨源毒株全部屬于基因I型，2個鴿源株(ND/05/028、ND/05/029)屬于基因VI型，另2個

8、鴿源株(ND/03/018、ND/03/044)屬于基因Ⅶ型，雞源株QH1和QH4屬于基因Ⅷ型。經(jīng)與GenBank檢索結果比較，本文首次測定了NDV基因Ⅷ型毒株的全基因序列，為開展不同基因型NDV全基因組的比較研究提供了條件。
　　5.鵝源新城疫病毒基因組密碼子使用分析。鵝的新城疫流行是最近10多年才出現(xiàn)的新情況。為了進一步研究鵝源NDV的分子生物學特性，本文以我國測序的兩株鵝源NDV全基因組序列(ZJ1、SF02)為材料，分析其

9、密碼子用法，并與鵝、雞的密碼子用法進行比較，在密碼子使用水平上作初步探索。結果表明，雖然鵝源NDV毒株與疫苗株、雞源株、鴿源株在基因組水平上存在較大差異，但在編碼蛋白的密碼子用法上沒有顯著差異，在密碼子的優(yōu)先使用和高頻密碼子數(shù)、有效密碼子數(shù)、密碼子使用指數(shù)等指標上基本一致。從與宿主密碼子使用頻率的比較上可以看出，鵝源NDV密碼子使用頻率與鵝、雞的密碼子使用頻率差異不大，僅有5-6個密碼子存在差異。
　　6.基于Z曲線的新城疫病毒基

10、因組生物信息學比較。利用從GenBank下載以及本室測序的NDV全基因組序列和本研究開發(fā)的Z曲線分析系統(tǒng)，分別繪制了NDV全基因組的X-n、Y-n、Z-n和三維Z曲線，比較了不同基因型NDV基因組Z曲線的特征，結果顯示，NDV全基因組中嘌呤堿基(A、G)、氨基堿基(A、C)占優(yōu)勢，在基因組序列的前1/3至前1/2區(qū)域，強氫鍵堿基(G、C)占優(yōu)勢。不同基因型的NDV基因組Z曲線，從ClassI到ClassⅡ的基因1-8型，有向Z軸靠攏的趨

11、勢，在基因組組成中表現(xiàn)為GC含量有降低的趨勢，這一進化規(guī)律的生物學意義值得進一步研究。通過對副粘病毒亞科小同屬病毒代表株全基因組序列Z曲線的比較，證明近年來將NDV與APMV2-9型劃為新的Avulavirus屬是正確的，從全基因組Z曲線上可以看出，它們與原副粘病毒屬的病毒有較大的區(qū)別。
　　7.不同基因型新城疫病毒全基因組的比較研究。
　　選擇53個NDV全基因組序列，包括ClassⅠ和ClassⅡ中基因Ⅰ-Ⅸ型毒株，利用

12、生物信息學分析方法在核酸水平和蛋白質水平開展比較研究，結果表明:NDV毒株中基因組長度有3種類型，即15186nt、15192nt、15198nt;基因組核酸序列中的堿基組成存在少許差異，且與病毒的基因型相關;3’端非編碼區(qū)RNA二級結構存在3種類型，其中有2個保守的莖.環(huán)結構，可能是與病毒復制與轉錄相關的重要功能區(qū);基因組編碼區(qū)序列和編碼蛋白序列的遺傳進化分析表明，NDV各基因的進化基本保持同步，核酸序列的同源性小于蛋白序列的同源性;

13、經(jīng)對F和HN蛋白功能位點預測和比較，發(fā)現(xiàn)不同基因型毒株存在一些特定的氨基酸位點;利用Insight II軟件模擬了雞源、鴿源和鵝源毒株的F、HN蛋白三級結構，對鵝源株三維結構中特征性位點及與周圍結構的關系進行了分析。
　　8.結論
　　(1)根據(jù)Z曲線理論，構建了基于Z曲線的基因組核酸序列可視化分析系統(tǒng)，提供了直觀地分析比較和研究基因組核酸序列的幾何學工具。
　　(2)根據(jù)新城疫病毒基因型分析的原理，利用JSP和數(shù)據(jù)庫

14、技術，構建了基于Web的新城疫病毒基因分型系統(tǒng)，實現(xiàn)了對NDV快速分型的目標。
　　(3)根據(jù)新城疫病毒基因組分析的需要，開發(fā)了基于Web的NDV生物信息分析系統(tǒng)，可以提供對NDV核酸和蛋白序列的提交、檢索查詢、核酸序列分析、蛋白質序列分析等功能。
　　(4)測定了9個NDV分離株的全基因組序列，其中基因Ⅷ型毒株為首次測序，并提交到GenBank核酸序列庫中，為開展不同基因型NDV全基因組的比較研究提供了條件。
　　(

15、5)對鵝源株NDV全基因組的密碼子使用情況進行了分析，確定了NDV基因組中優(yōu)先使用的密碼子情況，計算了密碼子使用的相關指數(shù)，并與鵝、雞的密碼子使用進行了比較。
　　(6)利用研究開發(fā)的Z曲線分析系統(tǒng)，對不同基因型NDV全基因組進行了幾何學的比較，探明了NDV基因組中堿基分布規(guī)律，發(fā)現(xiàn)在Z曲線表示上，從基因Ⅰ型到基因Ⅷ型有向Z軸靠攏的趨勢。
　　(7)對53株不同基因型的NDV全基因組序列進行了系統(tǒng)的比較研究，分別從基因組核酸