葛志琪 ，劉亞剛 *， 陳 舜 ，朱德康 ，劉馬峰 ，汪銘書 ，賈仁勇 ，孫昆峰 ，楊 喬 ，吳 英 ，趙新新 ，程安春 *

（1.西南民族大學(xué)生命科學(xué)與技術(shù)學(xué)院，成都 610041；2.四川農(nóng)業(yè)大學(xué)動(dòng)物醫(yī)學(xué)院，預(yù)防獸醫(yī)研究所，成都 611130；3.四川農(nóng)業(yè)大學(xué)動(dòng)物醫(yī)學(xué)院，禽病防治與研究中心，成都 611130；4.四川農(nóng)業(yè)大學(xué)動(dòng)物醫(yī)學(xué)院，動(dòng)物疾病與人類健康四川省重點(diǎn)實(shí)驗(yàn)室，成都 611130）

新型鵝細(xì)小病毒四川株的全基因組擴(kuò)增及遺傳進(jìn)化分析

葛志琪1，劉亞剛1*， 陳 舜2，3，4，朱德康2，3，4，劉馬峰2，3，4，汪銘書2，3，4，賈仁勇2，3，4，孫昆峰2，3，4，楊 喬2，3，4，吳 英2，3，4，趙新新2，3，4，程安春2，3，4*

（1.西南民族大學(xué)生命科學(xué)與技術(shù)學(xué)院，成都 610041；2.四川農(nóng)業(yè)大學(xué)動(dòng)物醫(yī)學(xué)院，預(yù)防獸醫(yī)研究所，成都 611130；3.四川農(nóng)業(yè)大學(xué)動(dòng)物醫(yī)學(xué)院，禽病防治與研究中心，成都 611130；4.四川農(nóng)業(yè)大學(xué)動(dòng)物醫(yī)學(xué)院，動(dòng)物疾病與人類健康四川省重點(diǎn)實(shí)驗(yàn)室，成都 611130）

【目的】了解四川地區(qū)新型鵝細(xì)小病毒（Novel goose parvovirus，NGPV）的分子生物學(xué)特征及遺傳變異情況?！痉椒ā扛鶕?jù)GenBank中發(fā)表的新型鵝細(xì)小病毒的全基因組序列，設(shè)計(jì)5對(duì)兼并引物采用PCR方法分段擴(kuò)增新型鵝細(xì)小病毒四川株的全基因組序列，并利用DNAstar等生物信息分析軟件對(duì)獲得的基因組序列進(jìn)行分子特性與遺傳演變分析?！窘Y(jié)果】獲得的新型鵝細(xì)小病毒四川株（命名為SC16）全長(zhǎng)5109 nt，5′端ITR長(zhǎng)408 nt，3′端ITR長(zhǎng)407 nt，NS基因長(zhǎng)1884 nt，VP基因長(zhǎng)2199 nt。序列比對(duì)和遺傳進(jìn)化分析顯示，SC16株與山東分離株SDLC01（KT343253.1）及QH15（KT751090.1）的基因組核苷酸同源性高達(dá)99%，與兩者的親緣關(guān)系最近，基于NS及VP蛋白的氨基酸序列系統(tǒng)進(jìn)化樹(shù)分析結(jié)果與此一致。但SC16株的5′端ITR還是發(fā)生了一定的變異：與SDLC01及QH15相比，多了幾個(gè)14 nt的插入片段。【結(jié)論】研究結(jié)果豐富了NGPV的分子生物學(xué)資料。

新型鵝細(xì)小病毒；鴨短喙侏儒綜合征；全基因組擴(kuò)增；遺傳進(jìn)化分析

傳統(tǒng)的水禽細(xì)小病毒主要包括鵝細(xì)小病毒（Goose parvovirus，GPV）和番鴨細(xì)小病毒（Muscovy duck parvovirus，MDPV）[1]。鵝細(xì)小病毒對(duì)雛鵝和雛番鴨有較高的致病性，10日齡以內(nèi)的雛鵝的發(fā)病率和死亡率可高達(dá)100%[2]，主要病理特征是纖維素性滲出性腸炎和腸道形成堵塞物。番鴨細(xì)小病毒主要侵害1～3周齡的雛番鴨，俗稱“三周病”，臨診特征為腹瀉、喘氣和腳軟，具有高度傳染性，其發(fā)病率和死亡率都較高，易感動(dòng)物除雛番鴨外，未見(jiàn)其他禽類或哺乳類動(dòng)物發(fā)生本病。

