鄭龍龍 ，朱玲云，宋春蓮，高利波，畢峻龍，徐 英

(1.云南農(nóng)業(yè)大學(xué)動物科學(xué)技術(shù)學(xué)院，云南昆明 650201；2.楚雄州動物疫病預(yù)防控制中心，云南楚雄 675000；3.云南農(nóng)業(yè)職業(yè)技術(shù)學(xué)院，云南昆明 650031 )

云南省犬細小病毒VP2基因檢測分析

鄭龍龍1，朱玲云1，宋春蓮1，高利波1，畢峻龍2，徐 英3*

(1.云南農(nóng)業(yè)大學(xué)動物科學(xué)技術(shù)學(xué)院，云南昆明 650201；2.楚雄州動物疫病預(yù)防控制中心，云南楚雄 675000；3.云南農(nóng)業(yè)職業(yè)技術(shù)學(xué)院，云南昆明 650031 )

為了解云南省犬細小病毒VP2基因分子流行變異情況，對來源于云南省7個地級市13個犬細小病毒樣品進行了VP2基因的擴增和序列分析，并與相關(guān)參考序列進行了遺傳進化比對分析。結(jié)果表明，13個樣品中9個基因型為CPV-2a，3個基因型為CPV-2b，1個基因型為CPV-2。氨基酸序列變異主要集中在297位-518位氨基酸處，其中324位和440位氨基酸變異在犬細小病毒進化過程中扮演著重要角色。遺傳進化分析顯示10個樣品序列形成了一個相對于國內(nèi)外參考序列獨立的遺傳進化分支，提示云南省犬細小病毒以CPV-2a 基因型流行為主。

犬細小病毒；VP2基因；基因型；遺傳進化

犬細小病毒病(Canine parvovirus disease，CPVD)是由犬細小病毒(Canine parvovirus disease，CPV)引起的一種急性、烈性傳染病，是當(dāng)前對犬類危害最大的傳染病之一。該病能夠引起犬類白細胞減少、劇烈嘔吐、出血性腸炎等多種癥狀。犬細小病毒屬于細小病毒科細小病毒屬成員，編碼2個結(jié)構(gòu)蛋白(VP1和VP2)，其中VP2蛋白可誘導(dǎo)機體產(chǎn)生中和抗體，是保護性顆?？乖P2蛋白中一些具有生物學(xué)意義的主要抗原位點的變異導(dǎo)致了病毒抗原特性和宿主范圍的改變。CPV-2于20世紀70年代后期被發(fā)現(xiàn)不久之后，CPV-2逐漸被新的抗原變異株CPV2-a和CPV2-b取代[1-3]。近年來，另外一種新的變異樣品序列CPV2-c被證實在歐美國家廣泛存在[4]，國內(nèi)學(xué)者于2010年發(fā)現(xiàn)并證實CPV-2C變異樣品序列也存在于我國[5-6]。隨著頻繁的貿(mào)易交流，犬細小病毒不斷出現(xiàn)新的變異樣品序列，流行病學(xué)呈現(xiàn)復(fù)雜態(tài)勢，導(dǎo)致免疫失敗現(xiàn)象時有發(fā)生，為犬細小病毒病的預(yù)防和控制帶來了很大困難。由于目前幾乎沒有關(guān)于犬細小病毒云南株流行樣品序列基因型的相關(guān)報道。本研究通過對細小病毒云南分離株VP2基因進行遺傳進化和氨基酸位點突變的分析，調(diào)查云南省犬細小病毒分子流行變異情況，為該病的防控提供一定的理論參考依據(jù)。

1 材料與方法

1.1 材料

1.1.1 病料來源 采自云南省不同地區(qū)發(fā)生疑似犬細小病毒犬場及寵物醫(yī)院病犬糞便,樣品相關(guān)信息見表1。

表1 樣品相關(guān)信息Table 1 Related information of sample gene sequence

