鼻咽癌和健康人群中EB病毒BARF0基因序列分析

王偉娜，高翔翔，沈智超，羅兵，王云

(1 青島大學醫學院微生物學教研室，山東 青島 266021；2 中國人民解放軍濟南軍區青島第一療養院檢驗科； 3 青島市城陽區流亭衛生院)

?

鼻咽癌和健康人群中EB病毒BARF0基因序列分析

王偉娜1，2，高翔翔3，沈智超1，羅兵1，王云1

(1 青島大學醫學院微生物學教研室，山東 青島 266021；2 中國人民解放軍濟南軍區青島第一療養院檢驗科； 3 青島市城陽區流亭衛生院)

目的 了解山東地區鼻咽癌(NPC)及健康人群中EB病毒(EBV)編碼BARF0基因變異特征，探討其變異意義。方法 采用PCR和DNA測序檢測NPC組織以及健康成年人咽漱液標本中EBV的BARF0基因序列，與EBV標準株進行比較，根據特征性突變對變異進行分類。結果 73例NPC病人和70例健康人群BARF0成功測序，共檢出4處共有突變。其中36例(49.3%)NPC和29例(41.4%)健康人群同時出現4處共有突變(T160473G、T160545C、A160701C和C160707G)，被分類為4M變異型；35例(48.0%)NPC和38例(54.3%)健康人群同時出現除T160473G外的3處共有突變，被分類為3M變異型，其余2例(2.7%)NPC和3例(4.3%)健康人群分類為Others。χ2檢驗精確概率法分析表明,BARF0變異型在2種人群中的分布差異無統計學意義(P=0.618)。結論 山東地區NPC和健康人群中EBV毒株BARF0序列變異一致，BARF0基因變異的地域性分布及與EBV相關腫瘤的關系值得進一步研究。

皰疹病毒4型,人；染色體結構變異；BARF0基因；鼻咽腫瘤

EB病毒(EBV)與多種人類淋巴細胞和上皮細胞性腫瘤的發生有關[1]。EBV BamHI A右向轉錄產物(BARTs)是EBV的一組選擇性剪接轉錄產物，在最后一個外顯子3′末端均含有一個共同的開放讀碼框BARF0。BARTs在所有EBV相關腫瘤組織中表達，有研究表明其可通過維持病毒潛伏感染或作為非編碼RNA參與EBV致癌過程[2-3]。雖然EBV廣泛分布于世界各地，但是鼻咽癌(NPC)等EBV相關腫瘤存在明顯的地域性，因此，尋找可能與腫瘤相關的EBV變異株一直是研究的熱點之一，但某一特定的EBV基因亞型代表地域相關性變異還是疾病相關性變異尚存在著爭議[4]。目前，有關BARTs基因多態性的研究報道較少[5-6]。本研究檢測山東地區NPC病人和健康人群中BARF0基因序列，探討其變異意義，旨在為深入分析BARF0基因變異與腫瘤發生的關系及其生物學功能提供依據。現將結果報告如下。

1 材料與方法

1.1 標本采集及DNA提取

NPC石蠟包埋組織取自青島大學醫學院附屬醫院和山東大學齊魯醫院病理科，通過原位雜交檢測組織標本中EBV編碼小RNA1鑒定EBV陽性NPC組織標本，具體步驟參見文獻的方法[7]。咽漱液(TW)標本取自健康查體人員，通過PCR檢測EBV BamHI W片段篩選和鑒定EBV陽性TW標本[8]。所有病例和健康人群來自相同地區，包括山東省多個地區，其祖籍均為山東且居住在山東的當地人群。采用FFPE石蠟包埋組織DNA提取試劑盒(德國QIAGEN公司)提取石蠟包埋組織DNA，使用酚-氯仿-異戊醇法常規提取新鮮組織和TW標本DNA。

1.2 PCR擴增

BARF0位于BARTs 3′末端，長525 bp(B85-8 coodinates 160470-160994)，推測編碼含174個氨基酸的蛋白質。根據EBV標準株B95-8(GenBank收錄號：V01555)序列，設計擴增BARF0基因的巢式PCR引物，引物序列、在基因組中的位置及擴增片段大小見表1。其中P1和P2為外巢PCR引物，P3和P4為內巢PCR引物。外巢PCR反應體系為20 μL，包括10×Buffer 2.0 μL，MgCl21.5 mmol/L，dNTPs 0.2 mmol/L，上下游引物各0.4 μmol/L，高保真Taq DNA聚合酶(大連TaKaRa公司)1.0 U，DNA 模板100 ng；取外巢PCR產物1 μL為內巢PCR模板，反應體積為30 μL。循環條件為94 ℃預變性5 min；然后，94 ℃、30 s，55 ℃、30 s，72 ℃、1 min，共擴增35個循環；最后72 ℃延伸10 min。取3 μL的 PCR產物在20 g/L的瓊脂糖凝膠上進行電泳并鑒定。

