黑龍江地區(qū)鯉皰疹病毒3型的流行病學調(diào)查及其主要囊膜蛋白的抗原表位預測

李 淵，徐黎明，趙景壯，劉 淼，任廣明，尹家勝，紀 鋒，盧彤巖

(1.中國水產(chǎn)科學研究院黑龍江水產(chǎn)研究所，哈爾濱 150070；2.上海海洋大學水產(chǎn)與生命學院，上海 201306)

黑龍江地區(qū)鯉皰疹病毒3型的流行病學調(diào)查及其主要囊膜蛋白的抗原表位預測

李 淵1,2，徐黎明1，趙景壯1，劉 淼1，任廣明1，尹家勝1，紀 鋒1，盧彤巖1

(1.中國水產(chǎn)科學研究院黑龍江水產(chǎn)研究所，哈爾濱 150070；2.上海海洋大學水產(chǎn)與生命學院，上海 201306)

2015-2016年對來自黑龍江地區(qū)46個養(yǎng)殖場送檢的鯉魚進行錦鯉皰疹病毒(Koi herpes virus，KHV)的分離鑒定，并應用DNAStar等軟件對檢出的陽性樣本的主要囊膜蛋白(Major envelop protein，MEP)基因進行生物信息學分析。結果顯示，KHV的檢出率約為6.5%，3株陽性樣本具有相同的MEP基因序列，稱之為KHV-HLJ1516。MEP基因開放閱讀框為771 bp，編碼256個氨基酸，相對分子質(zhì)量為28 280.58，等電點為8.40。對KHV-HLJ1516的MEP基因序列同源性分析表明，其與KHV-GZ11分離株僅有1個堿基不同，但并未導致氨基酸發(fā)生改變，堿基相似性為99%；而與HZ419分離株相比，KHV-HLJ1516的MEP基因有3個堿基發(fā)生突變，相應的氨基酸由Asn變?yōu)镮Ie。聚類分析結果顯示，KHV-HLJ1516與KHV-GZ11處于同一分支，而HZ419位于另一分支。B細胞抗原表位預測結果顯示，KHV-HLJ1516的MEP的抗原表位主要位于4～11，24～40，42～63，82～110，209～218，238～249氨基酸區(qū)段。

鯉皰疹病毒3型；進化分析；流行病學；B細胞抗原表位；預測

錦鯉皰疹病毒病(Koi herpes virus disease，KHVD)是由錦鯉皰疹病毒(Koi herpes virus，KHV)引起的錦鯉和鯉魚的一種高致病性疾病，出現(xiàn)癥狀后24～48 h死亡，死亡率高達80%～90%[1]。到目前為止，包括亞洲、歐洲和北美洲等多個國家均發(fā)現(xiàn)該病毒。所以快速準確的診斷已成為關注的重點，目前的診斷方法有細胞培養(yǎng)分離技術[2,3]、電鏡觀察[4,5]、PCR[6]、酶聯(lián)免疫吸附劑試驗(ELISA)[7,8]、原位雜交技術[9]及環(huán)介導等溫擴增技術(LAMP)[10]等。

錦鯉皰疹病毒(KHV)，又稱鯉皰疹病毒Ⅲ型(Cyprinid herpesvirus 3，CyHV-Ⅲ)，屬于皰疹病毒科，是雙鏈DNA病毒，有囊膜，與同科的鯉魚皰疹病毒(cyprinid herpesvirus 1，CyHV-1)、叉尾鮰病毒(channelcatfish virus，CCV)之間有免疫交叉反應[11]。主要囊膜蛋白(Major envelop protein，MEP)是錦鯉皰疹病毒囊膜的主要組成部分，其功能與病毒對細胞吸附、穿入、從細胞核膜出芽釋放及誘導細胞融合有關，并有誘導產(chǎn)生中和抗體和細胞毒作用[12]。本研究通過對2015～2016年間我國黑龍江地區(qū)送檢的鯉魚進行常規(guī)的錦鯉皰疹病毒(Koi herpes virus,KHV)病原監(jiān)測，采用了KHV國標檢測方法，得到3株具有相同MEP基因序列的陽性樣本，稱之為KHV-HLJ1516；通過對KHV-HLJ1516的MEP的基因序列的系統(tǒng)進化樹分析，確定其與其他區(qū)域KHV之間的遺傳進化關系，掌握該分離株的遺傳背景，為我國KHVD的防治提供參考；此外本研究應用相關生物信息學軟件對MEP序列進行了B細胞抗原表位綜合預測，旨在為進一步研究KHV的MEP的功能及其基因背景信息和分子流行病學的研究奠定基礎。

