《Antioxidants》：Convergent Lower Expression of Redox-Linked Stress-Adaptation and Synaptic-Plasticity Genes in Major Depressive Disorder Across Seven Postmortem dlPFC Cohorts
编辑推荐：
重性抑郁障碍（major depressive disorder，MDD）与氧化应激、线粒体功能障碍以及神经元可塑性受损有关，但相关转录组改变在不同死后脑组织队列中的可重复性仍不确定。研究人员对预先设定的14基因面板进行了靶向跨平台分析，该面板涵盖抗氧化防御、
重性抑郁障碍（major depressive disorder，MDD）与氧化应激、线粒体功能障碍以及神经元可塑性受损有关，但相关转录组改变在不同死后脑组织队列中的可重复性仍不确定。研究人员对预先设定的14基因面板进行了靶向跨平台分析，该面板涵盖抗氧化防御、线粒体氧化还原调控、细胞应激反应、神经营养信号传导、突触可塑性以及多胺代谢；分析对象为7个死后背外侧前额叶皮层（dorsolateral prefrontal cortex，dlPFC）队列，共包括146例MDD病例和179例对照。主要支持证据需同时满足以下条件：Fisher合并证据、在整个面板范围内进行Benjamini–Hochberg校正，以及在所有可用队列中MDD减对照（MDD-minus-control）方向一致。NPTX2、EGR1、VGF、BDNF和SAT1符合这些标准，且在MDD中均呈较低表达。同样的五基因模式还得到了加权有符号Stouffer分析、单阶段广义最小二乘（generalized least-squares，GLS）模型、随机效应Meta分析以及200,000次疾病标签置换分析的支持；在这些置换中，没有任何一次产生至少5个满足完整主要支持标准的基因（经验性p = 5.0 × 10?6）。最稳健的跨队列发现是：支持氧化还原相关应激适应、多胺稳态、神经营养信号传导、活动依赖性转录以及突触可塑性的基因，呈现一致的低表达模式。该模式提示，MDD中神经元应激韧性与适应性可塑性的分子能力受损。
本文发表于《Antioxidants》，围绕重性抑郁障碍（major depressive disorder，MDD）死后背外侧前额叶皮层（dorsolateral prefrontal cortex，dlPFC）中的可重复分子异常展开，重点检验与氧化还原调控、细胞应激适应、神经营养支持和突触可塑性相关的候选基因，是否能够在多个独立或半独立死后队列中表现出稳定一致的表达变化。研究背景在于，MDD的病理生物学已不再被视为单一神经递质异常，而被理解为涉及神经元适应受损、细胞应激调控紊乱、皮层环路可塑性障碍以及线粒体和氧化还原稳态失衡的复杂疾病。既往研究虽多次提示氧化应激、线粒体功能障碍、神经营养信号异常和活动依赖性转录受损与MDD有关，但死后脑转录组研究的一个核心难题是跨数据集可重复性不足。不同研究在样本构成、诊断流程、药物暴露、濒死状态、死后间隔、组织pH值、RNA完整性、检测平台以及预处理策略等方面存在明显异质性，使单一队列中观察到的表达信号难以确认为稳定的疾病相关改变。因此，开展这项研究的必要性在于，从“跨队列可重复性”而非“单队列显著性”出发，筛选真正稳定的MDD相关分子信号。
研究人员据此采用预先固定的14基因候选面板，覆盖抗氧化防御、线粒体-氧化还原生物学、应激反应、神经营养与突触可塑性、多胺代谢、细胞外基质重塑及转录调控等生物学领域，在7个dlPFC死后分析队列中进行靶向跨平台整合分析。结果显示，NPTX2、EGR1、VGF、BDNF和SAT1在MDD中表现出跨队列一致的低表达，且这一五基因模式经多种敏感性分析均得到支持。研究的重要意义在于，最稳定的异常并非经典抗氧化酶转录本的普遍抑制，而是一个连接氧化还原敏感性应激适应、神经营养支持、活动依赖性转录、突触组织和多胺稳态的收敛性低表达模式。这提示MDD的关键分子特征更可能是神经元应激韧性和适应性可塑性能力下降，而非单一路径的孤立缺陷。
本研究主要采用以下关键技术方法：研究对象来自6个独立公共项目形成的7个死后dlPFC分析队列，包括GSE54567、GSE54568、GSE102556、GSE101521、GSE53987、GSE208338和GSE213982，共325名供体。研究人员在各队列内提取14个候选基因表达值，采用配对t检验或Welch双样本t检验获得队列内效应；以Fisher法合并p值，并在14基因面板内进行Benjamini–Hochberg校正，以方向一致性作为主要支持标准。随后应用样本量加权有符号Stouffer分析、单阶段广义最小二乘（GLS）模型、随机效应Meta分析以及200,000次标签置换检验进行稳健性评估；并利用GSE213982的单核RNA测序（single-nucleus RNA-seq，snRNA-seq）数据进行候选基因细胞类型定位描述。
以下为研究结果解读。
3.1. Dataset Overview
主要分析纳入6个独立公共死后前额叶皮层研究形成的7个队列，共146例MDD和179例对照，合计325名供体。该部分结果说明本研究的证据基础来源于多个公开数据集的跨队列整合，而非单一样本来源，为后续评价可重复性提供了样本与平台多样性基础。
3.2. Primary Cross-Cohort Support
研究人员设定的主要支持标准要求同时满足：Fisher合并检验在14基因面板校正后达到显著，且在所有可用队列中MDD相对对照的效应方向完全一致。结果显示，仅有5个基因达到这一严格标准，即EGR1、NPTX2、VGF、BDNF和SAT1，并且它们在MDD中均为低表达。EGR1、NPTX2、VGF和BDNF在全部7个队列中方向一致下调；SAT1因GSE208338缺少处理后表达值，仅在6个可用队列中表现出一致下调。进一步的样本量加权有符号Stouffer分析、单阶段GLS分析同样支持这5个基因，且诊断-队列交互作用在多重校正后均不显著，说明其效应大小虽可能略有波动，但不存在某个队列单独驱动结果的证据。200,000次疾病标签置换中，没有一次出现至少5个基因同时满足完整主要支持标准，经验性p = 5.0 × 10
?6，表明该五基因模式并非随机标签重分配所致。
3.3. Genes with Partial or Discordant Support
该部分比较了未达到完整支持标准的候选基因。GPX4虽然在Fisher合并检验并经校正后仍显著，但不同队列间方向不一致，因此未被归类为主要支持基因。SKA2、ADAMTS16、PRDX2和SOD2具有名义显著或接近阈值的聚合证据，但同样缺乏完全的方向一致性；其余基因则表现为校正后不显著、方向不一致，或两者兼有。该结果表明，仅凭生物学合理性或合并显著性不足以认定为稳定的MDD相关信号，必须同时考虑跨队列方向收敛。
3.4. One-Stage Individual-Level Sensitivity Analysis
在对各基因、各队列内部标准化后，研究人员构建单阶段GLS模型，从个体层面估计共同诊断效应。结果表明，5个主要支持基因在MDD中均显著下调，标准化诊断系数范围为?0.435至?0.557，经Benjamini–Hochberg校正后均保留显著性。同时，诊断-队列交互作用在多重校正后不显著。这说明从个体层级、统一建模框架出发，也可重复得到同样的五基因低表达模式。
3.5. Random-Effects Meta-Analysis of Standardized Effects
随机效应Meta分析进一步验证了上述结果。NPTX2、EGR1、VGF、BDNF和SAT1的合并Hedges g分别为?0.551、?0.559、?0.490、?0.415和?0.423，95%置信区间均不跨0，提示MDD中这些基因表达稳定降低。异质性方面，NPTX2、EGR1、BDNF和SAT1的I
2为0–1.8%，几乎无明显队列间异质性；VGF的I
2为30.6%，属中等异质性。该部分结果说明五基因信号不仅在方向上收敛，且效应量层面也具有较好的跨队列稳定性。
3.6. Cell-Type Localization of Candidate Genes in the GSE213982 snRNA-Seq dlPFC Resource
研究人员进一步利用GSE213982的dlPFC单核RNA测序资源对5个主要支持基因进行细胞类型定位描述。结果显示，这5个基因中的大多数观测计数主要分布于兴奋性和抑制性神经元核，提示其表达具有明显神经元定位背景。不过作者明确指出，该分析仅用于提供细胞定位语境，不能被解释为细胞类型特异性差异表达的独立验证，也不能区分真实转录调控与组织细胞组成变化的贡献。
在讨论部分，研究人员首先强调，本研究最核心的发现是：MDD在dlPFC中的最可重复表达异常，并非经典抗氧化酶普遍下降，而是由BDNF、VGF、EGR1、NPTX2和SAT1构成的收敛性低表达模块。这一模块分别对应神经营养支持、活动依赖性转录、突触组织和多胺代谢，提示MDD中受损的是神经元长期应激适应与适应性可塑性的分子能力。BDNF、VGF、EGR1和NPTX2被视为一致的神经元可塑性相关模块：BDNF参与神经元存活、生长、突触发生和活动依赖性可塑性；VGF与神经营养和突触可塑性密切相关；EGR1代表即刻早期基因介导的活动依赖性转录调控；NPTX2参与AMPA受体（α-氨基-3-羟基-5-甲基-4-异噁唑丙酸受体）依赖的兴奋性突触组织以及对PV中间神经元（parvalbumin-positive interneurons，表达小清蛋白的中间神经元）的兴奋性募集。SAT1则将这一模式从经典突触和神经营养机制延伸至多胺周转与氧化还原稳态调节。研究同时指出，GPX4、PRDX2、PRDX3、SOD2等经典抗氧化酶相关基因未表现出稳定一致方向，因此当前数据不支持“MDD存在广泛统一的抗氧化防御转录抑制”这一简单模型，而更倾向于认为抗氧化相关反应具有情境依赖性和异质性。文章还讨论了探索性性别分析中GPX4和MTOR可能存在的男性相关模式，但由于诊断-性别交互在多重校正后未达显著，这些结果仅具有假设生成意义。
此外，作者强调本研究优势在于优先考察跨异质数据集的可重复性，并通过Fisher分析、Stouffer分析、GLS、Meta分析与置换检验形成方法学交叉验证。单核数据提供了神经元定位背景，但由于主体分析基于整体组织，低表达信号也可能受到细胞组成差异影响。研究局限包括：队列间临床与技术异质性较高，无法进行完全统一的协变量校正；整体组织数据不能区分细胞内在转录变化与细胞构成变化；转录本水平不能直接代表蛋白丰度、酶活性、突触功能或真实氧化还原状态；14基因面板范围有限，不能代表完整MDD转录组；此外，GSE54567与GSE54568虽为不同样本配对集，但来自同一母研究与平台，并非完全独立的技术复制；横断面死后设计亦无法进行因果推断。
研究结论部分可译为：经典抗氧化酶转录本未显示出一致的跨队列支持，这不支持MDD存在普遍性抗氧化缺陷或稳定的氧化应激转录组特征。它们不一致的表达模式更可能反映对慢性细胞应激的异质性、情境依赖性适应或代偿反应，尽管仅凭当前死后整体组织数据尚无法证实这一解释。在7个死后dlPFC队列中，MDD与一种可重复的低表达模式相关，该模式连接了氧化还原敏感性应激适应、神经营养信号传导、活动依赖性转录、突触组织和多胺代谢。NPTX2、EGR1、VGF、BDNF和SAT1的收敛性变化与MDD中神经元韧性和适应性可塑性的分子能力下降相一致。这些发现仍需进一步通过功能学和细胞类型特异性研究加以验证。
