IF50.5｜DNA病毒组随人类基因与环境变化

90万人WGS揭示常见病毒载量受宿主遗传、年龄、性别、季节和吸烟强烈影响，且EBV载量是霍奇金淋巴瘤的因果风险因素。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

许多病毒可在人体内终身持续感染，宿主对病毒载量的控制差异可导致清除或疾病。EBV感染是霍奇金淋巴瘤和多发性硬化的重要风险因素，但机制不明。大规模DNA测序已初步描绘血液DNA病毒组，常见类群包括疱疹病毒和指环病毒。病毒载量随时间变化，是免疫功能的潜在生物标志物。然而，病毒载量在个体间和生命周期中的变异程度、驱动因素及其对健康的影响尚不清楚。此前研究已发现MHC区域及少数其他位点影响HIV、丙型肝炎等病毒载量，但对健康人群中潜伏病毒的载量遗传基础了解甚少，非血液组织的病毒载量驱动因素也几乎未被探索。

科学问题

卡在哪里： 此前研究多聚焦于致病性病毒（如HIV、丙肝）的载量遗传学，而对健康人群中常见潜伏病毒（如EBV、HHV-7、指环病毒）的载量变异了解甚少；且病毒载量的遗传效应在血液和唾液等不同组织间的异质性，以及宿主遗传与病毒亚型间的相互作用，均缺乏大规模系统分析。

本文要回答： 利用大规模WGS数据系统描绘血液和唾液中31种常见DNA病毒的载量分布，识别影响病毒载量的宿主遗传位点、非遗传因素（年龄、性别、时间、吸烟等），并评估病毒载量与疾病的因果关系。

技术路线

作者主要用了：全基因组测序、病毒参考基因组比对、GWAS、孟德尔随机化、HLA等位基因分型。

作者从UK Biobank、All of Us和SPARK三个大型队列中提取未比对到人类基因组的测序读段，重新比对到31种常见病毒参考基因组上，定量病毒DNA载量（图1a）。通过严格的质量控制排除假阳性比对，获得血液和唾液中的病毒存在/丰度数据。随后，利用这些表型进行全基因组关联分析（GWAS），识别影响病毒载量的宿主遗传位点；结合HLA等位基因分型精细定位MHC区域信号；通过EBV亚型特异性比对探究宿主-病毒遗传互作；最后利用孟德尔随机化评估EBV载量与疾病的因果关系。

图 1　Fig. 1。图1展示了病毒DNA检测和定量流程及血液/唾液中的分布。a部分显示从UKB、AoU和SPARK队列的WGS数据中提取未比对读段并比对到31种病毒参考基因组。b部分显示UKB血液中EBV、HHV-7和HHV-6B的读段计数分布，大多数样本为0或1条读段，少数样本有大量读段（如eHHV-6整合）。c部分对比血液和唾液中不同病毒的检出率，唾液样本中疱疹病毒检出率远高于血液。d部分显示SPARK唾液中EBV和HHV-7读段计数的双峰分布。e部分显示SPARK唾液中各病毒检出率，MCPyV为5.7%。

核心亮点

亮点 1｜唾液病毒载量呈双峰分布

在SPARK唾液WGS中，EBV、HHV-7等疱疹病毒载量呈现明显的双峰分布：样本要么没有病毒读段，要么有大量读段（数十至数千条）。这提示唾液中的病毒DNA并非来自持续存在的潜伏感染，而是来自频繁的裂解性复制和释放事件，其强度差异很大。HHV-7的裂解事件似乎持续发生，几乎所有成人唾液样本均可检测到HHV-7 DNA；而EBV的唾液排毒常见但间歇性，约一半唾液样本无EBV DNA。相比之下，指环病毒DNA在唾液中仅比血液中略常见，且丰度低，提示其复制周期不同。

图 2　Fig. 2。图2展示病毒载量随年龄、性别、昼夜和季节的变化。a部分显示SPARK唾液中EBV、HHV-6B和HHV-7的检出率随年龄变化，EBV持续升高，HHV-6B从儿童期下降，HHV-7在中年后略降。b部分显示UKB血液中EBV和TTV随年龄升高，HHV-7和HHV-6B随年龄下降。c和d部分显示男性病毒载量高于女性。e部分显示EBV和HHV-7检出率在白天升高。f部分显示EBV冬季检出率高于夏季，HHV-7季节变化较弱。

亮点 2｜病毒载量随年龄、性别和时间变化

病毒DNA载量随年龄、性别、昼夜和季节显著变化。在UKB血液中，EBV和TTV载量随年龄增长而升高，而HHV-7和HHV-6B载量随年龄下降。在SPARK唾液中，多数病毒在生命最初几年感染率迅速上升，之后EBV持续升高，HHV-6B则从儿童期起持续下降。值得注意的是，7种病毒在男性和女性中均显示男性载量更高，且这种性别差异在青春期出现。EBV和HHV-7载量在白天升高（约1.2-1.3倍），EBV冬季比夏季高约1.3倍。这些模式在AoU队列中部分得到重复。

