睡眠障碍中的线粒体功能相关基因：一项多组学孟德尔随机化研究

《Journal of Molecular Neuroscience》：Mitochondrial Function-Related Genes in Sleep Disorders: A Multi-Omics Mendelian Randomization Study

　　

既往研究报道线粒体功能障碍与睡眠障碍相关，但特定基因的潜在作用仍不清楚。本研究旨在剖析不同亚型特异的遗传关联及其潜在机制。研究人员采用基于汇总数据的多组学孟德尔随机化（Summary-data-based Mendelian Randomization, SMR）方法，以识别线粒体功能相关基因与睡眠障碍之间的潜在因果关联。研究人员整合了来自FinnGen数据库的全基因组关联研究（Genome-Wide Association Study, GWAS）数据（发现集）、独立的GWAS数据集（涵盖不同睡眠障碍亚型并用于验证）以及顺式定量 trait 位点（cis-QTLs，包括mQTLs、eQTLs和pQTLs）进行系统性探索。特别地，研究人员利用基因型-组织表达（Genotype-Tissue Expression, GTEx）数据库内疾病相关脑区的基因表达数据，进行了组织特异性效应的靶向验证。SMR分析在多个生物学层面识别出可能调控睡眠障碍的线粒体功能相关基因，初步在甲基化水平鉴定出102个基因，在基因表达水平鉴定出48个基因，在蛋白质丰度水平鉴定出6个基因。随后的整合分析优先确定了DCXR和ACADVL，并揭示了它们不同的亚型特异性关联。DCXR在睡眠呼吸暂停中表现出保护作用，而ACADVL在日间嗜睡中表现出矛盾的风险促进作用。此外，分析还确定了DCXR的一种表观遗传调控机制，其表达和蛋白质水平受DNA甲基化调控。最后，在脑-下丘脑组织中的验证证实DCXR是一个显著的潜在保护因子（比值比OR = 0.929，95%置信区间CI: 0.887–0.973，P_HEIDI = 0.999，错误发现率FDR = 0.2449）。研究结果提示关键线粒体基因，尤其是DCXR和ACADVL，参与特定睡眠障碍亚型的病理生理过程，为精准医学提供了潜在途径。临床试验编号：不适用。

研究背景

睡眠障碍（包括失眠和阻塞性睡眠呼吸暂停[Obstructive Sleep Apnea, OSA]）表现为入睡或维持睡眠困难、日间疲劳及认知功能障碍等症状，在老年人和慢性病患者中尤为普遍，构成重大公共卫生挑战，显著降低患者生活质量并增加医疗负担。现有治疗策略包括失眠认知行为疗法（Cognitive Behavioral Therapy for Insomnia, CBT-I）和OSA的正压通气治疗等，但远期疗效仍不理想。线粒体作为通过氧化磷酸化产生能量的关键细胞器，在细胞代谢和稳态维持中发挥核心作用。近年来，线粒体功能障碍日益被发现在睡眠障碍的病理生理过程中扮演重要角色：线粒体活性氧（mitochondrial reactive oxygen species, mtROS）的积累可通过调节特定睡眠促进神经元的兴奋性影响睡眠稳态；孤立性快速眼动睡眠行为障碍中的皮质变薄也与线粒体功能相关基因表达模式显著关联。然而，特定线粒体功能相关基因在睡眠障碍中的潜在作用仍不清楚。传统的孟德尔随机化（Mendelian Randomization, MR）分析通常依赖单一组学数据（如基因表达），难以全面揭示疾病背后的复杂机制。为此，研究人员开展本研究，旨在通过多组学MR框架，首次整合DNA甲基化、基因表达和蛋白质丰度三个层面的数据，阐明线粒体功能相关基因在睡眠障碍发病中的潜在因果作用，探索睡眠-觉醒调控的潜在中枢机制，并揭示特定基因的新型亚型依赖性致病作用。

研究概述与意义

研究人员利用基于汇总数据的孟德尔随机化（Summary-data-based Mendelian Randomization, SMR）方法，整合FinnGen数据库的GWAS数据（发现集）、三个独立亚型验证队列以及多种顺式QTL数据，系统探索了线粒体功能相关基因与睡眠障碍的因果关联，发现DCXR和ACADVL具有亚型特异的差异性作用，其中DCXR在睡眠呼吸暂停中发挥保护作用，为精准医学干预提供了潜在靶点。该研究发表在《Journal of Molecular Neuroscience》。

主要技术方法

研究人员应用SMR方法，以线粒体基因顺式QTL（包括mQTL、eQTL和pQTL）作为工具变量，结合异质性检验（Heterogeneity in Dependent Instruments, HEIDI）排除连锁不平衡混杂，并采用Benjamini-Hochberg错误发现率（False Discovery Rate, FDR）校正控制假阳性，通过coloc共定位分析评估不同组学层面共享因果变异。数据来源：发现集为FinnGen数据库睡眠障碍表型（49,880例病例和358,194例对照）；验证集包括Lane队列（失眠）、DaytimeSleepiness数据集（日间嗜睡）和Campos队列（睡眠呼吸暂停）；QTL数据来自eQTLGen联盟、Pietzner等以及Brisbane/Lothian队列；组织特异性验证使用GTEx v8数据库脑-下丘脑组织数据。

