資海榮，李偉，周丹，李婕，宣楊，郭艷，衛(wèi)平民

(東南大學(xué) 公共衛(wèi)生學(xué)院，江蘇 南京 210009)

禽流感(Bird Flu或Avian Influenza)是禽流行性感冒的簡(jiǎn)稱，它是由甲型流感病毒的一種亞型(也稱禽流感病毒)引起的一種急性傳染病，也能感染人類，被國(guó)際獸疫局定為甲類傳染病，又稱真性雞瘟或歐洲雞瘟。到目前為止，已報(bào)道[1- 6]的可引起人類感染的禽流感病毒有H7N7、H9N2 和H5N1。自從1997年在香港暴發(fā)H5N1以來(lái)，此病癥引起世界各國(guó)衛(wèi)生組織的高度關(guān)注。以往僅在禽類發(fā)現(xiàn)甲型H7亞型病毒感染，2013年2月底開始，我國(guó)東部省市地區(qū)陸續(xù)首次發(fā)現(xiàn)了由甲型H7N9禽流感病毒引起人感染的病例[7]。截至2013年8月30日，新型甲型H7N9禽流感散發(fā)涉及12個(gè)城市(包括上海、安徽、浙江、江蘇、廣州、山東、湖南、河南、江西、福建、河北、北京)，確診感染病例達(dá)135例，其中死亡病例45例。進(jìn)一步研究[8]表明，H7N9流感病毒能夠發(fā)生有限的人傳人。

甲型流感病毒系單股負(fù)鏈RNA病毒，其基因組包括8個(gè)獨(dú)立片段：PB2、PB1、PA、HA、NP、NA、M、NS，它們依次編碼多聚酶成分(PB2、PB1、PA)、血凝素(HA)、核衣殼蛋白(NP)、神經(jīng)氨酸酶(NA)、基質(zhì)蛋白(M1)、膜蛋白(M2)、非結(jié)構(gòu)蛋白(NS1、NEP)[9- 10]。2001年，Chen等[11]在研究MHCⅠ類分子遞呈抗原肽時(shí)發(fā)現(xiàn)了一種不屬于流感病毒任何已知開放閱讀框(open reading frame，ORF)編碼的多肽——PB1- F2蛋白，該蛋白是導(dǎo)致甲型流感病毒致病性的關(guān)鍵蛋白之一。 研究[12]顯示，PB1- F2蛋白可以誘導(dǎo)宿主細(xì)胞線粒體途徑的細(xì)胞凋亡，部分流感病毒毒株的PB1- F2蛋白可以提高流感病毒聚合酶活性，從而促進(jìn)病毒在宿主細(xì)胞中復(fù)制。同時(shí)，完整的PB1- F2蛋白上存在特異的抗原表位，能誘導(dǎo)宿主體內(nèi)免疫T細(xì)胞的表達(dá)，引起其免疫水平的改變，這與二次細(xì)菌感染引發(fā)的肺損傷具有顯著關(guān)聯(lián)；還發(fā)現(xiàn)隨著基因的進(jìn)化，有部分甲型流感病毒的亞型在PB1- F2序列上提前出現(xiàn)了終止密碼子，導(dǎo)致PB1- F2蛋白的斷裂，形成PB1- F2長(zhǎng)度的多態(tài)性，而且從時(shí)間上呈現(xiàn)了一定的進(jìn)化特征。截短型PB1- F2蛋白仍保留部分功能，如A/Beijing/11/1956(H1N1)；已發(fā)現(xiàn)A/Hong Kong/1/1968(H3N2)的完整PB1- F2蛋白具有促炎癥反應(yīng)功能，H7N9作為新型甲型流感病毒的亞型之一，其PB1- F2蛋白的進(jìn)化現(xiàn)狀值得監(jiān)測(cè)并探討。本研究收集IRD數(shù)據(jù)庫(kù)上的自2013年2月19日到2013年10月31日的H7N9流感毒株63株，利用MEGA 5.2.1軟件進(jìn)行分析，旨在對(duì)新型甲型H7N9流感病毒的PB1基因進(jìn)化趨勢(shì)進(jìn)行分析，并對(duì)H7N9高致病性機(jī)制進(jìn)行探討。

1 資料與方法

1.1 資料來(lái)源