1 材料與方法

1.1 材料

取自2015-2016年間我國黑龍江地區(qū)送檢的鯉魚，采集肝、脾和頭腎為待檢樣品。采樣數(shù)量按SC/T7014-2016中的6.1的規(guī)定。DNA聚合酶rTaq premix(大連寶生物公司)；DNAzol ? Reagent(Invitrogen,USA)；DNA Marker 及PCR純化試劑盒為TaKaRa公司產(chǎn)品。引物均由哈爾濱博萊仕合生物公司合成。

1.2 DNA的提取

將100 μL組織懸液放入1 mL的離心管，再加入1 mL DNAzol，顛倒混勻5次，室溫下靜置5 min，然后10 000g離心10 min；吸取上清，加入0.5 mL無水乙醇，顛倒混勻5～8次，室溫下靜置1-3 min，然后5 000g離心5 min；吸棄上清，加入0.8～1 mL 75%乙醇漂洗兩次，然后離心；吸棄乙醇，干燥10-15 min；加入0.2～0.3 mL無菌水溶解，即為提取的DNA，作為PCR模板。

1.3 KHV的鑒定

以提取的DNA作為PCR模板。用于KHV檢測的引物為國標提供，引物序列：P3∶5′-GACACCACATCTGCAAGGAG-3′，P4∶5′-GACACATGTTACAATGGTCGC-3′。50 μL的體系：rTag 25 μL，引物2 μL，模板2 μL，水21 μL。PCR反應擴增特異性片段。反應參數(shù)為94 ℃ 30 s，隨后是94 ℃ 30 s，63 ℃ 30 s，72 ℃ 30 s，72 ℃ 7 min，共40個循環(huán)，最后保持在4 ℃。1%瓊脂糖凝膠電泳檢測PCR產(chǎn)物。

1.4 MEP基因擴增、序列分析及系統(tǒng)進化樹構建

以GenBank已登錄的錦鯉皰疹病毒MEP基因序列，設計合成引物：P1∶5′-ATGGCAGTCACCAAAGC-3′，P2∶5′-TCACCACATCTTGCCG-3′。以提取的DNA作為PCR模板，進行MEP基因的擴增，反應條件為：93 ℃ 2 min；隨后是93 ℃ 1 min，52 ℃ 1 min，72 ℃ 1 min，72 ℃ 10 min，35個循環(huán)，最后保持在4 ℃。結束后進行1%瓊脂糖凝膠電泳分析，對目的條帶進行膠回收，然后將回收產(chǎn)物送上海生工生物公司進行測序。利用DNAStar中的EditSeq軟件推導MEP基因的氨基酸序列。并用DNAman對MEP基因序列進行比對。將獲得的MEP基因序列BLAST進行序列同源性分析并用MEGA5.05軟件構建系統(tǒng)進化樹。選擇的參考序列來自Genbank中的日本株TUMST1(GenBank：AP008984.1)、美國株KHV-U(GenBank：DQ657948.1)和以色列株KHV-I(GenBank：DQ177346.1)，還有國內(nèi)的KHV-GZ11(GenBank：KJ627438.1[13]及HZ419 GenBank：KP004896.1)。本研究所分離的KHVHLJ-1516的MEP基因序列已提交Genbank，其accession no.為KY010497。

1.5 MEP的二級結構和抗原表位預測

