張艷艷，張靜遠(yuǎn)，楊金金，楊金梅，韓 桪，米立娟，張 菲，齊 宇，張守峰，王 穎，周鑫韜，岳慧賢，王述超，陳 騰，扈榮良

(軍事科學(xué)院 軍事醫(yī)學(xué)研究院 軍事獸醫(yī)研究所，吉林 長春 130122)

非洲豬瘟(African swine fever，ASF)是由非洲豬瘟病毒(African swine fever virus，ASFV)感染引起的一種主要以家豬急性出血性疾病為特征的高感染率和高死亡率的烈性傳染病，主要臨床癥狀為發(fā)熱、厭食、腹瀉或便秘、臥地、皮膚發(fā)紺和全身性出血等。ASF于1909年首先發(fā)生于東非地區(qū)[1]，隨后擴(kuò)散至歐洲、南美洲國家[2]，21世紀(jì)初擴(kuò)散至東歐國家和地區(qū)[3]，但大部分國家均通過撲殺政策和嚴(yán)格的生物安全措施實(shí)現(xiàn)了家豬ASF的消除[4-5]，除了非洲地區(qū)和意大利薩丁島呈地方性流行以外，歐洲地區(qū)的野豬群中ASF自流行以來持續(xù)存在[6]。2018年8月，我國首次報(bào)告家豬ASF[7-8]，隨后在各省區(qū)相繼暴發(fā)疫情，流行速度較快，死亡率較高，曾一度導(dǎo)致我國生豬、豬肉及相關(guān)產(chǎn)品出現(xiàn)供應(yīng)緊張局面。亞洲其他一些國家也陸續(xù)報(bào)道ASF流行情況[5]。目前，在我國ASF常態(tài)化防控政策下[9]，各養(yǎng)殖企業(yè)紛紛采取“拔牙”等方法控制疫情或采用高樓養(yǎng)豬、網(wǎng)格化管理等不同生物安全防范措施防止ASF發(fā)生。按照國家統(tǒng)計(jì)局和農(nóng)業(yè)農(nóng)村部發(fā)布的數(shù)據(jù)，至2020年9月，我國能繁殖母豬數(shù)量已達(dá)到3 822萬頭，能繁母豬逐漸恢復(fù)[10]。

經(jīng)過2年流行，ASF疫情在我國大多數(shù)地區(qū)已不再呈現(xiàn)暴發(fā)式發(fā)生，流行趨勢逐漸放緩，但在不少地區(qū)出現(xiàn)了臨床癥狀表現(xiàn)趨多，死亡率逐漸降低，發(fā)病不易察覺等情況。

ASFV是一種線性的有囊膜的雙鏈DNA病毒，是非洲豬瘟病毒科、非洲豬瘟病毒屬的唯一成員，也是唯一的媒傳DNA病毒[5]。病毒全基因組長170～194 kb[11-12]，變異程度較高。基因組編碼160多種蛋白。根據(jù)p72結(jié)構(gòu)蛋白基因3′端序列可將ASFV分為24個(gè)基因型[13]，不同基因型或同一基因型的不同分離株毒力差異較大[14]，感染豬表現(xiàn)超急性型、急性型、亞急性型、慢性型和亞臨床感染癥狀[15]。在ASF呈現(xiàn)地方性流行的非洲和意大利薩丁島等地區(qū)均發(fā)現(xiàn)了不同基因型或不同毒力株，ASFV自然變異株不斷有報(bào)道，如OUR T88/3[16](Ⅰ型，低毒力)、NH/P68[17](Ⅰ型，低毒力)、Lv17/WB/Rie1[18](Ⅱ型，低毒力)、Estonia 2014[19](Ⅱ型，低毒力)都是從田間分離得到的自然變異株，這些毒株均殘留某種程度的毒力，可引起如關(guān)節(jié)炎、皮膚潰瘍等?？梢灶A(yù)測，隨著我國ASF流行時(shí)間的延長，ASFV自然變異株也將不可避免地陸續(xù)出現(xiàn)。

