带状疱疹与带状疱疹后神经痛的血转录组免疫特征：基于干扰素通路模式的平行队列分析

《Frontiers in Cellular and Infection Microbiology》：Blood transcriptomic immune signatures in herpes zoster and postherpetic neuralgia: parallel cohort analyses of IFN-related pathway patterns

　　

摘要：带状疱疹（Herpes Zoster, HZ）由潜伏的水痘-带状疱疹病毒（Varicella-Zoster Virus, VZV）再激活引起，而带状疱疹后神经痛（Postherpetic Neuralgia, PHN）是其主要的慢性疼痛并发症。研究人员对两个独立的公共血转录组队列进行了平行分析：GSE242252（n=80；包括HZ发病期、1年恢复期及对照组）和PRJNA1142765（n=23；包括15例HZ和8例PHN）。由于两个队列在测序平台、文库制备、实验设计和采样结构上存在差异，表达矩阵未进行合并；跨队列的解读仅限于队列内部的结果及通路水平的方向一致性。HZ发病期表现出强烈的以I型干扰素（Interferon, IFN）为主的免疫激活（HALLMARK_INTERFERON_ALPHA_RESPONSE，标准化富集分数[NES]=2.21，校正后p=4.4×10-10），并且恢复期伴随样本水平的IFN评分下降（发病期 vs. 恢复期，p=0.005）。在较小的PHN队列中，仅有2个基因达到严格的差异表达阈值，但Hallmark基因集富集分析（GSEA）显示IFN-α和IFN-γ通路呈正向富集。排序敏感性分析证实，无论使用log2倍数变化还是DESeq2 Wald统计量作为排序依据，IFN-α和IFN-γ均为排名前两位的正向富集Hallmark通路，其前导边缘Jaccard相似度分别为0.848和0.897。样本水平IFN评分的精确患者标签置换检验未达到统计学显著性（双侧p=0.297），但在所有23次留一法迭代中，PHN减去HZ的差值方向始终保持为正。因此，新的生物学信息在于证据的不对称性：急性HZ的IFN激活及其恢复期的减弱得到了强有力的支持，而PHN相关的IFN模式在通路排序水平上是稳定的，但在患者水平上则具有统计不确定性。据此，先前关于PHN中存在持续性IFN激活的解读应被收窄而非强化。一项结构化的外部重复性检索未能找到包含PHN样本的合适的独立公开人类血液/外周血单个核细胞（PBMC）转录组队列；因此，外部证据仅被视为基于文献的一致性支持和辅助性背景，而非验证。受限的计算优先级筛选突出了IFN相关及CCL5/CCR5相关轴以供未来研究。

论文解读：带状疱疹与带状疱疹后神经痛的血转录组免疫特征——基于干扰素通路模式的平行队列分析

1. 研究背景、现存问题与研究目的

带状疱疹（Herpes Zoster, HZ）是由潜伏在感觉神经节中的水痘-带状疱疹病毒（Varicella-Zoster Virus, VZV）再激活引起的疾病，全球发病率约为每年每千人3-5例，且随年龄增长显著上升。虽然急性HZ通常具有自限性，但其最常见的并发症——带状疱疹后神经痛（Postherpetic Neuralgia, PHN），即皮疹愈合后持续超过三个月的疼痛，发生于10%-30%的HZ患者中，是最常见且最令人困扰的并发症。HZ和PHN的免疫生物学涉及全身性抗病毒免疫和外周感觉组织内的病理过程。尽管已有研究开始描绘HZ和PHN的外周血及免疫细胞特征，但仍有两个关键问题悬而未决：一是急性抗病毒程序在恢复期减弱的程度有多强；二是当面对小样本量的稳健性检验时，PHN相较于HZ的血液信号有多少能保持稳定。因此，本研究旨在通过平行分析两个独立的公共血转录组数据集，分别刻画HZ发病/恢复期以及PHN与HZ对比的免疫模式，并重点评估PHN相关IFN信号的稳健性，从而明确不同疾病背景下证据强度的不对称性。

2. 研究方法概述

研究人员使用了两个独立的公共数据集进行平行分析。第一个队列GSE242252（来源于Gene Expression Omnibus）包含80份样本，涵盖HZ发病期、一年恢复期及对照人群，采用QuantSeq 3‘ mRNA测序技术。第二个队列PRJNA1142765（来源于Sequence Read Archive）包含23份全血样本（15例HZ和8例PHN），来自已发表的研究（Wang et al., 2024），采用Illumina链特异性mRNA文库和双端测序。由于两个队列在测序平台、文库制备、实验设计等方面存在本质差异，研究人员未合并表达矩阵，也未进行跨研究的批次校正或荟萃分析，而是将跨队列解读严格限制在通路水平和调控因子方向的描述性一致上。主要分析方法包括：使用DESeq2进行差异表达分析，使用fgsea进行基因集富集分析（GSEA），使用GSVA包进行单样本基因集富集分析（ssGSEA）以估算免疫评分，使用xCell和MCP-counter进行细胞类型反卷积，使用WGCNA进行加权基因共表达网络分析，使用NicheNet从批量转录组特征推断配体-靶标调控关系，以及使用DoRothEA/decoupleR估计转录因子活性。针对PHN队列的小样本量，研究人员还进行了排序敏感性分析、精确患者标签置换检验、自助法置信区间估计和留一法分析等稳健性检验。