表1 BARF0基因序列分析所用引物

1.3 DNA測序

取25 μL的PCR擴增產物送上海邁浦生物科技有限公司，采用末端終止法進行雙向測序，測序引物為內巢PCR引物。測序結果用Chromas軟件查看峰圖文件，用DNAStar軟件(Lasergene，version 7.0)剪接序列和對排分析，以EBV標準株B95-8序列作為參比序列，同時與基因庫中相關序列進行對比分析。

1.4 統計學分析

以PPMS 1.5[9]軟件進行統計學處理，采用χ2檢驗精確概率法分析BARF0基因變異類型在NPC和健康人群之間分布的差異。以P<0.05為差異有顯著意義。

2 結 果

2.1 BARF0基因序列多態性

共有73例NPC和70例健康人群標本完成測序。與標準株B95-8比較，所有序列均出現變異，均為單核苷酸突變，其變異情況見圖1。共檢出22處突變，組成16種DNA序列，其中4處突變見于較多樣本，屬共有突變，其余18處突變均為見于1～3例樣本的散在突變。4處共有突變中，C160707G(2例為T160707G)見于所有的樣本中，T160545C和A160701C同時出現于幾乎所有的樣本(96.5%，138/143)，而T160473G突變則僅見于部分的樣本(45.5%，65/143)。根據4處共有突變對BARF0基因變異進行分類，將同時出現上述4處共有突變的序列分類為4M變異型，共同出現其中3處共有突變(T160545C、A160701C和C160707G)的序列分類為3M變異型。2種變異型樣本中，多數樣本僅具有共有突變，部分樣本還伴有其他位點的散在突變。有5例樣本除共有突變C160707G(或T160707G)外，未出現其他3處共有突變，但伴有其他位點的散在突變，分類為Others。

圖1中同時列出了GenBank中已經完整測序的19株EBV BARF0基因序列，來源于我國NPC高發區NPC的13株EBV毒株(GD1、GD2、M81、C666-1、HKNPC1-HKNPC9，GenBank基因收錄號分別為：AY961628、HQ020558、KF373730、KJ411974、JQ009376、KF992564-KF992571)[10-13]和另外來源于日本BL的Akata(KC207813)[14]均為4M變異型；而來源于非洲Burkitt淋巴瘤(BL)的AG876(DQ279927)[15]則為3M變異型，另外MUTU(KC207814)[14]只出現A160701C和C160707G的2處共有突變和1處散在突變，將其歸類為Others；來自美國健康成人的K4413Mi和K4123Mi(KC440852和KC440851)[16]與來自美國傳染性單核細胞增多癥的EBV標準株B95-8一致。

2.2 BARF0基因變異類型在NPC及健康人群中的分布

NPC和健康人群中BARF0基因變異類型的分布見表2。多數樣本為3M變異型和4M變異型，少數為Others。BARF0變異類型在2種人群中構成比的差異無統計學意義(P=0.618)。

表2 BARF0基因變異類型在NPC和健康人群中的分布

陰影中序列表示4M和3M變異型的共有序列，每一種序列以1例樣本名稱表示，括號中的數字表示來自NPC和健康人群與該樣本序列相同的樣本數量，N：NPC；T：健康人群。GenBank中已測序毒株以斜體表示，包括HKNPC1-HKNPC9、GD1、GD2、M81、C666-1、Akata、AG876、MUTU、B95-8、K4413Mi和K4213[11-17]。

3 討 論

本研究對山東地區73例NPC和73例健康人群中EBV BARF0基因進行序列多態性分析。與B95-8比較，所有標本BARF0基因均出現單核苷酸突變，無插入和缺失，BARF0基因出現的突變位點少，且不同樣本之間變異較為一致。共檢出4處共有突變，其中C160707G(或T160707G)、T160545C和A160701C同時出現于幾乎所有樣本(124/127)。目前僅有廣州地區20例NPC組織和19株完成測序的EBV毒株可獲知其BARF0序列，其中廣州地區20例NPC組織[6]、來源于廣東和香港NPC的13株EBV毒株[10-13]、來源于日本BL的Akata毒株和非洲BL的AG876[14-15]均同時出現該3處突變，提示該3處突變是BARF0基因的突變熱點。

