王 迪,姬洪衛(wèi),王澤東,李曉慧,張 力,魏 峰,劉 全 (1.吉林農(nóng)業(yè)大學(xué) 生命科學(xué)學(xué)院,吉林 長春 130118;.軍事醫(yī)學(xué)研究院 軍事獸醫(yī)研究所,吉林 長春 1301;3.佛山科學(xué)技術(shù)學(xué)院 生命工程與科學(xué)學(xué)院,廣東 佛山 58000)
蜱傳腦炎病毒(tick-borne encephalitis virus,TBEV)是一種臨床上以突發(fā)高熱、腦膜刺激癥、頭痛嘔吐、感覺過敏、意識障礙及頸和肢體癱瘓為特征的嗜神經(jīng)病毒[1]。該病毒分布廣泛,主要流行于歐洲和亞洲地區(qū)。在我國主要流行于東北、西北和西南等林區(qū)。蜱為TBEV的主要傳播媒介,全溝硬蜱及蓖子硬蜱是TBEV的主要媒介蜱種[2]。
TBEV的基因組為單股正鏈RNA,含有單一閱讀編碼框,編碼3個結(jié)構(gòu)蛋白(C、PrM、E)和7個非結(jié)構(gòu)蛋白(NS1、NS2A、NS2B、NS3、NS4A、NS4B、NS5)。其中E蛋白為TBEV的重要結(jié)構(gòu)蛋白,該蛋白內(nèi)單個氨基酸的改變可影響病毒的毒力和血凝活性等[3]。因此,對于TBEV的分子流行病學(xué)研究主要集中于E蛋白中。
近年來,隨著氣候變暖、人類活動范圍擴大,我國森林腦炎又稱蜱傳腦炎(TBE)發(fā)病率有明顯上升趨勢[4]。本研究為了解我國東北地區(qū)蜱源TBEV的流行及毒株進(jìn)化情況,通過RT-PCR方法對蜱攜帶TBEV的陽性率及毒株序列突變情況進(jìn)行分析,為我國東北地區(qū)森林腦炎的科學(xué)防控奠定基礎(chǔ)。
1.1 蜱樣品采集與鑒定2015年4—6月在吉林省吉林市和黑龍江省雙鴨山市使用布旗法采集游離成年蜱,按照采集地點和形態(tài)學(xué)進(jìn)行分類,平均15只/管。GenBank下載已發(fā)表的硬蜱科、革蜱屬和血蜱屬的16SrRNA序列,利用Vector NTI軟件設(shè)計蜱種鑒定引物,由吉林庫美生物有限公司合成。
1.2 TBEV檢測及統(tǒng)計學(xué)分析使用TIANamp Virus RNA Kit 試劑盒提取病毒組RNA,TaKaRa PrimeScript RT reagent Kit和gDNA Eraser合成cDNA鏈。利用軟件Primer Premier5.0在TBEV E蛋白、保守區(qū)設(shè)計半巢式檢測引物,上游引物F1:5′-GTBACCATAACAGCYGAGG-3′,上游引物F2:5′-GCCTDCACGCMAAGCTRTCG-3′,下游引物TB-EV-R1:5′-CCRTAGTHTCCCATGGTCA-AG-3′,擴增片段大小為361 bp。利用軟件Pooled-InfRate version 4.0中最大似然估算法計算[5],分析吉林省吉林市和黑龍江省雙鴨山市3種蜱種的TBEV陽性率。
1.3 TBEV分離及電鏡觀察使用DMEM將TBEV陽性樣本稀釋20倍,接種到非洲綠猴腎細(xì)胞(Vero細(xì)胞),盲傳3代,觀察細(xì)胞狀態(tài)[6]。將產(chǎn)生病變的Vero細(xì)胞上清液負(fù)染電子顯微鏡檢測。
1.4 TBEV全基因組擴增根據(jù)TBEV檢測結(jié)果和GenBank已發(fā)表的TBEV序列,利用Primer Premier 5.0、Vector NTI軟件設(shè)計TBEV全基因組擴增引物,使用RT-PCR方法陽性樣本進(jìn)行全基因組擴增(表1)。
表1 TBEV全基因組擴增引物序列
1.5 TBEV系統(tǒng)發(fā)育分析以GenBank登錄號、病毒株命名、年份、國家為順序?qū)ο螺d的病毒進(jìn)行命名。利用MegAlign軟件將所得到的序列進(jìn)行拼接和同源性分析,MEGA 5.0軟件的鄰位相連法對TBEV全基因組和E蛋白的核苷酸序列構(gòu)建系統(tǒng)進(jìn)化樹,選擇Bootstrap method作為測試方法,重復(fù)次數(shù)設(shè)置為1 000次。
2.1 蜱鑒定結(jié)果本試驗從我國吉林省吉林市和黑龍江省雙鴨山市采集蜱共1 155只,其中吉林省吉林市192只,黑龍江省雙鴨山市963只。通過形態(tài)學(xué)和分子生物學(xué)鑒定,共得到3種蜱:全溝硬蜱655只(56.7%)、森林革蜱90只(34.7%)、嗜群血蜱99只(8.6%);其中,吉林省吉林市全溝硬蜱102只(53.1%)、草原革蜱90只(46.9%)。黑龍江省雙鴨山市全溝硬蜱553只(57.4%)、草原革蜱311只(32.3%)、嗜群血蜱99只(10.3%)。全溝硬蜱為兩地優(yōu)勢蜱種(圖1、表2)。
表2 蜱種采樣信息統(tǒng)計 只/%
圖1 樣品采集地點分布
2.2 蜱中TBEV攜帶率在本研究中,TBEV僅在全溝硬蜱中檢出,陽性率為4.7%。吉林省吉林市全溝硬蜱陽性率為20.8%,黑龍江省雙鴨山市全溝硬蜱陽性率為1.8%。未在嗜群血蜱、草原革蜱檢出TBEV(表3)。
表3 東北地區(qū)蜱攜帶TBEV陽性率
2.3 TBEV分離及電鏡觀察應(yīng)用Vero細(xì)胞成功分離出1株TBEV,電子顯微鏡下觀察細(xì)胞出現(xiàn)脫落、變圓現(xiàn)象。對細(xì)胞上清液進(jìn)行負(fù)染,通過透射電子顯微鏡觀察發(fā)現(xiàn),其病毒形態(tài)呈現(xiàn)典型的黃病毒形態(tài):有包膜,病毒顆粒直徑為30~40 nm(圖2)。
圖2 電子顯微鏡下病變的Vero細(xì)胞(A)和負(fù)染狀態(tài)下的病毒形態(tài)(B)
2.4 TBEV全基因組拼接結(jié)果本研究成功獲得1株吉林省吉林市蜱源TBEV毒株,命名為JL-75,其全基因組全長為11 095 bp。其中5′-UTR包括131個核苷酸,開放閱讀框位于基因組的132~10 376 bp,包括10 245個堿基,3′-UTR包括727個核苷酸。JL-T75全基因組已上傳至GenBank,登錄號為MN615726。
2.5 TBEV編碼區(qū)序列同源性比較將JL-7株與下載的其他基因型代表株的編碼區(qū)序列進(jìn)行比對發(fā)現(xiàn):JL-T75與遠(yuǎn)東亞型TBEV的同源性最高,核苷酸同源性為94.1%~98.3%,氨基酸同源性為99.1%~99.4%。其中,與我國疫苗株森張株核苷酸同源性為98.3%,氨基酸同源性為99.1%。而與西伯利亞亞型、歐洲亞型、貝加爾湖亞型、喜馬拉雅亞型TBEV的同源性較低,核苷酸同源性為83.4%~87.6%,氨基酸同源性為92.9%~96.3%(表4)。
表4 TBEV分離株與參考株編碼區(qū)核苷酸和氨基酸同源性比較 %
2.6 TBEV全基因組系統(tǒng)進(jìn)化樹分析根據(jù)全基因組序列和編碼E蛋白片段核苷酸序列對JL-75株進(jìn)行系統(tǒng)進(jìn)化樹分析,發(fā)現(xiàn)JL-75株與我國森張株(JQ650523)親緣關(guān)系較近(圖3)。
圖3 蜱源TBEV全基因組和E蛋白核苷酸序列的系統(tǒng)進(jìn)化樹分析
2.7 TBEV編碼區(qū)氨基酸差異位點比較將JL-75株與我國疫苗株森株進(jìn)行氨基酸差異位點比較,發(fā)現(xiàn)共有29個氨基酸位點產(chǎn)生突變。其中,E蛋白有2個氨基酸位點產(chǎn)生突變,分別為Ile-292→Val和Val-351→Ala(表5)。
表5 分離株與森張株氨基酸差異位點比較
自我國1942年發(fā)現(xiàn)TBEV以來,TBEV在我國人群和野鼠、鳥類、野豬、牛、羊等動物中均有發(fā)現(xiàn),嚴(yán)重威脅了我國人民群眾的公共衛(wèi)生安全和農(nóng)業(yè)經(jīng)濟發(fā)展[7]。我國東北地區(qū)森林資源豐富,為蜱和其宿主提供了良好的生活環(huán)境[8]。有調(diào)查顯示,我國東北部分地區(qū)為長期穩(wěn)定的TBEV自然疫源地[9]。蔡增林等[10]在黑龍江省密山市采集的嗜群血蜱和森林革蜱以及吉林琿春采集的嗜群血蜱中檢測到TBEV,陽性率分別為0.78%、0.77%和0.17%。趙俊偉等[11]在吉林省琿春市采集的森林革蜱中檢測到TBEV,陽性率為1.16%。在本研究中,僅在吉林省吉林市和黑龍省雙鴨山市的全溝硬蜱中檢出TBEV,陽性率分別為20.8%和1.8%,在草原革蜱和嗜群血蜱中未檢出。分析是由于我國東北地區(qū)人類活動區(qū)域的擴大和退耕還林政策的實施等因素,導(dǎo)致東北地區(qū)的優(yōu)勢蜱種發(fā)生了一定的變化,使個別TBEV易感蜱種的數(shù)量下降,相對密度變低。
TBEV根據(jù)基因組特征可分為3個型,即遠(yuǎn)東型、西伯利亞型和歐洲型。不同基因型之間毒力也有所不同,其中遠(yuǎn)東型致病性最強,致死率為5%~20%[4]。我國發(fā)現(xiàn)或分離的TBEV流行株均為遠(yuǎn)東亞型[12]。張桂林等[13]于2017年在我國新疆中哈邊境地區(qū)自然保護區(qū)的全溝硬蜱和森林革蜱中分離出10株遠(yuǎn)東亞型和3株西伯利亞亞型TBEV,并表明同一自然疫源地可同時存在2種基因型。本研究中,在吉林省吉林市的全溝硬蜱中分離出1株TBEV。將獲得的蜱源TBEV毒株JL-T75與其他亞型代表株的編碼區(qū)序列和全基因組核苷酸序列進(jìn)行同源性比對和系統(tǒng)進(jìn)化樹分析發(fā)現(xiàn),分離株JL-T75屬于致病性最高的遠(yuǎn)東亞型TBEV,說明我國東北地區(qū)流行的主要TBEV亞型仍為遠(yuǎn)東亞型。
TBEV疫苗株森張株是我國50年代的分離株,與森張株相比,本研究新分離TBEV毒株JL-T75雖然存在一定程度的突變,但它們與森張株編碼區(qū)核苷酸和氨基酸同源性均高于98%。有研究表明發(fā)生在TBEV的E蛋白上關(guān)鍵位點(N52→S、D308→K、K311→E、T310→K、H392→Y)的氨基酸的突變可顯著影響病毒的毒力[14-15]。而新分離的JL-T75株與森張株相比,發(fā)現(xiàn)主要毒力位點未發(fā)生變異,并不會影響TBEV的毒力,推測接種TBEV疫苗后,依然可以產(chǎn)生對新分離的毒株JL-T75的保護作用。綜上所述,本研究為我國東北地區(qū)蜱傳TBEV的遺傳進(jìn)化和有效防控提供了理論基礎(chǔ)。