為了充分了解ASF在我國的流行特點(diǎn)，我們進(jìn)行了ASF臨床調(diào)研，并開展了病原生態(tài)學(xué)監(jiān)測，在疫情流行特征與早期流行表現(xiàn)不同的地區(qū)的自由菜市場采集了樣品。采用世界衛(wèi)生組織(OIE)推薦的實(shí)時(shí)熒光定量PCR (qPCR)方法對組織樣品進(jìn)行檢測，進(jìn)而采用p72(B646L)及CVR、p30(CP204L)、EP153R-EP402R和MGF505-3R等基因引物對陽性樣品做進(jìn)一步擴(kuò)增和測序分析，并進(jìn)行了病毒分離和全基因組測序。結(jié)果表明，我國可能已經(jīng)出現(xiàn)了不同于早期ASFV流行株的毒力自然變異株。

1 材料與方法

1.1 細(xì)胞原代豬肺泡巨噬細(xì)胞(PAM)，無菌采自無豬源性病毒污染的家豬肺臟組織。

1.2 樣品帶骨的豬肌肉組織，分別購買自吉林、河北、湖北3個(gè)省的自由菜市場。

1.3 主要試劑Premix Taq(Ex Taq Version 2.0 plus dye)和Premix Ex Taq(Probe qPCR)為TaKaRa公司產(chǎn)品；RPMI 1640 培養(yǎng)基為Gibco公司產(chǎn)品；特級胎牛血清為呼和浩特草原綠野生物公司產(chǎn)品；裂解液由本實(shí)驗(yàn)室自制；FITC-p30單克隆抗體由本實(shí)驗(yàn)室制備。

1.4 主要儀器組織研磨儀為上海凈信實(shí)業(yè)發(fā)展有限公司產(chǎn)品；Eco實(shí)時(shí)熒光定量PCR儀為Illumina公司產(chǎn)品。

1.5 樣品處理用滅菌的不銹鋼藥匙從帶骨豬肉的骨髓中刮取黃豆粒大小樣品，加入500 μL PBS(0.01 mol/L，pH7.2)，在組織研磨儀中研磨30 s(100次/s)，離心后取上清液備用。取50 μL上清液加入50 μL裂解液，沸水中5 min，12 000×g離心1 min，取上清進(jìn)行后續(xù)PCR擴(kuò)增。

1.6 引物qPCR和PCR擴(kuò)增各基因所對應(yīng)的引物信息見表1。

1.7 PCR擴(kuò)增

1.7.1反應(yīng)體系和反應(yīng)條件 針對ASFV診斷時(shí)，采用OIE推薦的qPCR方法；針對p72、p30、CVR、EP153R-EP402R、MGF505-3R基因測序時(shí)，采用普通PCR方法。qPCR檢測出的陽性樣品需進(jìn)一步進(jìn)行普通PCR擴(kuò)增和產(chǎn)物測序。qPCR和普通PCR反應(yīng)體系均為20 μL：Premix Taq(Ex Taq Version 2.0 plus dye)或Premix Ex Taq(Probe qPCR)10 μL，上下游引物(10 μmol/L)各1 μL，模板2 μL，滅菌蒸餾水補(bǔ)足至20 μL。qPCR反應(yīng)條件：95℃ 30 s；95℃ 5 s，60℃ 45 s，40個(gè)循環(huán)。PCR反應(yīng)條件：95℃ 5 min；95℃ 30 s，57℃ 30 s，72℃ 30～150 s，35個(gè)循環(huán)；72℃延伸10 min。

1.7.2結(jié)果判定 qPCR結(jié)果判定：Ct值≤38 判為ASFV陽性，無Ct值判為陰性。38

PCR結(jié)果判定：PCR產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳，OIE-p72引物擴(kuò)增產(chǎn)物出現(xiàn)278 bp、CVR引物擴(kuò)增產(chǎn)物出現(xiàn)763 bp、p30引物擴(kuò)增產(chǎn)物出現(xiàn)344 bp、MGF505-3R引物擴(kuò)增產(chǎn)物出現(xiàn)357 bp、EP153R-EP402R引物擴(kuò)增產(chǎn)物出現(xiàn)2 250 bp特異性條帶時(shí)，判定為陽性，否則，判為陰性。各引物擴(kuò)增產(chǎn)物大小見表1。

表1 ASFV基因檢測引物信息

1.8 病毒培養(yǎng)和病毒特性測定