1.1.2 主要試劑 DNA聚合酶、pMDTM18-T載體、限制性內(nèi)切酶BsrⅠ、HindⅢ等分子生物學(xué)試劑為寶生物工程(大連)有限公司產(chǎn)品；其他試劑由云南農(nóng)業(yè)大學(xué)動物科學(xué)技術(shù)學(xué)院動物傳染病實驗室提供。

1.1.3 儀器設(shè)備 移液器、超凈工作臺、高速冷凍離心機、恒溫水浴鍋、PCR儀、電泳槽、凝膠電泳成像儀等設(shè)備由云南農(nóng)業(yè)大學(xué)動物科學(xué)技術(shù)學(xué)院動物傳染病實驗室提供。

1.2 方法

1.2.1 引物設(shè)計與合成 VP2序列全長1 755 bp，編碼584個氨基酸，參考 GenBank中犬細小病毒VP2基因參考序列，即M19296、M24003、EU659121、GU569947和GU569948等序列，應(yīng)用Oligo6.0軟件，設(shè)計用于擴增CPV VP2基因的一對引物其中上游引物：5′-ACC ATG TAG ACT AAC ACA TAC-3′，下游引物：5′-ACA TAT TAG CAG CAT CAG-3′，擴增片段長度為741 bp。

1.2.2 VP2片段的擴增及克隆 使用DNA提取試劑盒提取組織樣品總DNA，根據(jù)PCR試劑說明書優(yōu)化反應(yīng)體系及條件，其中10×PCR buffer 2.5 μL，dNTP Mixture 1 μL，上游和下游引物分別0.5 μL，TaqDNA聚合酶0.5 μL，模板1 μL，加水補足25 μL。進行PCR擴增后，將擴增產(chǎn)物和上樣緩沖液混合后，用15 g/L的瓊脂糖凝膠進行電泳，20 min后觀察結(jié)果。 將陽性PCR產(chǎn)物經(jīng)切膠回收后與pMD18-T Vector進行連接反應(yīng)，連接產(chǎn)物轉(zhuǎn)化進入DH5α感受態(tài)細胞，挑取單克隆菌株后，擴大培養(yǎng)，提取質(zhì)粒，并使用限制性內(nèi)切酶進行酶切鑒定。

1.2.3 序列測定與分析 將鑒定為陽性的菌種進行測序分析，使用DNA Star和MEGA 4.0等生物學(xué)軟件分別對所獲得基因序列進行位點突變和遺傳進化分析，并與國內(nèi)外參考序列進行比較并繪制遺傳進化樹(表2)。

表2 分析所用的CPV參考序列及其來源Table 2 CPV strains used in this study and their origins

2 結(jié)果

2.1 VP2基因的擴增與分析

樣品DNA經(jīng)PCR擴增后，經(jīng)15 g/L的瓊脂糖凝膠進行電泳檢測，經(jīng)凝膠成像系統(tǒng)觀察后，發(fā)現(xiàn)陽性對照及樣品均在DNA Maker DL 700 bp附近處有一條帶，與預(yù)期741 bp相符(圖1)。

圖1 VP2基因PCR擴增結(jié)果Fig.1 PCR results of VP2 gene

2.2 VP2基因重組質(zhì)粒酶切鑒定

通過DNA Star 6.0分析，在pMDTM18-T載體上選擇HindⅢ，在目的片段第572位選擇BsrⅠ進行雙酶切，預(yù)期目的片段為2 834 bp和599 bp或3 237 bp和196 bp。將切膠回收純化的CPV VP2基因的PCR擴增產(chǎn)物(741 bp)，連接、轉(zhuǎn)化到感受態(tài)細胞DH5α中，經(jīng)過培養(yǎng)后提取重組質(zhì)粒進行酶切鑒定。重組質(zhì)粒經(jīng)過HindⅢ和BsrⅠ雙酶切后得到兩個大小約為2 834 bp和599 bp的片段，酶切結(jié)果與預(yù)期相符，表明成功構(gòu)建了重組質(zhì)粒(圖2)。

圖2 重組質(zhì)粒酶切鑒定Fig.2 Identification of recombinant plasmids digested with BsrⅠ and HindⅢ