自2015年以來(lái)，與鵝細(xì)小病毒密切相關(guān)的中等致病性的新型鵝細(xì)小病毒（Novel goose parvovirus，NGPV）在我國(guó)福建省[3]、山東省[4]、安徽省[5]等省份的半番鴨和櫻桃谷鴨中暴發(fā)流行，鴨群感染后臨床表現(xiàn)為發(fā)育遲緩、上下喙萎縮、舌頭外伸、腫脹、向下彎曲等，俗稱大（長(zhǎng)）舌病、鴨短喙長(zhǎng)舌綜合征、鴨短喙-侏儒綜合征（beak atrophy and dwarfism syndrome，BADS）等。該病多侵害10～30日齡的肉鴨,發(fā)病率10%～100%，死亡率2%～10%。雖然死亡率明顯低于鵝細(xì)小病毒感染和番鴨細(xì)小病毒感染，但感染鴨群料肉比顯著升高，患鴨出欄體重較健康鴨降低20%～30%[6]，嚴(yán)重者僅為正常體重的50%[7]。部分患鴨還出現(xiàn)單側(cè)行走困難、癱瘓等癥狀，在抓撲、屠宰過(guò)程中喙部、翅部骨骼和脛骨等易發(fā)生骨折。因此僵鴨淘汰率或殘鴨率高達(dá)20%～80%。目前，該病的發(fā)病區(qū)域不斷擴(kuò)大，給我國(guó)養(yǎng)鴨業(yè)造成了巨大的經(jīng)濟(jì)損失。

鴨短喙長(zhǎng)舌綜合征是近兩年來(lái)暴發(fā)流行的一種疾病，對(duì)其病原NGPV的研究目前還比較有限，特別是四川地區(qū)還未見(jiàn)報(bào)道。我們查閱GenBank后發(fā)現(xiàn)，目前僅獲得4個(gè)NGPV的全基因組序列，更多NGPV毒株的全基因序列測(cè)定將有助于理解NGPV的遺傳演化、病毒嗜性、致病機(jī)理及免疫應(yīng)答規(guī)律等相關(guān)信息，為鴨短喙長(zhǎng)舌綜合征的防治及疫苗研究提供科學(xué)依據(jù)。

1 材料和方法

1.1 病料來(lái)源

2016年11月，四川省峨眉山市某櫻桃谷鴨場(chǎng)暴發(fā)了疑似鴨短喙長(zhǎng)舌綜合征，該鴨場(chǎng)鴨子的發(fā)病率為10%左右。送檢的鴨子均具有明顯的短喙、長(zhǎng)舌、瘦小等癥狀。

1.2 試驗(yàn)材料

1.3 病理解剖

對(duì)送檢的鴨子進(jìn)行病理解剖，分別采集心、肝、脾、小腸等組織備用，隨時(shí)觀察并記錄患鴨及死鴨的臨床癥狀及病理變化。

1.4 PCR檢測(cè)

1.4.1 引物設(shè)計(jì)與合成

根據(jù)GenBank中發(fā)表的新型鵝細(xì)小病毒的全基因組序列，利用Primer Premier 5.0在NGPV VP3保守區(qū)設(shè)計(jì)一對(duì)特異性引物，引物序列如表1所示，由成都擎科梓熙生物技術(shù)有限公司合成。

1.4.2 DNA模板的提取

分別取病鴨及死鴨的肝臟和脾臟，按1∶5加入PBS，充分研磨，6000 r/min，離心 30 min，取上清，利用病毒基因組DNA/RNA提取試劑盒提取病料中的DNA。

表1 鑒定引物序列Table1 Identification primer sequences

1.4.3 PCR擴(kuò)增

取1.4.2中提取的DNA作為PCR擴(kuò)增的模板。PCR 反應(yīng)體系為：DNA 模板 1 μL，PrimeSTARMax DNA Polymerase 25 μL，上下游引物各 1 μL，ddH2O補(bǔ)足至50 μL。PCR 反應(yīng)條件如下：98℃ 2 min，98℃30 s，60℃ 15 s，72℃ 5 s，30 個(gè)循環(huán)，72℃延伸 5 min。PCR反應(yīng)結(jié)束后，產(chǎn)物經(jīng)1.5%瓊脂糖凝膠進(jìn)行電泳檢測(cè)，疑似陽(yáng)性者，用DNA純化回收試劑盒對(duì)PCR產(chǎn)物進(jìn)行液回收，然后分別連接pMD19-T載體，轉(zhuǎn)化至大腸桿菌SURE株，菌落PCR鑒定陽(yáng)性者送華大基因測(cè)序鑒定。

