多组学解析睡眠特征与癫痫之间共享遗传结构

《Biology》：Multi-Omics Dissection of the Shared Genetic Architecture Between Sleep Traits and Epilepsy

　　

癫痫是一种遗传性神经系统疾病，常与睡眠困难（包括短/长睡眠时间和失眠）共病。尽管流行病学研究一致报道了睡眠障碍与癫痫之间的共病性，但两者之间的共享遗传结构（shared genetic architecture）和分子机制仍不清楚，阻碍了治疗进展。在本研究中，研究人员整合了欧洲血统的大规模全基因组关联研究（GWAS）汇总统计，以解析睡眠特征与癫痫之间的遗传和分子联系。利用连锁不平衡分数回归（LDSC）和GWAS-pw方法，研究人员识别出睡眠行为与癫痫之间存在适度但具有统计学显著性（经Bonferroni校正）的全局和局部遗传相关性。随后的CPASSOC（Cross-Phenotype Association）跨性状元分析和转录组全关联研究（TWAS）精确定位了特定的多效性位点和共享候选基因，包括SPAG7、VRK2和LINC00925，这些基因在功能上与神经免疫信号相关。尽管对这些候选靶点进行的初步表型组全关联研究（PheWAS）分析未在当前数据库中识别出主要的不良关联，但研究人员强调，在将其纳入治疗策略之前，需要进行严格的体外和体内实验验证。最后，多效性稳健的双向孟德尔随机化（MR）分析提示，从癫痫到短睡眠时间存在单向因果效应。尽管估计的因果效应量很小，但它反映了终生的多基因结构，而非急性临床幅度。总之，研究人员采用的多组学方法揭示了睡眠-癫痫轴的共享遗传结构，并突出了未来功能研究的潜在生物标志物。

癫痫是一种常见且致残的神经系统疾病，给社会带来沉重健康和经济负担。睡眠特征（如短睡眠、长睡眠和失眠）在维持生理和神经稳态中起关键作用。既往流行病学及观察性基因组研究已建立癫痫与睡眠障碍的相关性，但二者共享的遗传结构和分子机制仍不明确。为填补这一空白，研究人员整合了欧洲血统的大规模全基因组关联研究（GWAS）汇总统计，从多组学角度系统解析睡眠行为与癫痫的遗传关系。研究发现了适度但显著的全局及局部遗传相关性，鉴定了多效性位点和共享候选基因，并通过双向孟德尔随机化（MR）提示癫痫对短睡眠的单向因果效应。该论文发表在《Biology》。

关键技术方法包括：利用连锁不平衡分数回归（LDSC）估算全基因组及功能分区遗传相关性；运用GWAS-pw评估1703个连锁不平衡（LD）独立基因组区域的局部遗传重叠；采用CPASSOC（Cross-Phenotype Association）算法进行跨性状元分析，以识别多效性位点；借助FUSION工具整合GTEx v8表达权重开展转录组全关联研究（TWAS），鉴定基因表达层面的共享基因；使用双向两样本MR（以逆方差加权为IVW主要方法）推断因果关系。样本来源：睡眠数据来自英国生物银行（UK Biobank）欧洲血统参与者（短睡眠106,192例、长睡眠34,184例、失眠109,548例，对照分别为305,742例、279,327例和277,440例）；癫痫数据来自欧洲血统GWAS（27,559例病例，42,436例对照）。

**3.1 全基因组遗传力和全局遗传相关**
通过LDSC估计单核苷酸多态性（SNP）遗传力：短睡眠0.0495、长睡眠0.029、失眠0.046、癫痫0.1049。全局遗传相关分析显示，无论采用无约束还是约束LDSC模型，短睡眠与癫痫均呈显著正相关（无约束：r_g=0.181，95% CI: [0.077, 0.285]，p=7.00×10^?4；约束：r_g=0.1098，95% CI: [0.067, 0.153]，p=5.76×10^?7）。长睡眠与癫痫无显著相关。失眠与癫痫仅在约束模型中显著（r_g=0.0654，95% CI: [0.026, 0.105]，p=1.10×10^?3）。

**3.2 功能类别分区遗传相关**
分区LDSC（经Bonferroni校正p<0.05/14）表明，短睡眠与癫痫在编码区（最强，约束模型r_g=0.3508）、转录区、保守区等多个功能类别中显著相关；失眠与癫痫在保守区、DGF、DHS等类别显著；长睡眠与癫痫在任何类别均无显著相关。