1.8.1病毒培養(yǎng) qPCR檢測為陽性和疑似陽性的樣品上清液，按5%比例接種原代PAM細(xì)胞，在37℃、5%CO2培養(yǎng)箱中培養(yǎng)96～120 h。培養(yǎng)過程中，每日顯微鏡下觀察。待細(xì)胞出現(xiàn)病變后，將細(xì)胞培養(yǎng)皿(瓶)在-20℃以下的低溫環(huán)境中凍融2次，分裝保存。

1.8.2病毒含量測定 將PAM細(xì)胞鋪于96孔板，2.5×106個(gè)/孔，在37℃、5%CO2培養(yǎng)箱培養(yǎng)。將0.1 mL 病毒液與0.9 mL細(xì)胞培養(yǎng)液混合，再取混合液0.1 mL與0.9 mL細(xì)胞培養(yǎng)液混合，依次獲得7個(gè)倍比稀釋液，每個(gè)稀釋度加入8個(gè)PAM細(xì)胞孔，每孔0.1 mL。37℃、5%CO2條件下培養(yǎng)3 d后，使用80%冷丙酮、4℃條件下固定細(xì)胞30 min，棄丙酮，晾干，加入0.1 mL 工作濃度的FITC熒光標(biāo)記的p30單克隆抗體稀釋液，在37℃避光孵育1 h；使用PBS(含0.5% Tween 20)洗滌3次，靜置3 min/次，熒光顯微鏡下觀察。病毒培養(yǎng)和含量檢測的操作均在BSL-3實(shí)驗(yàn)室進(jìn)行。

1.8.3血吸附特性 根據(jù)OIE推薦方法[20]，略作改良，采用原代PAM細(xì)胞進(jìn)行血吸附特性試驗(yàn)。簡言之，將100 μL PAM細(xì)胞鋪于96孔板，1×106個(gè)/孔，2 h后，將細(xì)胞培養(yǎng)物50 μL接種于PAM細(xì)胞(5 MOI)的96孔板上，同時(shí)設(shè)ASFV SY18毒株(本實(shí)驗(yàn)室2018年分離自國內(nèi)首次ASF疫情，為強(qiáng)毒株，具有血吸附特性)感染PAM細(xì)胞和正常PAM細(xì)胞對照。感染24 h后，將0.5%豬紅細(xì)胞按30 μL/孔加入96孔板，孵育0.5～3.0 d，顯微鏡下觀察。