1.5 全基因組擴(kuò)增

1.5 .1 引物設(shè)計(jì)與合成

根據(jù)GenBank中發(fā)表的新型鵝細(xì)小病毒的全基因組序列及結(jié)構(gòu)特征，并參照陳浩等人的引物設(shè)計(jì)方法[8]，利用Primer Premier5.0設(shè)計(jì)5對(duì)引物來(lái)擴(kuò)增全基因組，擴(kuò)增片段分別以P1～P5表示，引物序列如表2所示，由成都擎科梓熙生物技術(shù)有限公司合成。

表2 全基因組序列擴(kuò)增引物Table2 The primers of complete genome sequence amplification

1.5 .2 目的片段的分段擴(kuò)增

取1.4.2中提取的死鴨組織DNA作為PCR擴(kuò)增的模板。每個(gè)PCR反應(yīng)體系為50 μL：DNA模板1 μL，PrimeSTARMax DNA Polymerase 25 μL，上下游引物各 1 μL，ddH2O 補(bǔ)足至 50 μL。PCR 反應(yīng)條件參照PrimeSTAR高保真酶的產(chǎn)品說(shuō)明書進(jìn)行：98 ℃ 2 min，98 ℃ 30 s，60 ℃ 15 s，72 ℃ 5～15 s（根據(jù)擴(kuò)增片段長(zhǎng)度做相應(yīng)調(diào)整），35個(gè)循環(huán)，最后72℃延伸5 min。PCR反應(yīng)結(jié)束后，產(chǎn)物經(jīng)1.5%瓊脂糖凝膠進(jìn)行電泳檢測(cè)，疑似陽(yáng)性者，用DNA純化回收試劑盒對(duì)PCR產(chǎn)物進(jìn)行液回收，然后分別連接pMD19-T載體，轉(zhuǎn)化至大腸桿菌SURE株，涂布含氨芐的LB固體培養(yǎng)板中，37℃培養(yǎng)過(guò)夜，每個(gè)片段分別挑取3個(gè)菌落PCR鑒定為陽(yáng)性的單克隆進(jìn)一步送華大基因測(cè)序鑒定。

2 結(jié)果與分析

2.1 病理剖解

病鴨除了外觀上的短喙長(zhǎng)舌之外，剖檢發(fā)現(xiàn)沒(méi)有明顯的病變。死鴨體內(nèi)充血嚴(yán)重，肝臟、脾臟、腎臟、肺臟、小腸、胰臟等出現(xiàn)腫大，如圖1所示。

2.2 PCR檢測(cè)

根據(jù)在VP3保守區(qū)域設(shè)計(jì)的特異性引物，以從病鴨和死鴨肝臟和脾臟提取的DNA為模板來(lái)擴(kuò)增目的片段。結(jié)果顯示，病鴨和死鴨體內(nèi)均能擴(kuò)增出約268 bp的特異性片段，如圖2所示。對(duì)PCR產(chǎn)物進(jìn)行液回收，測(cè)序后發(fā)現(xiàn)擴(kuò)增片段大小為268 bp，通過(guò)BLAST比對(duì)發(fā)現(xiàn)，兩者與具有短喙長(zhǎng)舌癥狀的山東分離株SDLC01的相似性為100%。

圖1 死鴨臨床癥狀及病理剖解Figure1 The clinical symptoms and pathological anatomy of dead duck

圖2 268 bp的PCR擴(kuò)增圖Figure2 PCR amplification of 268 bp

2.3 全基因組擴(kuò)增

利用P1～P55對(duì)特異性引物，采用PCR方法擴(kuò)增得到各目的片段，如圖3所示。將各個(gè)目的片段回收純化后連接載體測(cè)序，測(cè)序后的各個(gè)片段經(jīng)DNAstar 5.0軟件分析、拼接，得到該病毒的全基因組序列。本株NGPV（命名為SC16）全長(zhǎng)5109 nt，5′端回文序列（ITR）全長(zhǎng) 408 nt，3′端 ITR 全長(zhǎng) 407 nt，NS基因長(zhǎng)1884 nt，推測(cè)編碼627個(gè)氨基酸。VP基因長(zhǎng)2199 nt，推測(cè)編碼732個(gè)氨基酸。SC16株的全基因組序列已上傳至GenBank，登錄號(hào)為KY679174。

