多组学解码重性抑郁障碍中核编码线粒体基因的因果作用

《Psychological Medicine》：Decoding nuclear-encoded mitochondrial genes in major depressive disorder: A multi-omics perspective

本研究针对重性抑郁障碍（MDD）中线粒体功能障碍的因果机制尚不明确的问题，通过整合GWAS与mQTL、eQTL、pQTL等多组学数据，采用SMR和贝叶斯共定位分析，发现TDRKH、METAP1D等五个核编码线粒体基因在MDD中具有潜在的因果作用。该研究为理解MDD的线粒体病理生理学提供了新视角，并为未来治疗靶点开发指明了方向。

　　

在全球范围内，有超过3亿人受到重性抑郁障碍（Major Depressive Disorder, MDD）的影响，这种复杂的精神疾病带来了巨大的社会经济负担。然而，其多因素病因至今仍未完全阐明，这给诊断和治疗带来了巨大挑战。遗传学研究显示MDD具有30%-50%的遗传力。近年来，线粒体——这个细胞内的“能量工厂”——在MDD中的作用日益受到关注，它不仅涉及能量代谢，还与氧化应激、细胞存活、炎症和免疫等过程密切相关。尽管先前的研究已经在线粒体与MDD之间建立了一些联系，例如蛋白质组学研究发现患者血浆神经元细胞外囊泡中线粒体组分存在异常，但这些证据往往存在局限：研究设计多局限于单一组学层面、样本量较小、人群覆盖范围有限；多使用外周组织数据，与大脑直接相关的整合分析不足；因果推断不够稳健，缺乏强有力的共定位验证；跨队列和跨组织的可重复性验证稀疏；并且对核编码线粒体基因缺乏系统性的评估。因此，将这些零散的发现整合成一个清晰且可重复的致病模型变得十分困难。

为了弥补这些空白，研究人员开展了一项大规模、方法学强化的研究。他们整合了两个独立的欧洲MDD全基因组关联研究（Genome-Wide Association Study, GWAS）资源（来自精神病学基因组学联盟-Psychiatric Genomics Consortium, PGC和FinnGen研究）以及三类分子数量性状基因座（Quantitative Trait Locus, QTL）数据——mQTL（DNA甲基化）、eQTL（基因表达）和pQTL（蛋白质丰度）——旨在追踪同一基因从DNA甲基化到基因表达再到蛋白质水平的调控级联。研究核心采用了基于摘要数据的孟德尔随机化（Summary-based Mendelian Randomization, SMR）方法，该方法将顺式QTL（cis-QTL）作为工具变量，以减轻环境混杂和反向因果关系的影响。同时，研究结合了HEIDI检验和贝叶斯共定位分析，以判断分子信号与疾病关联是否共享同一个因果遗传变异，从而减少连锁不平衡（Linkage Disequilibrium, LD）导致的假阳性。与传统的两样本孟德尔随机化不同，SMR非常适合在有效工具变量较少但效应集中的情况下，进行高通量的基因水平QTL-疾病关系筛选。为了确保结果的可重复性和生物学相关性，研究在PGC和FinnGen之间进行了交叉验证，并将分析从外周血扩展到多个脑区的组织特异性QTL，在保持临床可转化性的同时，直接探究中枢神经系统（Central Nervous System, CNS）的病理机制。与先前多关注线粒体DNA（mitochondrial DNA, mtDNA）的研究不同，本研究系统性地靶向了MitoCarta3.0数据库中注释的核编码线粒体基因，这些基因直接调控线粒体功能，拥有更强的遗传工具，并且是更具可操作性的治疗靶点。这一设计将先前零散的单层信号整合成了一条连贯、可重复且机制上可解释的证据链。

本研究主要应用了以下几种关键技术方法：研究利用了大样本的GWAS汇总统计数据（PGC和FinnGen，共超过80万个体）以及公开的分子QTL数据集（包括超过3.2万样本的mQTL、超过3.1万样本的血液eQTL和超过5.4万样本的血浆pQTL）。通过SMR分析评估线粒体相关分子性状（甲基化、表达、蛋白）对MDD风险的潜在因果效应，并采用HEIDI检验排除连锁不平衡造成的假关联。进一步使用贝叶斯共定位分析（通过R包coloc实现）判断分子性状与MDD是否共享同一因果变异（以PPH4/(PPH3+PPH4) > 0.7为标准）。此外，还利用基因型-组织表达（Genotype-Tissue Expression, GTEx）计划的12个脑区eQTL数据评估了大脑特异性效应。线粒体基因的定义则基于MitoCarta3.0数据库。

多组学证据整合

通过多层SMR分析，研究确定了五个与MDD存在潜在因果关联的线粒体基因，并根据跨甲基化、转录和蛋白质水平的证据一致性、独立数据集验证以及共定位证据对这些基因进行了分级优先排序。

TDRKH被列为Tier 1基因。该基因内部的一个CpG位点（cg24503712）的甲基化与基因表达和蛋白质丰度降低显著相关，而这二者又与较低的MDD风险相关。这些发现在FinnGen数据集中得到了重复。甲基化、表达和蛋白质水平的关联效应方向一致，并且通过了共定位阈值和HEIDI异质性检验，支持了稳健的因果关系。

METAP1D被列为Tier 2基因。虽然在甲基化水平未观察到显著关联，但其转录本和蛋白质丰度均显示出对MDD的一致保护性效应，这些关联在独立数据集中得到了验证。

