張 帥，顧文源，趙云環(huán)，劉 瑩，左玉柱，范京惠

（1. 河北農(nóng)業(yè)大學(xué) 動物醫(yī)學(xué)院，河北 保定 071001；2. 河北省動物疫病預(yù)防控制中心，河北 石家莊 050035；3. 河北省獸醫(yī)生物技術(shù)創(chuàng)新中心，河北 保定 071001）

非洲豬瘟（African swine fever，ASF）是由非洲豬瘟病毒（African swine fever virus，ASFV）感染引起的家豬和野豬的最具破壞性的出血性傳染病之一。早在撒哈拉以南非洲國家流行，2007 年傳入格魯吉亞地區(qū)，近年來，又傳播到西歐和東南亞各國，2018 年傳入我國后，迅速在各地豬場蔓延，給我國養(yǎng)豬業(yè)造成了毀滅性的打擊。根據(jù)毒力大小，ASFV可分為高毒力、中毒力、低毒力和無癥狀感染毒株。在非洲豬瘟的臨床病例中，高致病性病毒感染最為常見，發(fā)病率和死亡率可高達100%。因此，世界動物衛(wèi)生組織（OIE）將ASF 列為必須報告的動物疫病，在我國也將ASF 列為一級動物疫病［1］。

ASFV 是一種結(jié)構(gòu)復(fù)雜的大型、雙鏈、線性DNA 病毒，基因組龐大，大小約為170 ～190 kb，編碼超過200 種蛋白質(zhì)。由于其結(jié)構(gòu)的復(fù)雜性，盡管國內(nèi)外學(xué)者對其進行了大量研究，大多數(shù)蛋白質(zhì)的功能仍不清楚，阻礙了疫苗的研發(fā)和致病機制及免疫逃避機制的探明，導(dǎo)致目前尚無商品化疫苗及特效的藥物用于ASF 的有效防控［2］。在編碼的200多種蛋白質(zhì)中，超過50 種是構(gòu)成病毒結(jié)構(gòu)的主要蛋白，其中由EP402R 基因編碼的CD2v 蛋白是ASFV晚期表達蛋白［3］，可介導(dǎo)紅細胞對ASFV 感染細胞的吸附作用，抑制有絲分裂原誘導(dǎo)的淋巴細胞增殖，延長免疫系統(tǒng)消滅病毒的時間，同時還是ASFV 血清型分型的關(guān)鍵蛋白［4］，且在機體免疫調(diào)節(jié)、免疫逃逸和ASFV 病毒入侵宿主和傳播擴散中發(fā)揮重要作用。目前該蛋白的結(jié)構(gòu)尚未完全解析，因此深入了解CD2v 蛋白的結(jié)構(gòu)及功能，將有助于ASFV 的有效防控［5］。本研究以中國首次分離且在國內(nèi)流行的非洲豬瘟病毒Pig/HLJ/2018 株為研究對象，對ASFV CD2v 蛋白進行了遺傳進化分析、理化性質(zhì)分析、蛋白結(jié)構(gòu)和亞細胞定位預(yù)測、蛋白跨膜區(qū)和信號肽分析以及B 細胞和T 細胞抗原表位分析等生物信息學(xué)分析［6］，以期為ASFV 的進化分析和CD2v蛋白的功能研究、ASFV 與宿主細胞的作用機制研究以及疫苗研發(fā)提供一些基礎(chǔ)性數(shù)據(jù)，為ASF 的有效防控奠定基礎(chǔ)。

1 材料與方法

1.1 CD2v 蛋白遺傳進化分析

通過登錄NCBI 查詢中國流行株P(guān)ig/HLJ/2018株（GenBank：QBH90580）ASFV EP402R 基因并選取其他國內(nèi)、外具代表性的毒株20 條，獲取基因組和基因型等信息（見表1），用MEGA 7.0 軟件進行遺傳進化樹構(gòu)建，描述20 種不同基因組數(shù)據(jù)之間的關(guān)系；用Megalign 軟件對20 條氨基酸序列進行比較和分析。

表1 20 株ASFV EP402R 基因統(tǒng)計信息Table 1 Statistical information of EP402R gene of 20 ASFV strains

1.2 CD2v 蛋白理化性質(zhì)分析

利用ExPASy 在線分析工具ProtParam 對CD2v蛋白的分子量、等電點、氨基酸組成、原子數(shù)、疏水性等理化性質(zhì)進行分析。

1.3 CD2v 蛋白結(jié)構(gòu)預(yù)測