3. 研究结果

3.1 HZ发病期呈现以I型IFN为主导的免疫激活

通过对HZ发病期与所有对照组的差异表达分析，发现37个上调基因，包括IFN/抗病毒基因（如IFI27、ISG15、RSAD2）和细胞周期基因。基因集富集分析（GSEA）显示，IFN-α反应、IFN-γ反应、E2F靶点和G2/M检查点等Hallmark通路显著富集，而血红素代谢通路显著下调。细胞类型反卷积分析发现，HZ发病期浆细胞和Th2细胞比例显著增加。IFN评分在HZ发病期显著高于对照组和恢复期。

3.2 WGCNA鉴定出一个富含IFN的共表达模块

在GSE242252队列中进行的WGCNA分析鉴定出18个共表达模块，其中MEred模块与HZ发病期显著相关，并富含IFN/抗病毒基因。该模块的核心枢纽基因包括STAT1、RSAD2、GBP1、IFIH1等，其基因本体（GO）富集分析主要集中在I型IFN反应、抗病毒防御和免疫调节等生物过程。

3.3 PHN与HZ的差异主要体现在通路水平而非单基因水平

在PRJNA1142765队列中，PHN与HZ的对比在严格的差异表达阈值下仅发现极少数显著差异基因。然而，校正后的Hallmark GSEA显示，IFN-α反应和IFN-γ反应在PHN中呈显著正向富集。排序敏感性分析表明，无论使用log2倍数变化还是Wald统计量进行基因排序，IFN-α和IFN-γ始终是排名前两位的正向富集通路，且前导边缘重叠度高。但在样本水平的IFN评分分析中，精确的患者标签置换检验不显著，95%自助法置信区间包含零值。留一法分析显示，去除任一单个样本后，PHN与HZ的IFN评分差值方向均保持不变。综合来看，PHN相关的IFN信号被归类为部分稳健：在通路排序水平稳定且方向不受单个样本驱动，但在患者评分水平统计精度不足。

3.4 平行队列比较显示IFN相关方向的一致性

在两个队列各自的比较中，IFN评分在HZ发病期高于对照组，且在PHN中方向性地高于HZ。浆细胞和巨噬细胞的估计值也呈现方向性升高。转录因子活性分析发现，STAT2、IRF9和IRF1在两个比较中均显示出方向一致的活性变化趋势。

3.5 批量来源的配体-靶标推断凸显了PHN对比HZ中的候选调控轴

NicheNet分析在PHN对比HZ中识别出大量活跃的配体-靶标链接，预测的前列配体包括IL27、IFNG、IFITM1、IL1B等，突显了IFN相关、炎症和抗原呈递相关的调控轴。

3.6 文献一致性与外部重复性检索

一项结构化的多数据库检索未能找到适合直接重复PHN对比HZ结果的独立公开人类血液/PBMC转录组队列。因此，外部证据仅作为支持性的生物学背景，而非外部验证。

3.7 受限的计算优先级筛选确定了IFN和CCR5相关的候选轴

研究人员将药物靶向轴的优先级筛选范围缩小到ISG15、RSAD2、USP18、STAT2、IFNG/IL27和CCL5/CCR5这几个免疫学上合理且不过度外推的轴上。这些节点在多个分析层面（如HZ发病期差异表达基因、PHN探索性基因集、WGCNA枢纽基因）均有出现，并结合了外部文献支持。

4. 讨论与结论

讨论总结： 本研究的主要生物学贡献在于揭示了证据强度的不对称性：急性HZ的IFN激活及其在恢复期的减弱得到了独立纵向队列的有力支持，而PHN相关的IFN模式在通路排序水平上稳定，但在患者水平上统计不确定性高。因此，先前关于PHN中存在持续性IFN激活的解读应被收窄而非强化。研究强调，PHN是一种外周神经系统疼痛疾病，而所分析的数据来自全血，因此无法直接表征背根神经节、感觉神经元或局部神经-免疫相互作用，也不应推断任何神经机制或血液到神经元的因果通路。两个队列因平台、设计和临床采样结构不同而无法合并，跨队列比较仅限于描述性方向一致。PHN队列样本量小导致单基因显著性低和患者水平精度有限。此外，PHN对比HZ的横断面设计无法建立从急性HZ到PHN的时间持续性。本研究的优势在于提供了急性HZ及恢复期的有力队列内证据，对小PHN队列进行了压力测试，并明确区分了直接证据、探索性分析、支持性文献和不可用的验证，从而提高了可解释性和可重复性。

研究结论翻译： 本研究支持IFN主导的外周血免疫激活是急性HZ的一个主要特征，并在一个独立的纵向HZ队列中证明了其在恢复期的减弱。在单独的PHN对比HZ队列中，IFN-α和IFN-γ通路在敏感性分析中仍保持高排名，而患者水平的IFN评分置换检验不显著。因此，PHN的结果应被视为一个部分稳健但统计精度不足的通路水平模式，需要独立的重复验证。现有数据不能建立时间持续性、神经机制、向PHN的因果进展或经过验证的治疗靶点。受限的IFN相关和CCL5/CCR5相关节点最好被视为未来实验和前瞻性临床评估的假说。