2.3 VP2序列測定與分析

測序結(jié)果顯示，本試驗獲得的13個VP2基因包含著741個堿基，編碼247個氨基酸。對所獲得的基因序列使用生物學(xué)軟件DNA Star與其他參考序列進行核苷酸同源分析(表3)。結(jié)果表明，于2012年獲得的5株VP2基因的核苷酸同源性為99.5%～100%；于2013年獲得的8株VP2基因的核苷酸同源性為96.8%～100%；2012年序列和2013年樣品序列P2基因的核苷酸同源性98.8%～99.2%；2012序列與參考序列的核苷酸同源性為98.7%～99.5%，與疫苗序列EF112447核苷酸同源性為98.5%～98.8%；2013年序列與參考序列的核苷酸同源性為98.8%～99.75%，與疫苗序列EF112447核苷酸同源性為98.7%～99.6%。與2012年獲得的5株VP2基因的蛋白同源性為99.6%～100%；與2013年獲得的8株VP2基因的蛋白同源性為97.2%～100%；2012年序列和2013年序列VP2基因的蛋白同源性為98.8%～99.2%；2012年序列與參考序列的蛋白同源性為97.6%～99.2%，與疫苗序列EF112447蛋白同源性為97.2%～97.6%。2013年序列與參考序列的蛋白同源性為97.2%～100%，與疫苗EF112447蛋白同源性為96.8%～99.2%。

2.4 VP2基因氨基酸突變分析

將實驗獲得的13個VP2氨基酸序列和部分歐美、亞洲及我國樣品序列的相應(yīng)氨基酸序列進行比較(表4)。云南VP2氨基酸序列變異主要集中在297位～518位氨基酸處，其中297、300、305、324、375、386、426、427、440和518位點存在氨基酸替換。2012YN-1、2012YN-2、2012YN-3、2012YN-4、2012YN-5、2013YN-6、2013YN-10、2013YN-11、2013YN-12為CPV-2a型，2013YN-7株為CPV2型，2013YN-8和2013YN-9株為CPV-2b型，CPV-2a型是云南省2012年-2013年的主要流行序列。

2.5 VP2基因遺傳進化樹分析

將所得到的13個2012年-2013年云南流行序列VP2基因與參考序列進行遺傳分析后繪制進化樹，如圖3所示。本試驗所得病毒序列2013YN-7、2013YN-8和2013YN-9序列和其他云南毒株序列親源關(guān)系較遠；2012YN-1、2012YN-2、2012YN-3、2012YN-4、2012YN-5 2013YN-6、2013YN-10、2013YN-11、2013YN-12、2013YN-13相對于其他參考序列形成了一個獨立分支。

表3 不同序列VP2序列氨基酸同源性Table 3 Sequence similarities of vp2 gene from difficult strains on nucleotide acid and amino acid levels %

表4 VP2基因片段中發(fā)生置換的氨基酸Table 4 The substitutions of the amino acids in the VP2 gene

圖3 CPV VP2基因遺傳進化分析Fig.3 Phylogenetic relationship of CPV VP2 gene

3 討論

VP2基因全長1 755 bp，編碼584個氨基酸，它的關(guān)鍵位點的變異決定著犬細小病毒基因型。CPV-2a和CPV-2b 基因型序列VP2蛋白的在87、300、305位點氨基酸與CPV-2基因型不同[7]。CPV-2a和CPV-2b序列只在VP2基因的第426位氨基酸不同，其中426位點Asp是CPV-2b型序列特有的。在CPV-2a、CPV-2b基因型的基礎(chǔ)上，VP2蛋白的第300位點發(fā)生了G-D替換，由此產(chǎn)生了CPV-2c(a)和CPV-2c(b)兩個新的突變基因型[8]，因此，本研究擴增的741 bp VP2基因片段主要包括了決定犬細小病毒基因型的關(guān)鍵位點，可以根據(jù)位點的變異來判斷云南流行株的基因型。