1.9 擴(kuò)增產(chǎn)物測序、全基因組測序、驗(yàn)證和分析將普通PCR擴(kuò)增的陽性產(chǎn)物送至吉林省庫美生物科技有限公司測序。利用在線網(wǎng)站(https://blast.ncbi.nlm.nih.gov/Blast.cgi)對測序結(jié)果進(jìn)行同源性搜索和比對，分析核苷酸序列的同源性。對新分離毒株的細(xì)胞培養(yǎng)物利用二代測序技術(shù)(北京諾禾致源生物信息科技有限公司)進(jìn)行病毒全基因組測序，與SY18毒株全基因組序列(MH766894.2)比對，分析新分離毒株是否出現(xiàn)其他變異位點(diǎn)。對于出現(xiàn)變異的核苷酸位點(diǎn)，根據(jù)其所在基因或核苷酸片段位置兩側(cè)的核苷酸序列合成引物，同上方法PCR擴(kuò)增后送吉林省庫美生物科技有限公司測序進(jìn)行驗(yàn)證，以確定突變。

2 結(jié)果

2.1 PCR檢測結(jié)果自吉林、河北、湖北3個(gè)省共隨機(jī)購買帶骨豬肉樣品127份，通過qPCR方法檢測到4份陽性樣品。進(jìn)一步采用p72(B646L)、p30(CP204L)、EP153R-EP402R、MGF505-3R基因引物進(jìn)行擴(kuò)增，其中，2～4號樣品的p72(B646L)、p30(CP204L)、MGF505-3R、EP153R-EP402R基因的擴(kuò)增條帶大小和測序結(jié)果與SY18毒株的完全一致。而樣品1中其他基因大小和序列與SY18毒株完全一致，但EP153R-EP402R基因的擴(kuò)增條帶大小與SY18毒株明顯不同，擴(kuò)增產(chǎn)物經(jīng)過測序，其EP153R-EP402R基因均發(fā)生部分缺失，共缺失1 252 bp。

M.DL5000 DNA Marker;P.ASFV SY18毒株;1～4.樣品1～4

2.2 病毒分離及生長特性利用PAM細(xì)胞針對樣品1進(jìn)行了分離培養(yǎng)，結(jié)果接種的PAM細(xì)胞出現(xiàn)了細(xì)胞腫脹、空泡等細(xì)胞病變，獲得P1～P6代次細(xì)胞培養(yǎng)物，且第2，3代病毒含量分別為103.5，106.75TCID50/mL，隨后逐漸升高，第6代時(shí)病毒滴度達(dá)107.25TCID50/mL，與SY18毒株細(xì)胞培養(yǎng)的病毒含量相似。

2.3 全基因組測序及驗(yàn)證結(jié)果分析通過第2代測序技術(shù)，對第6代細(xì)胞培養(yǎng)物提取基因組進(jìn)行測序，獲得ASFV基因組總長為188 643 bp。與國內(nèi)提交GenBank的全部ASFV毒株相比，新分離毒株的基因組核苷酸變化主要表現(xiàn)在3個(gè)方面的差異：(1)發(fā)生EP153R基因和EP402R基因的部分大片段缺失，該區(qū)域與國內(nèi)ASFV毒株的比對結(jié)果見圖2，可以看出，EP153R基因?yàn)?′端部分缺失，CD2v為5′端部分缺失，缺失總長度為1 252 bp。(2)在多個(gè)基因編碼區(qū)和基因間隔區(qū)出現(xiàn)堿基突變、缺失或插入(表2)，這些突變和國內(nèi)已經(jīng)發(fā)表的10個(gè)ASFV毒株基因組中的堿基均不相同，具有本毒株特異性，其中部分突變導(dǎo)致編碼氨基酸的改變乃至開放閱讀框(ORF)的改變，后者包括多肽鏈的截短或延長，或失去編碼蛋白的功能。(3)基因組3′末端503 bp串聯(lián)重復(fù)片段的插入，與國內(nèi)毒株比對結(jié)果見圖3。

基于以上結(jié)果，故認(rèn)為該毒株為我國ASFV的1個(gè)新變異株，故將其命名為HuB20株。與ASFV SY18毒株相比，其基因組結(jié)構(gòu)簡圖如圖4所示。

2.4 血吸附特性HuB20株和ASFV SY18毒株分別感染PAM細(xì)胞 24 h后，在96孔板接種細(xì)胞的每孔中加入30 μL 0.5%豬紅細(xì)胞懸液，繼續(xù)放置37℃、5%CO2培養(yǎng)箱中培養(yǎng)。12 h后，感染ASFV SY18毒株的PAM細(xì)胞即出現(xiàn)了紅細(xì)胞吸附現(xiàn)象，感染HuB20株的PAM細(xì)胞和正常細(xì)胞未觀察到紅細(xì)胞吸附。至加入紅細(xì)胞72 h，感染HuB20株的PAM細(xì)胞和正常細(xì)胞仍無紅細(xì)胞吸附(圖5)。

3 討論

自2020年初以來，我國ASF不再呈大規(guī)模暴發(fā)流行，多表現(xiàn)疫情減緩，臨床癥狀不明顯，發(fā)病豬場死亡率30%～40%。根據(jù)農(nóng)業(yè)農(nóng)村部通報(bào)，2020年下半年以來，我國ASF疫情僅有零星發(fā)生，說明我國ASF的大流行可能已經(jīng)過去。由于一些疫情發(fā)生地ASF臨床癥狀大多不明顯，死亡率不高，與急性型ASF相比，在流行早期難以發(fā)現(xiàn)，采用之前的所謂“拔牙”等手段很難將ASF清除，場內(nèi)感染較難控制。本實(shí)驗(yàn)室根據(jù)承擔(dān)的課題任務(wù)開展了ASFV的主動監(jiān)測，在我國某些ASF疫情發(fā)生地區(qū)的自由菜市場購買了帶骨的豬肉樣品，經(jīng)檢測證明樣品具有一定比例的帶毒率。本研究中鑒定的自然變異毒株，其樣品來自湖北某地，但由于樣品來源于集市，故具體發(fā)病豬場不詳,不能完全排除其來源于進(jìn)口豬肉的可能。根據(jù)我們在當(dāng)?shù)卦儐柕那闆r，樣品來源地為ASF流行緩慢、病死率較低的地區(qū)。針對該毒株毒力的驗(yàn)證研究正在進(jìn)行中。