利 用 在 線 工 具GOR4、Swiss-Model、cNLS Mapper、Psort 等分析軟件對CD2v 蛋白NLS 序列和二、三級結(jié)構(gòu)和亞細胞定位情況進行預(yù)測［7］。

1.4 CD2v 蛋白跨膜區(qū)和信號肽預(yù)測

通過TMHMM-2.0 在線工具分析CD2v 蛋白跨膜區(qū)［8］，通過SignalP-5.0 在線工具分析CD2v 蛋白信號肽［9］。

1.5 CD2v 蛋白B 細胞和T 細胞抗原表位分析

采用在線分析工具BepiPred 1.0、IEDB 分析序列B 細胞抗原表位區(qū)域［10］，采用NetMHC 4.0、NetMHCpan 4.1 在線服務(wù)器分析序列T 細胞抗原表位區(qū)域［11］，并結(jié)合蛋白二級結(jié)構(gòu)預(yù)測中的α-螺旋、β-折疊、β-轉(zhuǎn)角等不易形成抗原表位的區(qū)域，篩選出優(yōu)勢表位區(qū)域。

2 結(jié)果與分析

2.1 基于Pig/HLJ/2018 株ASFV CD2v 蛋白序列遺傳進化分析結(jié)果

2.1.1 基于ASFV CD2v 蛋白基因EP402R 的進化樹 通 過 對Pig/HLJ/2018 株ASFV 和 從NCBI 查詢篩選出的ASFV CD2v 蛋白基因EP402R 序列組成20 條序列構(gòu)建進化樹（圖1）結(jié)果顯示，Pig/HLJ/2018 株 與 格 魯 吉 亞Georgia 2007、Georgia 2008 以及其他中國毒株處于同一分支，基因組分型同為Ⅱ型，同時顯示當(dāng)前格魯吉亞地區(qū)、中國、東南亞地區(qū)流行毒株均為基因Ⅱ型；西歐地區(qū)流行毒株主要為基因Ⅰ型；非洲地區(qū)當(dāng)前流行毒株多為Ⅸ型和Ⅹ型，從進化關(guān)系來看，Ⅸ型和Ⅹ型具有很高的相似性［12］。

圖1 20 株ASFV EP402R 基因進化樹Fig.1 20 ASFV EP402R gene evolutionary tree

2.1.2 ASFV CD2v 蛋白氨基酸序列分析 通過Megalign 軟件對20 條氨基酸序列進行比較和分析后結(jié)果（見圖2）顯示，不同基因型之間氨基酸序列存在著不同程度的突變、插入與缺失，同一基因型間氨基酸序列的一致性較高，Pig/HLJ/2018 株與格魯吉亞Georgia 2007、Georgia 2008 以及中國AnhuiXCGQ 等基因Ⅱ型毒株CD2v 蛋白氨基酸序列一致性高達100%，與其他基因型的氨基酸一致性在76.8%～81.7%之間。同時還可以看出，CD2v 蛋白氨基酸序列末端（C 端）出現(xiàn)了KPCPPP 的重復(fù)序列。

圖2 ASFV CD2v 氨基酸序列分析結(jié)果Fig.2 ASFV CD2v amino acid sequence analysis results

2.2 CD2v 蛋白理化性質(zhì)分析結(jié)果

根據(jù)ExPASy 在線工具中的ProtParam 分析結(jié)果顯示：CD2v 蛋白共360aa，分子量為41.008 89 kD，其中天冬氨酸（12.5%）、異亮氨酸（12.5%）、脯氨酸（11.4%）和蘇氨酸（9.7%）占比最高，其他氨基酸占比從0.6% ～8.9% 不等；理論等電點為6.21，帶負電荷的氨基酸總數(shù)為24，帶正電荷的氨基酸總數(shù)為22；原子總數(shù)為5 749，分子式為C1868H2858N460O548S15；預(yù)測280 nm 消光系數(shù)為51 800 M-1cm-1～52 550 M-1cm-1；在哺乳動物體外網(wǎng)織紅細胞中的半衰期為30 h，在酵母體內(nèi)半衰期大于20 h，在大腸桿菌體內(nèi)半衰期大于10 h；不穩(wěn)定指數(shù)為47.78，脂肪指數(shù)為88.81，親水性總平均值為-0.209，說明該蛋白屬于親水性蛋白。

2.3 CD2v 蛋白結(jié)構(gòu)和亞細胞定位預(yù)測

