山東地區水貂腸炎病毒VP2基因的克隆和測序分析

陳 童， 楊瑞梅， 單 虎

(青島農業大學山東省預防獸醫學重點實驗室， 山東青島266109)

山東地區水貂腸炎病毒VP2基因的克隆和測序分析

陳 童， 楊瑞梅， 單 虎

(青島農業大學山東省預防獸醫學重點實驗室， 山東青島266109)

本研究對濰坊、煙臺、威海和青島地區的疑似水貂病毒性腸炎病毒(MEV)糞便樣品進行檢測，并在9份樣品QD1309、QD1407、WH1407、WH1408、WH1508、WH1509、WF1310、WF1408、YT1408中擴增出VP2基因，和GenBank上已經公布的4株MEV參考株的VP2基因進行序列分析，結果顯示，13份MEV VP2基因的核苷酸和氨基酸均有較高的同源性，分別為97.5%～99.8%和96.9%～99.7%。系統發育進化樹表明，QD1309、YT1408、MEV-Beregovoj-Bincentr(MEV-BB)、MEV-e屬于同一個組群，不同毒株主要在第5、62、87、101、232、236、279、300、534位氨基酸發生了替換。本研究對水貂細小病毒流行株的主要衣殼蛋白VP2的編碼蛋白進行遺傳變異分析，為更好地預防和控制貂源細小病毒提供基礎。

水貂腸炎病毒； VP2基因； 遺傳變異分析

水貂病毒性腸炎病毒(minkenteritisvirus，MEV)引起的主要臨床癥狀是劇烈腹瀉，腸道出血，腸內容物呈煤焦油狀。MEV為單股負鏈DNA病毒，屬于細小病毒科細小病毒屬，基因組全長約為5 kb，其基因組含有2個開放閱讀框，其中一個調控病毒復制的非結構蛋白NS1和NS2，另一個組呈衣殼結構蛋白VP1和VP2[1]。VP2蛋白是細小病毒的主要衣殼蛋白，其對病毒的宿主范圍、組織嗜性、抗原性、血凝性等起重要作用。該蛋白中某幾個甚至一個氨基酸殘基的突變，都可能導致整個病毒的生物學特性發生改變[2-3]。本研究旨在初步調查山東省主要地區的MEV流行情況,系統地分析MEV VP2基因的變異狀態并進行分析，以探索山東省MEV分子流行病學趨勢[4]。

1 材料與方法

1.1 病料采集 病料取自濰坊、威海、煙臺、青島地區的某些水貂養殖場，臨床排血便，剖檢為腸內容物呈煤焦油狀，小腸黏膜充血。

1.2 主要試劑 ExTaqDNA聚合酶、dNTP、DL-2 000 DNA Marker、DL-5000 DNA Marker、pMD18-T、DH5α、EcoRⅠ、SalⅠ限制性內切酶，均購自TaKaRa公司；膠回收試劑盒和質粒提取試劑盒，購自OMEGA公司。

1.3 樣品處理 用棉棒蘸取腸內容物加入PBS緩沖液中，顛倒混勻5～8次，7 000r/min離心10 min，取上清，-70 ℃反復凍融3次。12 000r/min離心10 min，取上清。

1.4 DNA的提取 取500 μL病毒上清液，加入500 μL DNAiso Reagant，顛倒混勻6～8次，靜置10 min。10 000 r/min離心10 min。取上清800 μL，加入400 μL無水乙醇，顛倒混勻6～8次，4 000 r/min離心2 min，棄上清。加入1 mL 75%乙醇，顛倒混勻6～8次，12 000 r/min離心5 min，棄上清。用吸水紙把底部的75%乙醇吸干，用20 μL ddH2O吹打混勻制備成PCR擴增的模板，-20 ℃保存備用。

1.5 MEV的PCR鑒定 參考GenBank中公布的MEV NS基因序列，利用Primer 5.0軟件設計一對引物，上游引物為：5′-GTAAGCTTCCAGGAGACTTT-3′,下游引物為：5′-GTAAGCTTCGTCGTGTTCTT-3′，擴增片段長度671 bp。引物由北京賽百盛基因技術有限公司合成。PCR反應體系為25 μL：模板1.5 μL，10× ExTaqBuffer 2.5 μL,上下游引物各1 μL，dNTP 2 μL，ExTaq酶0.2 μL，滅菌純水補足至25 μL。PCR擴增條件為：94 ℃預變性4 min；94 ℃ 變形40 s；52 ℃退火30 s。72 ℃延伸1 min，30個循環；72 ℃延伸10 min。經PCR擴增檢測樣品，同時設定陽、陰性對照，瓊脂糖凝膠電泳，驗證為陽性的樣品擴增VP2基因全長。

