戴 銀，張丹俊，趙瑞宏，沈?qū)W懷，胡曉苗，侯宏艷，潘孝成，周學(xué)利

(安徽省農(nóng)業(yè)科學(xué)院畜牧獸醫(yī)研究所，安徽合肥 230031)

?

番鴨細(xì)小病毒安徽分離株結(jié)構(gòu)蛋白基因的克隆及序列分析

戴 銀*，張丹俊*，趙瑞宏，沈?qū)W懷，胡曉苗，侯宏艷，潘孝成，周學(xué)利

(安徽省農(nóng)業(yè)科學(xué)院畜牧獸醫(yī)研究所，安徽合肥 230031)

為研究番鴨細(xì)小病毒(MDPV)安徽分離株的遺傳變異特征，通過(guò)PCR擴(kuò)增獲得了MDPV結(jié)構(gòu)蛋白(VP)基因全長(zhǎng)序列AH-MDPV-VP，并將該序列與GenBank中登錄的12條MDPV和鵝細(xì)小病毒(GPV)VP基因序列進(jìn)行比對(duì)。結(jié)果顯示，AH-MDPV-VP基因全長(zhǎng)2 199 bp，包括完整的VP1、VP2和VP3蛋白編碼區(qū)。MDPV與GPV的VP基因部分序列一致，但具有明顯的差異。進(jìn)化分析顯示，MDPV安徽分離株與基因重組型MDPV上海分離株SAAS-SHNH為同一分支，親緣關(guān)系較近。同源性分析顯示二者核苷酸序列同源性最高，為99.9%，且AH-MDPV-VP與GPV毒株SHFX1201的序列同源性也有89.5%。此外安徽分離株與其他MDPV的VP1、VP2和VP3基因同源性有逐步下降的趨勢(shì)，而與GPV則相反呈上升趨勢(shì)。進(jìn)一步顯示MDPV安徽分離株與基因重組型MDPV上海分離株SAAS-SHNH相似，可能為MDPV和GPV基因重組型水禽細(xì)小病毒。

番鴨細(xì)小病毒；結(jié)構(gòu)蛋白基因；克??；序列分析

番鴨細(xì)小病毒(Muscovy duck parvovirus，MDPV) 屬于細(xì)小病毒科細(xì)小病毒亞科依賴病毒屬成員[1]。MDPV可引起番鴨的細(xì)小病毒病，該病是一種嚴(yán)重危害番鴨養(yǎng)殖業(yè)的急性傳染性病，可通過(guò)多種方式傳播，無(wú)明顯的季節(jié)性，但冬春兩季多發(fā)，1周齡～3周齡雛番鴨較為易感，一旦感染，發(fā)病率和病死率均很高。耐過(guò)病愈鴨大多生長(zhǎng)發(fā)育緩慢，其經(jīng)濟(jì)價(jià)值受到很大影響[2-5]。MDPV基因組均為線性、單鏈DNA，大小約為5 kb，包括兩個(gè)開放閱讀框(open reading frame，ORF)，左側(cè)編碼非結(jié)構(gòu)蛋白(NS)，右側(cè)編碼結(jié)構(gòu)蛋白(VP)。依據(jù)不同的起始密碼子和蛋白裂解位置，VP可生成VP1、VP2、VP3三種結(jié)構(gòu)蛋白。編碼VP1、VP2和VP3的核苷酸分別為2 199、1 764、1 605 bp，三者共用同一終止密碼，但起始密碼子的位置不同。VP2基因序列中包含VP3的全部基因序列，而VPl基因序列中又包含VP2、VP3的全部基因序列，形成套式結(jié)構(gòu)[6-7]。番鴨細(xì)小病毒病在世界各地均可發(fā)生，近年來(lái)，我國(guó)福建、浙江、安徽等地均有不同程度發(fā)生，給番鴨養(yǎng)殖業(yè)帶來(lái)巨大的經(jīng)濟(jì)損失[8-11]。本研究對(duì)安徽番鴨細(xì)小病毒分離株的結(jié)構(gòu)蛋白全長(zhǎng)基因進(jìn)行了克隆和序列分析，旨在研究該病毒的基因遺傳變異特征，為進(jìn)一步研究該病毒的基因功能和致病機(jī)理奠定基礎(chǔ)。

1 材料與方法

1.1 材料