Pig/HLJ/2018 株ASFV CD2v 蛋白一級結(jié)構(gòu)表明該基因全長1 080 bp，G+C 含量為27.88%，A+T含量為72.12%；表達的蛋白為41.008 89 kD，通過在線工具GOR4 預(yù)測的CD2v 蛋白二級結(jié)構(gòu)結(jié)果（圖3）顯示，α-螺旋占5.56%，擴展鏈占31.11%、無規(guī)卷曲占63.33%；將CD2v 蛋白序列提交在線工具Swiss-Model 預(yù)測該蛋白的三級結(jié)構(gòu)模型（見圖4），證明了該蛋白主要以無規(guī)卷曲為主。

圖3 CD2v 蛋白二級結(jié)構(gòu)預(yù)測結(jié)果Fig.3 Rediction results of secondary structure of CD2v protein

圖4 CD2v 蛋白三級結(jié)構(gòu)預(yù)測結(jié)果Fig.4 Prediction of tertiary structure of CD2v protein

根據(jù)在線工具cNLS Mapper 和Psort 預(yù)測結(jié)果顯示，該蛋白可能包括一段NLS 單分型序列（表2），無雙分型序列，其cut-off 分數(shù)為7.5；CD2v 蛋白可能定位于高爾基體的概率為22.2%，定位于內(nèi)質(zhì)網(wǎng)和內(nèi)膜的概率均為11.1%。

表2 ASFV CD2v 蛋白序列單分型NLS 預(yù)測Table 2 Single NLS prediction of ASFV CD2v protein sequence

2.4 CD2v 蛋白跨膜區(qū)和信號肽分析結(jié)果

根據(jù)TMHMM-2.0 在線工具分析結(jié)果顯示，CD2v 蛋白存在一個跨膜區(qū)（見圖5-A），信號肽分析結(jié)果的可能性為0.925，說明存在信號肽（見圖5-B）。

圖5 CD2v 蛋白跨膜區(qū)和信號肽預(yù)測結(jié)果Fig.5 Prediction results of CD2v protein transmembrane region and signal peptide

2.5 CD2v 蛋白B 細胞和T 細胞抗原表位分析結(jié)果

根據(jù)在線分析工具BepiPred 1.0、IEDB 分析序列B 細胞抗原表位區(qū)域結(jié)果顯示，存在9 個B 細胞抗原表位；根據(jù)NetMHC 4.0、NetMHCpan 4.1 在線服務(wù)器分析序列T 細胞抗原表位區(qū)域結(jié)果顯示，CD2v 蛋白存在14 個T 細胞抗原表位；結(jié)合2.3 中預(yù)測的二級結(jié)構(gòu)，去掉其中的α-螺旋等篩選出4 個優(yōu)勢B 細胞抗原表位和5 個優(yōu)勢T 細胞抗原表位（見表3）。

表3 CD2v 蛋白優(yōu)勢B、T 細胞抗原表位分析統(tǒng)計結(jié)果Table 3 Statistical results of CD2v protein dominant B cell and T cell epitope analysis

3 討論與結(jié)論

ASFV 是一種基因組龐大而復(fù)雜的DNA 病毒，可導(dǎo)致豬的急性死亡，給全球養(yǎng)豬業(yè)造成了巨大的經(jīng)濟損失。我國于2018 年在遼寧首次爆發(fā)ASF［13］，之后迅速蔓延至多個省市和地區(qū)，并呈區(qū)域性流行，嚴重威脅著我國養(yǎng)豬業(yè)的健康發(fā)展。盡管政府部門按照一類傳染病的防控措施，對疫區(qū)及受威脅地區(qū)進行了嚴格的管控［14］，但防控形勢依然嚴峻，急需安全有效的疫苗用于本病的防控。然而，由于對ASFV 與宿主之間的作用機制、病毒在宿主體內(nèi)的復(fù)制調(diào)控機制尚不清楚，對與病毒免疫逃避相關(guān)蛋白的研究相對較少，目前尚無商品化的非洲豬瘟疫苗。因此，深入研究病毒蛋白的結(jié)構(gòu)和功能，將有助于病毒致病機制的闡明以及非洲豬瘟疫苗的研發(fā)。

3.1 ASFV CD2v 蛋白的遺傳變異分析