1.6 VP2基因的克隆和變異分析 根據GenBank上公布的MEV VP2基因序列設計引物，上游引物為：5′-CGGGATCCATGAGTGATGGAGCAGTTCAA-3′，下游引物為：5′-CCGCTCGAGATATAATTTTCTAGGTGCT-3′ ，擴增片段長度為1 755 bp。PCR擴增條件為：94 ℃預變性2 min；94 ℃變性30 s；55 ℃退火30 s。72 ℃延伸40 s，35個循環；72 ℃延伸3 min。將陽性產物連接到pMD18-T載體上，轉化到DH5α感受態，經PCR和雙酶切(EcoRⅠ和SalⅠ)鑒定后送到北京華大基因科技服務有限公司測序。使用DNAStar和Mega 5.0軟件進行同源性分析及構建進化樹。

2 結果與分析

2.1 PCR鑒定 取5 μL PCR擴增產物經1%瓊脂糖凝膠電泳檢測，目的片段約為671 bp，其大小與預期相符表明檢測的病毒為MEV(圖1)。

圖1 PCR鑒定結果

M：DL-2 000 DNA Marker；1：QD1309；2：QD1407；3：WH1407；4：WH1408；5：WF1408 ；6：YT1408

2.2 VP2基因擴增和雙酶切鑒定 擴增VP2基因，瓊脂糖凝膠電泳出現與預期大小一致目的條帶，即1 755 bp。PCR產物經膠回收、連接、轉化后，菌液PCR鑒定為陽性的質粒，用EcoRⅠ和SalⅠ雙酶切鑒定大小正確，即2 694 bp 的載體帶和1 755 bp的目的條帶(圖2)。

圖2 陽性質粒雙酶切鑒定結果

M：DL-5 000 DNA Marker；1：陽性質粒雙酶切鑒定結果

2.3 MEV VP2基因核苷酸和氨基酸同源性分析 測序結果表明，9份樣品VP2序列全長為1 755 bp，編碼548個氨基酸。應用DNAStar分析軟件，將其與GenBank上公布的4株國內外MEV VP2基因進行序列比較，發現13份樣品MEV VP2基因的核苷酸和氨基酸均有較高的同源性，分別為97.5%～99.8%和96.9%～99.7%，遺傳變異率較低。見圖3和圖4。

2.4 MEV VP2基因進化分析 用DNAStar軟件對13株MEV VP2基因繪制了系統進化樹，見圖5。從圖上可以看出，根據系統進化樹將13份MEV VP2基因分為兩個組群，第1個族群以樣品WH1407為代表，還包括樣品WF1310、WH1508、WH1408、WH1407、WH1509、QD1407、WF1408、JL2010和MEV-vaccine(MEV-vac)。第二個族群以QD1309為代表，還包括YT1408、MEV-BB、MEV-e。其中，組群一中的9株MEV又分為兩個亞群，WF1310、WH1508、WH1408、WH1407、Jl2010屬于第一亞群，WH1509、QD1407、WF1408和MEV-vac屬于第二亞群.

圖3 MEV VP2基因核苷酸序列分析

圖4 MEV VP2基因氨基酸序列分析

圖5 MEV VP2基因進化樹分析

2.5 9株MEV VP2蛋白的氨基酸變異分析 本研究通過對比13份MEV VP2蛋白的氨基酸序列表明，MEV VP2的第80、93、103、323、564、568位氨基酸是決定宿主范圍的關鍵氨基酸，相對比較保守；不同毒株主要在第5、62、87、101、232、236、279、300、534位氨基酸發生了替換，QD1309在第5、87、279、534位分別發生了A→T、M→T、W→R、V→I替換,QD1407分別在第236、300、534位發生了S→T、I→A、V→I替換，WH1408分別在第62、101、232、236、300位發生了T→A、T→I、V→I、S→T、V→I替換，WH509和WF1408在第5位點發生了A→T替換。 WF1310分別在第62、232、236、300位點發生了T→A、V→I、S→T、I→V替換。YT1408分別在第5、87、279位發生了A→T、M→T、W→R替換。