LONP1、FIS1和SCP2被归类为Tier 3基因。每个基因均在蛋白质水平显示出与抑郁风险的关联，并伴有单个数据集内的甲基化或基因表达关联。例如，FIS1的多个CpG位点（如cg17825709, cg25385322等）与表达呈负相关，与其保护效应一致。SCP2和LONP1则表现出更复杂的双向调控模式，不同CpG位点的甲基化对MDD风险的影响方向相反，提示存在情境特异性的调控效应。

10(p) 值。'>

线粒体基因甲基化与MDD

在排除P-HEIDI < 0.01的关联后，研究人员观察到614个CpG位点（涉及210个基因）与MDD存在名义显著关联（p < 0.05）。经过错误发现率（False Discovery Rate, FDR）校正后，保留了位于20个独特基因区域的46个CpG位点，其中43个信号显示出强共定位证据（PPH4 > 0.70），支持DNA甲基化与MDD之间存在共享的遗传调控。多个基因（如MSRA、SCP2、LONP1）的不同CpG位点表现出对MDD风险方向相反的影响，凸显了甲基化调控的位点特异性。TDRKH中cg24503712的保护性关联在FinnGen队列中得到了验证。

10(p)值代表了通过mQTL并与eQTL和pQTL数据整合后定位到的位点的关联强度。染色体交替着色以便视觉清晰。'>

线粒体基因表达与MDD

在eQTL-SMR分析中，排除P-HEIDI < 0.01的关联后，共有87个关联达到边际显著性。TDRKH和FIS1的表达水平与MDD风险呈正相关，而METAP1D、SCP2和LONP1的表达则与风险降低相关。COQ8A、NT5DC2和METAP1D的关联在FinnGen数据库中得到了重复。利用GTEx脑区eQTL数据的分析显示，TDRKH在尾状核中与MDD存在显著正相关，METAP1D在皮层中显示出保护性关联（仅见于FinnGen GWAS），而LONP1在多个脑区（如杏仁核、前扣带回皮层、海马体）显示出保护性关联（主要见于FinnGen GWAS）。这些大脑特异性关联的方向与血液分析中观察到的趋势基本一致。

线粒体蛋白与MDD

在蛋白质水平，经过筛选和共定位分析后，在PGC发现数据集中鉴定出11个线粒体蛋白与MDD风险存在名义显著关联（p < 0.05）。基因预测的METAP1D和SCP2蛋白丰度与MDD风险呈负相关，其中METAP1D的效应方向在FinnGen队列中得到了重复。在FinnGen中，LONP1蛋白水平与MDD风险负相关，而TDRKH和FIS1则显示正相关。

研究结论与讨论

这项研究通过多组学孟德尔随机化框架，为核编码线粒体基因在MDD中的潜在因果作用提供了证据。TDRKH作为一个突出的候选基因，其特定CpG位点（cg24503712）的高甲基化可能通过降低基因表达和蛋白质丰度，从而起到降低MDD风险的作用，这一调控级联在血液和大脑（如尾状核）组织中均得到部分验证。TDRKH编码一种位于线粒体膜的含有Tudor和KH（K Homology）结构域的蛋白，已知参与piRNA（PIWI-interacting RNA）生物生成，本研究提示了其在MDD中可能具有新的表观遗传作用。METAP1D在转录和蛋白水平均表现出对MDD的一致保护效应，该基因编码线粒体基质氨肽酶，负责新生蛋白的N端甲硫氨酸切除，对线粒体蛋白成熟和功能至关重要，尽管其大脑复制结果不一致，但仍是一个有前景的靶点。LONP1、FIS1和SCP2则显示出更复杂的调控模式。LONP1是线粒体基质中的ATP依赖性蛋白酶，负责维持蛋白质稳态，其不同CpG位点的甲基化对MDD风险的影响方向相反，但在FinnGen队列中，其表达和蛋白水平显示出保护性关联，并在多个脑区得到验证，提示其可能通过维持神经元线粒体功能完整性来对抗MDD。FIS1是位于线粒体外膜的裂变蛋白，其甲基化介导的下调可能具有保护作用，而较高的蛋白丰度则与风险增加相关，但缺乏大脑组织的显著关联，表明其作用可能具有组织特异性。SCP2参与细胞内脂质转运和代谢，其不同甲基化位点的效应存在异质性，且其保护性关联未能在FinnGen中重复，也未见于脑组织，因此其在MDD中的作用尚不确定，需要进一步验证。

本研究也存在一些局限性。首先，eQTL等数据主要来源于血液组织，限制了对大脑特异性线粒体调控的直接推断，尽管整合了脑QTL数据，但样本量有限和组织异质性可能导致复制不一致。其次，尽管使用了HEIDI过滤和贝叶斯共定位来减少假阳性，但不能完全排除水平多效性的可能性。第三，分析主要在欧洲人群中进行，限制了其向其他祖先人群的普适性。最后，缺乏实验验证，且研究焦点仅限于MitoCarta3.0注释的核编码线粒体基因，未涵盖mtDNA变异。因此，从这些统计关联中得出的机制解释应被视为假设，而非已确立的功能性后果。未来的实验研究对于证实这些发现至关重要，整合mtDNA变异、拷贝数和单倍群数据与核编码基因的研究将能更全面地揭示线粒体对MDD的贡献。

总之，这项研究为理解MDD的线粒体病理机制提供了新的多组学证据，突出了TDRKH、METAP1D等关键基因的重要性，为未来的机制探索和治疗干预提供了有价值的靶点。