《European Journal of Neurology》：Distinct UNC13A Haplotype Blocks Define Disease Severity and Survival in Chinese Amyotrophic Lateral Sclerosis
编辑推荐：
背景：UNC13A是欧洲人群肌萎缩侧索硬化症(amyotrophic lateral sclerosis, ALS)的已知遗传修饰基因，但其在中国人群中的作用尚不清楚。研究人员探讨了中国ALS队列中UNC13A变异谱及其对疾病易感性与进展的影响。方法：对153
背景：UNC13A是欧洲人群肌萎缩侧索硬化症(amyotrophic lateral sclerosis, ALS)的已知遗传修饰基因，但其在中国人群中的作用尚不清楚。研究人员探讨了中国ALS队列中UNC13A变异谱及其对疾病易感性与进展的影响。方法：对1533例中国散发性ALS患者和1405例对照进行整合遗传学分析，包括罕见变异负荷检验、全基因组生存分析、单倍型作图及条件分析，并建立并验证临床-遗传联合预后评分。结果：罕见有害UNC13A变异与ALS易感性无关。研究人员鉴定出两个具有不同临床影响的独立单倍型区块(Block)：Block 1(以rs75421007为标签)与基线肌力降低相关(p = 0.030)；Block 2(以rs78549703为标签，为脑特异性剪接数量性状位点splicing QTL即sQTL)是生存异质性的主要驱动因素。欧洲人群报道的变异rs12608932在单标记分析中显示生存关联(p = 0.024)，但条件分析揭示其效应不独立于Block 2。结合临床因素与Block 2单倍型的整合预后评分将患者分层为低、中、高风险组(中位生存期分别为52.6、37.1和32.0个月；p < 0.001)，决策曲线分析证实其临床效用。结论：本研究阐明了中国ALS人群中UNC13A的遗传架构，鉴定出分别影响疾病严重程度与生存的两个独立单倍型区块。含脑sQTL的Block 2单倍型是生存异质性的主要决定因素，可为未来研究的患者分层提供依据。
论文解读——《Distinct UNC13A Haplotype Blocks Define Disease Severity and Survival in Chinese Amyotrophic Lateral Sclerosis》
该研究发表于《European Journal of Neurology》。
一、研究背景与立项依据
肌萎缩侧索硬化症（amyotrophic lateral sclerosis, ALS）是一种致死性神经退行性疾病，以上、下运动神经元进行性退变为特征，最终导致肌无力、瘫痪及呼吸衰竭。虽多为散发性，但遗传因素显著影响疾病易感性及表型异质性。目前已发现包括C9orf72、SOD1、TARDBP和FUS在内的40余个致病基因，而超过95%的ALS病例存在TAR DNA结合蛋白43（TDP-43）胞浆错位聚集及核功能缺失，导致RNA加工异常（如UNC13A前体mRNA中隐秘外显子cryptic exon的异常剪接）。UNC13A编码突触前膜释放必需蛋白，是TDP-43调控的关键靶基因；欧洲人群全基因组关联研究（GWAS）已确立rs12608932等多态性与ALS易感性及生存相关。然而，UNC13A在中国等非欧人群中的遗传架构及表型影响尚未明确，已报道的欧洲风险变异在不同族裔中效应不一致，且缺乏对中国人群UNC13A罕见与常见变异综合影响的分析。因此，研究人员开展此项研究以界定中国ALS人群UNC13A群体特异性遗传结构，识别调控疾病严重程度与生存的关键变异。
二、主要关键技术方法
研究人员纳入中南大学湘雅医院确诊的1533例散发性ALS患者（符合Gold Coast标准）及同地区1405例健康对照。对患者行全基因组测序（WGS, GRCh38/hg38），对WGS覆盖不全的rs12608932行Kompetitive Allele-Specific PCR（KASP）补充分型。罕见有害变异定义为MAF<0.01且经≥2种软件预测为有害（REVEL>0.7等），并经Sanger测序验证。常见变异进行Hardy–Weinberg平衡检验，采用PLINK和R进行连锁不平衡（LD）及单倍型块（haplotype block）结构分析。生存定义为症状发作至死亡或气管切开时间，采用Cox比例风险回归（调整性别、发病年龄、延髓起病、基线MRC及ALS FRS-R评分）和Kaplan–Meier法（log-rank检验）评估变异-生存关联，显著性阈值p<5×10??，多重检验用Benjamini-Hochberg校正。通过条件分析判断变异间独立性，以4个标签SNP定义Block 2风险单倍型（AGTA），建立临床评分、遗传评分及二者整合评分，并用时间依赖ROC、校准图及决策曲线分析（DCA）验证模型效能。排除携带SOD1/TARDBP/FUS/C9orf72致病突变患者进行敏感性分析。
三、研究结果
3.1 研究人群与临床特征（Study Population and Clinical Characteristics）
最终纳入1533例患者（男62.5%，女37.5%），中位发病年龄59.0岁，延髓起病占24.3%，首发症状起中位生存期38.3个月，基线ALS FRS-R中位40.0分，MRC总分中位114.0分。人口学与临床数据汇总于表1。
3.2 罕见UNC13A编码变异负荷与ALS风险或严重程度无关联（Burden of Rare UNC13A Coding Variants Shows no Association With ALS Risk or Severity）
罕见错义变异总负荷（OR=1.39, p=0.080）及预测为有害的罕见变异子集（OR=1.96, p=0.112）在患者中均未显著富集；未检出功能缺失（LoF）变异。按频率与致病性分层后，各类罕见变异与发病年龄、起病部位及生存均无显著关联。表明UNC13A罕见编码变异不是本散发性ALS队列的主要致病或表型异质性机制。
3.3 已知变异验证与新型修饰位点的全基因组发现（Validation of Known Variants and Genome-Wide Discovery of Novel Modifiers）
欧洲风险变异rs12608932（基因型和p=0.024，隐性p=0.042）及隐秘外显子内rs12973192（基因型和p=0.029，隐性p=0.016）均与生存显著相关。条件分析显示二者互相条件后关联均不显著（rs12608932条件p=0.19；rs12973192条件p=0.12），提示信号源于同一功能性遗传单元。全基因组生存分析发现两个新常见变异：rs78549703（β=-0.190, p=0.003）和rs75421007（β=-0.208, p=0.004），均为UNC13A脑特异性剪接数量性状位点（sQTL）。rs78549703 AA基因型中位生存期32.1月短于GG型40.5月（加性p=0.0018，隐性p=0.0014），多变量Cox示AA型死亡风险增加67%（HR=1.67, 95%CI:1.20–2.34, p=0.003）；rs75421007 AA型HR=1.69（p=0.034）。
3.4 具不同表型影响的两个独立单倍型块的界定（Delineation of Two Independent Haplotype Blocks With Distinct Phenotypic Impacts）
LD分析明确两个独立单倍型块（块间平均r2<0.1）：Block 1以新变异rs75421007为标签；Block 2以顶级sQTL rs78549703为标签并包含rs12973192，rs12608932位于块外但与Block 2呈中度LD（r2=0.34）。多变量Cox确认两块的领衔SNP独立关联生存——rs78549703（Block 2）隐性模型AA vs GA+GG: HR=1.61, p=0.005；rs75421007（Block 1）AA vs GG: HR=1.69, p=0.034。表型特异性分析显示功能分化：Block 1（rs75421007）显性模型与基线肌力降低相关（回归系数=-2.056, p=0.030），提示影响初始疾病严重程度；Block 2变异与任何基线表型（发病年龄、延髓起病、ALS FRS-R、MRC）无关，仅与发病后进展速率关联。
3.5 条件分析确立Block 2为生存异质性的主要遗传驱动因素（Conditional Analysis Establishes Block 2 as the Primary Genetic Driver of Survival Heterogeneity）
固定rs75421007(GG)背景，rs78549703(Block 2)生存关联仍显著（加性HR=1.61, p=0.011；隐性HR=1.59, p=0.013）；反之固定rs78549703(GG)背景，rs75421007整体加性模型不显著（p=0.102）。Kaplan–Meier图显示rs75421007在rs78549703 GG携带者中效应微弱，而rs78549703在rs75421007 GG中维持强梯度效应。证实Block 2是更原发且独立的生存异质性驱动因子；条件于Block 2后rs12608932生存关联完全衰减（条件p=0.28），确认欧洲变异信号由Block 2介导。
3.6 单倍型风险评分揭示生存上的累积遗传负荷（A Haplotype-Based Risk Score Reveals a Cumulative Genetic Burden on Survival）
用rs78549703-rs12973192-rs12975883-rs36108613定义Block 2 AGTA风险单倍型，将患者分为4组：Group 1（AGTA纯合，n=36）、Group 2（AGTA杂合）、Group 3（无AGTA但≥4风险等位基因）、Group 4（参照）。Log-rank检验p=0.0069，呈基因剂量效应；Group 1较Group 4死亡风险升高89%（HR=1.89, 95%CI:1.27–2.80, p=0.002），调整临床因素后HR=1.69（p=0.011）；每升一级风险类别调整后死亡风险增9.3%（HR=1.093, p=0.017）。
3.7 排除已知ALS致病基因突变携带者的敏感性分析（Sensitivity Analysis Excluding Carriers of Known ALS-Associated Gene Mutations）
剔除112例SOD1/TARDBP/FUS/C9orf72致病变异携带者后，rs78549703(AA vs GG: HR=1.80, p<0.001)、rs75421007(AA vs GG: HR=1.69, p=0.040)及Block 2风险单倍型组（Group 1 vs 4: HR=1.81, p=0.003）关联均保持稳健，与全队列一致，证实UNC13A的独立预后价值。
3.8 整合临床-遗传预后模型的开发与验证（Development and Validation of an Integrated Clinical-Genetic Prognostic Model）
建立三种评分：临床评分（发病年龄>55岁、延髓起病、MRC<110、ALS FRS-R<38、男性各1分，范围0–5）、遗传评分（Block 2单倍型四分类）及整合评分（遗传评分数值化加临床评分，范围0–8）。整合评分有效分层为低（0–2分，中位生存52.6月）、中（3–4分，37.1月）、高（5–8分，32.0月）风险组（p<0.001）；决策曲线分析示整合模型在宽泛阈值范围内净获益高于单纯临床模型，时间依赖ROC及校准分析验证其判别能力与泛化性。
四、讨论总结与结论翻译
讨论指出本研究首次系统描绘中国ALS人群UNC13A罕见与常见变异谱，发现罕见编码变异不关联疾病或表型；确认rs12608932生存关联但非独立，被Block 2解释；鉴定出两独立单倍型块——Block 1影响基线疾病严重程度（肌力），Block 2（含脑sQTL，与TDP-43调控的隐秘外显子区域相连）主导生存/进展速率，且具单倍型剂量效应；整合预后模型可识别初诊表现轻但具快速进展遗传负荷的患者，有助于临床试验富集设计及未来TDP-43靶向治疗（如反义寡核苷酸）的患者分层；研究受限于风险单倍型纯合子样本量小、回顾性设计及尚需外部验证与功能实验证实。
结论（Conclusion）翻译：
综上，本研究阐明了中国ALS人群UNC13A的群体特异性遗传架构。研究人员确定UNC13A内两个独立单倍型区块（分别以rs75421007〈Block 1〉和rs78549703〈Block 2〉为标签）的常见变异差异性地影响疾病严重程度与进展。研究人员阐明rs12608932的生存关联并非独立效应，而可能由Block 2介导。整合预后模型为风险分层提供了实用工具，Block 2风险单倍型的累积效应进一步支持UNC13A作为极具潜力的治疗靶点。未来研究应优先对鉴定出的变异进行功能表征，并在不同人群中评估其预后效用。