1.1.1 毒株基因組DNA 番鴨細(xì)小病毒安徽分離株基因組DNA，安徽省農(nóng)業(yè)科學(xué)院畜牧獸醫(yī)研究所實(shí)驗(yàn)室分離、鑒定并保存。

1.1.2 菌株、質(zhì)粒、工具酶和試劑 大腸埃希菌DH5α菌株，本實(shí)驗(yàn)室提供;載體pMD18-T、DNA Marker DL 2000、T4連接酶等，寶生物工程(大連)有限公產(chǎn)品司；dNTP、Taq酶、瓊脂糖(電泳純)等，上海生工生物工程技術(shù)服務(wù)有限公司產(chǎn)品；DNA清潔/PCR產(chǎn)物清潔試劑盒、DNA凝膠回收試劑盒、質(zhì)粒小量抽提試劑盒，杭州維特潔生化技術(shù)有限公司產(chǎn)品。

1.2 方法

1.2.1 引物設(shè)計(jì) 參照GenBank中登錄的MDPV基因組序列(登錄號(hào)：KC171936)，用Primer Premier 5.0軟件設(shè)計(jì)兩對(duì)特異性擴(kuò)增引物。P-1：5′-TGCTATTGATGATATGGAGA-3′；P-2：5′-CAGGTTCGGAGCCTTCGGTG-3′；P-3：5′- CCGAACCTGTGGCAGCATCTAAAATGGCAGAG-3′；P-4：5′- TGATTGGCTGGTTCGAACGAACGAA-3′。預(yù)期目的擴(kuò)增片段大小分別為890 bp和1 733 bp，拼接后基因片段包括完整的VP基因的編碼區(qū)，命名為AH-MDPV-VP。引物由寶生物工程(大連)有限公司合成。

1.2.2 AH-MDPV-VP基因的分段擴(kuò)增 以保存的基因組DNA為模板，分別采用P1 、P2和P3 、P4為上、下游引物擴(kuò)增。PCR反應(yīng)體系：10×PCR buffer 5.0 μL，dNTP Mixture 3.0 μL, 上、下游引物各1.0 μL，DNA 2.5 μL，TaqDNA 聚合酶0.5 μL，加雙蒸水補(bǔ)足60 μL。分別按以下程序進(jìn)行擴(kuò)增：94℃ 5 min;94℃ 1 min，51℃ 30 s，72℃ 2 min，32個(gè)循環(huán);72℃ 10 min。94℃ 5 min；94℃ 30 s，56℃ 30 s，72℃ 3 min，32個(gè)循環(huán)；72℃延伸10 min。擴(kuò)增產(chǎn)物經(jīng)8 g/L瓊脂糖凝膠電泳檢測(cè)，用DNA凝膠回收試劑盒，然后連入pMDl8-T載體，16℃水浴過(guò)夜，轉(zhuǎn)化感受態(tài)大腸埃希菌DH5α，涂布于含氨芐青霉素(Amp)的LB瓊脂平板上，37℃溫箱培養(yǎng)過(guò)夜，初步篩選陽(yáng)性克隆，并用小量堿法提取質(zhì)粒，送上海生工生物工程技術(shù)服務(wù)有限公司測(cè)序。