2.4 核苷酸序列同源性比對(duì)

為了探討本研究獲得的NGPV核苷酸序列與其他水禽細(xì)小病毒分離株的核苷酸序列的同源性，利用Sequence Demarcation Tool-version 1.2軟件（http://web.cbio.uct.ac.za/SDT）進(jìn)行了全基因組比對(duì)，根據(jù)核苷酸序列間的相似性構(gòu)建了二維彩色編碼矩陣圖，如圖4所示。同時(shí)該軟件分析結(jié)果顯示SC16株核苷酸序列與具有短喙長(zhǎng)舌癥狀的山東分離株SDLC01、QH15的同源性分別達(dá)到了99.6%和99.3%，而與同樣具有短喙長(zhǎng)舌癥狀的福建株M15的同源性僅為96.9%，與GPV分離株的同源性在95.4%～97%之間，而與MDPV分離株的同源性比較低，在81.9%～87%之間。涉及的水禽細(xì)小病毒分離株的詳細(xì)信息如表3所示。

圖3 SC16株的全基因組擴(kuò)增Figure3 Amplification of complete genome of SC16 strain

2.5 5′端回文序列（ITR）分析

利用Sequence Demarcation Tool-version 1.2軟件進(jìn)行了SC16株與其他水禽細(xì)小病毒分離株5′ITR核苷酸序列比對(duì)，如圖5所示。從圖可以看出，SC16株5′ITR的核苷酸序列仍然是與山東分離株SDLC01和QH15最接近。利用DNAMAN軟件分析發(fā)現(xiàn)，SC16株與SDLC01分離株的不同之處在于NGPV比SDLC01株多了3個(gè)長(zhǎng)14 nt的插入片段（TCCGGTGAC GCACA、GGTGACGTGTTTCC、GGATGTGCGTCACC），除此之外僅有一個(gè)堿基發(fā)生了替換。其中插入片段TCCGGTGACGCACA、GGATGTGCGTCACC在QH15的5′ITR核苷酸序列中也未曾發(fā)現(xiàn)。

表3 水禽細(xì)小病毒分離株信息Table3 Description of waterfowl parvovirus isolates involved in this study

圖4 全基因組比對(duì)的2維編碼矩陣圖Figure4 The 2D color coded matrix of complete genome alignment

圖5 5′ITR序列比對(duì)的2維編碼矩陣圖Figure5 The 2D color coded matrix of 5′ITR alignment

2.6 NS及VP蛋白的遺傳進(jìn)化分析

使用MEGA5.0軟件 Boot-stapped Neighbor-Joining方法，對(duì)SC16株和其他水禽細(xì)小病毒的NS蛋白、VP蛋白的氨基酸序列分別構(gòu)建系統(tǒng)進(jìn)化樹(shù)并進(jìn)行遺傳進(jìn)化分析，結(jié)果見(jiàn)圖6、圖7所示。由圖可知，不論是從NS蛋白的氨基酸序列還是從VP蛋白的氨基酸序列來(lái)看，SC16株都與山東分離株SDLC01和QH15在同一分支上，而與同樣具有短喙長(zhǎng)舌癥狀的福建分離株M15的親緣關(guān)系較遠(yuǎn)。并且SC16株同其他NGPV一樣，與GPV分離株在同一亞支，而與MDPV分離株相距較遠(yuǎn)。同時(shí)通過(guò)Sequence Demarcation Tool-version 1.2軟件進(jìn)行NS蛋白的氨基酸序列比對(duì)發(fā)現(xiàn)，SC16株的NS蛋白高度保守，其氨基酸序列與山東分離株QH15完全一致，與山東分離株SDLC01的同源性為99.5%。而VP蛋白的氨基酸序列與山東分離株SDLC01的同源性最高，為99.5%，與山東分離株QH15的同源性為98.9%。

3 討論

3.1 NGPV SC16株的全基因組擴(kuò)增