**3.3 局部遗传相关**
GWAS-pw评估1703个LD独立区域时发现，长睡眠与癫痫存在超过10个共享因果变异区域（PPA3>0.5），如chr2:113,929,896-116,771,011等；短睡眠和失眠与癫痫未检测到共享区域（PPA3>0.5），但存在多个独立关联区域（PPA4>0.5）。

**3.4 跨性状元分析鉴定多效性位点**
CPASSOC（P_CPASSOC<5×10^?8且P_单性状<0.05）在剪枝后分别鉴定出：短睡眠-癫痫9个独立多效位点（如rs7556815近LOC101927400、rs2717076近SCN1A、rs60055328近TTC21B等）；长睡眠-癫痫6个位点（如rs2678901近VRK2、rs7534398近CAMTA1等）；失眠-癫痫4个位点（如rs6561715近CNIH2和YIF1A、rs57272144近LINC00925等）。关键区域如2p16.1、2q13、2q24.3在多个性状对中共有。

**3.5 通路富集**
Metascape富集分析显示：短睡眠-癫痫显著富集于心肌细胞动作电位、上皮管形成、氧化应激反应、微管运输和小分子生物合成过程；长睡眠-癫痫富集于昼夜节律、行为、细胞器生物发生、高尔基体囊泡运输、蛋白质稳定性调控和区域化；失眠-癫痫仅富集于膜运输（R-HSA-199991）。

**3.6 转录组全关联研究（TWAS）**
TWAS（FDR校正p<0.05）在脑组织中鉴定出短睡眠与癫痫的共享基因SPAG7（小脑半球、小脑、前额叶皮层BA9、海马）和VRK2（黑质）；失眠与癫痫的共享基因LINC00925（全血）；长睡眠与癫痫未发现共享基因。

**3.7 遗传药物靶点**
通过DGIdb和DSigDB查询，SPAG7对应伏立诺他和丙戊酸；VRK2对应达沙替尼、博舒替尼、氯唑沙宗和伏立诺他；LINC00925作为长链非编码RNA（lncRNA）无候选药物。

**3.8 表型组全关联研究（PheWAS）**
基于UK Biobank的AstraZeneca PheWAS门户分析显示，SPAG7和VRK2与感染、肿瘤、血液/免疫、内分泌/代谢、精神、神经、心血管等广泛表型类别无显著关联。

**3.9 孟德尔随机化分析**
双向两样本MR使用F统计量>10的仪器变量，IVW为主要方法。短睡眠不显著导致癫痫（IVW: BETA= -0.0005, p=0.0856），但癫痫显著增加短睡眠风险（IVW: BETA=0.0006, p=1.1×10^?7；加权中位数BETA=0.00054, p=1.1×10^?5）。长睡眠与失眠均与癫痫无相互因果效应（所有IVW p>0.05）。MR-Egger截距检验未提示水平多效性（所有p>0.05），Cochran's Q检验在长睡眠分析中显示显著异质性时采用随机效应IVW。

**讨论与结论**
讨论部分指出，本研究通过整合多组学方法提供了睡眠特征与癫痫共享遗传结构的统计证据。全局及分区遗传相关确认了短睡眠与癫痫的显著重叠，局部分析发现长睡眠的多区域共享，CPASSOC和TWAS共同鉴定了SPAG7、VRK2、LINC00925等候选基因。这些基因在神经免疫、突触调控和细胞内运输中发挥功能，与海马、黑质等脑区相关。MR结果支持癫痫对短睡眠的单向因果效应，但效应量极小（BETA=0.0006），反映终生的多基因易感性而非急性临床效应。同时，研究局限性包括仅限于欧洲血统、失眠表型基于单一问卷条目、缺乏纵向数据和环境协变量、未进行多队列重复等。结论部分翻译如下：总之，研究人员采用的多组学分析为睡眠特征与癫痫之间的共享遗传架构提供了稳健的统计证据。通过整合跨性状元分析和TWAS框架，识别出显著的遗传相关性和特定共享位点。值得注意的是，SPAG7和VRK2作为涉及神经免疫和突触调控的多效性候选基因，值得进一步功能研究。此外，双向MR提示从癫痫到短睡眠时间存在单向因果效应；然而，极小的效应量（BETA=0.0006）反映的是微妙的多基因影响而非急性临床驱动因素。必须结合若干关键局限性来解释这些发现，包括严格限制于欧洲血统数据集、依赖粗略的自我报告睡眠表型、使用汇总级统计数据以及缺乏独立队列重复。因此，虽然本研究突出了共享的遗传通路，但强调这些统计关联并不能直接转化为临床应用。在任何候选基因被考虑用于未来治疗策略之前，必须进行严格的体外和体内实验验证。