研究结果

多组学SMR鉴定候选基因

在FinnGen发现队列中，多组学SMR分析依据预设标准（P_SMR-multi < 0.05且P_SMR < 0.05，P_HEIDI > 0.01）鉴定出大量与睡眠障碍表型相关的基因。甲基化水平鉴定出102个基因共187个显著位点，其中ABAT、C5orf63、COX6B2等为风险因素，ACADVL、CAT、CPT1C等为保护因素；表达水平发现48个显著基因，NAGS、BCL2L11、FDPS为风险因素，ACADVL、DCXR、ISCU为保护因素；蛋白质水平发现6个相关蛋白，BOLA3和PC为风险因素，ACADVL、CAT、DCXR、NFU1为保护因素。整合分析揭示了DCXR的表观遗传调控机制：cg07073120位点高甲基化导致基因表达降低，进而降低蛋白质丰度，从而增加睡眠障碍风险；ACADVL则呈现相反的调控模式。

亚型特异队列中的候选基因验证

在失眠队列中，甲基化和表达水平的关联得到部分验证，NAGS基因高表达被证实与失眠风险增加相关（OR = 1.047，95% CI: 1.001–1.095，P_HEIDI = 0.433，P_SMR = 0.046）。在日间嗜睡队列中，NDUFB8表达和BCL2、MRRF甲基化位点的保护关联与发现队列一致，但ACADVL在mQTL、eQTL和pQTL三个层面的效应方向均发生逆转，高水平ACADVL（高甲基化OR = 1.017，高表达OR = 1.01，高蛋白丰度OR = 1.034）均与日间嗜睡风险增加相关。在睡眠呼吸暂停队列中，DCXR的保护作用得到强力验证：高基因表达（OR = 0.978，95% CI: 0.964–0.992，FDR = 0.041）和高蛋白水平（OR = 0.927，95% CI: 0.878–0.979，FDR = 0.018）均与风险降低相关，同时ISCU、NDUFB8和NFU1的保护作用也得到验证。

下丘脑组织中的关键发现验证

利用脑-下丘脑组织eQTL数据进行的组织特异性验证中，DCXR、MRPL43、MTRF1L和SETD9被验证具有显著eQTL效应，其中DCXR被确认为重要的潜在保护因子（OR = 0.929，95% CI: 0.887–0.973，P_HEIDI = 0.999，FDR = 0.2449）。

讨论与结论

讨论部分指出，DCXR作为参与碳水化合物代谢和精子成熟的多功能蛋白，在发现队列中显示甲基化水平与睡眠障碍风险正相关、表达和蛋白质水平与风险负相关，整合分析证实cg07073120高甲基化降低DCXR表达和蛋白质水平，睡眠呼吸暂停队列的验证强烈支持这一调控级联机制。ACADVL作为线粒体β-氧化途径中催化极长链脂肪酸（Very-Long-Chain Fatty Acids, VLCFAs）分解初始步骤的关键酶，对维持线粒体能量代谢至关重要，在总体睡眠障碍表型中表现为保护因素，但在日间嗜睡亚型中效应发生逆转。这一矛盾可能源于广义睡眠障碍表型（通常反映夜间睡眠紊乱）与其特定后果日间嗜睡之间的概念区分：针对夜间睡眠紊乱根本病因具有保护作用的遗传变异，反而可能表现为其最突出症状（日间嗜睡）的风险因素。此外，NDUFB8（呼吸链复合物I亚基）和ISCU（铁硫簇生物合成所必需）对睡眠呼吸暂停具有保护作用，COX19（复合物IV组装因子）与失眠特异性关联，提示不同睡眠障碍亚型可能与线粒体呼吸链不同复合物的易损性相关。研究局限性包括：发现表型异质性可能掩盖亚型特异效应、部分验证队列基于自我报告数据易受回忆偏倚影响、主要采用血液QTL数据可能无法完全反映中枢神经系统的生物学过程、样本均来自欧洲人群限制外推性、以及基于汇总数据无法进行年龄性别等分层分析。

结论：本研究强调了线粒体相关基因在睡眠障碍发展中的重要作用，并揭示其贡献具有高度亚型特异性。研究为DCXR在睡眠呼吸暂停中的保护作用提供了稳健的遗传学证据，同时揭示了ACADVL的复杂情境依赖作用——其对广义睡眠障碍表型具有保护作用，却是日间嗜睡的风险因素。基于这些发现，未来研究应超越对这些疾病的广义理解，深入阐明线粒体功能障碍的亚型特异性生物学机制，为开发新型、个体化的预防和治疗策略奠定基础。