p72基因末端的測序結(jié)果表明，4個(gè)陽性樣品均為ASFV基因Ⅱ型，p72、p30等基因序列未發(fā)生改變，但其中樣品1(命名為HuB20株)CD2v和EP153R基因處發(fā)生了缺失。從缺失的位點(diǎn)來看，CD2v為5′端部分缺失，其上游的EP153R基因也發(fā)生了3′端部分缺失，缺失總長度為1 252 bp。這和目前已知實(shí)驗(yàn)室在研疫苗的CD2v缺失區(qū)域不同，在該處無熒光蛋白基因插入或人為設(shè)計(jì)的核苷酸序列。為了對該流行株進(jìn)行更為全面的分析，我們開展了全基因組二代測序。與國內(nèi)提交的不同株ASFV的全基因組比對后發(fā)現(xiàn)，該毒株有多處獨(dú)有的堿基變化。除了在EP153R、CD2v處發(fā)生缺失外，第2個(gè)明顯的變化是位于基因組3′末端的503 bp串聯(lián)重復(fù)序列的插入，該插入片段與國內(nèi)毒株Wuhan2019-1基因組5′末端的部分序列高達(dá)77%同源性，與國內(nèi)其他毒株基因組3′末端達(dá)28%同源性，與很多國外ASFV基因組的兩端序列高度同源，有的高達(dá)96%；第3個(gè)明顯的變化是在A238L、I9R等多個(gè)基因的編碼區(qū)出現(xiàn)堿基突變、缺失或插入，導(dǎo)致編碼氨基酸的改變乃至開放閱讀框變化。這些缺失、插入、突變是否和國內(nèi)報(bào)道的疫情狀況相關(guān)聯(lián)，需要進(jìn)一步通過動物試驗(yàn)進(jìn)行驗(yàn)證。

灰色區(qū)域?yàn)橥葱?00%，藍(lán)色區(qū)域?yàn)闊o同源性；虛線區(qū)域?yàn)镠uB20株獨(dú)有序列，與其他序列無同源性

該自然變異株體外表型的變化導(dǎo)致其血吸附特性喪失。ASFV的血吸附特性據(jù)報(bào)道與CD2v有關(guān)。CD2v由EP402R基因編碼，與T細(xì)胞CD2黏附受體同源，直接參與了ASFV在感染細(xì)胞中誘導(dǎo)的血吸附現(xiàn)象[22]。EP153R基因的表達(dá)產(chǎn)物C型凝集素由153個(gè)氨基酸組成，包含1個(gè)中央跨膜區(qū)、1個(gè)C型動物凝集素樣結(jié)構(gòu)域和1個(gè)細(xì)胞附著(RGD)序列，與人及動物CD44分子氨基酸序列具有較高同源性。有報(bào)道認(rèn)為，CD2v的血吸附特性的穩(wěn)定性似乎與C型凝集素有關(guān)聯(lián)[23]。文獻(xiàn)中報(bào)道的血吸附特性喪失的自然變異株包括NH/P68、OUR T88/3、Lv/17/WB/Rie1等。其中，OUR T88/3是2004年從西班牙豬場軟蜱中分離得到的，為基因Ⅰ型，其缺失核苷酸片段較大，與本自然變異株不屬于同一基因型[24]。NH/P68為慢性感染豬場中分離到的1株弱毒，也是基因Ⅰ型，其血吸附特性的喪失是因CD2v基因中堿基突變導(dǎo)致讀框發(fā)生移位所致[17]。Lv/17/WB/Rie1株是2017年從野豬體內(nèi)分離到的一個(gè)基因Ⅱ型的自然弱毒株，無血吸附特性，則是由1個(gè)核苷酸缺失導(dǎo)致讀框移位提前終止所致[18]。本變異株血吸附特性的喪失與國外報(bào)道的毒株均不一樣，是因基因片段的缺失所致。因此，據(jù)我們所知，本變異株(HuB20株)應(yīng)為我國境內(nèi)首次出現(xiàn)和鑒定的無血吸附特性及毒力可能致弱的自然變異株。