水禽細(xì)小病毒的基因組均為線性、單鏈DNA，大小約為5 kb。雖然基因組較小，但兩端的末端倒置重復(fù)序列（ITR）的結(jié)構(gòu)非常復(fù)雜，是擴(kuò)增全基因組的難點(diǎn)。Zádori Z等[9]在測(cè)定GPV和MDPV的全基因組時(shí)發(fā)現(xiàn)，GPV基因組的兩端存在著由444個(gè)核苷酸形成的高GC含量的、相同的ITR結(jié)構(gòu)：ITR的前407個(gè)核苷酸可形成“U”型發(fā)夾結(jié)構(gòu)，發(fā)夾結(jié)構(gòu)中有181個(gè)核苷酸形成無(wú)錯(cuò)配“莖”部和一個(gè)44 nt的“泡”區(qū)，外形如“箭”狀，在“箭頭”有一個(gè)SphI限制性內(nèi)切酶位點(diǎn)，形成2分對(duì)稱中心。MDPV的基因組結(jié)構(gòu)類似于GPV。ITR是發(fā)生高頻重組的熱點(diǎn)區(qū)，在細(xì)菌細(xì)胞基因組缺乏主要的重組基因時(shí)，“泡”區(qū)都可能發(fā)生“翻”“轉(zhuǎn)”互變現(xiàn)象，因此ITR完整的“泡”區(qū)克隆到E.coli中，會(huì)發(fā)生“翻”“轉(zhuǎn)”方向的互變，造成測(cè)序困難。并且ITR形成的發(fā)夾結(jié)構(gòu)使得聚合酶難以通過(guò)，PCR方法一次性擴(kuò)增基因組難以實(shí)現(xiàn)[10]。因此，在設(shè)計(jì)引物時(shí)需要考慮將ITR結(jié)構(gòu)剖開(kāi)，而不是一次性擴(kuò)增完整的ITR序列。本實(shí)驗(yàn)參考GenBank上發(fā)表的具有短喙長(zhǎng)舌癥狀的NGPV的全基因組序列及基因組結(jié)構(gòu)特征，通過(guò)核苷酸比對(duì)后，選擇基因組中相對(duì)較保守的區(qū)域設(shè)計(jì)了5對(duì)兼并引物來(lái)分段擴(kuò)增，并且使擴(kuò)增相鄰片段有重疊，PCR擴(kuò)增后分別得到了長(zhǎng)204 bp、1397 bp、1489 bp、2419 bp、1461 bp的片段，經(jīng)過(guò)軟件拼接后獲得NGPV四川株的全基因組，全長(zhǎng)為5109 nt,豐富了GenBank數(shù)據(jù)庫(kù)，為NGPV的分子生物學(xué)研究奠定一定的基礎(chǔ)。

圖6 基于NGPV NS蛋白的系統(tǒng)進(jìn)化樹(shù)Fig.6 Phylogenetic tree of NGPV NS protein

圖7 基于NGPV VP蛋白的系統(tǒng)進(jìn)化樹(shù)Figure7 Phylogenetic tree of NGPV VP protein

3.2 NGPV SC16株的遺傳變異情況

NGPV SC16株的基因組全長(zhǎng)為5109 nt，其5′端ITR由408個(gè)核苷酸組成，比NGPV的山東分離株SDLC01株多了3個(gè)長(zhǎng)14 nt的插入片段。細(xì)小病毒基因組兩端的ITR在病毒DNA復(fù)制和包裝過(guò)程中起重要作用[11]。有推測(cè)稱ITR的長(zhǎng)短和病毒毒力有關(guān)，且ITR較短的多為低致病力毒株或疫苗株[12]。但此推測(cè)卻與萬(wàn)春和等[13]的研究結(jié)果不盡一致。目前，關(guān)于水禽細(xì)小病毒全基因組特征的研究特別是NGPV全基因組的研究有限，所以此推測(cè)仍需進(jìn)一步研究驗(yàn)證。