CD2v 蛋白位于ASFV 的外層囊膜，可誘導(dǎo)機體產(chǎn)生中和抗體，是病毒毒力及免疫逃避的相關(guān)蛋白之一［15-16］。本研究對ASFV CD2v 蛋白氨基酸序列分析結(jié)果顯示，不同基因型之間氨基酸序列存在較大差異，以中國流行株P(guān)ig/HLJ/2018 株ASFV CD2v 蛋白氨基酸序列為例，在C 末端出現(xiàn)了連續(xù)的重復(fù)序列（KPCPPP），研究表明，該重復(fù)序列是一種遺傳標記或抗原表位［17-18］，Yang 等［18］確認了該重復(fù)序列（KPCPPP）為細胞穿透肽(CPP)，CPP 是多種多肽家族，通常由5 ～30 個氨基酸組成，可通過能量依賴或能量獨立機制通過組織和細胞膜，不與特定受體相互作用［19］，它可以攜帶EGFP 通過網(wǎng)格蛋白和微胞飲介導(dǎo)的內(nèi)吞進入CHO 細胞，這為ASFV 顆粒進入細胞的方式提供了新的見解，同時也提示ASFV研究可以發(fā)展成為藥物傳遞的工具。另外，連續(xù)的脯氨酸序列與宿主SH3 結(jié)構(gòu)域蛋白SH3P7/mAbp1 結(jié)合，并在囊泡運輸中發(fā)揮作用［20］。對ASFV CD2v 蛋白遺傳進化與氨基酸序列分析表明該蛋白可作為ASF 治療藥物和預(yù)防疫苗的研究靶點，為ASF 的有效防控奠定基礎(chǔ)。

3.2 ASFV CD2v 蛋白結(jié)構(gòu)和抗原表位

對ASFV CD2v 蛋白的結(jié)構(gòu)分析結(jié)果顯示，CD2v 是一種跨膜蛋白，蛋白二級結(jié)構(gòu)顯示α-螺旋結(jié)構(gòu)占5.56%，存在跨膜區(qū)和信號肽，具有由2 個免疫球蛋白樣結(jié)構(gòu)域組成的細胞外N 端區(qū)域和胞質(zhì)C 端區(qū)域，與Rodríguez 等［21］的研究結(jié)果相一致。Sereda 等［22］研究表明CD2v 蛋白能介導(dǎo)紅細胞對ASFV 感染細胞的吸附作用，并指出ASF 保護性免疫可能是血清型特異性的表現(xiàn)，血液吸附(HAI)抑制同一血清型內(nèi)的病毒互相交叉保護，而同一血清型外的病毒則沒有交叉保護；同時Malogolovkin 等［23］根據(jù)CD2v 蛋白HAI 特性進行了血清學(xué)分型，發(fā)現(xiàn)同一基因型下的毒株表現(xiàn)出不同的血清型分型，這可能更好地解釋ASFV 分離株之間的同源交叉保護問題和影響疾病出現(xiàn)的病毒決定因素，表明ASFV CD2v 蛋白的血清分型對疫苗的設(shè)計和開發(fā)非常重要。T 細胞介導(dǎo)的細胞免疫和B 細胞介導(dǎo)的體液免疫，在機體抵抗病原微生物感染的過程中，發(fā)揮著重要作用，B、T 細胞的優(yōu)勢抗原表位常用于靶向藥物或表位疫苗的研制，Burmakina 等［24］在免疫動物攻毒后7 ～10 d，觀察到T 細胞對病毒的強大反應(yīng)，并鑒定出了在ASF 保護性免疫環(huán)境中誘導(dǎo)T 細胞應(yīng)答的T 細胞表位區(qū)域，這表明ASFV 的T 細胞表位對保護性宿主反應(yīng)具有重要意義。同時，Burmakina等［24］提示T 細胞表位對給定的病毒血清型也有特異性，T 細胞宿主反應(yīng)可能與觀察到的血清型特異性保護有關(guān)。本研究通過在線工具篩選出4 個優(yōu)勢B 細胞抗原表位和5 個優(yōu)勢T 細胞抗原表位，可作為靶向藥物、表位疫苗篩選的對象［6］，有助于ASFV 亞單位疫苗的設(shè)計和開發(fā)。

通過對ASFV 的CD2v 蛋白及其基因進行生物信息學(xué)分析，深入解析了中國流行毒株P(guān)ig/HLJ/2018株CD2v 蛋白基因的遺傳進化情況以及CD2v 蛋白的結(jié)構(gòu)與功能關(guān)系，為進一步闡明病毒-宿主相互作用機制以及ASFV 疫苗的研發(fā)奠定了基礎(chǔ)。