《Influenza and Other Respiratory Viruses》：A Mechanistic Viral Host Perturbation Score (VHPS) Reveals Shared and Virus-Specific Transcriptional Programs Across Human Respiratory Infections
编辑推荐：
背景：急性病毒感染会触发重叠但具有病毒特异性的血液转录反应；然而，当前的分析很少量化各通路、调控因子和免疫细胞组成中激活的机制深度。
方法：对来自GSE157240的RNA-seq数据进行重新分析，使用DoRothEA/decoupleR推断转录因子（TF）
背景：急性病毒感染会触发重叠但具有病毒特异性的血液转录反应；然而，当前的分析很少量化各通路、调控因子和免疫细胞组成中激活的机制深度。
方法：对来自GSE157240的RNA-seq数据进行重新分析，使用DoRothEA/decoupleR推断转录因子（TF）活性，使用PROGENy推断信号通路活性，并评估应激通路和转录扰动广度指标。将这四个组成部分整合为病毒宿主扰动评分（VHPS），这是一个旨在量化宿主转录扰动程度的概念验证框架。采用差异表达分析、基因本体论（GO）富集、模块评分和基于CIBERSORTx的免疫去卷积，对呼吸道病毒感染与系统性病毒感染进行比较。
结果：VHPS将病毒沿宿主转录扰动轴进行分层，其中DNA病毒显示出最高的中位扰动值，其次是登革病毒和肠道病毒/鼻病毒，而甲型流感病毒和其他呼吸道RNA病毒处于中间水平，健康对照组最低。上调基因揭示了病毒特异性和共享的特征，免疫/应激GO条目显示呼吸道病毒具有更强的T细胞和干扰素富集。一个包含26个基因的“共享呼吸道抗病毒”模块与VHPS及TF/应激组分相关，但与PROGENy活性无关。免疫去卷积揭示了病毒特异性的白细胞重塑，其中呼吸道病毒表现为显著的调节性T细胞（Treg）耗竭和中性粒细胞升高，登革病毒表现为单核细胞偏移，而DNA病毒的Treg比例最低。
结论：综合来看，VHPS和免疫去卷积描绘了急性病毒感染中共享和病毒特异性的调控程序，但VHPS的普适性仍需在独立数据集中进行评估。
**论文解读：基于VHPS的急性病毒感染宿主转录扰动机制分析**
**一、研究背景与科学问题**
急性呼吸道病毒感染仍是全球发病率和死亡率的主要原因。尽管病原体各异，不同病毒感染患者常呈现重叠的临床综合征，在诊疗现场难以区分。大规模血液转录组研究已表明，许多急性病毒感染共享保守的干扰素刺激基因（ISG）特征，并叠加反映不同免疫和组织嗜性的病毒特异性转录程序。然而，现有分析大多集中在差异表达基因（DEG）列表和通路过表达上，对上游调控因子（UR）的协调活性以及决定宿主激活程度的特定信号回路知之甚少。因此，亟需一个能够整合转录因子（TF）活性、通路参与、细胞应激和全网络扰动为统一严重度指标的正式定量框架。
GSE157240数据集由Tsalik等人建立，包含多种急性病毒感染（甲型流感病毒、肠道病毒/鼻病毒、偏肺病毒、登革病毒及其他呼吸道RNA和DNA病毒）成人全血RNA-seq数据及健康对照。虽然原始分析确认了以I型干扰素信号为主导的常见抗病毒反应，并指出个别病毒激活部分不同的免疫模块，但仍缺乏将机制性活动推断方法应用于多病毒队列、以区分共享与病毒特异性调控架构的系统性研究。此外，现有严重度概念多依赖临床评分或住院等替代指标，缺乏机制可解释的宿主扰动评分来跨病毒比较。
**二、研究目的与意义**
本研究通过重新分析GSE157240，引入病毒宿主扰动评分（VHPS），旨在量化每个样本中宿主转录扰动的深度，实现以下三个目标：（i）量化各病毒组和样本在TF、通路、应激和网络组分上的差异；（ii）通过整合DEG重叠和GO富集区分共享与病毒特异性调控程序；（iii）确定保守的“共享呼吸道抗病毒”基因模块是否与VHPS在全队列中相关。VHPS提供了以机制而非纯基因中心术语衡量严重度的综合模板，为未来感染研究中系统级严重度评估奠定基础。
**三、主要技术方法**
研究人员使用GSE157240队列的全血RNA-seq数据（包含22例健康对照及六种病毒感染组，共160例样本）。核心方法包括：（1）使用limma-voom流程进行差异表达分析；（2）应用DoRothEA/decoupleR框架（置信水平A–C）经多元线性模型（MLM）推断TF活性；（3）使用PROGENy足迹估计信号通路活性（包括JAK–STAT、NF-κB、MAPK和TNF）；（4）计算预定义应激基因集的模块评分和基于蛋白质-蛋白质相互作用数据库的转录扰动广度；（5）将四个标准化组分（UR_norm、CP_norm、Stress_norm、Network_norm）加权（0.3/0.3/0.2/0.2）求和并min–max归一化生成VHPS；（6）使用Enrichr/gseapy进行GO生物过程富集分析；（7）采用CIBERSORTx（绝对模式，100次置换）推断22种白细胞亚群的相对丰度。统计分析使用Kruskal–Wallis检验和Benjamini–Hochberg校正的Dunn事后检验，Spearman秩相关用于模块与VHPS组分间的关联分析。
**四、研究结果**
**4.1 VHPS的构建**
VHPS揭示了不同病毒病因间的宿主扰动层级（表1；图1）。DNA病毒感染表现出最高的中位转录扰动（p<0.001 vs. 对照），其次是登革病毒和肠道病毒/鼻病毒（p<0.01）。Levene方差检验（p=0.012）确认DNA病毒组内宿主激活存在显著异质性。甲型流感病毒和其他呼吸道RNA病毒处于中间位置，而偏肺病毒聚类接近健康对照。
**4.2 不同病毒严重度下的调控架构**
将VHPS分解为组成成分后发现，总严重度评分并非由单一反应层驱动，而是由病毒特异性调控策略决定（图3）。DNA病毒、登革病毒和肠道病毒/鼻病毒感染的UR评分较高（p<0.01 vs. 对照），表明其扰动主要由TF网络的广泛重连驱动；而“其他呼吸道RNA病毒”组呈现独特轮廓——CP和转录扰动广度评分显著升高而TF活性较低，提示其反应以广泛信号级联激活为主。应激通路评分在所有病毒组中均升高，肠道病毒/鼻病毒和甲型流感病毒中位数最高。
**4.3 VHPS主要受病毒驱动且独立于生物学性别**
无论总体还是按病毒分层，男性和女性的VHPS分布均高度重叠，Mann–Whitney U检验无显著差异（p>0.05），表明生物学性别在本队列中不对转录扰动深度产生系统性调节。
**4.4 呼吸道与全身性病毒激活不同转录程序**
维恩图（图5）显示，甲型流感病毒、肠道病毒/鼻病毒和偏肺病毒共享26个上调DEG的核心集，但甲型流感病毒诱导255个独特DEG，反映广泛转录重塑。DNA病毒和登革病毒与呼吸道RNA病毒共享7个“全身性核心”基因，但DNA病毒拥有71个独特DEG。GO富集（图6）显示，经典呼吸道病毒在T细胞活化、趋化性和干扰素应答方面富集强而凝聚；而全身性感染（DNA病毒、登革病毒）呈现异质性模式，干扰素和免疫调节过程富集较弱。
**4.5 共享呼吸道抗病毒模块连接TF活性与严重度**
26基因“共享呼吸道抗病毒”模块评分与总VHPS呈中等正相关（Spearman r=0.319，p=8.9×10??），且与UR_norm（r=0.349，p<0.001）、Stress_norm（r=0.272，p=0.003）和Network_norm（r=0.179，p=0.042）显著相关，但与CP_norm无显著相关性（r=0.095，p=0.281）。这提示保守的呼吸道抗病毒程序是TF驱动的特异性应激反应，而非广泛信号通路激活的替代指标。
**4.6 免疫细胞景观的病毒特异性重塑**
CIBERSORTx去卷积发现，除偏肺病毒外的所有感染组均显示调节性T细胞（Tregs）显著耗竭（p<10??），DNA病毒、甲型流感病毒和肠道病毒/鼻病毒中的降幅最大。呼吸道感染的特征为Treg耗竭伴中性粒细胞比例升高；全身性感染（登革病毒和DNA病毒）以单核细胞偏移为主；偏肺病毒呈中间表型。
**五、讨论与结论**
**讨论部分：** VHPS框架将研究焦点从静态基因列表转向整合调控活性。DNA病毒和登革病毒驱动最深刻的系统级扰动，呼吸道RNA病毒居中，偏肺病毒反应相对温和，这一层级与原始GSE157240分析一致。研究人员提出，全血转录组扰动不仅反映倍数变化幅度，更反映病原体重塑宿主调控架构的程度。关键概念是TF活性（UR）与转录扰动广度的解耦：DNA病毒因核内复制而强占宿主TF；胞质RNA病毒则触发广泛PRR信号而不一定进行集中的核重编程。应激组分与保守抗病毒模块的相关性提示“抗病毒状态”与细胞应激机制（如ISR和UPR）内在关联。Treg在严重感染中的耗竭与高VHPS相关，提示VHPS能有效捕捉外周免疫调节的丧失。26基因模块与TF和应激组分相关而与PROGENy通路不相关，支持“抗病毒骨架”假说。登革病毒结果与部分文献不一致，可能反映队列临床异质性。由于所有分析限于单一GSE157240队列，VHPS主要作为可重复应用的分析框架，需在外部队列中验证。
**结论：** 总之，VHPS框架证明急性病毒感染可沿连续的机制负担轴进行分层，该轴整合了调控活性、应激和网络拓扑。研究显示“严重度”并非单一概念：DNA病毒的重负担由TF重连和Treg耗竭驱动，而呼吸道RNA病毒往往依赖于应激通路和信号网络扩展。通过将保守的呼吸道基因模块连接到特定调控轴，本研究为超越基因列表、走向机制知情严重度谱系分析提供了模板，对监测新发病原体和识别宿主反应治疗靶点具有直接意义。
