溶酶体基因罕见变异负荷提示唾液酸化与神经节苷脂代谢参与帕金森病

Rare-Variant Burden across Lysosomal Genes Implicates Sialylation and Ganglioside Metabolism in Parkinson's Disease

作者信息Konstantin Senkevich, Sitki Cem Parlar, Cloe Chantereault, Lang Liu, Eric Yu, Uladzislau Rudakou, Jamil Ahmad, Jennifer A Ruskey, Farnaz Asayesh, Dan Spiegelman, Cheryl Waters, Oury Monchi, Yves Dauvilliers, Nicolas Dupré, Lior Greenbaum, Sharon Hassin-Baer, Irina Miliukhina, Alla Timofeeva, Anton Emelyanov, Sofya Pchelina, Roy N Alcalay, Ziv Gan-Or
PMID42272075
发布时间2026-06-10
DOI10.1002/mds.70391

实验完整度

中

多队列靶向/全基因组测序及罕见变异负荷分析,包含FDR校正和结构域分析,但无功能验证。

主要模型

PD患者与对照人群队列 McGill大学四个队列(靶向测序) UK Biobank全基因组测序队列 AMP-PD全基因组测序队列

重点核对

罕见变异定义:MAF≤0.01 变异功能分类:非同义、LOF、CADD≥20 统计方法:SKAT-O检验及MetaSKAT荟萃分析 多重检验校正:Benjamini-Hochberg FDR,Pfdr<0.05 协变量调整:年龄、性别、种族

摘要

背景:溶酶体功能障碍是帕金森病(PD)发病机制的核心,其中GBA1是最强的已确立遗传风险因素。众多参与溶酶体鞘脂、糖鞘脂和神经酰胺代谢的其他基因也被认为与PD相关,凸显了对这些通路进行遗传分析的必要性。目的:评估溶酶体基因罕见变异对PD风险的贡献。方法:我们分析了8267例PD个体和68208例对照中36个溶酶体基因的罕见变异(次要等位基因频率≤0.01),包括793例早发性PD(≤50岁)病例。在McGill大学的四个队列(3456例病例和2664例对照)中进行靶向测序,并与来自英国(UK)Biobank(2848例病例,62451例对照)和Accelerating Medicines Partnership-PD队列(1963例病例,3093例对照)的全基因组测序数据合并。使用Sequence Kernel Association Test-Optimal across variant classes(罕见变异、非同义、功能丧失和预测有害变异,combined annotation-dependent depletion score >20)检验关联,随后进行跨队列荟萃分析。对位于蛋白质结构域内的变异进行结构域水平分析。应用错误发现率(FDR)校正。结果:荟萃分析发现ST3GAL3罕见变异与帕金森病显著相关(Pfdr=0.04)。基于结构域的分析显示HEXA的β-乙酰氨基己糖苷酶样结构域内非同义变异富集(P=8.0×10−4),但该信号未通过校正(Pfdr=0.154)。在早发性PD中,基于结构域的分析发现NAGLU(Pfdr=7.3×10−6)和ST3GAL5(Pfdr=0.03)存在显著关联。结论:多个溶酶体通路中的罕见变异,特别是与唾液酸化、神经节苷脂代谢、神经酰胺生物学和溶酶体蛋白水解相关的变异,可能超出GBA1范围促进PD易感性,突出了未来复制和研究的通路。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
评估溶酶体基因罕见变异对帕金森病风险的贡献,并识别除GBA1外的其他溶酶体风险基因。
核心机制
唾液酸化、神经节苷脂代谢、神经酰胺生物学和溶酶体蛋白水解通路中的罕见变异可能促进PD易感性。
主要证据
跨队列荟萃分析显示ST3GAL3罕见变异与PD显著相关(Pfdr=0.04);早发性PD中NAGLU和ST3GAL5结构域关联显著(Pfdr分别为7.3×10−6和0.03);HEXA结构域富集信号未通过校正。
研究意义
研究结果扩展了除GBA1外的溶酶体遗传结构,提示这些通路可能作为未来复制和研究的方向,但需要更大规模和祖先多样性队列及功能研究验证。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

研究人群与队列组建

构建PD病例和对照队列以评估罕见变异负荷。

纳入McGill大学四个队列(靶向测序)、UK Biobank和AMP-PD全基因组测序数据,共8267例PD和68208例对照,包括793例早发性PD。

2

靶向测序与变异检测

对溶酶体基因进行深度测序以识别罕见变异。

设计包含36个溶酶体基因的靶向panel,使用MIPs捕获,在Illumina NovaSeq 6000平台测序(paired-end 100 bp)。

3

质量控制与变异注释

确保变异检测可靠性并注释功能类别。

对变异和样本进行QC,包括call rate、GQ、测序深度过滤,使用ANNOVAR注释,定义罕见变异(MAF≤0.01)并按功能分组。

4

罕见变异负荷关联分析

检验溶酶体基因罕见变异与PD的关联。

在每个队列中使用SKAT-O检验,调整年龄、性别和种族,随后使用MetaSKAT进行跨队列荟萃分析。

5

结构域水平分析

评估蛋白质结构域内变异的富集关联。

将变异映射至UniProt/Pfam定义的功能区域,重复负荷分析。

6

早发性PD亚组分析

评估早发性PD中溶酶体基因罕见变异负荷。

对793例早发性PD进行基因水平和结构域水平负荷分析。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
分子倒置探针----
Illumina NovaSeq 6000测序仪IlluminaNovaSeq 6000

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
研究人群与队列组建
病例和对照的样本量、性别比例、平均年龄;早发性PD的样本量
阅读提示:Methods中Study Population小节及Table 1
靶向测序与变异检测
靶向panel基因列表(36个基因)、MIPs方案、测序平台和读长
阅读提示:Methods中Gene Selection, Sequencing, and Processing小节及Table S1
质量控制与变异注释
各队列特异的QC阈值(call rate、GQ、深度)、HWE过滤、罕见变异MAF定义
阅读提示:Methods中Quality Control和Annotation, Grouping, and Domains小节
罕见变异负荷关联分析
SKAT-O检验的协变量(年龄、性别、种族)、荟萃分析方法、FDR校正方法
阅读提示:Methods中Statistical Analysis小节
结构域水平分析
结构域的定义来源(UniProt/Pfam)、变异分组方式
阅读提示:Methods中Annotation, Grouping, and Domains小节