李本強(qiáng),劉惠莉,劉佳佳,陶 潔,程靖華,石 迎,喬長(zhǎng)濤,沈曉暉*

(1 上海市農(nóng)業(yè)科學(xué)院畜牧獸醫(yī)研究所∕上海市農(nóng)業(yè)遺傳育種重點(diǎn)實(shí)驗(yàn)室,上海 201106;2 上海市種豬工程技術(shù)中心,上海 201106;3 上海海洋大學(xué),上海 201306)

豬薩佩羅病毒(Porcine sapelovirus,PSV)屬小RNA 病毒科,是微球形、無(wú)囊膜、單股正鏈RNA 病毒的一種,為腸道病毒屬II 群、豬腸道病毒8 型[1]。 通過對(duì)腸道病毒基因組及進(jìn)化關(guān)系進(jìn)一步研究,TSENG等[2]將該病毒劃為薩佩羅病毒屬。 PSV 基因組為單股正鏈RNA,全長(zhǎng)約7 400 bp,包含1 個(gè)大的開放閱讀框架,編碼1 個(gè)多聚蛋白前體,在5’端有一個(gè)Vpg 區(qū),3’末端連有Poly A 結(jié)構(gòu)[3-4]。

PSV 于20 世紀(jì)60 年代在英國(guó)首先發(fā)現(xiàn)[5],隨后在北美、日本、澳大利亞等地先后出現(xiàn)[6-10]。 PSV在飼養(yǎng)豬群中感染率很高,多呈散發(fā),不同年齡的豬均有易感性[4,11],幼齡豬易感性較強(qiáng),病豬、康復(fù)豬和隱性感染豬為主要傳染源。 感染PSV 后可導(dǎo)致豬腸道、呼吸道、生殖器官疾病以及腦脊髓炎等多系統(tǒng)綜合征[12-13],還可能呈現(xiàn)亞臨床感染。 但無(wú)論其以何種形式爆發(fā),都會(huì)給養(yǎng)豬業(yè)帶來(lái)潛在威脅。 目前對(duì)PSV 的研究主要集中于病原分離與鑒定、基因組與結(jié)構(gòu)蛋白分析方面。 本研究通過監(jiān)測(cè)上海地區(qū)豬場(chǎng)PSV 感染情況,并對(duì)其中2 個(gè)PSV 分離株進(jìn)行全基因組測(cè)定與分析,旨在進(jìn)一步明確PSV 在腹瀉疫情中的作用。

1 材料與方法

1.1 材料

1.1.1 樣品采集

于2016—2020 年從上海地區(qū)11 個(gè)豬場(chǎng)采集具有腹瀉癥狀的豬糞便樣品1 446 份,用PBS 液將樣品按照1∶10 的比例稀釋,-80 ℃冰箱凍存。

1.1.2 主要試劑及儀器

Trizol、M-MLV RNA 聚合酶、PCR 試劑及pMD-18T 載體等均購(gòu)自TaKaRa 公司;5’RACE kit 試劑盒購(gòu)自北京百泰克生物技術(shù)有限公司;pEASY-Blunt Zero Cloning Kit、DNA 凝膠回收試劑盒、Trans1-T1 Phage Resistant Chemically Competent Cell 購(gòu)自北京全式金生物公司;Phanta Max Super-Fidelity DNA Polymerase購(gòu)自諾唯贊生物公司。

PCR 儀和核酸電泳儀為美國(guó)Bio-Rad 公司生產(chǎn)。

1.2 方法

1.2.1 樣品的處理及總RNA 的提取

將-80 ℃保存的樣品漩渦震蕩5 min,取250 μL 溶液,采用Trizol 法提取總RNA,具體操作方法按照TaKaRa 產(chǎn)品說明書進(jìn)行。

1.2.2 引物的設(shè)計(jì)及合成

根據(jù)Genbank 中HQ875059、KJ821021、KJ821020 和KX574284 序列,在3D 保守區(qū)域設(shè)計(jì)1 對(duì)PSV 特異檢測(cè)引物(PSV-F:5’-GATGTGGCGCATGCTCTT-3’,PSV-R:5’-TGCTGCCTCCTGTGTTGTTAT-3’),擴(kuò)增片段長(zhǎng)度約為624 bp[14]。 采用實(shí)驗(yàn)室已公布的引物用于擴(kuò)增PSV 全基因組序列,引物序列由鉑尚生物技術(shù)(上海)有限公司合成。

1.2.3 臨床樣品的RT-PCR 檢測(cè)

