靶向性爱泼斯坦-巴尔病毒（EBV）捕获测序鉴定出BBLF4-L322M是鼻咽癌中一个独立的预后变异体

《Microbiology Spectrum》：Targeted Epstein–Barr virus capture sequencing identifies BBLF4-L322M as an independent prognostic variant in nasopharyngeal carcinoma

　　

爱泼斯坦-巴尔病毒（Epstein–Barr virus, EBV）是鼻咽癌（nasopharyngeal carcinoma, NPC）的致病因素，并且表现出地域异质性。为了确定病毒基因型是否影响预后，研究人员对114例广西鼻咽癌活检组织的EBV进行了基于捕获的全基因组测序（平均覆盖率约93%，深度约979×），并在BNRF1-G696R位点进行了Sanger测序验证。研究人员共检测到5,683个高置信度变异体，平均每个基因组有1,112个，这些变异体在潜伏基因（LMP1、LMP2A/B和EBNA家族）中富集。次要等位基因频率≥0.15且缺失率≤0.05的变异体被用于与5年总生存期（overall survival, OS）进行关联性分析。单变量Cox模型中共有46个位点显著相关，随后对19个错义变异体进行了基于基因型相关性的无监督聚类，得到了两个模块。从这些模块中，研究人员选择BZLF1-A205S和BBLF4-L322M作为代表性变异体。在调整了年龄和二分化的美国癌症分期联合委员会（American Joint Committee on Cancer, AJCC）分期后进行的多变量Cox模型中，BBLF4-L322M仍然与较差的总生存期（OS）独立相关，而BZLF1-A205S则不显著。整合了本研究序列与公共数据的720个基因组的系统发育分析显示，存在多个具有部分地域结构性的中国亚谱系。广西的EBV毒株与广东和香港的毒株关系密切，但又部分不同。这些结果支持将EBV基因型，特别是BBLF4-L322M，纳入预后模型，并推动功能研究以明确其机制并验证其临床实用性。 重要性 目前，与鼻咽癌（NPC）预后相关的爱泼斯坦-巴尔病毒（EBV）变异体在很大程度上仍是未知的。本研究首次在中国广西的鼻咽癌患者中进行了EBV序列多样性的全基因组调查。识别出的突变，包括BZLF1-A205S和BBLF4-L322M，是与患者生存相关的潜在病毒生物标志物。潜伏期基因中丰富的非同义变异体表明病毒基因型在疾病进程中可能扮演着重要角色。这些见解为基于EBV的预后生物标志物和为鼻咽癌流行人群量身定制的精准医疗策略奠定了基础。

研究背景、问题与目标

鼻咽癌（Nasopharyngeal carcinoma, NPC）具有显著的地理和种族异质性，超过70%的病例发生在东亚和东南亚。在流行区，爱泼斯坦-巴尔病毒（Epstein–Barr virus, EBV）在恶性上皮细胞中被一致地检测到，凸显了其病因学作用。EBV具有明显的遗传地理结构，其遗传特征（如BALF2_CCT单倍型和EBER非编码RNA区多态性）已被证明与NPC风险相关，血浆EBV DNA载量也与肿瘤分期和预后相关。然而，直接将特定的EBV序列变异体与临床预后（如总生存期）联系起来的研究仍然有限。此外，中国广西壮族自治区是高危人群，其NPC年发病率高达20-40例/10万人，但广西地区EBV的种群结构尚未得到系统分析。为了填补这一空白，本研究旨在通过全基因组测序鉴定能预测NPC预后的EBV基因组标志物，并分析广西地区EBV亚型的遗传多样性及其与预后的关联。

研究方法与技术路径

研究人员在广西招募了114名经病理确诊、既往未经治疗、有EBV感染证据且可获得至少5年随访数据的鼻咽癌（NPC）患者。从诊断活检样本中提取DNA，使用针对EBV基因组设计的安捷伦SureSelectXT定制探针进行靶向捕获富集，随后在Illumina HiSeq/NovaSeq平台上进行双端测序。数据分析方面，将原始测序数据比对到EBV参考基因组（NC_007605.1），遵循GATK最佳实践流程进行变异检测，并注释变异的功能类别。利用公共数据库中的EBV基因组序列构建了包含720个毒株的系统发育树。预后关联分析采用Cox比例风险回归模型，筛选出与总生存期（OS）相关的EBV变异体，并通过调整年龄和AJCC分期等临床协变量进行多变量分析，以识别独立的预后因素。

研究结果

1. 探针捕获测序与注释

研究人员对114株从NPC活检组织中获得的EBV毒株进行了捕获测序，平均测序覆盖度约为93%，平均测序深度为979×。对与NPC风险相关的BNRF1-G696R变异体位点的Sanger测序验证结果与捕获测序结果完全一致。过滤后，共检测到5,683个SNP/InDel变异体，平均每个样本1,112个变异。变异体在基因组中的分布主要集中在编码区（73.2%），其中非同义替换占47.4%，同义替换占46.9%。短插入缺失（1-5 bp）构成了大多数插入缺失事件。

2. EBV基因多态性

研究人员分析了包含至少10个独立变异体的EBV基因，发现10个变异最丰富的基因中，有7个是潜伏期相关基因（LMP-2A, LMP-2B, LMP-1, EBNA-3A, EBNA-3B/3C, EBNA-2, EBNA-1），它们占所有检测到变异体的44.6%。长度归一化的对数比率分析进一步显示，LMP-1, LMP-2B, EBNA-3B/3C, EBNA-2和EBNA-3A的变异体丰富度显著高于预期。这表明潜伏基因是EBV基因组中变异最密集的区域。