應用DNAStar軟件中的Protean模塊分析KHV-HLJ1516的MEP的B細胞抗原表位；Chou-Fasman法和Gamier-Robson法預測其二級結構；Karplus-Schulz法預測其柔性區(qū)段；Kyte-Doolittle法分析其親水性；Emini法預測其表面可能性；Jameson-Wolf 法和Kolaskar-Tongaonkar法(http://tools.immuneepitope.org/tools/bcell/iedb_input)分別預測其抗原指數(shù)和抗原表位指數(shù)。

2 結果與分析

2.1 病魚剖檢

感染KHV后的魚皮膚上出現(xiàn)蒼白的塊斑與血泡，體表分泌大量黏液，眼睛凹陷，腹部充血和出血(圖1)。

圖1 病鯉的癥狀Fig.1 Symptoms of infected common carp

2.2 KHV的鑒定

對送檢的46份鯉魚采樣(肝、脾、頭腎)，提取核酸，利用國標提供的引物進行以Sph為目標基因的PCR擴增。凝膠電泳結果顯示在預計大小位置出現(xiàn)特異性目的片段(292 bp)(圖2)。通過測序對比三株陽性樣本Sph基因片段完全相同，且與Genbank收錄的鯉KHV的Sph基因具有100%的相似性，說明三個樣本均為KHV。本次調(diào)研結果顯示黑龍江地區(qū)2015～2016年間的KHV檢出率約為6.5%。

圖2 KHV PCR擴增產(chǎn)物的凝膠電泳分析Fig.2 Gel electrophoresis of PCR product of KHV1，2，3陽性樣本PCR產(chǎn)物；N，陰性對照； M，DL2000 DNA marker

2.3 KHV-HLJ1516的MEP基因的擴增及序列測定

以檢出的陽性樣本為模板，P1、P2為引物，進行PCR基因擴增，然后通過瓊脂糖凝膠電泳分析，結果顯示得到了與目的條帶大小一致的特異性條帶(圖3)。序列分析結果表明MEP基因開放閱讀框長為771 bp，編碼含有256個氨基酸的主要囊膜蛋白。MEP預測的相對分子質(zhì)量為28 280.58，等電點為8.40。值得注意的是這三株陽性樣本的MEP具有相同的核酸序列。因此，本研究將這三株KHV分離株稱之為KHV-HLJ1516。

圖3 KHV-HLJ1516主要囊膜蛋白基因PCR 擴增產(chǎn)物的凝膠電泳分析Fig.3 Gel electrophoresis of PCR product of MEP gene of KHV-HLJ15161，2，3陽性樣本PCR產(chǎn)物；M，DL2000 DNA marker

2.4 KHV-HLJ1516的MEP基因序列比對及遺傳進化分析

將測序結果與GenBank中收錄的MEP基因序列進行比較，結果表明，廣州某地區(qū)的KHV-GZ11、KHV-I、KHV-U和TUMST1分離株的MEP核酸序列相似性均為100%，而KHV-HLJ1516的MEP基因與KHV-GZ11的MEP基因相似性為99%，其MEP基因有1個堿基發(fā)生了突變，但并未導致該位點的氨基酸發(fā)生改變；而與國內(nèi)的HZ419分離株相比，KHV-HLJ1516的MEP基因有3個堿基發(fā)生了突變，并影響了相應的氨基酸由Asn變成了IIe。系統(tǒng)進化樹結果顯示，KHV-HLJ1516與KHV-GZ11與KHV-I、KHV-U和TUMST1分離株均處于同一個分支，而國內(nèi)的HZ419分離株卻位于另一分支(圖4)。

圖4 基于KHV-HLJ1516主要囊膜蛋白基因序列構建的系統(tǒng)進化樹Fig.4 Phylogenetic tree based on the MEP gene sequence of KHV-HLJ1516

2.5 KHV-HLJ1516的MEP二級結構和抗原表位分析

應用Garnier-Robson法和Chou-Fasman法預測MEP的二級結構發(fā)現(xiàn)，其預測的α-螺旋、β-折疊和β-轉角存在較多的重合區(qū)域(圖5)。預測的二級結構顯示：KHV-HLJ1516的MEP是以α-螺旋和β-折疊為主有一定的β-轉角和無規(guī)則卷曲結構的混合型蛋白，α-螺旋和β-折疊一般是起支架作用，蛋白質(zhì)的功能部位和酶的活動中心一般在無規(guī)則卷曲部位。

圖5 不同方法預測的KHV-HLJ1516的主要囊膜蛋白的二級結構Fig.5 Secondary structure of MEP of KHV-HLJ1516 predicted by various methods A: α-螺旋；B: β-折疊；T: β-轉角；C: 無規(guī)則卷曲。

Karplus-Schulz法預測的MEP柔性區(qū)段位置為：5，10～16，22～30，40～41，71～73，108～110，149～152，177～180，192～196，201～205，

218～239，243～253，這些區(qū)域有一定的柔韌性，形成抗原表位的可能性較大(圖6)。

圖6 Karplus-Schulz法預測KHV-HLJ1516的主要囊膜蛋白的柔性區(qū)段Fig.6 Flexible regions of MEP of KHV-HLJ1516 predicted by Karplus-Schulz method

Kyke-Doolittle法氨基酸親水性預測結果顯示，KHV-HLJ1516的MEP的親水性較低，親水性區(qū)段主要在氨基酸序列的后段，提示該區(qū)段位于蛋白表面的可能性最大，作為抗原表位的幾率也相對較高。Emini法分析MEP的表面可能性結果顯示，可能性較大的區(qū)域為：7～13，68～75，147～151，198～214，221～237，245～253，因此，這些區(qū)段很有可能位于蛋白的表面，可能是B細胞的優(yōu)勢區(qū)域(圖7)。

圖7 KHV-HLJ1516的主要囊膜蛋白的表面可能性分析Fig.7 Surface probability of MEP of KHV-HLJ1516

2.6 KHV-HLJ1516的MEP抗原指數(shù)和抗原表位指數(shù)

Jameson-Wolf法對MEP的抗原指數(shù)進行預測，可看出MEP中存在許多抗原指數(shù)較高的區(qū)域，其中區(qū)段最大的位置在218～239，因此該段有可能是優(yōu)勢抗原表位；還有一些其他抗原指數(shù)較高的區(qū)段，如7～16，22～31，36～43，69～77，104～111，184～196，198～211，243～255形成抗原表位的可能性也較大(圖8)。

圖8 KHV-HLJ1516的主要囊膜蛋白的抗原指數(shù)分析Fig.8 Antigenic index of MEP of KHV-HLJ1516

Kolaskar-Tongaonkar法[14]預測的MEP抗原表位平均抗原指數(shù)為1.051，最大值為1.265，最小值為0.917。

2.8 KHV-HLJ1516的MEP B細胞表位預測結果綜合分析

如果某一氨基酸區(qū)段中：抗原指數(shù)≥1.033，親水性指數(shù)≥0，氨基酸的抗原表位可能性指數(shù)≥1，氨基酸個數(shù)大于等于7個，并且其內(nèi)部或附近又含有柔性結構區(qū)域，則這一區(qū)段有可能是抗原表位。Kolaskar-Tongaonkar方法預測大于7個氨基酸的KHV-HLJ1516的MEP抗原表位分析結果可以直接給出抗原表位的序列，見表1。表位大小至多不超過30個氨基酸殘基。綜合以上考慮，本研究預測的KHV-HLJ1516的MEP可能的B細胞抗原表位區(qū)段主要位于4～11，24～40，42～63，82～110，209～218，238～249氨基酸區(qū)段。

表1 KHV-HLJ1516的主要囊膜蛋白的B細胞抗原表位的氨基酸序列Tab 1 Amino acid sequence of B cell antigenic epitope of MEP of KHV-HLJ1516

3 討論

KHVD自1997年首次報道以來，傳播、流行于亞洲、歐洲、北美及非洲多個國家[15]。2002年，劉葒等[16]首次報道KHVD傳到中國廣東，時隔三年后，周永燦等[17]報道已從中國自然感染鯉魚中檢測到KHV，2010年閆春梅等[18]又從吉林省某網(wǎng)箱養(yǎng)殖場檢測出KHV，至此表明，KHVD已在中國成蔓延趨勢，應引起重視。中國南北地理、氣候、環(huán)境等差異顯著，所以針對北方的KHVD監(jiān)測就很重要了。2015-2016年，在KHVD易發(fā)的季節(jié)，課題組對黑龍江省多家鯉魚養(yǎng)殖場送來的監(jiān)測鯉魚樣品進行病原檢測，從46家養(yǎng)殖場共檢出3個陽性樣本，檢出率約為6.5%，該結果說明KHV已經(jīng)存在并且流行于黑龍江省的鯉魚養(yǎng)殖場。KHV 主要感染鯉和錦鯉，其魚苗、幼魚、成魚階段均可感染[19]，其他鯉科魚、錦鯉和鯽的雜交后代、錦鯉和金魚的雜交后代也對其易感[20]。魚體感染KNV后，在其腦、眼、脾、鰓、腎、腸和腸內(nèi)容物都能檢測出該病毒[21]。所以認為該病是通過水傳播的。接種疫苗可能是一種控制KHVD發(fā)生的容易想到的辦法，該病的減毒苗已經(jīng)應用于臨床[22]，但目前還沒有理想的疫苗來防治KHVD[23]，主要靠加強監(jiān)測，切斷傳播途徑等措施防控。總之，為更好的防制該病暴發(fā)流行，需要進行區(qū)域的監(jiān)測管理[24]，加強全球的合作交流，建立有效的檢測方法，研制有效的KHV疫苗，保障我國錦鯉和鯉魚養(yǎng)殖業(yè)的健康發(fā)展。

在本研究中分離得到的3株KHV陽性樣本的MEP基因序列完全一致。在NCBI中對測序的結果直接BLAST分析，結果顯示KHV-HLJ1516的MEP核苷酸與參考病毒株的MEP核苷酸相似性均為99%，而KHV-I(以色列)、KHV-U(美國)、TUMST1(日本)和KHV-GZ11(中國廣州)的MEP核酸序列100%相同。該結果說明本研究分離的KHV毒株與我國其他地區(qū)分離株具有基因序列差異性，可能為新毒株。跟廣州KHV-GZ11分離株相比，黑龍江分離株KHV-HLJ1516的MEP基因有1個堿基發(fā)生了突變，但并未導致該位點的氨基酸發(fā)生改變。而與國內(nèi)另一參考毒株HZ419相比，其MEP基因有3個堿基發(fā)生了突變，并影響了相應的氨基酸由Asn變成了IIe。系統(tǒng)進化樹結果顯示，KHV-HLJ1516仍與這四株參考毒株處于同一個分支，說明KHV-HLJ1516與廣州KHV-GZ11及KHV-I(以色列)、KHV-U(美國)、TUMST1分離株具有相對較近的親緣關系，而與中國另一KHV分離株HZ419親緣關系相對較遠。

囊膜蛋白通常與病毒的感染甚至病毒的毒力相關[25]。本實驗成功獲得了KHV長771 bp的MEP基因，編碼256個氨基酸，預測的相對分子質(zhì)量為28 280.58，等電點為8.40。預測的二級結構顯示：KHV-HLJ1516的MEP是以α-螺旋和β-折疊為主有一定的β-轉角和無規(guī)則卷曲結構的混合型蛋白，α-螺旋和β-折疊一般是起支架作用，蛋白質(zhì)的功能部位和酶的活動中心一般在無規(guī)則卷曲部位。抗原表位即抗原決定簇，決定抗原性的特殊化學基團，有效的分析抗原表位有助于疾病的診斷以及疫苗分子的設計。抗原表位篩選方法主要集中在研究B細胞抗原表位，作為B細胞抗原表位首先應位于或易于移動于蛋白質(zhì)抗原表面，有利于與抗體結合，這樣就傾向于參與構成β轉角結構。另外，具有一定柔韌性，親水性強、表面極性大的氨基酸殘基通常參與表位構成[26]。表位大小一般由5～7個氨基酸和單糖殘基組成，至多也不超過30個氨基酸殘基。綜合以上考慮，本研究預測的KHV-HLJ1516的MEP可能的B細胞抗原表位區(qū)段主要位于4～11，24～40，42～63，82～110，209～218，238～249氨基酸區(qū)段，為進一步探索KHV分子流行病學及表位疫苗的構建研究奠定了理論基礎。

[1]何賢臣,李慶東,康安娜.鯉魚病毒病的診斷方法與預防措施[J].黑龍江水產(chǎn),2004,(4):21-26.

[2]李 茂,肖 丹,劉天強,等.鯉科皰疹病毒Ⅱ型射陽株的分離與鑒定[J].淡水漁業(yè),2015,(3):93-96.

[3]Imajoh M,Fujioka H,Furusawa K,et al.Establishment of a new cell line susceptible to Cyprinid herpesvirus 3 (CyHV-3) and possible latency of CyHV-3 by temperature shift in the cells [J].J Fish Dis,2015,38(6): 507-514.

[4]劉秋鳳,曾令兵,周 勇,等.微載體規(guī)模化培養(yǎng)草魚細胞與病毒的工藝及優(yōu)化[J].淡水漁業(yè),2013,43(4):69-75.

[5] Miyazaki T,Kuzuya Y,Yasumoto S,et al.Histopathological and ultrastructural features of Koi herpesvirus (KHV)-infected carp Cyprinus carpio,and the morphology and morphogenesis of KHV [J].Dis Aqua Org,2008,80(1): 1-8.

[6] Bergmann S M,Riechardt M,Fichtner D,et al.Investigation on the diagnostic sensitivity of molecular tools used for detection of koi herpesvirus[J].J Virol Methods,2010,163(2): 229-233.

[7] Ronen A,Perelberg A,Abramowitz J,et al.Efficient vaccine against the virus causing a lethal disease in cultured Cyprinus carpio [J].Vaccine,2003,21(32): 4677-4684.

[8] Adkison M A,Hedrick R P.An nzyme linked immunosorbent assay (ELISA) for detection of antibodies to the Koi herpesvirus (KHV) in the serum of Koi Cyprinus carpio [J].Fish Pathol,2009,40(2): 53-62.

[9] Bergmann S M,Kempter J,Sadowskim J,et al.First detection,confirmation and isolation of koi herpesvirus (KHV) in cultured common carp (Cyprinus carpio L.) in Poland [J].Bull Eur Assoc Fish Pathol,2006,26(2): 97-104.

[10]Soliman H,Elmatbouli M.Immunocapture and direct binding loop mediated isothermal amplification simplify molecular diagnosis of Cyprinid herpesvirus-3[J].J Virol Methods,2009,162(1-2): 91-95.

[11] Calle P P,Mcnamara T,Kress Y.Herpesvirus-associated papillomas in koi carp (Cyprinus carpio) [J].J Zoo Wildlife Med,1999,30(1): 165-169.

[12]孟祥俊,賀 冰.單純皰疹病毒糖蛋白的特性及其臨床意義[J].國際病毒學雜志,2002,9(6):184-188.

[13]董 剛,李雪竹,翁少萍,等.近三年華南地區(qū)錦鯉皰疹病毒的分子流行病學調(diào)查[A].中國水產(chǎn)學會魚病專業(yè)委員會2013年學術研討會論文摘要匯編[C],2013:1.

[14] Kolaskar A S,Tongaonkar P C.A semi-empirical method for prediction of antigenic determinants on protein antigen [J].Febs Letters,1990,276(1): 172-174.

[15]鄭樹城,王 慶,李瑩瑩,等.鯉皰疹病毒3型研究進展[J].病毒學報,2016,(1):108-120.

[16]劉 葒,史秀杰,高隆英,等.進口錦鯉暴發(fā)病病原的nested-PCR鑒定[J].華中農(nóng)業(yè)大學學報,2002,21(5):414-418.

[17]周永燦,袁軍法,任武澤,等.錦鯉皰疹病毒的檢測與人工感染試驗[J].武漢大學學報(理學版),2005,(S2):249-252.

[18]閆春梅,張雅斌,鄭 偉,等.應用聚合酶鏈式反應(PCR)檢測錦鯉皰疹病毒[J].河北漁業(yè),2010,(2):37-39.

[19]李玉平,毛洪順,唐紹林,等.錦鯉皰疹病毒病的診斷與防控[J].畜牧獸醫(yī)科技信息,2016,(10):32-33.

[20] Dishon A,Perelberg A,Bishara-Shieban J,et al.Detection of carp interstitial nephritis and gill necrosis virus in fish droppings [J].Appl Environ Microbiol,2005,71(11): 7285-7291.

[21]羅 丹,梁利國,謝 駿,等.鯉皰疹病毒Ⅰ、Ⅱ、Ⅲ型研究進展[J].水生態(tài)學雜志,2014,35(3):94-100.

[22]王 好,劉麗凡,于 丹,等.錦鯉皰疹病毒減毒株的篩選及免疫原性研究[J].黑龍江畜牧獸醫(yī),2016,(19):171-174.

[23]李俊超,馬 杰,曾令兵,等.抗錦鯉皰疹病毒的植物藥物篩選及其效果比較[J].淡水漁業(yè),2014,(3):62-67.

[24]賈敬德.淡水漁業(yè)環(huán)境現(xiàn)狀及保護對策[J].淡水漁業(yè),2004,34(5):59-61.

[25]王友超,周井祥,王 好,等.錦鯉皰疹病毒CJ株ORF27基因克隆及生物信息學分析[J].水產(chǎn)學雜志,2014,(2):20-24.

[26]李海俠,毛旭虎.蛋白質(zhì)抗原表位研究進展[J].微生物學免疫學進展,2007,35(1):54-58.

(責任編輯：張瀟峮)

Epidemiologic investigation of Koi herpes virus type 3 in Heilongjiang province and prediction of antigenic epitopes of major envelope protein

LI Yuan1,2,XU Li-ming1,ZHAO Jing-zhuang1,LIU Miao1,REN Guang-ming1,YIN Jia-sheng1,JI Feng1， LU Tong-yan1

(1.HeilongjiangRiverFisheryResearchInstitute,ChineseAcademyofFisherySciences,Harbin150070,China;2.CollegeofFisheryandLifeScience,ShanghaiOceanUniversity,Shanghai201306,China)

Carps from 46 fish farms of Heilongjiang province were subjected to isolation of Koi herpes virus (KHV) during 2015～2016.The major envelope protein (MEP) gene of the KHV was subjected to bioinformatics analyses by using the DNAStar and other software.The results showed that the detection rate of KHV in Heilongjiang province during 2015～2016 was about 6.5%.The 3 isolated KHV had the same MEP gene sequence,and which were called KHV-HLJ1516.The open reading frame (ORF) of the MEP was 771 bp,encoding 256 amino acids.The predicted relative molecular mass of the MEP was 28280.58 and isoelectric point was 8.40.The homology analysis showed that compared with KHV-GZ11,the MEP gene of the KHV-HLJ1516 had a single site nucleotide mutation which did not cause the change of amino acid.The MEP of KHV-HLJ1516 was 99% identical to KHV-GZ11 at the nucleotide level.Compared with HZ419,KHV-HLJ1516 had 3 nucleotide mutations which caused a mutation from Asn to IIe.The results of cluster analysis showed that KHV-HLJ1516 was clustered in the same branch with KHV-GZ11,while another isolate of HZ419 was clustered in another branch.B cell epitopes of the MEP showed that the likely B cell epitopes of KHV-HLJ1516 MEP were mainly located at 4～11,24～40,42～63,82～110,209～218,238～249 amino acids.

Cyprinid herpesvirus 3;phylogenetic analysis;epidemiology;B cell antigenic epitope;prediction

2016-10-26；

2017-03-18資助項目：黑龍江省應用技術研究與開發(fā)計劃項目(GA13B401)

李 淵(1989-)，女，碩士研究生，專業(yè)方向為魚類病害防治。E-mail: 1154662843@qq.com

盧彤巖。E-mail: lutongyan@hotmail.com*

S941

A

1000-6907-(2017)04-0037-07