以提取的樣品RNA 為模板,按照Reverse Transcriptase M-MLV 的說明書進(jìn)行反轉(zhuǎn)錄,合成cDNA。RT-PCR 反應(yīng)體系為:cDNA 1 μL,10 × PCR mix 2 μL,PSV-F(0. 5 μmol∕L)和PSV-R(0. 5 μmol∕L)各0.5 μL,用水補(bǔ)至20 μL;反應(yīng)程序?yàn)?94 ℃4 min;94 ℃30 s;52 ℃30 s;72 ℃35 s,30 個(gè)循環(huán);72 ℃8 min;4 ℃結(jié)束反應(yīng)。

1.2.4 病毒的增殖及鑒定

選取生長(zhǎng)狀態(tài)良好的PK15 細(xì)胞,胰酶消化后鋪于6 孔板中。 第2 天,棄去細(xì)胞上清,用PBS 洗滌2遍,加入含0.5 μg∕mL 胰酶的DMEM 病毒培養(yǎng)液;吸附2 h 后,棄去上清,加入含胰酶的DMEM 培養(yǎng)液。置于37 ℃培養(yǎng)箱,觀察細(xì)胞病變情況,若無(wú)病變,則5 d 后收毒;凍融2 次后,重復(fù)上述過程,進(jìn)行病毒的連續(xù)傳代培養(yǎng)。 不同代次的細(xì)胞培養(yǎng)物各取250 μL,進(jìn)行RT-PCR 檢測(cè),驗(yàn)證病毒傳代成功與否。

1.2.5 間接免疫熒光檢測(cè)

將傳至第5 代的PSV 接種PK15 細(xì)胞,吸附1 h 后,吸掉病毒液,用DMEM 洗滌2 次。 加入含有0.5 μg∕mL TPCK 的DMEM 培養(yǎng)基于37 ℃培養(yǎng)72 h 后,棄去細(xì)胞上清,PBST 洗1 遍,70%乙醇冰上固定15 min;PBST 洗3 遍并甩干,加入1∶1 000 稀釋的PSV 單抗,37 ℃濕盒中孵育1 h;重復(fù)洗滌,加入1∶3 000稀釋的FITC 標(biāo)記羊抗鼠IgG,37 ℃孵育30 min;洗滌后,置于倒置熒光顯微鏡(Axio Observer)下觀察。

1.2.6 PSV 全基因組序列擴(kuò)增及序列分析

取第5 代PSV 的病毒液,擴(kuò)增全基因組序列。 反應(yīng)體系參照Phanta Max Super-Fidelity DNA Polymerase 說明書進(jìn)行PCR 反應(yīng)。 反應(yīng)程序:94 ℃4 min;94 ℃30 s,52 ℃30 s,72 ℃60 s,35 個(gè)循環(huán);72 ℃8 min;4 ℃結(jié)束反應(yīng)。

PSV 的5’末端序列擴(kuò)增采用5’RACE 方法進(jìn)行[15]。 將目的片段與pEASY-Blunt Zero Cloning Kit 載體進(jìn)行連接,隨后用Trans1-T1 Phage Resistant Chemically Competent Cell 進(jìn)行克隆,挑取陽(yáng)性重組質(zhì)粒進(jìn)行測(cè)序。

采用Lasergene 分析軟件對(duì)各測(cè)定片段進(jìn)行拼接,獲得PSV 全基因組序列并進(jìn)行基因分析,利用MEGA 6.0軟件制作進(jìn)化樹。

2 結(jié)果與分析

2.1 PSV 流行情況調(diào)查

對(duì)2016—2020 年采集的1 416 份糞便樣品進(jìn)行RTPCR 檢測(cè),結(jié)果顯示,共有256 份PSV 陽(yáng)性樣品,總陽(yáng)性率 為 17.7%。 各 年 PSV 陽(yáng) 性 率 分 別 為 18.26%(65∕356)、3.03% (3∕99)、18.72% (44∕235)、21.58%(126∕584)和10.47%(18∕172)(圖1),表明PSV 在豬群中普遍存在。

圖1 豬群中PSV 的年度流行情況Fig.1 Prevalence trends of PSV infection in pig herds

2.2 PSV 病毒分離

將PCR 檢測(cè)為陽(yáng)性的PSV 樣品接種PK15 細(xì)胞,連續(xù)傳5 代后,2 株P(guān)SV 上出現(xiàn)細(xì)胞病變。 PCR 檢測(cè)顯示其第1 代次、第2 代次的擴(kuò)增條帶都較弱,從第3 代開始,病毒增殖比較穩(wěn)定。 將成功分離的2 株P(guān)SV 毒株分別命名為CM2019008 和CM2019337(圖2)。

圖2 分離株的PCR 鑒定Fig.2 PCR identification of the isolated strain

2.3 病毒間接免疫熒光鑒定

將第5 代的CM2019008 和CM2019337 病毒分別接種于PK15 細(xì)胞,間接免疫熒光檢測(cè)顯示,2 株分離病毒均能與PSV 自制單抗發(fā)生免疫反應(yīng),出現(xiàn)特異性綠色熒光(圖3),表明這兩個(gè)毒株在PK15 細(xì)胞上增殖成功。