盡管CPV廣泛存在于世界，它的基因型在不同的國家和地區(qū)具有很大差異。近年來，在歐洲的許多國家，如英國、德國和意大利，CPV-2b型 和 CPV-2c型逐漸超過CPV-2a型，成為CPV優(yōu)勢流行基因型；CPV-2b型和 CPV-2c型在北美地區(qū)同許多歐洲國家一樣是CPV的優(yōu)勢基因型；而在南美，以CPV-2c型為CPV主要流行基因型[9-10]。CPV-2a是亞洲和澳大利亞的CPV主要優(yōu)勢基因型[11-13]。Zhang R等[14-15]通過分析2006年-2009年的CPV序列，證明我國CPV主要流行CPV-2a基因型。將云南13個CPV VP2的氨基酸序列使用生物學(xué)軟件進行氨基酸位點變異分析后，結(jié)果表明2012YN-1、2012YN-2、2012YN-3、2012YN-4、2012YN-5、2013YN-6、2013YN-10、2013YN-11、2013YN-12基因型為CPV2a，2013YN-7基因型為CPV2，2013YN-8和2013YN-9基因型是CPV2b，CPV-2a

流行程度超過CPV-2b，是當(dāng)前云南省CPV主要流行基因型。本試驗所獲得的VP2氨基酸序列變異主要集中在VP2基因297位～518位氨基酸處，其中97、300、305、324、375、386、426、427、440、518位點存在氨基酸替換，2012YN-1、2012YN-2、2012YN-3、2012YN-4、2012YN-5、2013YN-6、2013YN-10、2013YN-11、2013YN-12、2013YN-13的324位氨基酸由Tyr變成Ile位點，同324位點臨近的323位點對病毒具有重要的生物意義，該位點影響著病毒和犬運載蛋白受體的結(jié)合及感染宿主的范圍[16]，2012YN-1、2012YN-3、2013YN-6、2013YN-11、2013YN-12、2013YN-13的440氨基酸位點由Thr變成Ala，相關(guān)研究表明了該位點在細小病毒進化過程中扮演著重要角色。所得到的13個2012年-2013年云南流行VP2基因序列與國內(nèi)外VP2參考序列繪制遺傳進化樹，云南流行毒2013YN-9和2013YN-8樣品序列和其他云南毒親源關(guān)系較遠，2012YN-1、2012YN-2、2012YN-3、2012YN-4、2012YN-5、2013YN-6、2013YN-10、2013YN-11、2013YN-12、2013YN-13形成了一個獨立分支，提示：犬細小病毒云南流行基因型在云南本地的流行過程中形成了相對獨立于其他參考序列的云南分支。

本研究通過對云南省犬細小病毒序列進行基因型和遺傳進行分析，表明了當(dāng)前云南省流行的病毒基因型主要有 CPV-2a、CPV-2b和CPV-2，其中以CPV-2a占主要地位，初步了解云南省犬細小病毒VP2基因序列遺傳進化趨勢，為該病的預(yù)防和控制提供了一定的理論基礎(chǔ)和參考依據(jù)。

[1] Marina G C, Maximiliano W, Lorena B, et al. Study of canine parvovirus evolution: comparative analysis of full-length VP2 gene sequences from Argentina and international field strains[J].Virus Genes,2012，44:32-39.

[2] Decaro N, Desario C, Amorisco F, et al. Detection of a canine parvovirus type 2c with a non-coding mutation and its implications for molecular characterisation[J].Vet J,2013,196:555-557.

[3] 鞠厚斌，周錦萍，劉 健，等.上海地區(qū)犬細小病毒的分離及VP2基因序列分析[J].中國動物傳染病學(xué)報,2013,21(1):37-41.

[4] Calderon M G, Romanutti C, Antuono D, et al. Evolution of canine parvovirus in Argentina between years 2003 and 2010: CPV2c has become the predominant variant affecting the domestic dog population[J].Virus Res,2011，157:106-110.