图 3　Fig. 3。图3展示病毒载量的全基因组关联结果。a-c部分分别为EBV、HHV-7和指环病毒TUS01的曼哈顿图，MHC区域信号最强。d部分显示MHC区域对EBV、HHV-7和指环病毒载量遗传力的贡献比例。e和f部分放大MHC区域，显示EBV和HHV-7的关联模式不同，EBV最强关联为HLA-DRB104:04，HHV-7为其他HLA等位基因。g部分显示ERAP1-ERAP2-LNPEP位点与EBV和HHV-7载量的关联。*

亮点 3｜MHC区域主导病毒载量遗传效应

GWAS在血液中鉴定出82个与病毒载量相关的位点（P<5×10-8），涉及EBV（45个位点）、HHV-7（37个位点）、HHV-6B、MCPyV和三种指环病毒。除HHV-6B外，MHC区域均产生最强关联信号，且各病毒的关联模式具有病毒特异性。例如，EBV载量最强关联由HLA-DRB1*04:04驱动，该等位基因携带者EBV阳性率几乎翻倍（23.0% vs 13.6%）。MHC区域贡献了EBV载量56%的遗传力，HHV-7为32%，指环病毒为14-30%。此外，ERAP1/ERAP2位点与EBV和HHV-7载量强相关，涉及抗原加工肽酶。

图 4　Fig. 4。图4展示EBV亚型特异性分析。a部分显示SPARK唾液中EBV 1型和2型参考基因组上读段覆盖度，灰色区域为用于区分亚型的EBNA2和EBNA3A-3C区域。b部分显示不同出生地人群中EBV 2型相对1型的比例，非洲出生者2型比例更高。c部分显示英国出生者中2型比例的区域差异。d和e部分分别为EBV 1型和2型载量的MHC区域关联图。f部分显示病例-病例关联分析中HLA-B08:01是最强信号。g部分显示条件分析后MHC区域剩余关联信号。*

亮点 4｜HLA-B*08:01与EBV亚型特异性互作

通过分别定量EBV 1型和2型DNA，发现HLA-B08:01等位基因对EBV亚型具有相反效应：对1型EBV具有保护作用（OR=0.62，P=7.1×10-72），但对2型EBV增加风险（OR=1.63，P=4.1×10-24）。病例-病例关联分析中，B08:01是最强信号（P=7.4×10-70）。该等位基因可能对1型EBNA-2或EBNA-3蛋白衍生的肽段具有高亲和力，从而介导1型特异性细胞毒性T淋巴细胞反应。条件分析显示B*08:01可解释MHC区域大部分亚型特异性关联信号。

图 5　Fig. 5。图5展示病毒载量与临床表型的关联及孟德尔随机化结果。a部分显示UKB中病毒载量与多种疾病表型的关联，包括吸烟相关疾病、免疫抑制状态等。b和c部分显示吸烟包年数和每日吸烟支数与EBV和HHV-7阳性率的关系，EBV随吸烟增加，HHV-7随吸烟下降。d部分显示EBV载量对多发性硬化的孟德尔随机化结果不显著。e和f部分显示EBV载量对霍奇金淋巴瘤的因果效应显著。g部分显示基线EBV阳性者未来霍奇金淋巴瘤风险增加。

亮点 5｜吸烟对EBV和HHV-7载量相反效应

吸烟暴露与EBV DNA载量增加和HHV-7 DNA载量降低显著相关。在UKB中，重度吸烟者EBV阳性率几乎是非吸烟者的两倍（P=3.3×10-195），而HHV-7阳性率随吸烟量增加而下降（P=9.6×10-22）。每日吸烟支数也显示类似趋势。这些关联在AoU血液和唾液样本中得到重复，且吸烟不仅增加EBV阳性率，还增加EBV阳性个体中的病毒丰度。

亮点 6｜EBV载量是霍奇金淋巴瘤的因果风险因素

孟德尔随机化分析支持EBV DNA载量对霍奇金淋巴瘤风险的因果效应（P=1.2×10-3至0.014），但对多发性硬化无因果效应（MR-Egger P=0.52）。此外，基线无霍奇金淋巴瘤的个体中，EBV DNA阳性者未来发生霍奇金淋巴瘤的风险是阴性者的2.01倍（OR=2.01，95%CI 1.40-2.84）。这些结果提示，较高的慢性EBV载量增加淋巴瘤风险，而EBV感染与自身免疫病的关联可能反映宿主对特定病毒表位的免疫反应。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 全基因组测序、病毒参考基因组比对、GWAS、孟德尔随机化、HLA等位基因分型
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事病毒基因组学、感染与免疫遗传学、生物信息学及流行病学的研究人员。

带走一句

利用大规模WGS数据可系统解析宿主遗传、环境因素与病毒载量的关系，并通过孟德尔随机化区分因果与关联，为病毒相关疾病机制研究提供范式。