不論從全基因組比對(duì)情況還是基于NGPV NS和VP蛋白的遺傳進(jìn)化分析來(lái)看，SC16株都與山東分離株SDLC01和QH15親緣關(guān)系最近，相似性都高達(dá)99%左右。通過(guò)調(diào)查發(fā)現(xiàn)，送檢的病死鴨的鴨苗均購(gòu)買自山東省。Chen H.等[14]發(fā)現(xiàn)NGPV在自然狀況下有垂直傳播的可能，可以從種鴨傳至鴨胚和雛鴨。因此，我們推測(cè)，本研究分離的NGPV四川株很可能來(lái)自山東省，但隨著時(shí)間、空間等的改變，病毒不斷發(fā)生遺傳演化，且VP基因上的遺傳進(jìn)化速度較NS基因更快，這與GPV分離株GoPV[13]的遺傳演變情況相似，有報(bào)道稱這可能與宿主免疫系統(tǒng)的壓力有關(guān)[15]。NS蛋白主要參與病毒復(fù)制以及基因表達(dá)，這部分結(jié)構(gòu)的穩(wěn)定可以很好地維持DNA復(fù)制、轉(zhuǎn)錄等活動(dòng)的穩(wěn)定進(jìn)行。而VP蛋白的作用不僅在于病毒顆粒衣殼蛋白的合成及子代病毒的裝配釋放，而更主要的作用是作為配體與相關(guān)的宿主細(xì)胞表面的受體結(jié)合，是病毒感染過(guò)程中最為重要的一個(gè)因素。因此，病毒為了成功感染宿主細(xì)胞，VP基因可能發(fā)生了一些利于自身生存的變異，從而逃避宿主免疫系統(tǒng)的攻擊。

本研究成功獲得了NGPV四川株的全基因組序列，比較了NGPV四川株與山東分離株及與其他水禽細(xì)小病毒分離株的共性與特性，解釋了引起鴨短喙長(zhǎng)舌綜合征的病原在四川的發(fā)生及變異情況，豐富了新型鵝細(xì)小病毒的分子生物學(xué)資料，為該病的流行病學(xué)調(diào)查及防控奠定了理論基礎(chǔ)。

［1］竇硯國(guó)，陳浩，鄭肖強(qiáng)，等.鴨細(xì)小病毒套式PCR檢測(cè)方法的建立與應(yīng)用［J］.中國(guó)獸醫(yī)學(xué)報(bào)，2016，36（10）：1706-1709.

［2］黃宇翔，劉力威，李洪彬，等.鵝細(xì)小病毒病、副黏病毒病的流行病學(xué)調(diào)查及病原分離鑒定［J］.中國(guó)家禽，2013，35（2）：22-25.

［3］CHEN S，WANG S，CHENG X，et al.Isolation and characterization of a distinct duck-origin goose parvovirus causing an outbreak of duckling short beak and dwarfism syndrome in China［J］.Archives of Virology，2016，161（9）：1-10.

［4］YU K，MA X，SHENG Z，et al.Identification of goose-origin parvovirus as a cause of newly emerging beak atrophy and dwarfism syndrome in ducklings［J］.Journal of Clinical Microbiology，2016，54（8）：JCM.03244-15.

［5］LI C，LI Q，CHEN Z，et al.Novel duck parvovirus identified in Cherry Valley ducks(Anas platyrhynchos domesticus)，China［J］.Infection Genetics&Evolution，2016（44）：278-280.

［6］殷冬冬，唐井玉，王瑞，等.新型鴨源細(xì)小病毒安徽株AHD15株的分離鑒定［J］.中國(guó)預(yù)防獸醫(yī)學(xué)報(bào)，2016，38（7）：528-531.

［7］袁萬(wàn)哲，崔元，經(jīng)美，等.鴨短喙——侏儒綜合征研究初報(bào)［J］.中國(guó)動(dòng)物保健，2016，18（5）：73-74.

［8］CHEN H，DOU Y，TANG Y，et al.Isolation and Genomic Characterization of a Duck-Origin GPV-Related Parvovirus from Cherry Valley Ducklings in China［J］.Plos One，2015，10（10）：e0140284.

［9］ZA'DORI Z，STEFANCSIK R，RAUCH T，et al.Analysis of the complete nucleotide sequences of goose and muscovy duck pervoviruses indicates common ancestral origin with adenoassociated virus 2［J］.Virology，1995，212（2）：562-573.

［10］王建業(yè)，黃鈺，段進(jìn)坤，等.鵝細(xì)小病毒強(qiáng)毒株YZ99-6感染性質(zhì)粒的構(gòu)建及其在鵝胚中的拯救［J］.中國(guó)預(yù)防獸醫(yī)學(xué)報(bào)，2016，38（3）：185-189.

［11］QIU J，CHENG F，YOTO Y，et al.The expression strategy of goose parvovirus exhibits features of both the Dependovirus and Parvovirus genera［J］.Journal of Virology，2005，79（17）：11035.