[5] 張仁舟，楊松濤，馮 昊，等.中國國內(nèi)首次檢測到犬細小病毒CPV-2c[J].中國病原生物學(xué)雜志,2010,54(2):246-249.

[6] 姜燕霞，張淑玲，王金良，等.犬細小病毒JN株的分離鑒定及VP2基因序列分析[J].動物醫(yī)學(xué)進展,2013,34(5):83-86.

[7] Hoelzer K, Shackelton L A, Parrish C R, et al. Phylogenetic analysis reveals the emergence, evolution and dispersal of carnivore parvoviruses[J]. J Gen Virol,2008,89(9):2280-2289.

[8] Martella V, Cavalli A, Pratelli A, et al. A canine parvovirus mutant is spreading in Italy[J].J Clin Microbiol,2004，42:1333-1336.

[9] Hong C, Decaro N, Desario C, et al. Occurrence of canine parvovirus type 2c in the United States[J].J Vet Diagn Invest,2007，19:535-539.

[10] Perez R, Francia L, Romero V, et al. First detection of canine parvovirus type 2c in South America[J].Vet Microbiol,2007，124:147-152.

[11] Kang B K, Song D S, Lee C S, et al. Prevalence and genetic characterization of canine parvoviruses in Korea[J].Virus Genes,2008，36:127-133.

[12] Meers J, Kyaw-Tanner M, Bensink Z, et al. Genetic analysis of canine parvovirus from dogs in Australia[J].Aust Vet J,2007，85:392-396.

[13] Nandi S, Chidri S, Kumar M, et al. Occurrence of canine parvovirus type 2c in the dogs with haemorrhagic enteritis in India[J].Res Vet Sci,2010，88:169-171.

[14] 李世靜，嵇辛勤，主 性，等.犬細小病毒VP2基因測序分析[J].動物醫(yī)學(xué)進展,2013,34(6):96-100.

[15] Zhang R, Yang S, Zhang W, et al. Phylogenetic analysis of the VP2 gene of canine parvoviruses circulating in China[J].Virus Genes,2010，40:397-402.

[16] Jaime Aldaz, Juan García-Díaz, Lucía Calleros, et al. High local genetic diversity of canine parvovirus from Ecuador[J].Vet Microbiol,2013,166(1-2):214-219.

Detection and Analysis of VP2 Gene of Canine Parvovirus in Yunnan Province

ZHENG Long-long1,ZHU Ling-yun1,SONG Chun-lian1,GAO Li-bo1,BI Jun-long2,XU Ying3

(1.CollegeofAnimalScienceandTechnology,YunnanAgriculturalUniversity,Kunming,Yunnan,650201,China; 2.ChuxiongProvinceCenterforDiseaseControlandPrevention,Chuxiong,Yunnan,675000,China; 3.YunnanVocationalandTechnicalCollegeofAgriculture,Kunming,Yunnan,650031,China)

In order to investigate the prevalence and variation of VP2 gene of CPV in Yunnan province, the gene amplification and sequence analysis were performed for thirteen strains of canine parvovirus, comprising of 9 strains of CPV-2a, 3 strains of CPV-2b and 1 strain of CPV-2. It was shown that most of the amino acid variations occurred within the region of position 297-518. A high level of substitution at 324 and 440 sites was associated with the evolution of antigenic variants in circulating parvovirus types. Based on the phylogenetic analysis on the amino acid level of the VP2 proteins, the thirteen strains of CPVs were grouped into 3 clusters. Nine of them belonged to CPV-2a, three belonged to CPV-2b, and one belonged to CPV-2. The result showed that CPV-2a was the main prevalent strain in Yunnan province and formed the independent branch.

Canine parvovirus;VP2 gene;typing;genetic evolution

2014-09-28

云南省公益性關(guān)鍵技術(shù)研究開發(fā)計劃子項目(2013CH001)

鄭龍龍(1988-)，男，山西運城人，碩士研究生，主要從事動物傳染病學(xué)與分子流行病學(xué)研究。*

S852.659.2

:A

:1007-5038(2015)01-0066-05