1.2.3 基因序列分析 測(cè)序后的基因序列在NCBI (http://blast.ncbi.nlm.nih.gov/blast.cgi)數(shù)據(jù)庫(kù)中進(jìn)行同源性比對(duì)搜索, 利用SeqMan進(jìn)行序列拼接，獲得VP基因全長(zhǎng)核苷酸序列，并推導(dǎo)其氨基酸序列。選擇番鴨和鵝細(xì)小病毒株共12條參比序列作為分析樣本，參考毒株信息見表1。使用Clustal X 1.83軟件比對(duì)序列，采用DNA Star、MegAlign等軟件分析各序列特征。

表1 參考毒株信息

2 結(jié)果

2.1 番鴨AH-MDPV-VP基因的擴(kuò)增

以 MDPV的基因組DNA為模板，P1、 P2和P3、 P4為引物PCR分別擴(kuò)增獲得了一條約 890 bp 和 1 733 bp 大小的特異性DNA片段(圖1)，將擴(kuò)增的基因產(chǎn)物純化、測(cè)序。測(cè)序結(jié)果與GenBank中登錄的相應(yīng)基因序列進(jìn)行比對(duì)，均為目的基因序列。

A.P1、 P2引物的PCR產(chǎn)物;B.P3、 P4引物的PCR產(chǎn)物；M.DNA 標(biāo)準(zhǔn)DL 2 000;1、2.基因組DNA的PCR產(chǎn)物

2.2 基因序列特征分析

序列分析可見，AH-MDPV-VP基因全長(zhǎng)2 199 bp，包含完整MDPV VP1、VP2和VP3蛋白質(zhì)閱讀框，分別為2 199、1 764、1 605 bp，編碼733、588、534個(gè)氨基酸。MDPV和GPV的VP基因核苷酸均未見缺失，VP2的起始密碼子ACG均位于VPl基因序列的436 bp處，VP3的起始密碼子ATG位于VPl基因序列的595 bp處，三者共用同一終止密碼子。同時(shí)可見MDPV與GPV的VP基因序列存在較多一致的核苷酸位點(diǎn)，但部分序列差異較為明顯，尤其是MDPV位于VP3區(qū)域的核苷酸替代較為顯著。

2.3 序列進(jìn)化分析

采用MEGA分析軟件，將獲得的13條番鴨和鵝的細(xì)小病毒VP基因序列構(gòu)建系統(tǒng)發(fā)育樹(圖2)。結(jié)果可見，所比對(duì)序列分為番鴨和鵝細(xì)小病毒兩個(gè)類群。AH-MDPV-VP與國(guó)內(nèi)MDPV上海分離株SAAS-SHNH的基因序列KC171936，以及廣西分離株MDPV-GX5的基因序列KM093740聚集為同一組，尤其是與分離株SAAS-SHNH的遺傳距離較近；與國(guó)內(nèi)較早分離株的基因序列JF926698、AY510603、JF926697，以及歐洲分離株基因序列X75093和Z68272遺傳距離較遠(yuǎn)。

遺傳距離在圖的下方；聚類分析的可靠性用bootstrap置信區(qū)限檢測(cè)，同時(shí)每1 000 bootstrap測(cè)試結(jié)果大于50%的在節(jié)點(diǎn)上顯示

2.4 同源性分析

將安徽分離株全長(zhǎng)基因序列AH-MDPV-VP與其他在GenBank獲得的番鴨和鵝的VP基因序列進(jìn)行同源性比較(表2)。結(jié)果表明，AH-MDPV-VP與MDPV分離株SAAS-SHNH的基因序列同源性最高，為99.9%，與MDPV福建P1毒株同源性較低，為88.9%，而與GPV上海SHFX1201毒株序列的同源性達(dá)89.5%。進(jìn)一步分析可見，安徽分離株與其他MDPV的VP1、VP2、VP3同源性分別介于99.9%～88.9%、99.9%～86.6%、99.9%～85.4%之間，從VP1到VP3同源性有逐步下降的趨勢(shì)；與鵝細(xì)小病毒VP1、VP2、VP3同源性分別介于89.5%～88.6%、91.0%～89.0%、93.2%～91.3%之間，從VP1到VP3的同源性則有逐步上升的趨勢(shì)。

3 討論

本研究通過(guò)分段擴(kuò)增獲得了MDPV安徽分離株的全長(zhǎng)結(jié)構(gòu)蛋白基因序列AH-MDPV-VP，該基因全長(zhǎng)2 199 bp，包括完整MDPV VP1、VP2和VP3蛋白質(zhì)編碼區(qū)，分析可見3個(gè)蛋白質(zhì)具有不同的起始密碼子，但共用同一終止密碼子。這與目前已有相關(guān)報(bào)道基本一致[6,12]。此外，序列比對(duì)可見，雖然MDPV與GPV的VP基因部分序列一致，但仍具有明顯的差異，尤其是MDPV位于VP3區(qū)域的核苷酸替代較為明顯，變異程度較高。系統(tǒng)發(fā)育樹進(jìn)一步顯示，AH-MDPV-VP與近年在國(guó)內(nèi)分離的MDPV上海毒株SAAS-SHNH和廣西MDPV-GX5分離株的基因序列KC171936和KM093740為同一小分支，尤其是與分離株SAAS-SHNH親緣關(guān)系較近，二者來(lái)源于共同祖先；與國(guó)內(nèi)較早分離的4個(gè)毒株，以及歐洲分離的2毒株基因序列遺傳距離較遠(yuǎn)。該研究結(jié)果表明MDPV的遺傳變異規(guī)律與分離株發(fā)生的地理和時(shí)間呈現(xiàn)一定程度的相關(guān)性，但是否確切仍需要研究更多不同地點(diǎn)和時(shí)間分離的MDPV毒株。

表2 VP基因核苷酸序列同源性

研究發(fā)現(xiàn)MDPV與GPV存在一定交叉反應(yīng)， MDPV只能感染番鴨，但GPV既可感染鵝，也可感染番鴨，在生產(chǎn)中小鵝瘟活疫苗和抗血清也能夠提高雛番鴨對(duì)MDPV的抵抗能力[13-15]。依據(jù)13個(gè)VP基因核苷酸序列的同源性比對(duì)結(jié)果，可見MDPV和GPV的全長(zhǎng)VP基因序列同源性介于89.5%～79.6%之間，安徽分離株與上海分離株SAAS-SHNH的同源性最高，達(dá)99.9%，與MDPV福建P1毒株同源性相對(duì)較低，為88.9%，但與GPV上海SHFX1201毒株也有89.5%的同源性。MDPV和GPV基因核苷酸序列如此高的同源性，也正說(shuō)明了二者具有交叉反應(yīng)特性的根本原因。此外，有研究表明，某些病毒在自然遺傳進(jìn)化過(guò)程中會(huì)發(fā)生基因變異、交叉重組現(xiàn)象，從而產(chǎn)生一種新的致病病毒。這種現(xiàn)象可能更有利病毒逃避寄主的免疫應(yīng)答，提高本身適應(yīng)環(huán)境的能力。一些動(dòng)物細(xì)小病毒就存在這種基因重組現(xiàn)象，其中MDPV分離株SAAS-SHNHB被認(rèn)為是由MDPV和GPV基因重組而來(lái)的水禽細(xì)小病毒[1,16]。VP1、VP2、VP3三者同源性分析表明，安徽分離株與其他番鴨細(xì)小病毒株，從VP1到VP2再到VP3基因同源性有逐步下降的趨勢(shì)，而與鵝細(xì)小病毒從VP1到VP3的同源性則有逐步上升的趨勢(shì)，表明安徽分離株與MDPV和GPV不同節(jié)段的VP基因各有較高的同源性。這些特點(diǎn)顯示了安徽分離株與親緣關(guān)系較近分離株SAAS-SHNH相似，可能為基因重組型水禽細(xì)小病毒。本試驗(yàn)初步研究了安徽MDPV分離株結(jié)構(gòu)蛋白VP基因的分子結(jié)構(gòu)特征，將為深入研究該基因型番鴨細(xì)小病毒的遺傳變異特點(diǎn)，以及獲得有效的預(yù)防控制方法奠定了基礎(chǔ)。

[1] 萬(wàn)春和，陳紅梅，傅秋玲，等.番鴨細(xì)小病毒浙江分離株VP基因的克隆與序列分析[J].中國(guó)畜牧獸醫(yī)，2015,42(10):2600-2605.

[2] 謝麗基，謝芝勛，劉加波，等.番鴨細(xì)小病毒熒光定量PCR方法的建立及初步應(yīng)用[J].中國(guó)預(yù)防獸醫(yī)學(xué)報(bào),2013,35(3):218-221.

[3] 陳 瑤，關(guān) 淼.雛番鴨細(xì)小病毒病的診斷與防制[J].現(xiàn)代畜牧獸醫(yī),2010(4):42-44.

[4] 程曉霞，陳仕龍，陳少鶯，等.番鴨細(xì)小病毒和鵝細(xì)小病毒的抗原相關(guān)性研究[J].福建農(nóng)業(yè)學(xué)報(bào),2013,28(9):869-871.

[5] 崔國(guó)杰，陳少鶯.番鴨細(xì)小病毒和小鵝瘟病毒混合感染的診治[J].畜牧與獸醫(yī)，2015,47(3):147-148.

[6] 王建業(yè)，黃 鈺，龔建森，等.番鴨細(xì)小病毒鵝胚化弱毒株FZ91-30的全基因組克隆及序列分析[J].中國(guó)家禽，2015,37(22):46-48.

[7] 季 芳，張毓金，楊增岐，等.番鴨細(xì)小病毒和鵝細(xì)小病毒廣東株VPl基因的克隆與序列分析[J].中國(guó)預(yù)防獸醫(yī)學(xué)報(bào),2004,26(4):245-247.

[8] Woolcock P R,Jestin V,ShivaPrasad H L,et al.Evidence of muscovy duck parvovirus in muscovy duckling in California[J].Vet Rec,2000,1:68-72.

[9] 宋永峰，吳發(fā)興，宋延華，等.2008年—2009年廣東番鴨細(xì)小病毒分子流行病學(xué)研究[J].中國(guó)動(dòng)物檢疫,2010(1):39-40.

[10] 王文秀，張松林，付 強(qiáng)，等.番鴨細(xì)小病毒ZQ分離株非結(jié)構(gòu)蛋白基因的克隆及序列分析[J].動(dòng)物醫(yī)學(xué)進(jìn)展,2013,34(12):79-83.

[11] 張丹俊，戴 銀，沈?qū)W懷，等.番鴨細(xì)小病毒的VP3基因克隆和序列分析[J].動(dòng)物醫(yī)學(xué)進(jìn)展,2015,36(11):7-10.