圖3 PSV 分離株的間接免疫熒光檢測(cè)(200 ×)Fig.3 Indirect immunofluscent test of PSV isolates(200 ×)

2.4 PSV 基因組遺傳分析

2.4.1 PSV 全基因組的遺傳進(jìn)化分析

利用Mega 6.0 軟件對(duì)2 株P(guān)SV 和國(guó)內(nèi)外其他代表性毒株進(jìn)行全基因組遺傳進(jìn)化分析。 結(jié)果表明,分離株CM2019337 與上海PSV 株MN685785.1 進(jìn)化關(guān)系最近,分離株CM2019008 與SHPD202150 在中國(guó)PSV 的一個(gè)大的分支上。 中國(guó)PSV 分離株與日本2 株P(guān)SV 分離株同處一個(gè)大的分支上,遺傳進(jìn)化關(guān)系最為接近(圖4)。

圖4 PSV 全基因組序列構(gòu)建的系統(tǒng)進(jìn)化樹Fig.4 The phylogenetic tree constructed based on the whole gene sequence of PSV

進(jìn)一步利用DNAStar 軟件對(duì)上述毒株的核苷酸相似性進(jìn)行分析,發(fā)現(xiàn)CM2019008 和CM2019337 的核苷酸相似性為94.8%,其次是MN685785.1,兩者核苷酸相似性為94.4%;CM2019337 和MN685785.1的核苷酸相似性為99.5%,其次是JX286666.1,核苷酸相似性為89.8%。 中國(guó)的PSV 毒株與韓國(guó)PSV 毒株核苷酸相似性為87%—88.5%,與歐洲PSV 毒株核苷酸相似性為84.9%—87.3%(圖5)。

圖5 蛋白編碼區(qū)序列相似性比較Fig.5 Homology comparison of sequence in protein coding region

2.4.2 分離株蛋白編碼區(qū)的重組分析

在SimPlot 3.5.1 軟件中CM2019008 和CM2019337 作為獨(dú)立的重組事件進(jìn)行標(biāo)準(zhǔn)相似節(jié)點(diǎn)分析以查找有無(wú)重組事件的發(fā)生,結(jié)果所示,CM2019008 和CM2019337 分離株在VP1、VP2 和VP3 區(qū)域與HeBO4核苷酸序列高度相似,而在3D 區(qū)域與ISU-SHIC 高度相似(圖6A)。 利用RDP、Chimaera、Bootscan、3Seq、Maxchi、Phylpro、GENECONV、LARD 及SISCAN 九種算法對(duì)重組毒株進(jìn)行分析以尋找重組毒株的潛在斷裂點(diǎn)。 Bootscan 分析結(jié)果所示,CM2019008 和CM2019337 分離株在3D 區(qū)域存在重組斷點(diǎn)(6B)。 利用鄰接法構(gòu)建重組毒株的進(jìn)化樹,結(jié)果如圖6C 所示,CM2019008 和CM2019337 分離株的重組區(qū)域位于4 696—5 945 bp 處。

圖6 潛在重組毒株多聚蛋白重組事件Fig.6 Recombinant events identified according to the polyprotein gene of potential recombinant strains

2.4.3 分離株VP1 基因的氨基酸序列分析

采用MegAlign 和GeneDoc 2.7 軟件對(duì)VP1 基因?qū)?yīng)的氨基酸序列進(jìn)行分析,結(jié)果如圖7 所示。 與PSV 經(jīng)典株相比,CM2019008 共存在6 個(gè)氨基酸位點(diǎn)的突變,分別為第28、88、98、223、229、271 位,且突變位點(diǎn)多集中在VP1 氨基酸的兩端序列;CM2019337 共存在5 個(gè)氨基酸位點(diǎn)的突變,分別為第28、94、95、218、236 位。 CM2019008 和CM2019337 中間氨基酸99—217 位點(diǎn)相對(duì)保守。

圖7 PSV VP1 基因氨基酸突變位點(diǎn)的分析Fig.7 Analysis of amino acid mutation sites of PSV VP1 gene

2.4.4 VP1 蛋白的潛在糖基化位點(diǎn)和抗原表位的分析

運(yùn)用在線分析軟件NetNGlyc 1.0 Server 和GeneDoc 2.7 對(duì)VP1 蛋白的潛在糖基化位點(diǎn)進(jìn)行分析,結(jié)果表明所有PSV 分離株和參考株都含有潛在NSTS 糖基化位點(diǎn),且都含173 位潛在糖基化位點(diǎn)。 其中分離株CM2019008 和CM2019337 分別含有1 個(gè)和2 個(gè)潛在糖基化位點(diǎn);日本PSV 分離株LC326556.1 含有2 個(gè)糖基化位點(diǎn),德國(guó)PSV 分離株LT900497.1 和法國(guó)PSV 分離株MH513612.1 分別含有2 個(gè)和1 個(gè)糖基化位點(diǎn)。