與上述3處突變不同，T160473G僅見于部分樣本，根據其出現與否，將所檢測樣本中BARF0基因分為4M和3M兩種主要變異類型。兩種變異型在NPC和健康人群中的分布無顯著性差異，提示BARF0變異可能與山東地區NPC發生的易感性無相關性。但值得關注的是，來自NPC高發區的所有20例NPC和13株來源于NPC的EBV毒株均為4M變異型[6,10-13]，而4M變異型在山東地區NPC病人和健康人群中檢出率分別為49.3%(36/73)和41.4%(29/70)，提示4M變異型在高發區NPC具有更高的檢出率。以往有研究表明，EBV其他潛伏期基因如潛伏膜蛋白1的China 1亞型或核抗原1的V-val亞型在高發區NPC的檢出率高于同一地區健康人群，認為它們與NPC的發生相關[4,17]。但有研究顯示，China 1或V-val在NPC和健康人群中檢出率類似，認為它們與NPC發生無關[4]。張琳琳等[5]研究顯示，BARTs A73基因A157154 C突變在廣東地區NPC中的檢出率高于健康人群。由于目前尚缺乏來自NPC高發區正常人群及其他地域NPC和非NPC人群BARF0基因多態性的報道，BARF0 4M變異型毒株是否在高發區NPC中流行及與NPC的相關性，需要進一步采集不同地區樣本及擴大檢測樣本例數證實。

BARTs的功能尚未明確，但由于其在EBV相關腫瘤組織中普遍表達，尤其在上皮性腫瘤中呈高表達，推測其在腫瘤發生過程中發揮作用[1-2]，其可能反義調節互補鏈裂解期基因的表達參與維持病毒的潛伏感染[1-2]，可能作為長鏈非編碼RNA而起作用[3]，BARTs體外表達的某些蛋白如RK-BARF0可以通過與Notch相互作用誘導EBV致癌基因潛伏膜蛋白1的表達[18]。本研究4處共有突變在推測的BARF0蛋白中均引起氨基酸突變。由于對該基因功能、作用機制、基因結構等了解較少，因而現在還難以預測這些變異對其功能的影響。BARF0變異的研究，對明確BARF0基因變異的地域分布、與EBV相關腫瘤的關系及生物學意義具有指導意義，為進一步明確EBV基因變異與腫瘤的相關性提供了實驗依據。

[1] TSAO S W, TSANG C M, TO K F, et al. The role of Epstein-Barr virus in epithelial malignancies[J]. The Journal of Pathology, 2015,235(2):323-333.

[2] MARQUITZ A R, RAAB-TRAUB N. The role of miRNAs and EBV BARTs in NPC[J]. Semin Cancer Biol, 2012,22(2):166-172.

[3] STRONG M J, XU G R, COCO J, et al. Differences in gastric carcinoma microenvironment stratify according to EBV infection intensity: implications for possible immune adjuvant the-rapy[J]. PLoS Pathogens, 2013,9(5): e1003341.

[4] CHANG C M, YU K J, MBULAITEYE S M, et al. The extent of genetic diversity of Epstein-Barr virus and its geogra-phic and disease patterns: a need forreappraisal[J]. Virus Research, 2009,143(2):209-221.

[5] 張琳琳,李大疆,李之珩,等. EB病毒A73基因多態性及其與鼻咽癌易感性的關系[J]. 癌癥, 2007,26(10):1047-1051.

[6] 宋立兵,曾益新,馬英紅,等. 鼻咽癌組織中 EBV-BARF0基因序列的測定[J]. 癌癥, 2001,20(5):464-467.

[7] 王海娟,李慧,王笑峰,等. 鼻型結外NK/T細胞淋巴瘤與EB病毒感染的相關性[J]. 齊魯醫學雜志, 2014,29(6):471-474.

[8] IKUTA K, SATOH Y, HOSHIKAWA Y, et al. Detection of Epstein-Barr virus in salivas and throat washings in healthyadults and children[J]. Microbes and Infection / Institut Pasteur, 2000,2(2):115-120.