[12] 張 云，耿宏偉，郭東春，等.鵝和番鴨細(xì)小病毒全基因克隆與序列分析[J].中國(guó)預(yù)防獸醫(yī)學(xué)報(bào),2008,30(6):415-419.

[13] Zadori Z，Stefancsik R,Rauch T,et al.Analysis of complete nueleotide sequence of goose and muscovy duck parvovirus indicates common ancestral origin with adeno-associated virus [J].Virol,1995,212:562-573．

[14] 劉 偉，程曉霞，陳少鶯.番鴨源小鵝瘟病毒和番鴨細(xì)小病毒的結(jié)構(gòu)蛋白及其抗原性[J].福建農(nóng)林大學(xué)學(xué)報(bào)：自然科學(xué)版,2015,44(1):64-68.

[15] 范文勝,吳陽(yáng)開，方文棋，等.番鴨細(xì)小病毒與鵝細(xì)小病毒鑒別診斷方法的建立及初步應(yīng)用[J].中國(guó)家禽,2014,36(22):20-23.

[16] Zhu Y,Zhou Z,Huang Y,et al.Identification of a recombinant Muscovy duck parvovirus (MDPV) in Shanghai,China [J].Vet Microbiol,2014,174(3-4):560-564.

Cloning and Sequence Analysis of Structural Protein Gene of Muscovy Duck Parvovirus Anhui Isolate