采用在線分析軟件BepiPred 1.0 Server 和DNAStar Protean 對(duì)分離株VP1 蛋白的抗原表位進(jìn)行分析,發(fā)現(xiàn)中國(guó)分離株CM2019008 和CM2019337 與中國(guó)已有的分離株及國(guó)外分離株的抗原表位差異不大,抗原表位均勻分散在整個(gè)VP1 蛋白分子表面。

3 結(jié)論與討論

本研究利用RT-PCR 方法對(duì)上海地區(qū)腹瀉豬群中PSV 的攜帶情況進(jìn)行檢測(cè),發(fā)現(xiàn)PSV 總體陽(yáng)性率為17.7%,與Lan 等[16]在中國(guó)華東地區(qū)豬糞便中17.2%的PSV 檢出率、Chen 等[17]在中國(guó)華南地區(qū)豬糞便中18.21%的PSV 檢出率、Li 等[18]在四川地區(qū)豬糞便中20.4%的PSV 檢出率相近。 表明PSV 在上海地區(qū)豬場(chǎng)普遍存在。

本研究從陽(yáng)性樣品中分離得到兩株P(guān)SV。 對(duì)兩株P(guān)SV 全基因組序列進(jìn)行分析,發(fā)現(xiàn)其與日本部分流行PSV 毒株之間親緣關(guān)系最近,與韓國(guó)的分離株親緣關(guān)系較遠(yuǎn)。 我國(guó)的PSV 分離株呈簇聚集在同一分支中,說明我國(guó)的流行毒株可能存在獨(dú)立進(jìn)化。 核苷酸相似性分析發(fā)現(xiàn)PSV 分離株與韓國(guó)PSV 毒株核苷酸相似性為87%—88.5%,與歐洲PSV 毒株核苷酸相似性為84.9%—87.3%。 病毒突變和重組是病毒進(jìn)化的主要方式,在進(jìn)化過程中扮演著至關(guān)重要的角色。 本研究發(fā)現(xiàn),CM2019008 和CM2019337 分離株的重組區(qū)域位于4 696—5 945 bp 處,此區(qū)域編碼病毒特異性RNA 聚合酶,具有合成病毒RNA 能力。

VP1 衣殼蛋白是PSV 最外部的蛋白,是誘導(dǎo)機(jī)體產(chǎn)生中和抗體的主要蛋白質(zhì),也是誘導(dǎo)豬場(chǎng)產(chǎn)生體液免疫應(yīng)答的主要成分。 目前,關(guān)于VP1 基因的氨基酸突變對(duì)生物學(xué)特性的意義還未見報(bào)道。 本研究對(duì)PSV 分離株和已報(bào)道的PSV 毒株進(jìn)行了序列分析,發(fā)現(xiàn)PSV 分離株VP1 基因的氨基酸較少發(fā)生突變,同時(shí)整個(gè)VP1 基因的氨基酸序列未發(fā)現(xiàn)插入和缺失現(xiàn)象,該結(jié)果對(duì)病毒防治具有指導(dǎo)意義。

蛋白質(zhì)糖基化修飾與很多重要的功能蛋白或酶相關(guān),特別是與蛋白質(zhì)的空間構(gòu)想、生物活性、運(yùn)輸和定位密切相關(guān),參與很多復(fù)雜的生理生化反應(yīng),如:分子識(shí)別、細(xì)胞通信、信號(hào)轉(zhuǎn)導(dǎo)等,因此糖基化位點(diǎn)分析對(duì)研究病毒機(jī)制具有十分重要的意義[19]。 本研究發(fā)現(xiàn)分離株P(guān)SV 和參考PSV 都含有潛在NSTS 糖基化位點(diǎn),其中分離株CM2019008 和法國(guó)分離株MH513612.1 含有1 個(gè)潛在糖基化位點(diǎn),CM2019337 和日本LC326556.1 及德國(guó)分離株LT900497.1 含有2 個(gè)潛在糖基化位點(diǎn)。 而所有的PSV 都含有173 NSTS 糖基化位點(diǎn)。 通過對(duì)分離株VP1 蛋白的抗原表位分析,發(fā)現(xiàn)分離株CM2019008 和CM2019337 與中國(guó)目前已有的分離株及國(guó)外分離株的抗原表位差異不大,這可能是VP1 蛋白保持良好抗原性的重要原因之一,研究結(jié)果對(duì)進(jìn)一步研究PSV 高感染率及分子機(jī)制具有指導(dǎo)意義。