3. EBV潜伏基因的dN/dS分析

对包含编码区替换的6个潜伏基因（LMP-2A, LMP-1, EBNA-3A, EBNA-3B/3C, EBNA-2, EBNA-1）的分析显示，EBNA-3A, EBNA-3B/3C, LMP-1, EBNA-2和EBNA-1的dN/dS比值均大于1，表明存在强烈的正选择，这与免疫介导的压力推动病毒在长期宿主持续感染中发生氨基酸改变以适应宿主的选择压力一致。

4. 识别与NPC预后相关的风险SNP

为了研究EBV基因组变异对NPC预后的影响，研究人员对114名具有至少5年随访数据的NPC患者的EBV基因组进行了单变量和多变量Cox回归分析。经过质量控制后，保留了796个SNP。单变量Cox回归分析发现，有46个SNP与总生存期显著相关。其中，BZLF1基因中的六个高频错义变异体（A205S, Q163L, V152A, T124P, S76P, T68A）显示出与较好预后相关的保护性效应。而BBLF4-L322M变异体则与死亡风险增加显著相关。对19个错义候选变异体进行无监督聚类，得到了两个不同的基因型模块。从每个模块中选择代表性变异体（BZLF1-A205S和BBLF4-L322M）进行多变量Cox模型分析。在调整了年龄和二分化的AJCC分期后，晚期AJCC分期是预后不良的强预测因子，而BBLF4-L322M仍与更高的死亡率独立相关，其保护性效应在调整后不再显著。这些结果将BBLF4-L322M确定为NPC早期死亡的一个独立风险变异体。

5. 地理和表型中EBV基因组的系统发育结构

研究人员重建了一个包含720个EBV毒株全基因组的系统发育树，并根据地理来源和疾病表型对每个毒株进行注释。系统发育分析揭示了EBV的多个病毒谱系。中国的毒株（广西、广东、香港）构成了序列的大部分，并形成了几个大的亚谱系，显示出部分地理隔离性。进一步分析BBLF4-L322M在不同地理和疾病表型中的分布发现，该变异体在香港、广西和广东的数据集中被发现，但未在非洲或欧洲基因组中检测到。当按疾病类别分组时，在现有公开数据集中，BBLF4-L322M在健康携带者中的频率高于NPC患者。这表明BBLF4-L322M并非NPC特异性的致病性突变，而是表现出明显的地理结构特征，并且在东亚EBV谱系中富集。

讨论与结论总结

本研究通过广西NPC患者的EBV全基因组测序，揭示了具有潜在临床相关性的病毒变异谱系特异性模式。遗传变异在潜伏期基因中高度富集，dN/dS分析表明多个潜伏期抗原基因处于正选择之下，这与免疫驱动的EBV多样化一致。短框移插入缺失（1-5 bp）构成了观察到的大多数突变，表明插入缺失变异是EBV在NPC中进化的一个共同特征。本研究的主要创新在于发现了EBV基因变异与NPC患者预后的关联。特别是在调整临床分期后，BBLF4基因中第322位的亮氨酸-甲硫氨酸替换（L322M）仍是生存期的强有力预测因子。相比之下，BBLF4是一个裂解周期基因，编码EBV解旋酶-引发酶复合物的解旋酶组分。BBLF4-L322M的预后关联可能反映了改变的病毒复制生物学，而非潜伏基因的机制。系统发育分析显示了明显的地理结构，中国NPC来源的EBV基因组形成了几个亚谱系，广西的毒株与广东和香港的毒株部分分离。这种模式表明，EBV毒株的多样性受到宿主种群历史和地理因素的影响。BBLF4-L322M在东亚数据集中被观察到，但在非洲或欧洲数据集中未检测到。因此，BBLF4-L322M不应被视为一种全球通用的预后标志物。相反，目前的证据更支持它是一种在EBV相关的中国南方流行区NPC人群中可能具有预后价值的谱系相关变异体。值得注意的是，汇总的公开比较数据显示，BBLF4-L322M在NPC中的频率并未高于健康携带者，这提示它不太可能是NPC特异性的易感突变。然而，预后关联和疾病特异性是两个不同的问题。一个病毒变异体可能在背景人群中很常见，但仍能在发展为EBV相关NPC的患者中对临床结果进行分层。因此，研究人员将BBLF4-L322M更狭义地解释为广西NPC队列中的一个候选预后标志物，而非定义疾病的EBV改变。

本研究存在一定局限性，如样本来自单一流行中心，需要进行独立、前瞻性、多中心队列的验证；未评估宿主遗传学相互作用或肿瘤体细胞改变；缺乏对L322M生物学效应的直接功能数据；比较的公共数据组在地理上不完全匹配，缺乏来自广西的健康对照等。因此，BBLF4-L322M应被视为候选的地区相关预后标志物，而非NPC特异性的EBV突变。

结论翻译

总而言之，本研究确定BBLF4-L322M是一种候选的地区相关性EBV变异体，在来自广西的EBV相关鼻咽癌患者中，它与较差的生存期独立相关。这些数据扩展了病毒遗传学在疾病易感性之外的作用，并指明了有潜力的研究方向。将病毒基因分型与临床协变量结合可能优化风险分层，并且有必要开展对BBLF4的机制研究，以阐明其在鼻咽癌中的病毒-宿主相互作用。EBV序列变异可以影响NPC的临床进程，而我们将BBLF4-L322M确定为生存期较差的标志物，这为理解病毒在癌症进展中的作用提供了新的视角。