鵝細小病毒遼寧分離株全基因測序與遺傳演化分析

李閣錦，郭維軍，曹祁峰，周鐵忠，李 冰*

（1. 錦州醫科大學畜牧獸醫學院，遼寧錦州121000；2. 錦州市動物疫病預防控制中心，遼寧錦州121001）

小鵝瘟又稱鵝細小病毒病，由鵝細小病毒（goose parvovirus，GPV）引起4～20日齡雛鵝的一種急性或亞急性敗血癥，以發生滲出性腸炎為主要特征，傳播快、病死率高。本病最早由方定一于1956年發現，在我國江蘇揚州地區流行，并首次報道該病的流行特點和臨診特征[1]。1961年，方定一團隊對該病病原進行系統研究，確定該病原是一種和新城疫、鴨病毒性肝炎無關的病毒，建議命名為小鵝瘟病毒[2]。1967年，匈牙利學者Derzsy利用13日齡鵝胚成功分離出小鵝瘟病毒[3]。1973年，Schettler 利用鵝胚成纖維細胞也成功分離到該病毒，病毒的理化特性與其他已知的禽類病毒不同，具備細小病毒的特征，建議將該病毒歸類為細小病毒[4]。1974年，世界家禽學會（World's Poultry Science Association，WPSA）為紀念Derzsy 對小鵝瘟研究所做出的貢獻，將該病命名為Derzsy's 病[5]。1979年，國際病毒分類委員會（The International Committee on Taxonomy of Viruses，ICTV）第3次報告將該病毒正式命名為鵝細小病毒[6]。

GPV是細小病毒科（Parvoviridae）成員，歸屬于細小病毒亞科（Parvovirinae）、依賴病毒屬（Dependovirus）、雁形目依 賴 細 小 病 毒1 型（Anseriform dependoparvovius 1）[7]。GPV為單股線狀DNA病毒，基因組全長約5 kb，編碼區依次編碼非結構蛋白NS1、NS2 和結構蛋白VP1、VP2、VP3[8]，非編碼區由回文序列（inverted terminal repeats，ITR）組成位于編碼區兩側，依次為5'-ITP-VS1-VS2-VP1-VP2-VP3-ITR-3'。非結構蛋白NS1 和NS2，主要是參與病毒的轉錄與復制[9]；結構蛋白VP1、VP2和VP3共同組成衣殼蛋白。Tullis等[10]報道，VP1蛋白在病毒感染中發揮重要作用，可以引導衣殼蛋白進入細胞核，與宿主細胞受體結合。VP3是衣殼蛋白的主要組成部分，可占衣殼蛋白總量的80%，內含GPV 主要抗原決定簇，是誘導機體產生中和抗體的主要蛋白[11]。ITR 為倒置重復序列，位于基因組的兩端，頭部可以折疊形成U型發夾結構，ITR主要參與病毒復制與組裝的過程[12]。本研究對遼寧地區分離的鵝細小病毒進行全基因測序，并與國內外代表毒株進行序列比對分析，為遼寧地區鵝細小病毒的分子流行病學研究提供參考。

1 材料與方法

1.1 試驗材料

鵝細小病毒由遼寧省錦州市疾病預防控制中心分離保存；病毒基因組DNA提取試劑盒、質粒小提試劑盒購自天根生化科技有限公司；2×Taq Master Mix（含染料）購自北京康為世紀生物科技有限公司；快速瓊脂糖凝膠DNA回收購自南京維諾贊生物科技有限公司；DL2000 DNA Marker、pMD-18T Vector 購自寶生物工程有限公司；DH5α感受態細胞購自北京索萊寶科技有限公司。

1.2 試驗方法

1.2.1 病毒總DNA的提取

取200 μL 細胞培養物于滅菌1.5 mL 離心管中，用病毒基因組DNA 提取試劑盒，按照使用說明書從細胞培養物中提取病毒總DNA后，置于-20 ℃保存備用。

1.2.2 引物的設計與合成

根據GenBank 中GPV 毒株DY16（登錄號：MH209633）設計8對引物，對分離的病毒全基因序列進行分段擴增，引物信息詳見表1，由上海生工生物工程股份有限公司合成。

表1 全基因引物序列Tab.1 Whole gene primer sequence

1.2.3 全基因序列擴增及克隆測序

