吳 青,李 銀,李祥瑞,趙冬敏,劉青濤,劉宇卓,韓凱凱,黃欣梅
(1.江蘇省農(nóng)業(yè)科學(xué)院,江蘇 南京 210014;2.南京農(nóng)業(yè)大學(xué),江蘇 南京 210095)
江蘇H9N2亞型禽流感病毒HA和NA的遺傳進(jìn)化分析
吳 青1,2,李 銀1*,李祥瑞2*,趙冬敏1,劉青濤1,劉宇卓1,韓凱凱1,黃欣梅1
(1.江蘇省農(nóng)業(yè)科學(xué)院,江蘇 南京 210014;2.南京農(nóng)業(yè)大學(xué),江蘇 南京 210095)
為了解江蘇H9N2亞型禽流感病毒(AIV)的變異情況,采用RT-PCR技術(shù)對10株從江蘇不同地區(qū)分離的H9N2亞型禽流感病毒的HA和NA基因進(jìn)行了擴(kuò)增、克隆、測序和遺傳進(jìn)化分析。結(jié)果顯示:10株H9N2病毒間HA基因核苷酸和氨基酸的同源性分別是94.3%~99.0%、93.8%~99.6%,NA基因核苷酸和氨基酸的同源性分別是94.9%~99.9%、95.2%~100.0%,它們均屬于歐亞分支中的A/Duck/Hong Kong/Y280/97亞群;10株病毒的HA裂解位點(diǎn)均為RSSR-GLF,受體結(jié)合位點(diǎn)第198位有8株為蘇氨酸,2株為丙氨酸,第191位有1株為組氨酸;10株病毒NA基因在62~64位均存在缺失的現(xiàn)象,從而導(dǎo)致61位的糖基化位點(diǎn)缺失,NA基因第402位均出現(xiàn)了由天冬酰胺(N)到天冬氨酸(D)的變異,8株病毒只具有6個(gè)糖基化位點(diǎn)。
江蘇;H9N2亞型;禽流感病毒;HA;NA;遺傳進(jìn)化
禽流感(Avian Influenza)被世界動物衛(wèi)生組織(OIE)列入A類傳染病,被中國列為一類傳染病。H9N2亞型禽流感病毒(Avian Influenza Virus, AIV)屬于低致病性禽流感病毒(LPAI),家禽感染此病毒后發(fā)病率和死亡率較低,臨床表現(xiàn)為無癥狀帶毒(隱性感染)、產(chǎn)蛋量下降、呼吸系統(tǒng)疾病等。但是,因?yàn)镠9N2亞型禽流感病毒具有傳播性強(qiáng)、變異快、逐步突破哺乳動物屏障等特點(diǎn)[1],因此其不僅給家禽養(yǎng)殖業(yè)帶來了損失,而且對人和哺乳動物的健康也構(gòu)成了潛在威脅。我國1994年首次報(bào)道從雞群分離到H9N2亞型AIV[2],目前,全國各省都有H9亞型AIV的流行。自1999年以來,我國不斷有人感染H9N2病毒的報(bào)道[3],其中2013年暴發(fā)的人感染H7N9亞型AIV與H9N2亞型AIV關(guān)系密切,已證實(shí)其內(nèi)部基因是通過基因片段重組從H9N2亞型AIV獲得的[4]。
血凝素(Hemagglutinin, HA)是構(gòu)成病毒囊膜纖突的主要成分之一,在病毒吸附及穿膜過程中起關(guān)鍵作用,是所有流感病毒蛋白中最重要的1個(gè)。在AIV感染時(shí),只有在HA蛋白裂解的情況下才能使病毒囊膜和小泡膜融合以及破裂。因此,HA蛋白的可裂解性對病毒的毒力有著至關(guān)重要的作用[5]。HA的第234位氨基酸非常關(guān)鍵[6],若是谷氨酰胺(Q),則易與唾液酸(sialic acid, SA)α-2,3-Gal受體結(jié)合;若為亮氨酸(L),則與SA α-2,6-Gal受體結(jié)合;人流感病毒優(yōu)先結(jié)合SA α-2,6-Gal受體,AIV則易與SA α-2,3-Gal受體結(jié)合[7]。HA糖基化位點(diǎn)不僅在介導(dǎo)病毒傳播和受體結(jié)合方面起作用,而且能夠保護(hù)病毒免受機(jī)體免疫系統(tǒng)的作用[8]。有關(guān)研究表明HA糖基化的改變不僅影響病毒對宿主細(xì)胞的識別,而且影響自然殺傷細(xì)胞對感染細(xì)胞的識別。神經(jīng)氨酸酶(Neuraminidase, NA)也是流感病毒粒子表面的1個(gè)重要抗原,NA可以將位于宿主細(xì)胞表面的流感病毒受體末端的唾液酸殘基切斷,從而使病毒成功侵染細(xì)胞[9];同時(shí),NA可以切斷子代病毒受體末端的唾液酸殘基,使得病毒釋放和成熟[10]。HA的受體結(jié)合特性和NA的受體裂解特性相互適應(yīng),共同發(fā)揮作用。HA和NA的變異會引起AIV的抗原漂移以及抗原轉(zhuǎn)變,對于禽流感病毒的毒力具有很大的影響[11]。因此,對H9N2病毒的HA和NA進(jìn)行遺傳進(jìn)化分析有著十分重要的意義。本研究中,我們對2016年從江蘇各地區(qū)分離的H9N2亞型禽流感病毒HA和NA基因進(jìn)行了遺傳進(jìn)化分析,以期為深入研究H9N2亞型禽流感病毒提供分子依據(jù)。
1.1 病毒
采集江蘇省內(nèi)部分地區(qū)家禽的肛拭子和咽拭子,將拭子置于運(yùn)輸液(加入四抗的PBS)中,將棉拭子擠壓后棄去,以5000 r/min離心10 min后取上清;吸取上清液0.2 mL,接種9~10日齡的SPF雞胚(購于南京天邦生物科技有限公司);收集接種后24~96 h死亡的雞胚尿囊液,經(jīng)HA和HI鑒定,將H9N2陽性尿囊液在-70 ℃下保存?zhèn)溆谩?/p>
1.2 主要試劑
液體病毒RNA/DNA抽提試劑盒、凝膠回收試劑盒購自Axygen生物科技有限公司;first-strand反應(yīng)緩沖液、High Pure dNTPs (10 mmol/L)、Rnase Inhibitor(40 U/μL)、反轉(zhuǎn)錄酶AMV、dNTPs (2.5 mmol/L)、Mg2+(25 mmol/L)、Ex Taq反應(yīng)緩沖液(10×)、Ex Taq酶均購自TAKARA公司;通用引物Uni12由南京金斯瑞生物科技有限公司合成。
1.3 引物合成
參照GenBank中H9N2 AIV的基因序列,應(yīng)用Primer 5.0設(shè)計(jì)引物(見表1),交由南京金斯瑞生物科技有限公司合成。
表1 HA和NA的PCR擴(kuò)增引物序列
1.4 提取病毒RNA并進(jìn)行反轉(zhuǎn)錄
用RNA提取試劑盒提取雞胚尿囊液中H9N2流感病毒的總RNA。反轉(zhuǎn)錄采用20 μL體系,包括dNTPs 2 μL、5×AMV Buffer 4 μL、Reverse Transcriptase XL(AMV) 0.5 μL、Recombinant Rnase Inhibitor 0.5 μL、Unit12 1 μL、RNA 12 μL。反應(yīng)程序?yàn)?5 ℃ 5 min,42 ℃ 60 min,98 ℃ 5 min。
1.5 擴(kuò)增目的片段并測序
以反轉(zhuǎn)錄產(chǎn)物為模板進(jìn)行PCR擴(kuò)增。PCR反應(yīng)體系為25 μL,其中雙蒸水14.75 μL,10×PCR Buffer 2.5 μL,25 mmol/L MgCl21.5 μL,dNTP (2.5 mmol/L)2.0 μL,上、下游引物(25 pmol/μL)各1.0 μL,模板DNA 2.0 μL,Ex Taq (5 U/μL)0.25 μL。PCR反應(yīng)程序: 94 ℃預(yù)變性5 min;94 ℃變性40 s,52 ℃退火40 s,72 ℃延伸3 min,共35個(gè)循環(huán);72 ℃延伸10 min。在循環(huán)反應(yīng)結(jié)束后,對擴(kuò)增產(chǎn)物進(jìn)行1.5%瓊脂糖凝膠電泳檢測,對PCR產(chǎn)物進(jìn)行膠回收。將回收產(chǎn)物送至南京金斯瑞生物科技有限公司進(jìn)行測序。
1.6 HA、NA基因序列分析
用Lasergene 7.1 DNAStar軟件對HA、NA基因進(jìn)行序列拼接、同源性分析,并對潛在糖基化位點(diǎn)、唾液酸結(jié)合位點(diǎn)進(jìn)行分析。用MEGA 5.0對10株病毒的HA、NA基因序列與GenBank中已發(fā)表的H9N2亞型AIV經(jīng)典毒株進(jìn)行遺傳進(jìn)化分析,繪制進(jìn)化樹。參考毒株的名稱如下: TU/WS/1/66、CK/HK/Y280/1997、CK/BJ/1/1994、CK/SH/F/98、CK/HK/G9/1997、 CK/GX/10/1999、Quai/HK/G1/97、DK/HK/Y439/1997、CK/JS/HA/2012、CK/HB/31/2000。
2.1 HA、NA基因擴(kuò)增及序列測定
分別提取10株H9N2亞型禽流感病毒的RNA,用RT-PCR方法擴(kuò)增HA、NA基因片段。10株病毒的代號及簡稱如下:A/Chicken/Jiangsu/01/2016(簡稱A/CK/JS/01/16)、A/Chicken/Jiangsu/02/2016(簡稱A/CK/JS/02/16)、A/Chicken/Jiangsu/03/2016(簡稱A/CK/JS/03/16)、A/Chicken/Jiangsu/04/2016(簡稱A/CK/JS/04/16)、A/Chicken/Jiangsu/05/2016(簡稱A/CK/JS/05/16)、A/Chicken/Jiangsu/06/2016(簡稱A/CK/JS/06/16)、A/Chicken/Jiangsu/07/2016(簡稱A/CK/JS/07/16)、A/Chicken/Jiangsu/08/2016(簡稱A/CK/JS/08/16)、A/Chicken/Jiangsu/09/2016(簡稱A/CK/JS/09/16)、A/Chicken/Jiangsu/10/2016(簡稱A/CK/JS/10/16)。
2.2 HA、NA基因同源性分析
通過DNAStar軟件比較分析,得出10株H9N2病毒間HA基因核苷酸的同源性為94.3%~99.0%,推導(dǎo)氨基酸同源性為93.8%~99.6%;NA基因核苷酸的同源性為94.9%~99.9%,推導(dǎo)氨基酸同源性為95.2%~100.0%。10株病毒與H9N2亞型禽流感病毒歐亞分支中Y280-like亞系代表株CK/HK/Y280/1997的親緣性較近,HA基因核苷酸同源性為89.5%~90.4%,推導(dǎo)氨基酸同源性為91.0%~92.5%;NA基因核苷酸同源性為91.1%~94.1%,推導(dǎo)氨基酸同源性為91.3%~94.3%。
2.3 HA、NA蛋白氨基酸序列分析
將10株病毒與歐亞分支Y280-like亞系代表株CK/HK/Y280/1997進(jìn)行比較分析。HA蛋白的基因裂解位點(diǎn)第335~341位為RSSR-GLF,只有1個(gè)堿性氨基酸,符合H9N2亞型AIV典型的低致病性特征;第198位受體結(jié)合位點(diǎn)有8株病毒由纈氨酸(V)變異為蘇氨酸(T),2株變異為丙氨酸(A),可能會導(dǎo)致H9N2亞型AIV與細(xì)胞受體的結(jié)合力發(fā)生變異,對AIV的流行性可能起重要作用;第234位氨基酸與受體唾液酸α-2,3-半乳糖和唾液酸α-2,6-半乳糖的結(jié)合特性有關(guān),禽類的H9N2亞型 AIV的第234位大多為谷氨酰胺(Q),與唾液酸α-2,3-半乳糖受體結(jié)合。而在本研究中10株病毒第234位均突變?yōu)榱涟彼?L)(表2),易與唾液酸α-2,6-半乳糖受體結(jié)合,可能引發(fā)人感染。上述結(jié)果表明這10株病毒存在引發(fā)人和其他哺乳動物流感的潛在可能。
N-糖基化位點(diǎn)(擁有N-X-T/S基序,X為除P以外的氨基酸)被證實(shí)對于維持蛋白分子的結(jié)構(gòu)和正常生物學(xué)功能具有重要意義。10株病毒HA蛋白的第313位均由脯氨酸(P)變異為絲氨酸(S),因而在第313位新增加了1個(gè)糖基化位點(diǎn);9株病毒的第218位糖基化位點(diǎn)由蘇氨酸(T)變異為異亮氨酸(I),1株病毒由蘇氨酸(T)變異為纈氨酸(V)(表2),導(dǎo)致第218位糖基化位點(diǎn)缺失;10株病毒NA蛋白第63位缺失,導(dǎo)致第61位的糖基化位點(diǎn)缺失;第402位由天冬酰胺(N)變異為天冬氨酸(D),導(dǎo)致10株病毒第402位糖基化位點(diǎn)缺失;8株病毒第264位由絲氨酸(S)變異為精氨酸(R),1株由異亮氨酸(I)變異為纈氨酸(V),導(dǎo)致9株病毒缺失了第264位糖基化位點(diǎn)(表3)。
唾液酸結(jié)合位點(diǎn),即紅細(xì)胞結(jié)合位點(diǎn)(hemadsorbing site, HB)位于NA蛋白的表面,遠(yuǎn)離神經(jīng)氨酸酶酶活性區(qū)域,對病毒的復(fù)制具有重要作用。對10株病毒NA蛋白HB位點(diǎn)氨基酸變異情況的分析結(jié)果如表3所示:10株病毒在431~433位氨基酸比較保守,均為PQE;而在366~373區(qū)域的368、369 位氨基酸出現(xiàn)較大變異。
表2 H9N2亞型AIV HA基因的特殊位點(diǎn)
2.4 HA、NA遺傳進(jìn)化分析
將10株病毒的HA、NA基因核苷酸序列與參考毒株進(jìn)行比較,并繪制系統(tǒng)進(jìn)化樹。如圖1、圖2所示:10株病毒的HA、NA基因均屬于歐亞分支中以DK/HK/Y280/97為代表株的Y280-like亞分支。
H9N2亞型AIV是低致病性流感病毒,主要表現(xiàn)為高感染率和低死亡率,能引起家禽生產(chǎn)性能下降和免疫抑制,繼發(fā)細(xì)菌感染后能造成家禽死亡,嚴(yán)重影響家禽業(yè)生產(chǎn)。禽流感病毒變異后能在人、家禽、豬體之間循環(huán),有的能直接感染人,有的為感染人的病毒提供基因,如:2003年香港發(fā)生H9N2亞型禽流感病毒感染人的疫情,2013年感染人并引起較高死亡率的H7N9亞型禽流感病毒的6個(gè)內(nèi)部基因也均來自H9N2病毒[12]。H9N2病毒的 8個(gè)基因片段為適應(yīng)新的宿主環(huán)境在不斷發(fā)生變化,以表面基因HA和NA的變化最為顯著。本研究中的10株H9N2亞型流感病毒與參考毒株HA基因核苷酸及推導(dǎo)的氨基酸序列的同源率分別在89.5%~90.4%、91.0%~92.5%,NA基因核苷酸及推導(dǎo)的氨基酸序列的同源率分別在91.1%~94.1%、91.3%~94.3%。這就表明HA基因和NA基因隨著時(shí)間的變化,其突變在積累。
表3 H9N2亞型AIV NA基因的特殊位點(diǎn)
HA蛋白是流感病毒重要的表面糖蛋白,含有與病毒生物學(xué)功能相關(guān)的多個(gè)功能域,這些區(qū)域的一些氨基酸變化可能導(dǎo)致病毒的一些生物學(xué)功能的改變[13]。在一般情況下,高致病性禽流感病毒的裂解位點(diǎn)附近有多個(gè)堿性氨基酸插入,而低致病力毒株的裂解位點(diǎn)附近只有1個(gè)或沒有堿性氨基酸。本研究中的10株病毒裂解位點(diǎn)均為RSSR-GLF,只有1個(gè)堿性氨基酸,具有低致病性禽流感病毒的分子特點(diǎn)。HA基因第234位受體結(jié)合位點(diǎn)變化顯著,均由谷氨酰胺(Q)突變?yōu)榱涟彼?L),受體結(jié)合位點(diǎn)的突變可能引起流感病毒在自然進(jìn)化過程中與細(xì)胞受體的親和力發(fā)生變異。
糖基化位點(diǎn)能防止抗體結(jié)合對病毒的中和作用[14],是改變禽流感病毒毒力的重要因素之一。在本研究中,9株病毒HA蛋白第218位由蘇氨酸(T)變異為異亮氨酸(I),1株病毒由蘇氨酸(T)變異為纈氨酸(V)(表2),導(dǎo)致第218位糖基化位點(diǎn)缺失;第 313位由脯氨酸(P)變異為絲氨酸(S),因而第313位新增加了1個(gè)糖基化位點(diǎn)。較早的研究表明,中國大陸地區(qū)H9N2亞型AIV的NA蛋白中第63、264、402位潛在糖基化位點(diǎn)高度保守[15]。本研究中10株病毒NA蛋白第63位缺失,從而導(dǎo)致61位的糖基化位點(diǎn)缺失;第402位由天冬酰胺(N)變異為天冬氨酸(D),導(dǎo)致10株病毒第402位糖基化位點(diǎn)缺失;8株病毒第264位由絲氨酸(S)變異為精氨酸(R),1株由異亮氨酸(I)變異為纈氨酸(V),導(dǎo)致9株病毒缺失了第264位糖基化位點(diǎn)。與早期報(bào)道相比,本研究中病毒NA蛋白潛在的糖基化位點(diǎn)出現(xiàn)了較多的變異,這是否會影響病毒的生物學(xué)特性還有待進(jìn)一步研究。
在NA基因的血凝素結(jié)合位點(diǎn)367~372位、399~404位和431~433位共同形成了3個(gè)環(huán)狀結(jié)構(gòu),而環(huán)狀結(jié)構(gòu)中的氨基酸與唾液酸直接結(jié)合[16]。在本研究中,10株病毒的NA基因的血凝素結(jié)合位點(diǎn)HB的367~372位的氨基酸發(fā)生了變異,這可能會降低與唾液酸結(jié)合的活性,但是否會引起H9N2禽流感病毒的致病力增強(qiáng)尚未可知。
本研究對江蘇各地區(qū)2016年分離的10株H9N2亞型AIV的HA、NA基因進(jìn)行了序列分析。結(jié)果表明,HA蛋白第218位潛在糖基化位點(diǎn)缺失,在第313位新增潛在糖基化位點(diǎn),此變異對病毒致病力的影響有待進(jìn)一步研究。此外,NA蛋白的糖基化位點(diǎn)、HB位點(diǎn)發(fā)生了較大變異;值得注意的是,HA蛋白的第234位氨基酸均為L,說明這10株病毒具有潛在感染人的特性。H9N2亞型AIV能引起家禽的免疫抑制,若繼發(fā)細(xì)菌感染則可造成家禽的大量死亡,會給養(yǎng)禽業(yè)帶來巨大的經(jīng)濟(jì)損失。因此加強(qiáng)H9N2亞型AIV流行病的監(jiān)測與研究,對畜禽生產(chǎn)和公共衛(wèi)生都有非常重要的意義。
[1] Xu K M, Smith G J, BahI J, et al. The genesis and evolution of H9N2 inflluenza viruses in poultry from southern China, 2000 to 2005 [J]. J Virol, 2007, 81(19): 10389-10401.
[2] 陳伯倫,張澤紀(jì),陳偉斌.禽流感研究:I雞A型流感病毒分離與血清學(xué)初步鑒定[J].中國獸醫(yī)雜志,1994,22(10):3-5.
[3] Sun Y, Pu J, Jiang Z, et al. Genotypic evolution and antigenic drift of H9N2 influenza viruses in China from 1994 to 2008 [J]. Vet Microbiol, 2010, 146: 215-225.
[4] Lam T T, Wang J, Shen Y, et al. The genesis and source of H7N9 influenza viruses causing human infection in China [J]. Nature, 2013, 502(7470): 241-244.
[5] Graaf D M, Fouchier R A. Role of receptor binding specificity in influenza A virus transmission and pathogenesis [J]. EMBO, 2010, 33(8): 823-841.
[6] Matrosovich M, Tuzikov A, Bovin N, et al. Early alterations of the receptor-binding properties of H1, H2, and H3 avian influenza virus hemagglutinins after their introduction into mammals [J]. Virol, 2000, 74(18): 8502-8512.
[7] Rogers G N, Paulson J C. Receptor determinants of human and animal influenza virus isolates: differences in receptor specificity of the H3 hemagglutinin based on species of origin [J]. Virology, 1983, 127(2): 361-373.
[8] Perdue M L, Latimer J, Greene C, et al. Consistent occurrence of hemagglutinin variants among avian influenza virus isolates of the H7 subtype [J]. Virus Res, 1994, 34(1): 15-29.
[9] Hay A J. The virus genome and its replication [C]//Nicholson K G, Webster R G, Hay A J. Textbook of influenza. London: Blackwell Science Ltd, 1998: 65-73.
[10] Colman P M. Structure and function of neuraminidase [C]//Nicholson K G, Webster R G, Hay A J. Textbook of influenza. London: Blackwell Science Ltd, 1998: 65-73.
[11] 金元昌,李景鵬,張龍,等.禽流感病毒分子生物學(xué)研究進(jìn)展[J].動物醫(yī)學(xué)進(jìn)展,2003(1):1215.
[12] lius J K, Chen J. Phylogenetic diversity and distribution of type A influenza virus [J]. PLoS One, 2009, 4(3): e5022.
[13] Guan Y, Shortridge K F, Krauss S, et al. H9N2 influenza viruses possessing H5N1 like internal genomes continue to circulate in poultry in Southeastern China [J]. Virology, 2000, 74(20): 9372-9380.
[14] Baigent S J, McCauley J W. Glycosylation of hacmalutinin and stalk-length of neuraminidase combine to regulate the growth of avian influenza virus in tissue culture [J]. Virus Res, 2001, 79(1/2): 177-185.
[15] 郭霄峰,廖明,辛朝安. H9N2亞型禽流感病毒纖突蛋白基因的分析[J].中國人獸共患病雜志,2002,18(4) :17-19.
[16] Matrosovich M N, Krauss S, Webster R G. H9N2 influenza A viruses from poultry in Asia have human virus-like receptor specificity [J]. Virology, 2001, 281(2): 156-162.
[17] 鄧顯文,謝芝勛,謝志勤,等.H9亞型禽流感病毒血凝素單克隆抗體制備及鑒定[J].南方農(nóng)業(yè)學(xué)報(bào),2016,47(4):679-683.
(責(zé)任編輯:黃榮華)
Genetic Evolution Analysis of H9N2-subtype Avian Influenza Viruses HA and NA in Jiangsu Province
WU Qing1,2, LI Yin1*, LI Xiang-rui2*, ZHAO Dong-min1,LIU Qing-tao1, LIU Yu-zhuo1, HAN Kai-kai1, HUANG Xin-mei1
(1. Jiangsu Academy of Agricultural Sciences, Nanjing 210014, China;2. Nanjing Agricultural University, Nanjing 210095, China)
In order to understand the variation situation of H9N2-subtype avian influenza virus (AIV) in Jiangsu area, we cloned the HA gene and NA gene of 10 H9N2-subtype AIV strains isolated from different regions of Jiangsu by RT-PCR amplification, and carried out the genetic evolution analysis of the HA and NA sequences. The results showed that the homology of nucleotide and amino acid of HA gene among 10 strains was 94.3%~99.0% and 93.8~99.6%, respectively; the homology of nucleotide and amino acid of NA gene among 10 strains was 94.9%~99.9% and 95.2%~100.0%, respectively; ten strains all belonged to A/Duck/Hong Kong/Y280/97 subgroup in Eurasian clade; the HA cleavage sites of all 10 isolates were RSSR-GLF; at the 198th receptor binding site, eight of ten strains had threonine, and two of ten strains had alanine; at the 191st receptor binding site, one of ten strains had histidine; the deletion of the 62~64th sites in NA gene of all 10 strains resulted in the deletion of the 61st glycosylation site; the variation from asparagine (N) to aspartic acid (D) occurred at the 402nd site of NA gene, and eight of ten H9N2-subtype AIV strains had only six glycosylation sites.
Jiangsu; H9N2 subtype; Avian influenza virus; HA; NA; Genetic evolution
2017-03-28
江蘇省農(nóng)業(yè)科技自主創(chuàng)新資金項(xiàng)目[cx(14)2091、cx(15)1058]。
吳青,女,碩士,主要從事動物病毒病分子生物學(xué)研究。*通訊作者:李銀、李祥瑞。
S858.3
A
1001-8581(2017)07-0104-06