灰色區(qū)域?yàn)樗卸局晖葱?00%，藍(lán)色區(qū)域?yàn)镠uB20株獨(dú)有序列;虛線區(qū)域?yàn)槠渌局晖磪^(qū)域

表2 HuB20株單(或多重)核苷酸變異統(tǒng)計(jì)表

續(xù)表2

圖4 新分離ASFV HuB20株與國內(nèi)毒株基因組差異

圖5 ASFV HuB20變異株和SY18毒株血吸附結(jié)果

HuB20株在體外的復(fù)制能力未受到影響。GLADUE等[25]在ASFV-G-Δ9GL基礎(chǔ)上對CD2v和EP153R進(jìn)行缺失，缺失毒的生長曲線顯示，ASFV-G-Δ9GL/ΔCD2v/ΔEP153R在PAM細(xì)胞中復(fù)制能力低于ASFV-G-Δ9GL和ASFV-G-Δ9GL/ΔCD2v，后兩者的復(fù)制能力無差異，三者的復(fù)制能力均低于親本毒株Georgia 2007；說明9GL和EP153R基因?qū)Σ《緩?fù)制可能具有重要影響。但本實(shí)驗(yàn)室在對HuB20株培養(yǎng)過程中發(fā)現(xiàn)，最初3代病毒繁殖滴度較低，后期病毒增殖滴度和SY18毒株沒有明顯差別，在沒有9GL基因缺失的情況下，EP402R、EP153R的聯(lián)合缺失對于病毒的復(fù)制并未產(chǎn)生明顯影響。但由于本變異株中存在多處堿基變化，不能排除其他基因的改變對本病毒的體外復(fù)制產(chǎn)生了有利影響。

自然變異株的產(chǎn)生，可能系多方面原因促成。首先，我國生豬養(yǎng)殖基數(shù)大，品種多，ASFV在不同種群中長期存在，為病毒在不同易感性豬品種內(nèi)的適應(yīng)和變異提供了機(jī)會；第二，隨著ASF大流行的過去，耐過豬數(shù)量增多，在抗體和免疫選擇的壓力下，病毒出現(xiàn)變異是可能的；第三，ASFV本身的特性決定了其易變的可能。其基因組長度變異較大，不同毒株為170～194 kb，一般認(rèn)為其中間125 kb為恒定區(qū)，左側(cè)38～47 kb和右側(cè)13～15 kb為可變區(qū)，病毒自然發(fā)生變異的頻率較高[26]。本研究證明的自然變異株僅為我國某一地區(qū)的分離株，可能是眾多變異株中的1株，不能排除在我國其他地區(qū)出現(xiàn)了其他變異株的可能。因此，有必要加強(qiáng)監(jiān)測，根據(jù)流行和變異株的特征來調(diào)整和確定防控策略。

另外，根據(jù)ASFV的遺傳特性，在不實(shí)行撲殺的常態(tài)化疫情防控措施下，尤其是目前很多養(yǎng)殖場通過“拔牙”的方法，或使用免疫增強(qiáng)藥物，或使用基因工程改造的微生態(tài)制劑以及試圖使用治療藥物控制疫情的情況下，都將會隨著流行時(shí)間的延長以及對野毒株流行的“縱容”，加速變異株的產(chǎn)生。因此，加強(qiáng)對ASF流行病學(xué)和分子流行病學(xué)的監(jiān)測，應(yīng)當(dāng)是我們每個(gè)獸醫(yī)從業(yè)者尤其是研究人員的責(zé)任。

自然變異株的出現(xiàn)，都是病毒為了自身更好地生存而作出的繁殖策略的改變，與體內(nèi)復(fù)制有限、不能水平傳播的人工缺失致弱株相比，自然變異株在豬體內(nèi)更容易生存和繁殖，也更容易在豬群中散播。自然變異株由于致病力發(fā)生一定程度減弱，其傳播不易引起注意，更容易造成大面積擴(kuò)散，對其根除也將變得困難。因此，其危害性與引起最急性和急性型癥狀的強(qiáng)毒株相比更加隱蔽，因而更加廣泛和嚴(yán)重，值得我國養(yǎng)豬行業(yè)的關(guān)注。應(yīng)當(dāng)毫不猶豫地清除自然變異株感染豬群，只有這樣，通過生物安全清除ASF才不是空談。