［12］SHIEN J H，WANG Y S，CHEN C H，et al.Identification of sequence changes in live attenuated goose parvovirus vaccine strains developed in Asia and Europe ［J］.Avian Pathology，2008，37（5）：499-505.

［13］萬(wàn)春和，朱海俠，黃瑜，等.一株鵝細(xì)小病毒全基因特征分析［J］.中國(guó)動(dòng)物傳染病學(xué)報(bào)，2011，19（4）：19-24.

［14］CHEN H，TANG Y，DOU Y，et al.Evidence for vertical transmission of novel duck-origin goose parvovirus-related parvovirus［J］.Transboundary&Emerging Diseases，2016，63（3）：243-247.

［15］LE G R G，JESTIN V，CHAGNAUD P，et al.Expression of muscovy duck parvovirus capsid proteins(VP2 and VP3)in a baculovirus expression system and demonstration of immunity induced by the recombinant proteins ［J］.Journal of General Virology，1996，77（9）：2159.

Amplification and Phylogenetic Analysis of Novel Goose Parvovirus Sichuan Strain

GE Zhi-qi1，LIU Ya-gang1*，CHEN Shun2，3，4，ZHU De-kang2，3，4，LIU Ma-feng2，3，4，WANG Ming-shu2，3，4，JIA Ren-yong2，3，4，SUN Kun-feng2，3，4，YANG Qiao2，3，4，WU Ying2，3，4，ZHAO Xin-xin2，3，4，CHENG An-chun2，3，4*
（ 1.College of Life Science and Technology，Southwest University for Nationalities，Chengdu 610041，China；2.Institute of Preventive Veterinary Medicine，Sichuan Agricultural University，Chengdu 611130，China；3.Avian Disease Research Center，College of Veterinary Medicine of Sichuan Agricultural University，Chengdu 611130，China；4.Key Laboratory of Animal Disease and Human Health of Sichuan Province，Sichuan Agricultural University，Chengdu 611130，China)

【Objective】The aim of this study was to learn the molecular biological characteristics and genetic variation of novel goose parvovirus（NGPV）in Sichuan area.【Methods】5 pairs of degenerate primers were designed to amplify the whole genome sequence of NGPV Sichuan strain using PCR method，according to the genome sequences of the NGPV deposited in GenBank.Then，the study usd DNAstar and other bioinformatics softwares to analyze the molecular characteristics and genetic evolution of the obtained sequences.【Results】The genomic sequence of Novel goose parvovirus from Sichuan(named SC16)was 5109 nucleotides(nt)in length.The inverted terminal repeats(ITRs)at the 5′and 3′end of the genome were 408 nt and 407 nt respectively.The NS protein was encoded by 1884 nt and the VP protein was encoded by 2199 nt.Sequence alignment and phylogenetic analysis demonstrated that the nucleotide identities of SC16 genome and Shandong isolates SDLC01(KT343253.1)andQH15(KT751090.1)was up to 99%，showing a very close genetic relationship with them.The phylogenetic tree analyses of amino acids of NS and VP protein were consistent with this result.However，SC16 still had some genetic variations：compared with SDLC01 and QH15，SC16 had two or three 14-nt extra fragments at the 5′end of the genome.【Conclusion】These findings will enrich the molecular biology information of NGPV.

novel goose parvovirus；beak atrophy and dwarfism syndrome of duck；whole genome amplification；phylogenetic analysis

S834+8

A

1000-2650（2017）04-0574-07

10.16036/j.issn.1000-2650.2017.04.018

2017-05-02

“十二五”國(guó)家科技支撐計(jì)劃（2015BAD12B05）；國(guó)家現(xiàn)代農(nóng)業(yè)（水禽）產(chǎn)業(yè)技術(shù)體系專項(xiàng)（CARS-43-8）；國(guó)家重點(diǎn)研發(fā)計(jì)劃“家禽重要疫病診斷與檢測(cè)新技術(shù)研究”（2016YFD0500800）；四川省重點(diǎn)實(shí)驗(yàn)室專項(xiàng)（2016JPT0004）。

葛志琪，碩士研究生。*責(zé)任作者：劉亞剛，教授，博士，主要從事獸醫(yī)免疫學(xué)的研究，E-mail：LYG_068@163.com；程安春，教授，博士，主要從事禽病學(xué)的研究，E-mail：chenganchun@vip.163.com。

（本文審稿：陳宗艷；責(zé)任編輯：秦碧雯；英文編輯：劉益平）