采用PCR 方法對全基因序列進行分段擴增。PCR反應體系（30 μL）如 下：2×Taq Master Mix 15 μL、ddH2O 10 μL、引物2 μL、模板DNA 3 μL，按照表1 中退火溫度進行擴增，擴增產物進行膠回收，純化后產物連接到pMD-18T 載體上并轉入DH5α 大腸桿菌感受態細胞中，涂布含氨芐西林（100 μg/mL）的LB 固體培養基上進行篩選，過夜培養12～16 h，挑取單個菌落，經PCR 驗證為陽性，使用質粒小提試劑盒提取質粒，放于-20 ℃保存，將所提陽性質粒送至上海生工生物科技有限公司進行測序，每個樣品重復擴增及克隆測序3 次。

1.2.4 全基因序列比對分析（見表2）

測序結果經NCBI-BLAST 比對正確后，采用DNAMAN 進行各段序列拼接，獲得完整的全基因序列。采用MEGA 7.0、DNAStar 7.0將遼寧分離株鵝細小病毒全基因序列與GenBank 中公布的國內外代表毒株進行序列 比對分析，并制作進化樹。

表2 全基因比對代表毒株信息Tab.2 Full gene comparison representative strain information

2 結果與分析

2.1 GPV 遼寧分離株全基因各片段PCR 擴增結果（見圖1）

由圖1 可知，利用PCR 方法，成功擴增出GPV 遼寧分離株全基因各片段，PCR 產物經瓊脂糖凝膠電泳鑒定，目的片段與預期片段大小基本一致。

圖1 GPV遼寧分離株全基因各片段PCR產物Fig.1 PCR products of all gene fragments of GPV Liaoning isolate

2.2 GPV 遼寧地區分離株全基因序列拼接與分析（見表3）

表3 LNGPV20毒株與代表毒株各段基因序列比對Tab.3 Alignment of gene sequences between LNGPV20 strain and representative strain

采用DNAMAN 軟件，將測序結果進行拼接、校正，經BLAST 驗證獲得GPV 遼寧分離株全基因序列，命名為LNGPV20。分離毒株基因全長5106 bp，共有兩個非結構蛋白NS1、NS2和3個結構蛋白VP1、VP2、VP3，ITR序列長為444 bp。由表3可知，LNGPV20毒株各段基因序列長度與歐洲標準株B 毒株各段基因序列長度一致，LNGPV20毒株的結構蛋白與非結構蛋白序列長度與DY16、RC16、HB-2019、SYG61v毒株相同，與82-0321、SQ0412、E、FJ01毒株各蛋白序列長度均存在一定差異。其中ITR 基因序列長度變化最大，序列長度在414～447 bp之間，相差33個堿基。14個毒株各段基因序列同源性為91.13%～98.00%，各基因片段核苷酸同源性較高，其中非結構蛋白NS1核苷酸同源性為98.00%，同源性最高。

2.3 基于ITR序列比對分析（見圖2）

由圖2可知，GPV全基因組5'端與3'端ITR序列相同，呈末端倒置重復序列。LNGPV20 毒株ITR 由444 個堿基組成，與標準毒株B 株相比ITR 的5'端缺少一個堿基C，ITR外側405個堿基形成一個U型的發夾結構，由181個堿基對形成的“莖部”和43個堿基形成的“泡區”構成，“泡區”197～202 陰影區為SphⅠ內切酶的酶切位點。ITR 內側有39 個堿基，不參與U 型發夾結構的形成，但全基因組5'端與3'端ITR內側的39個堿基呈反向互補，不影響基因組形成回文結構。LNGPV20毒株與B毒株相比，在181個堿基對區域存在20個堿基的差異，但不影響互補配對。

圖2 LNGPV20毒株5'端ITR序列和二級結構Fig.25'ITR sequence and secondary structure of LNGPV20 strain

2.4 LNGPV20毒株全基因序列比對分析（見圖3）

由圖3 可知，采用DNAstar 軟件Jotun Hein 方法對LNGPV20 毒株與已公布于GenBank 的14 個代表毒株進行基因序列的同源性分析，LNGPV20 毒株與14 個毒株全基因序列同源性在93.2%～98.9%之間，同源性均在90%以上。由此可見，鵝細小病毒全基因序列變化差異不大，比較保守。LNGPV20毒株與DY16、RC16毒株同源性最高，均為98.9%；與標準毒株B 株同源性為96.6%，而與GDaGPV、GPV-98E、SYG61v 疫苗毒株同源性較低，分別為93.2%、93.4%、93.5%。

