《Neurology Genetics》：Identification of Nuclear Genetic Loci Linked to Clinical Features of the m.3243A>G Mitochondrial DNA Variant
编辑推荐：
背景与目的：线粒体DNA（mtDNA）疾病表现出显著的临床变异性，而已知因素如变异体异质性、年龄或性别对此解释不足。核遗传修饰因子可能在这种异质性中起重要作用。研究人员的目的是针对常见致病变异体m.3243A>G的2种常见综合征表现，即线粒体脑肌病伴乳酸酸中毒
背景与目的：线粒体DNA（mtDNA）疾病表现出显著的临床变异性，而已知因素如变异体异质性、年龄或性别对此解释不足。核遗传修饰因子可能在这种异质性中起重要作用。研究人员的目的是针对常见致病变异体m.3243A>G的2种常见综合征表现，即线粒体脑肌病伴乳酸酸中毒和卒中样发作（MELAS）及母系遗传性糖尿病和耳聋（MIDD），阐明核遗传参与的性质。
方法：研究人员组建了一个多中心临床确诊的m.3243A>G携带者队列（总计n = 488），在76个家系中鉴定出198名适合进行遗传连锁分析的个体。研究人员调查了MELAS和MIDD的4个临床特征：糖尿病、听力障碍、卒中样发作和脑病。采用基于Haseman-Elston回归的遗传连锁分析来鉴定与这些特征共分离的核基因组区域。使用逻辑回归考虑了m.3243A>G异质性、年龄和性别的影响；通过特征特定的基因滴落模拟确定了经验显著性阈值。在247名个体中使用单变异体（SAIGE）和基于基因的方法（SAIGE-GENE+和MAGMA）进行关联分析，以细化显著连锁区域内的候选位点。
结果：研究人员发现了与脑病显著遗传连锁的区域（7号染色体q22；LOD = 3.72），以及在1、5、6、11和13号染色体上与脑病和卒中样发作提示遗传连锁的区域。未发现糖尿病或听力障碍的连锁。7号染色体区域内的关联分析鉴定出p值最低的变异体rs62500792（SDHAF3和TAC1之间的基因间区域）（3.7 × 10?5），但没有变异体达到比例显著性阈值（5.3 × 10?6）。基于基因的分析强调PLOD3（p = 3.9 × 10?3）和IMMP2L（p = 6.4 × 10?3）为候选基因，因为在互补负荷测试方法中，它们在连锁区域内显示了最强的基因级信号，尽管两者均未达到校正后的显著性阈值。
讨论：修饰m.3243A>G的核遗传结构因临床特征而异。严重的神经特征（脑病和卒中样发作）可能受少数效应量相对较大的核基因影响，而核遗传对糖尿病和听力障碍的贡献似乎更具多基因性。本研究突出了大型且表征良好的患者队列在鉴定修饰位点以及推动理解mtDNA疾病表型变异机制方面的价值。
线粒体DNA（mtDNA）疾病具有极端的临床异质性，即使在携带相同致病变异体的患者群体中也存在显著差异。m.3243A>G是最常见的致病性mtDNA变异体，可导致线粒体脑肌病伴乳酸酸中毒和卒中样发作（MELAS）及母系遗传性糖尿病和耳聋（MIDD）等综合征。已知风险因素如mtDNA异质性水平、年龄和性别仅能解释不到20%的表型变异，这给患者的疾病诊断和预后遗传咨询带来了巨大挑战。由于核基因组编码了约1,150种维持线粒体功能所需的蛋白质，推测核遗传变异在m.3243A>G相关疾病的表型多样性中扮演着关键的修饰角色。为了探究这一未知的核遗传因素，研究人员开展了一项多中心队列研究，对m.3243A>G携带者及其母系亲属进行了非参数遗传连锁分析。研究得出结论，严重的神经学特征受少数效应量较大的核基因影响，而糖尿病和听力损伤则呈现出多基因遗传特征。该研究揭示了mtDNA疾病表型变异的潜在机制，为未来的个性化医疗和靶向治疗提供了理论依据。本研究发表在《Neurology Genetics》期刊上。
在关键技术方法方面，研究人员从5个临床中心招募了488名m.3243A>G携带者（涵盖英国、德国和意大利的队列），经Newcastle Mitochondrial Disease Adult Scale（NMDAS）评估并转化为4项二分类临床特征表型。使用逻辑回归排除年龄、性别和异质性影响后提取残差作为代理表型。随后利用Haseman-Elston回归进行多点遗传连锁分析，通过1,000次基因滴落模拟确定经验显著性阈值。最后，在显著连锁区域内利用SAIGE、SAIGE-GENE+和MAGMA工具进行单变异体和基因负荷关联分析，以细化候选位点。
研究结果如下：
m.3243A>G异质性与年龄不能解释临床变异性
通过在扩大队列中使用二分类特征数据进行分析，证实了m.3243A>G异质性、年龄在脑病、卒中样发作、糖尿病和听力障碍的严重程度变异中所占比例不足20%。在脑病分析中观察到性别效应，但加入性别变量后可解释的表型变异比例仍低于10%。这表明需要寻找其他如核基因组的因素来解释剩余的变异性。
m.3243A>G相关疾病神经学特征的遗传连锁
针对脑病，在第1、5、11和13号染色体上发现了提示性遗传连锁区域，并在第7号染色体上（7q21.13–31.32）鉴定出显著遗传连锁区域。该区域包含200多个基因，其中14个具有潜在的线粒体功能。对于卒中样发作，在第5和第6号染色体上发现了提示性遗传连锁区域。由于脑病与卒中样发作在临床上存在因果相关性，两者在第5号染色体上的重叠区域暗示可能存在影响这两种特征的共同核因子。
m.3243A>G相关疾病非神经学特征表现为多基因性
针对MIDD的临床特征进行分析，证实了糖尿病与听力障碍之间存在相关性，但未发现任何连锁区域。这表明这些非神经学特征的潜在遗传结构可能更具多基因性，与普通人群中的风险模式相似。
提示连锁区域的进一步特征分析
通过基因滴落模拟数据进一步评估发现，脑病分析中识别出的6个超过提示阈值的区域具有高度统计学显著性（仅在千分之一模拟中出现），支持对这些区域进行后续研究。而卒中样发作中识别出的提示区域在此类模拟检验中不具备同等显著性。
关联分析
为进一步研究第7号染色体上与脑病显著连锁的区域，研究人员测试了跨30 Mb区域的54,895个单核苷酸多态性。尽管未达到比例显著性阈值，但发现位于SDHAF3和TAC1之间的基因间变异体rs62500792具有最低的p值。基于基因的分析表明，PLOD3和IMMP2L分别在外显子和全基因负荷测试中信号最强，但均未通过多重检验校正。
讨论部分表明，本研究通过大型多中心家系队列，表征了m.3243A>G相关疾病不同临床特征的核遗传修饰图谱。影响严重神经学特征（脑病和卒中样发作）的核基因具有相对较大的效应量，而影响糖尿病和听力损伤的核基因则呈现多基因模式。在连锁区域中，编码线粒体转录终止因子1的MTERF1被认为是特别值得关注的候选基因，因m.3243A>G变异会降低其结合效率。然而，由于样本量相对较小以及遗传异质性，未能在连锁区域内发现达到全基因组显著性的关联变异。研究强调了临床表现定义、年龄相关特征对错误分类的影响，并指出需要扩大国际队列以覆盖不同遗传背景的个体。本研究为未来识别mtDNA疾病的核修饰因子及进一步了解线粒体疾病的分子机制奠定了基础。