DAI Yin,ZHANG Dan-jun,ZHAO Rui-hong,SHEN Xue-huai,HU Xiao-miao,HOU Hong-yan,PAN Xiao-cheng,ZHOU Xue-li

(InstituteofAnimalHusbandryandVeterinaryScience,AnhuiAcademyofAgriculturalScience,Hefei,Anhui,230031,China)

To investigate the variation status of muscovy duck parvovirus(MDPV) Anhui isolate,the structural protein (VP) gene was cloned by PCR,which was named as AH-MDPV-VP. And we compared AH-MDPV-VP with the sequence of twelve VP genes of MDPV and goose parvovirus (GPV) in the GenBank database.The results showed that AH-MDPV-VP gene was 2 199 bp,coding full-length open reading frames of VP1,VP2,and VP3.The partial sequences were consistent,but obvious difference was observed between MDPV and GPV genes.Evolution analysis revealed that Anhui isolate and the genetic recombinant MDPV of Shanghai isolate SAAS-SHNH gathered together in the same branch,they had a close genetic relationship.Similarity analysis showed that the nucleotide sequence homology between two isolates was the highest,99.9%,and the sequence homology of Anhui isolate and GPV isolate SHFX1201 was 89.5%.Moreover,the VP1,VP2 and VP3 homology of Anhui isolate and other MDPV had a tendency to decline gradually.On the contrast,the rising trend of the homology with GPV was found.The result further showed that MDPV Anhui isolate and genetic recombinant Shanghai SAAS-SHNH isolate were similar,they may be waterfowl parvoviruses,recombined with MDPV and GPV.

Muscovy duck parvovirus；structural protein gene；cloning；sequence analysis

2016-07-19

安徽省農(nóng)業(yè)科學(xué)院人才發(fā)展專項(xiàng)資金項(xiàng)目(17F0404);安徽省家禽產(chǎn)業(yè)技術(shù)體系資金項(xiàng)目

戴 銀(1980-)，女，安徽蒙城人，副研究員，博士，主要從事獸醫(yī)微生物與免疫學(xué)研究。

S852.659.4;Q789

A

1007-5038(2017)04-0024-04

*通訊作者