3 討論

通過對山東地區的9份MEV樣品和國內外4份樣品MEV VP2基因進行同源性分析表明，核苷酸和氨基酸具有很高的同源性，分別為97.5%～99.8%和96.9%～99.7%。說明山東地區的MEV VP2基因遺傳變異率較低。山東地區的9份MEV樣品和目前國內使用的疫苗毒株MEV-vac的核苷酸同源性為99.0%～99.5%，氨基酸同源性為98.6%～99.5%；通過進化樹結果表明，9份MEV 樣品屬于兩個組群，彼此都有一定的親緣性。山東地區大部分MEV和與疫苗毒株屬于同一個群組，說明分離馴化于20年以前的疫苗毒株仍然適用于目前山東地區水貂病毒性腸炎的防治。通過 13株 MEV 氨基酸序列比較分析，不同樣品主要在第5、62、87、101、232、236、279、300、534位氨基酸發生了替換，其中62、87、101、279、534位是在山東地區發現的新的變異位點。其中，QD1309、WH1509、WF1408、YT1408在第5位發生了A→T的替換，QD1407、WH1408、WF1310都在第236位發生了S→T的替換，是否發生生物學特性變化需要進一步研究。針對山東地區的MEV分子流行病學調查研究，對于控制該病在山東地區的流行就顯得非常有必要。本研究對發現的9份MEV樣和GenBank上公布的4株MEV參考株VP2基因序列進行分析對比，對山東地區MEV的變異情況有了初步了解，可為山東地區MEV的診斷及防控提供參考。

[1] Reed A P,Jones E V,Miller T J.Nucleotied sequence and genome organization of canine parvovirus[J].J Virol，1988，62(1):266-276.

[2] Schofield F W.virus enteritis in mink[J].North Am Vet,1949,30;651-654.

[3] Truyen U,Agbandje M,Parrish C R.Characterization of the feline host range and a specific epitope of feline panleukopenia virus[J].Virology,1994,200(2)：494-503.

[4] 楊玲，徐向明，殷俊,等.犬細小病毒病毒分離株 VP2基因的克隆與序列分析[J].揚州大學學報，2002，2(3)：12-13.

[5] 閆喜軍，張海玲，柴秀麗,等.水貂腸炎細小病毒結構蛋白VP2基因的克隆和序列分析[J].特產研究，2011(4):1-3.

Cloning and Sequencing of VP2Gene of Mink Enteritis Parvovirus in SHAN DONG Province

CHEN Tong ， YANG Rui-mei , SHAN Hu

(Qingdao Agricultural University,Key laboratory Preventive veterinary medicine in Shandong province,Qingdao 266109,China)

We tested samples of suspected parvovirus feces in Weifang,Weihai,Yantai,and Qingdao and got 9 VP2genes from QD1309,QD1407,WH1407,WH1408,WH1508,WH1509,WF1310,WF1408,and YT1408.We compared them with 4 MEV VP2genes out or in our country from GenBank.Results showed that the 13 MEV VP2genes had high homology of nucleotide and amino acids,which were 97.5%～99.8% and 96.9%～99.7%.Phylogenetic analysis showed that QD1309,YT1408,MEV-BB,and MEV-e belonged to the same group,and the difference for different strains were mainly at the site of 5,62,87,101,232,236,279,300,and 534 amino acid.

mink enteritis virus(MEV);VP2gene; genetic variation analysis

s:SHAN Hu;YANG Rui-mei

2015-12-25

科技部科技基礎性工作專項(2012FY111000)；公益性行業(農業)科研專項(201303042)；山東省自主創新及成果轉化專項(2014ZZCX07105);山東省農業重大應用技術創新課題;山東省農業科技成果轉化資金項目

陳童(1983-)，男，碩士，從事預防獸醫學領域研究，E-mail:1141549063@qq.com

單虎，E-mail: shanhu67@163.com; 楊瑞梅，E-mail: yrm.cc@163.com

Q78

A

0529-6005(2016)12-0028-03