《Cell Biochemistry and Biophysics》：Consensus transcriptional signature of dengue virus infection reveals coordinated interferon suppression and stress-response activation
编辑推荐：
登革病毒感染仍然是全球重大公共卫生威胁，重症登革病可进展为危及生命的出血与休克综合征。尽管其临床表现已被充分描述，但与登革病严重程度相关的转录组学图谱仍未被完全阐明。研究人员对公开数据集开展了多队列转录组分析，以建立登革病毒感染的一致性基因表达特征。研究中采用
登革病毒感染仍然是全球重大公共卫生威胁，重症登革病可进展为危及生命的出血与休克综合征。尽管其临床表现已被充分描述，但与登革病严重程度相关的转录组学图谱仍未被完全阐明。研究人员对公开数据集开展了多队列转录组分析，以建立登革病毒感染的一致性基因表达特征。研究中采用功能富集分析和蛋白质—蛋白质相互作用（PPI）网络鉴定相关生物学过程及高连接节点，并通过计算机模拟（in silico）药物—基因相互作用分析筛选候选化合物。随后，研究人员利用一个独立队列进行了临床验证，以评估不同登革病严重程度分组中的基因表达模式。差异表达分析显示，干扰素相关基因下调，包括STAT1相关通路，同时应激反应基因如CDKN1A（p21）、FOS和JUNB上调。网络分析进一步将这些基因识别为高连接节点，提示抗病毒信号与内皮应激反应之间存在协调性相互作用。基于整理的药物—基因相互作用信息，研究确定了Garcinol、Paclitaxel和Romidepsin等候选化合物，作为这一转录特征的潜在调节因子。本研究界定了与登革病毒感染及登革病严重程度相关的一致性转录特征，强调了宿主反应通路中高连接基因的协调性变化。这些发现为后续实验验证及宿主导向治疗策略（host-directed therapeutic strategies）的开发提供了框架。
该论文发表于《Cell Biochemistry and Biophysics》，聚焦于登革病毒（DENV）感染过程中宿主转录反应的系统性重构，尤其关注与疾病严重程度相关的分子特征。登革热是热带和亚热带地区重要的虫媒传染病，临床表型从自限性发热到伴随血浆渗漏、出血以及休克的重症登革病不等。虽然既往研究已较充分揭示临床分型与部分免疫学现象，但不同实验系统、不同样本来源与不同检测平台之间存在显著异质性，使得能够稳定反映登革病毒感染核心病理过程的一致性转录特征仍不明确。特别是在血管内皮损伤、抗病毒干扰素反应失衡以及炎症/应激通路激活之间究竟呈现何种协调关系，仍缺乏跨队列、跨模型的综合证据。因此，开展本研究的必要性在于，通过整合公开转录组数据，提炼出具有可重复性的核心基因网络，从而为理解宿主反应机制及寻找宿主导向干预靶点提供依据。
研究人员采用多队列系统生物学整合策略，综合分析多个GEO转录组数据集，分别覆盖NS1刺激的内皮模型、登革病毒直接感染的内皮细胞模型以及包含多类细胞反应的宿主数据集。在此基础上，研究人员先进行差异表达分析，再通过跨数据集交集筛选方向一致的差异基因，以建立“共识性转录特征”。随后结合功能富集分析、蛋白质—蛋白质相互作用（PPI）网络分析及药物—基因相互作用数据库挖掘，识别关键生物学过程、网络枢纽节点和潜在调节化合物。最后，研究使用独立临床队列验证这些关键基因在不同疾病严重程度中的表达模式。研究结果表明，登革病毒感染相关的核心转录改变以干扰素相关通路受抑和细胞应激/炎症相关通路被激活为主要特征，STAT1、CDKN1A、FOS和JUNB构成该特征中的高连接核心节点。这一发现的重要意义在于，它并未局限于单一实验模型，而是从跨情境数据中提炼出宿主反应的保守模式，为后续机制实验、分层生物标志物探索和宿主导向治疗策略提供了可操作的分子框架。
在技术方法方面，研究人员主要采用以下几类关键手段：首先，从GEO数据库获取GSE9378、GSE34628、GSE139603三个发现队列，以及独立验证队列GSE17924；前者分别代表不同生物学背景下的登革病毒相关宿主反应，后者为按登革热（DF）、登革出血热（DHF）和登革休克综合征（DSS）分层的临床转录组队列。其次，使用R语言中的limma进行背景校正、分位数标准化与差异表达分析，并以log
2 fold change＞0.5且校正后p值＜0.05为阈值。再次，通过交集分析建立一致性差异基因集，结合Enrichr/MSigDB、GO分析和STRING-db构建PPI网络并进行拓扑分析；同时利用DGIdb进行药物—基因相互作用筛选。最后，采用单因素方差分析、Tukey事后检验、受试者工作特征（ROC）曲线、曲线下面积（AUC）及留一交叉验证（LOOCV）评估临床区分能力。
在研究结果部分，论文首先以“Suppression of Interferon Signaling and Activation of Cellular Stress Pathways”为题，说明研究人员通过多步骤整合分析描绘了登革感染相关的转录图景。差异表达分析在各数据集中独立实施，结果显示不同数据来源间存在相对一致的模式：抗病毒相关基因整体下调，而应激与炎症相关基因整体上调。进一步的交集分析识别出一组在全部数据集中方向一致的核心差异基因，其中干扰素刺激基因（ISGs）如ISG15、OAS3和MX1持续下调，支持登革病毒感染过程中抗病毒反应被协调性抑制这一结论；相对地，CDKN1A、FOS和JUNB等与细胞应激和炎症反应相关的基因上调，提示宿主细胞处于显著的应激激活状态。
其次，在“STAT1 and CDKN1A as Central Regulatory Hubs of Dengue Virus Infection”部分，研究人员借助功能富集和PPI网络分析进一步明确了这些差异基因的生物学组织结构。下调基因显著富集于Interferon Alpha Response和Interferon Gamma Response通路，上调基因则显著富集于经NF-κB介导的TNF-α信号及缺氧相关通路。这一结果表明，登革病毒感染并非单纯诱导炎症，而是伴随抗病毒干扰素轴被压制、应激与炎症轴增强的双向重塑。网络分析显示，下调模块以STAT1、IRF7和MX1等抗病毒调控因子为中心，上调模块则包含CDKN1A、FOS和JUNB等与细胞周期调控及应激反应相关的节点。基于连接度中心性，STAT1和CDKN1A被识别为最关键的高连接节点，说明它们在整个转录响应网络中处于核心位置。
第三，在“Mapping Candidate Modulators of the Core Signature”部分，研究人员将上述枢纽基因纳入药物—基因相互作用框架中，筛选潜在的转录特征调节因子。结果显示，Paclitaxel可同时关联FOS和CDKN1A，Garcinol与STAT1相关调控网络相关，而Romidepsin及4-Phenylbutyric acid（4-PBA）则涉及CDKN1A相关机制。论文明确指出，这些结果基于整理数据库中的已知或预测相互作用，仅构成假说生成层面的候选信息，并不能直接等同于治疗有效性证据。尽管如此，这一部分仍然为宿主导向药物重定位研究提供了有价值的起点。
第四，在“Hub Gene Expression Stratifies Dengue Disease Severity”部分，研究人员使用独立临床队列GSE17924验证关键基因与疾病严重程度的关系。结果显示，核心特征中的枢纽基因能够区分DF、DHF和DSS三个临床分组。值得注意的是，STAT1在临床重症组DSS中高于DF，这与体外内皮模型中观察到的下调模式并不一致；而CDKN1A和FOS则从DF到DHF再到DSS呈渐进性升高，提示应激相关转录程序会随着疾病加重而不断增强；JUNB在DSS组中相较DF和DHF也显著升高。研究人员据此认为，不同组织背景下的宿主反应存在明显情境依赖性：内皮模型主要反映局部细胞对病毒成分的反应，而全血转录组则更多反映系统性免疫激活状态。
在该部分中，研究人员还通过ROC分析评估这些枢纽基因的临床判别效能。STAT1、CDKN1A、FOS和JUNB在区分DF与DSS、DHF与DSS时均获得AUC=1.00，且灵敏度和特异度均为1.00，说明在该数据集中这些基因呈现出极强的分组区分能力。然而论文同时强调，由于每组样本量仅n=5，且组内变异较低，因此这种近乎完美的判别性能很可能被高估；尽管留一交叉验证（LOOCV）结果一致，但仍必须在更大规模独立队列中验证其稳健性和泛化能力。这种谨慎表述是论文讨论部分的重要特点，体现出研究人员对生物信息学结论边界的严格把握。
在讨论部分，论文系统总结了上述发现的生物学意义。首先，研究支持这样一种观点：登革病毒感染相关宿主反应的核心，不仅是炎症激活，更是干扰素信号与应激反应之间的协调失衡。研究人员认为，STAT1相关通路下调与NS1刺激模型、感染时间点差异及细胞类型差异均有关，也提示干扰素反应具有高度动态性，可能经历早期抗病毒激活与后期补偿性或失调性改变。其次，CDKN1A作为上调枢纽节点的识别，提示细胞应激、细胞周期停滞与内皮功能障碍之间可能存在机制联系；FOS和JUNB等早期反应转录因子上调，则提示炎症与应激适应程序被快速启动。再次，论文特别讨论了STAT1在发现队列与临床队列中的表达分歧，指出这并不构成结论冲突，而是揭示了局部内皮反应与系统免疫反应之间的生物学差异。最后，研究人员强调，本研究完全基于计算分析，网络枢纽并不等于因果调控因子，药物候选也不代表临床可用药物，因此所有发现仍需进一步实验和临床验证。
研究结论部分可译为：本研究通过整合性转录组分析识别出与登革病毒感染严重程度相关的一致性分子特征。通过结合多个实验数据集与临床验证，研究突出显示了以干扰素相关通路抑制以及应激反应和炎症机制激活为特征的保守基因网络。STAT1、CDKN1A、FOS和JUNB被一致识别为该特征的核心组成部分。它们在不同临床严重程度分组中的差异表达支持其与疾病进展的潜在相关性，尽管情境依赖性差异也凸显了登革病毒感染宿主反应的复杂性。最后，候选化合物的识别为宿主导向治疗策略中的假说生成提供了初步框架。