2.5 基于全基因序列的系統進化分析（見圖4）

由圖4可知，采用MEGA7.0軟件，根據LNGPV20毒株與已公布于GenBank的20個參考毒株基因序列，構建系統進化樹。進化樹顯示鵝細小病毒的分離株與番鴨細小病毒的分離株分別處于兩大遺傳分支上，鵝細小病毒的疫苗株與野毒株親緣關系較遠，處于兩個小分支內。LNGPV20 毒株與DY16、RC16、HB-2019、SQ0412 毒株親緣關系較近，與標準毒株B株親緣關系較遠處于不同分支內。由此可見，國內流行的鵝細小病毒基因組變異不明顯，同源性較高，親緣關系較近。

圖3 LNGPV20毒株核苷酸比對分析Fig.3 Nucleotide comparison analysis of LNGPV20 strain

圖4 20株GPV全基因進化樹Fig.420 GPV whole gene evolutionary trees

3 討論

在細小病毒科所有成員中，GPV、MDPV 與人細小病毒AAV-2 親緣關系最近。但AAV-2 只有與腺病毒或皰疹病毒等輔助病毒共同感染時，才能進行病毒復制[11]。而GPV和MDPV不需要輔助病毒，主要通過內吞途徑將衣殼蛋白帶入細胞，進入細胞質，將病毒DNA 轉運到細胞核中，進行DNA 復制[13]。GPV 的基因組比MDPV 的基因組短，GPV 的ITR 在414～447 之間，MDPV 的ITR 在421～457之間。由此可知，GPV、MDPV 的ITR 長度差異較大[14]。LNGPV20 毒株基因組全長5106 bp，ITR 序列長444 bp，ITR內側有39個堿基不參與ITR的U型發夾結構的形成，而B 株ITR 內側有38 個堿基不參與ITR 的U 型發夾結構形成，可能與LNGPV20毒株ITR的5'端缺少一個堿基C有關。有研究表明，ITR 的長度，是影響毒株毒力的關鍵[15]，但宋娜[16]對2 株江蘇省鵝細小病毒分離株，ITR 分別為416 bp 和444 bp，且毒力均較強，因此推斷GPV 的毒力強弱并不能依據ITR 的長度來判斷。各毒株結構蛋白與非結構蛋白同源性在97.64%～98.00%之間，相對于ITR 差異不大。近年來，學者對GPV非結構蛋白與結構蛋白抗原表位區進行定位分析[17]，經郭鷺等[18]、黎明等[19]對GPV 結構蛋白、非結構蛋白進行抗原表位的精準定位，抗原表位為VP3 蛋白430～435 aa 和643～647 aa、NS1 蛋白503～509 aa，為疫苗制備、特異性抗體制備等奠定理論基礎。

遼寧分離毒株與14 個毒株全基因序列同源性比對分析，與GDaGPV、GPV-98E、SYG61v 疫苗毒同源性均在93%以上，除去ITR 序列，疫苗株與國內11 個野毒株的NS、NP 蛋白同源性均在97%以上，由于疫苗株與野毒株NS、NP同源性較高，且大部分抗原決定簇集中分布于NS、NP 蛋白[20]，可見國內小鵝瘟疫苗對小鵝瘟應有較高的保護率，但由于養鵝業飼養環境的差異，導致母源抗體參差不齊，大多數種鵝未進行免疫接種且呈隱性感染，雖然無發病癥狀，但仍可以散播病毒，導致小鵝瘟對養殖業的危害越來越嚴重[21]。另外，遼寧分離株與國內野毒株同源性在96.6%～98.9%之間，其中有4 個病毒株同源性在98%以上，分別源于安徽、重慶、河北、山東。由此可見，國內流行的鵝細小病毒基因組變異不明顯，同源性較高，可能是近年來養鵝業迅速發展，鵝苗流通增加，交叉感染所致。

4 結論

本研究采用PCR 方法對鵝細小病毒全基因分段進行擴增，首次獲得遼寧分離株全基因序列，該研究結果進一步豐富我國鵝細小病毒基因數據庫，同時也為遼寧省鵝細小病毒的防控提供參考。