流感病毒基因序列通過(guò) The NIAID Influenza Research Database(IRD) online(http://www.fludb.org)獲得，下載自2013年2月19日到2013年10月31日的新型甲型H7N9流感病毒PB1基因序列63條，其中上海22株、江蘇13株、安徽1株、山東7株、廣東2株、浙江8株、湖南2株、江西2株、河南2株、福建2株、臺(tái)灣2株；環(huán)境樣本毒株13株，禽流感毒株33株(其中雞源26株，鴨源2株，鴿子源5株)，人流感17株。

1.2 構(gòu)建進(jìn)化樹

利用Molecular Evolutionary Genetics Analysis(MEGA)Version 5.2.1 軟件進(jìn)行核苷酸序列比對(duì)，采用鄰接(neighbor- joining，NJ)法(bootstraps=1 000)繪制H7N9流感病毒PB1基因的系統(tǒng)進(jìn)化樹，分析2013年新型甲型H7N9流感病毒PB1基因的進(jìn)化狀況。

1.3 PB1- F2氨基酸位點(diǎn)改變

利用MEGA 5.2.1 軟件進(jìn)行氨基酸序列比對(duì)，對(duì)PB1- F2氨基酸位點(diǎn)改變進(jìn)行分析。

1.4 各型氨基酸構(gòu)成比例

利用MEGA 5.2.1軟件進(jìn)行氨基酸序列比對(duì)，對(duì)PB1- F2蛋白中各類氨基酸構(gòu)成比進(jìn)行運(yùn)算、輸出。Excell表格對(duì)酸性氨基酸、堿性氨基酸、疏水性氨基酸、親水性氨基酸的構(gòu)成比進(jìn)行計(jì)算，利用PASW 18.0 進(jìn)行統(tǒng)計(jì)學(xué)分析。

2 結(jié)果與分析

2.1 H7N9流感病毒PB1系統(tǒng)進(jìn)化樹

見(jiàn)圖1。

從圖1 PB1系統(tǒng)進(jìn)化樹可以發(fā)現(xiàn)，63株新型甲型H7N9流感毒株主要進(jìn)化為兩大系別，18株病毒株屬于系別Ⅰ，45株病毒株屬于系別Ⅱ。系別Ⅰ與系別Ⅱ之間的核苷酸同源性以及氨基酸同源性分別為 96.3%～100%和98.6%～100%，系別Ⅰ核苷酸同源性以及氨基酸同源性分別為96.8%～100%和98.6%～100%，系別Ⅱ核苷酸同源性以及氨基酸同源性分別為98.5%～100%和99.2%～100%。

2.2 H7N9流感病毒中發(fā)現(xiàn)PB1- F2斷裂株存在

將下載的PB1序列通過(guò)MEGA 5.2.1軟件進(jìn)行比對(duì)，在PB1的氨基酸第95～367從起始密碼子ATG的位點(diǎn)到TAA或者TGA的終止密碼子中截選出來(lái)，利用MEGA 5.2.1軟件進(jìn)行氨基酸序列比對(duì)，在63株毒株中發(fā)現(xiàn)有17株序列同時(shí)在氨基酸第26位提前出現(xiàn)了終止密碼子，導(dǎo)致PB1- F2蛋白在第25位氨基酸后出現(xiàn)斷裂。出現(xiàn)截短型PB1- F2蛋白毒株中分別有環(huán)境樣本6株，其中3株來(lái)自山東，2株來(lái)自河南，1株來(lái)自上海；禽流感8株，其中7株為雞源性，1株為鴿子源性；人流感3株，分別為A/Nanchang/1/2013、A/Changsha/1/2013、A/Changsha/2/2013。具體見(jiàn)圖2。

2.3 H7N9流感病毒的氨基酸位點(diǎn)改變

通過(guò)對(duì)H7N9流感病毒PB1- F2病毒的氨基酸序列進(jìn)行分析，發(fā)現(xiàn)本研究中的63株毒株中共有23處位點(diǎn)發(fā)生了改變。由圖2可以看出，PB1- F2蛋白截短型毒株基因發(fā)生突變的位點(diǎn)明顯多于PB1- F2蛋白完整型毒株，且PB1- F2蛋白截短型毒株基因的突變位置相對(duì)較固定。然而PB1- F2蛋白截短型毒株中A/chicken/Rizhao/875/2013的發(fā)生突變的位點(diǎn)相對(duì)于其他毒株要多。完整型PB1- F2蛋白毒株相對(duì)較少，卻也同樣存在發(fā)生突變位點(diǎn)較高的毒株，比如A/chicken/Shanghai/S1080/2013發(fā)生突變的位點(diǎn)卻相對(duì)較多，此外還有A/chicken/Zhejiang/SD019/2013、A/Guangdong/1/2013、A/duck/Zhejiang/SC410/2013、A/duck/Zhejiang/SC410/2013。

2.4 H7N9流感病毒PB1- F2各類氨基酸構(gòu)成

見(jiàn)表1。

氨基酸根據(jù)其結(jié)構(gòu)與特性可分為疏水性氨基酸、親水性氨基酸、酸性氨基酸、堿性氨基酸[13]，63株流感毒株中構(gòu)成PB1- F2蛋白的各類氨基酸比例見(jiàn)表1。經(jīng)方差分析，完整型PB1- F2蛋白毒株與截短型PB1- F2蛋白毒株的各類氨基酸構(gòu)成比差異有統(tǒng)計(jì)學(xué)意義(F1=44.25，P1<0.01；F2=17.83，P2<0.01；F3=719.02，P3<0.01；F4=1 016.19，P4<0.01)(1=疏水性氨基酸，2=親水性氨基酸，3=酸性氨基酸，4=堿性氨基酸)。數(shù)據(jù)顯示，由于PB1- F2蛋白的斷裂導(dǎo)致構(gòu)成各類氨基酸的比例發(fā)生改變。

表1 2013年H7N9流感毒株P(guān)B1-F2蛋白各類氨基酸比例

毒株類型毒株數(shù)PB1-F2氨基酸長(zhǎng)度PB1-F2蛋白各類氨基酸構(gòu)成疏水性氨基酸比例/%親水性氨基酸比例/%酸性氨基酸比例/%堿性氨基酸比例/%完整型PB1-F24690aa25.33±0.9743.52±0.948.92±0.2922.23±0.25截短型PB1-F21752aa27.15±0.9342.31±1.185.43±0.7525.11±0.47

3 討 論

通過(guò)對(duì)甲型H7N9流感病毒PB1系統(tǒng)進(jìn)化樹分析可以發(fā)現(xiàn)，系別 Ⅰ、Ⅱ的劃分主要是依據(jù)PB1- F2蛋白斷裂與否。系別 Ⅰ主要是完整型PB1- F2蛋白毒株，除A/chicken/Rizhao/875/2013外。系別Ⅱ主要是截短型PB1- F2蛋白毒株，A/chicken/Shanghai/S1358/2013、A/chicken/Jiangsu/SC099/2013除外。本研究結(jié)果表明，截短型PB1- F2蛋白已經(jīng)成為新型甲型H7N9進(jìn)化的重要特征。

● 代表人H7N9毒株，◇代表環(huán)境樣本毒株；分支處數(shù)據(jù)是Bootstrap檢驗(yàn)的可信度；0.005標(biāo)尺代表單位長(zhǎng)度內(nèi)核苷酸差異水平

圖1 2013年期間63株甲型H7N9流感毒株P(guān)B1系統(tǒng)進(jìn)化樹

方框內(nèi)*為終止密碼子部位，即發(fā)生斷裂部位

圖2 2013年甲型H7N9流感毒株P(guān)B1-F2氨基酸序列特征

PB1- F2基因有40個(gè)左右的位置可以通過(guò)單個(gè)核甘酸的突變提前引入終止密碼子，從而編碼不同長(zhǎng)度的C端截短型PB1- F2蛋白[14]。McAuley等[15- 16]采用一系列方法，對(duì)來(lái)源于不同毒株、完整長(zhǎng)度的PB1- F2蛋白的相關(guān)功能分別進(jìn)行研究，通過(guò)反向遺傳方法，將PR8病毒的7個(gè)基因片段(PB2、PA、HA、NP、NA、M、NS)和不同來(lái)源的PB1片段[A/Brevig Mission/1/1918(H1N1)、A/Vietnam/1203/2004(H5N1)]分別克隆至pHW2000質(zhì)粒，將每一種PB1片段和PR8病毒的7個(gè)基因片段進(jìn)行重組轉(zhuǎn)化，合成多種來(lái)源PB1片段的重組病毒；再對(duì)PR8和新病毒進(jìn)行細(xì)胞和小鼠致病性實(shí)驗(yàn)，比較不同毒株P(guān)B1- F2蛋白的促細(xì)胞凋亡作用、致小鼠肺部炎癥及小鼠二次肺炎鏈球菌感染的嚴(yán)重程度，并監(jiān)測(cè)病毒感染細(xì)胞時(shí)聚合酶活性的變化；又合成既定長(zhǎng)度PB1- F2的C端肽鏈(第67至90位aa)。研究結(jié)果發(fā)現(xiàn)，來(lái)源于不同毒株的完整PB1- F2蛋白的功能存在一定差異；截短型PB1- F2蛋白仍保留部分功能,如A/Beijing/11/1956(H1N1)，PB1- F2蛋白長(zhǎng)度為57aa，截短型PB1- F2蛋白能增加病毒在293T細(xì)胞中的聚合酶活性。以上差異同樣被其他相關(guān)研究所證實(shí)[17- 19 ]。

在本研究中的63株新型甲型H7N9流感毒株中發(fā)現(xiàn)17株出現(xiàn)斷裂，且均生在PB1- F2的第25位氨基酸后，提示H7N9流感病毒在25位氨基酸后發(fā)生斷裂是比較固定的斷裂位點(diǎn)，而H7N9流感PB1- F2截短型毒株的出現(xiàn)對(duì)于H7N9致病性的影響尚需要進(jìn)一步研究。

已有研究[20]發(fā)現(xiàn)，Brevig Mission(1918年，H1N1)與香港(1997年，H5N1)均有發(fā)現(xiàn)流感毒株中PB1- F2蛋白的第66位堿基發(fā)生點(diǎn)突變，產(chǎn)生氨基酸N66S，并證實(shí)其具可提升流感病毒的致病性。本研究中63株新型甲型H7N9流感毒株中第66位均為N，未發(fā)生堿基突變。有研究[21- 22]發(fā)現(xiàn)，H7N9人流感病毒是一個(gè)來(lái)自于H7N3(A/duck/Zhejiang/12/2011，ZJ12亞型)、H7N9(A/wild bird/Korea/A14/2011，KO14亞型)、H9N2(A/brambling/Beijing/16/2012- like viruses)的禽流感病毒重配株。其中6個(gè)內(nèi)部基因(M、NS、NP、PB1、PB2、PA)則來(lái)自H9N2?，F(xiàn)有關(guān)于H9N2基因進(jìn)化的研究[14]發(fā)現(xiàn)，H9N2進(jìn)化中存在66N與N66S毒株，本研究中未發(fā)現(xiàn)N66S毒株，提示在進(jìn)行新型甲型H7N9的重配過(guò)程中，參與重配的H9N2屬于66N毒株。

本研究通過(guò)對(duì)自2013年2月19日到2013年10月31日所有的甲型H7N9流感病毒毒株P(guān)B1- F2蛋白核酸以及氨基酸序列的分析，發(fā)現(xiàn)在H7N9流感毒株中存在著截短型PB1- F2蛋白毒株，并且H7N9流感病毒毒株P(guān)B1- F2蛋白斷裂位點(diǎn)均在25位氨基酸后，這將成為H7N9流感毒株進(jìn)化過(guò)程中重要的特征。本研究結(jié)果發(fā)現(xiàn)，甲型H7N9流感病毒截短型PB1- F2蛋白表現(xiàn)為更加明顯的突變?；谘芯拷Y(jié)果猜想，截短型流感毒株易于導(dǎo)致基因結(jié)構(gòu)的不穩(wěn)定，從而使流感毒株易于發(fā)生基因突變，進(jìn)而利于流感病毒的復(fù)制與傳播。這一猜想有待進(jìn)一步研究和驗(yàn)證。此外，完整型PB1- F2蛋白與截短型PB1- F2蛋白氨基酸構(gòu)成存在明顯的差異，這種差異對(duì)于H7N9流感病毒的致病性機(jī)制的影響也有待進(jìn)一步的研究。

[1] LUPIANI B，REDDY S M.The history of avian influenza[J].Comp Immunol Microbiol Infect Dis，2009，32(4)：311- 323．

[2] FOUCHIER R A，SCHNEEBERGER P M，Rozendaal F W，et al.Avian influenza A virus(H7N7)associated with human conjunctivitis and a fatal case of acute respiratory distress syndrome[J].Proc Natl Acad Sci USA，2004，101(5)：1356- 1361．

[3] KOOPMANS M，WILBRINK B，CONYN M，et al.Transmission of H7N7 avian influenza A virusto human beings during a large outbreak in commercial poultry farms in the Netherlands[J].Lancet，2004，363(9409)：587- 593．

[4] LIN Y P，SHAW M，Gregory V，et al.Avian- to- human transmission of H9N2 subtype influenza A viruses：relationship between H9N2 and H5N1 human isolates[J].Proc Natl Acad Sci USA，2000，97(17)：9654- 9658．

[5] PEIRIS J S，de JONG M D，GUAN Y.Avian influenza virus(H5N1)：a threat to human health[J].Clin Microbiol Rev，2007，20(2)：243- 267．

[6] GAMBOTTO A，BARRATT- BOYES S M，de JONG M D，et al.Human infection with highly pathogenic H5N1 influenza virus[J].Lancet，2008，371(9622)：1464- 1475．

[7] 劉春艷，艾軍紅.甲型H7N9禽流感病毒的病毒學(xué)特征[J].中國(guó)當(dāng)代兒科雜志，2013，15(6)，405- 408.

[8] Avian influenza A H7N9 in zhejiang，China[DB/OL]published online. April.25,2013,http://dx.doi.org/10.1016/s0140- 6736(13)60936- 8.

[9] WEBSTER R G，BEAN W J，GORMAN OT，et al.Evolution and ecology of influenza A viruses[J].Micro Rev，1992，56(1):152- 179．

[10] McHARDY A C，ADAMS B.The role of genomics in tracking the evolution of influenza A virus[J].PLoS Pathog，2009，5(10)：e1000566．

[11] CHEN W，CALVO P A，MALIDE D，et al.A novel influenza A virus mitochondrial protein that induces cell death[J].Nature Medicine，2001，7(12):1306- 1312．

[12] 資海榮，李偉，李婕，等.甲型流感病毒PB1- F2蛋白功能的研究進(jìn)展[J].東南大學(xué)學(xué)報(bào)：醫(yī)學(xué)版，2013，32(5):631- 635．

[13] 莊榮輝.生物化學(xué)基礎(chǔ)[M/OL].臺(tái)灣：網(wǎng)頁(yè)系統(tǒng)之生物化學(xué)基礎(chǔ)，2001，1.

[14] 黃艷艷，胡北俠，文心田.A型流感病毒中國(guó)分離株P(guān)B1- F2基因進(jìn)化分析[J].中國(guó)農(nóng)業(yè)科學(xué)，2009，42(6):2156- 2163.

[15] McAULEY J L，CHIPUK J E，BOYD K L，et al.PB1- F2 proteins from H5N1 and 20 century pandemic influenza viruses cause immunopathology[J].PLoS Pathog，2010，6：e1001014．

[16] McAULEY J L，ZHANG K，MCCULLERS J A.The effects of influenza A virus PB1- F2 protein on polymerase activity are strain specific and do not impact pathogenesis[J].J Virol，2010，84(1)：558- 564．

[17] CHEN C J，CHEN G W，WANG C H，et al.Differential localization and function of PB1- F2 derived from different strains of influenza a virus[J].J Virol，2010，84(19)：10051- 10062．

[18] HAI R，SCHMOLKE M，VARGA Z T，et al.PB1- F2 expression by the 2009 pandemic H1N1 influenza virus has minimal impact on virulence in animal models[J].J Virol，2010，86(9)：4442- 4450．

[19] KOSIK I，KREJNUSOVA I，BYSRRICKA M，et al.N- terminal region of the PB1- F2 protein is responsible for increased expression of influenza A viral protein PB1[J].Acta Virologica，2011，55：45- 53．

[20] CONENLL G M，TISONCIK J R，ROSENZWEIG E，et al.A single N66S mutation in the PB1- F2 protein of influenza avirus increases virulence by inhibiting the early interferon responseinvivo[J].J Virol，2011，85(2):652- 662．

[21] GAO R，CAO B，HU Y，et al.Human infection with a novel avian- origin influenza A(H7N9) virus[J].N Engl J Med，2013,368:1888- 1897.

[22] CHENY，LIANG W，YANG S，et al.Human infections with the emerging avian influenza A H7N9 virus from wet market poultry：clinical analysis and characterization of viral genome[J].Lancet，2013，6736(13)：60903- 60904．