[9] 周曉彬,紀新強,徐莉. PPMS 1.5統計軟件的功能及其應用[J]. 青島大學醫學院學報, 2009,45(1):91-93.

[10]ZENG M S, LI D J, LIU Q L, et al. Genomic sequence analysis of Epstein-Barr virus strain GD1 from a nasopharyngeal carcinoma patient[J]. Journal of Virology, 2005,79(24):15323-15330.

[11]LIU P, FANG X D, FENG Z Z, et al. Direct sequencing and characterization of a clinical isolate of Epstein-Barr virus from nasopharyngeal carcinoma tissue by using next-generation sequencing technology[J]. Journal of Virology, 2011,85(21):11291-11299.

[12]TSAI M H, RAYKOVA A, KLINKE O, et al. Spontaneous lytic replication and epitheliotropism define an Epstein-Barr virus strain found in carcinomas[J]. Cell Reports, 2013,5(2):458-470.

[13]KWOK H, WU C W, PALSER A L, et al. Genomic diversity of Epstein-Barr virus genomes isolated from primary nasopharyngeal carcinoma biopsy samples[J]. Journal of Virology, 2014,88(18):10662-10672.

[14]LIN Z, WANG X, STRONG M J, et al. Whole-genome sequencing of the Akata and Mutu Epstein-Barr virus strains[J]. Journal of Virology, 2013,87(2):1172-1182.

[15]DOLAN A, ADDISON C, GATHERER D, et al. The genome of Epstein-Barr virus type 2 strain AG876[J]. Virology, 2006,350(1):164-170.

[16]LEI H Y, LI T W, HUNG G C, et al. Identification and cha-racterization of EBV genomes in spontaneously immortalized human peripheral blood B lymphocytes by NGS technology[J]. BMC Genomics, 2013,14:804.

[17]LI D J, BEI J X, MAI S J, et al. The dominance of China 1 in the spectrum of Epstein-Barr virus strains from Cantonese patients with nasopharyngeal carcinoma[J]. Journal of Medical Virology, 2009,81(7):1253-1260.

[18]KUSANO S, RAAB-TRAUB N. An Epstein-Barr virus protein interacts with notch[J]. Journal of Virology, 2001,75(1):384-395.

(本文編輯 張蓓 于國藝)

SEQUENCE ANALYSIS OF BARF0 GENE OF EPSTEIN-BARR VIRUS IN NASOPHARYNGEAL CARCINOMA AND HEALTHY CROWD

WANGWeina,GAOXiangxiang,SHENZhichao,LUOBing,WANGYun

(Department of Microbiology, Qing-dao University Medical College, Qingdao 266021, China)

ObjectiveTo investigate the characteristics of polymorphisms of BARF0 gene of Epstein-Barr virus (EBV) in nasopharyngeal carcinoma (NPC) and healthy crowd of Shandong province, and explore the significance of BARF0 gene variations.MethodsEmploying PCR and DNA sequencing, the BARF0 sequences of isolates from NPC biopsies and throat washing samples of healthy persons were detected. The sequences were compared with the EBV type strain, and the variations were classified according to distinctive mutation.ResultsA total of 143 isolates from 73 NPC and 70 healthy donors were successfully sequenced and four common mutations were detected. Thirty-six (49.3%) NPC and 29 (41.4%) healthy donor isolates with four common mutations (T160473G, T160545C, A160701C and C160707G) were arranged into one group and named variant 4M, and 35 (48.0%) NPC and 38 (54.3%) healthy donor isolates with three common mutations (T160545C, A160701C and C160707G) were characterized into one group and named variant 3M. The remaining two (2.7%) NPC and three (4.3%) healthy donor isolates were classified as “Others” group. Chi-square test showed that the distribution of the BARF0 variants among two population groups was not significantly different (P=0.618).ConclusionThe isolates in NPC and healthy donors in Shandong province demonstrates consistent variation pattern. The geographical distribution of the BARF0 variants and the relationship between BARF0 variation and EBV-associated tumors are worth further study.

herpesvirus 4, human; genomic structural variation; BARF0 gene; nasopharyngeal neoplasms

2015-02-09;

2015-04-23

國家自然科學基金(81171571)和青島市科技局科研基金(13-1-3-50-nsh)資助項目

王偉娜(1984-)，女，碩士研究生。

王云(1970-)，女，博士，教授，碩士生導師。

R373.9;R739.63

A

1008-0341(2015)04-0383-04