《Nucleic Acids Research》：MitImpact 4: a genome-wide resource for the interpretation of coding and non-coding mitochondrial DNA variants
编辑推荐：
线粒体DNA变异是遗传性疾病的常见病因，然而其计算解读几乎仅集中于13个蛋白编码基因。tRNA与rRNA基因中的变异承载了已知致病等位基因中不成比例的份额，却长期处于诊断盲区，其原因在于此类变异的效应取决于RNA结构而非蛋白质序列，且专用的预测工具极少。Mi
线粒体DNA变异是遗传性疾病的常见病因，然而其计算解读几乎仅集中于13个蛋白编码基因。tRNA与rRNA基因中的变异承载了已知致病等位基因中不成比例的份额，却长期处于诊断盲区，其原因在于此类变异的效应取决于RNA结构而非蛋白质序列，且专用的预测工具极少。MitImpact 4消除了这一不对称性。该数据库现已覆盖线粒体基因组中可能发生的全部49 704个单核苷酸变异，并报告其功能后果、系统发育保守性、局部突变约束以及在gnomAD、HelixMTdb与ToMMo中的群体频率。错义变异由APOGEE 2以及新增的预测工具进行注释，包括AlphaMissense、ESM1b、VARITY与LIST-S2；而tRNA与rRNA变异则接受RNA专用模型nAPOGEE以及MitoTIP与PON-mt-tRNA的预测。对两类变异，评分均伴随致病性后验概率，可用于线粒体DNA的ACMG/AMP框架。来自ClinVar与MITOMAP的临床注释，包括ClinGen专家组的分类，被整合于重新设计的界面中，该界面同时可对用户提交的VCF文件进行注释。MitImpact 4可在https://mitimpact.mcb2lab.org免费获取。
# 四、论文主体解读
**研究背景与存在问题**
线粒体是生成细胞大部分三磷酸腺苷（ATP）的膜结合细胞器，通过氧化磷酸化（OXPHOS）发挥核心能量代谢作用。线粒体病由核DNA（nDNA）或线粒体DNA（mtDNA）序列变异所致，其中mtDNA变异在成人（>16岁）确诊的线粒体病病例中占多数，尽管mtDNA仅编码线粒体蛋白质组1136个蛋白中的13个。这一比例失衡部分归因于mtDNA相对nDNA具有更高的突变率，因其修复机制有限。截至统计时，MITOMAP共报告1091个与疾病相关的mtDNA变异，其中仅107个被确认为致病；这些变异在基因组上分布不均：mt-tRNA基因仅占mtDNA约10%，却承载了当前评估的致病与可能致病变异中的大部分（45%，107个中48个，依据遵循ACMG/AMP指南的ClinGen专家组）。此前的计算解读几乎仅聚焦13个蛋白编码基因，tRNA与rRNA变异的效应依赖RNA结构而非蛋白序列，专用预测工具稀缺，形成诊断盲区。MitImpact自2015年首次发布以来为mtDNA全部可能的非同义变异提供功能注释，本研究将其扩展至覆盖全部可能的线粒体单核苷酸变异（SNV），旨在建立编码与非编码变异的全景功能注释框架，支撑临床解读与科研应用。
**主要关键技术方法**
研究以修订版剑桥参考序列（rCRS；NC_012920.1）为参考基因组，对每个SNV按HGVS基因组命名法赋值，并以SnpEff注释功能后果。错义变异致病性预测主要来自AlphaMissense、ESM1b、VARITY、LIST-S2及线粒体专用预测器APOGEE 2；同义与无义变异仅注释保守性、约束与群体频率；tRNA变异由MitoTIP、PON-mt-tRNA与nAPOGEE注释，rRNA变异由nAPOGEE预测。系统发育保守性以470路多重比对计算的PhyloP与PhastCons评分评估，全基因组约束采用线粒体局部约束（MLC）评分。群体频率来自gnomAD v3.1、ToMMo 61KJPN与HelixMTdb；临床注释取自ClinVar与MITOMAP（含ClinGen专家组分类），更新至2026年8月。数据库以Python（pandas、NumPy、openpyxl）构建，网页端由Flask、gunicorn、nginx与MongoDB支撑，前端采用Bootstrap与NGL.js结构渲染。
**研究结果**
*MitImpact 4*：本版本为16 569个mtDNA位点上全部49 704个可能SNV提供功能注释（不含rCRS占位位点m.3107处的三个替换）。对编码与非编码基因重叠区（58个核苷酸、174个SNV），按用户所选目标基因提供不同功能注释。错义变异为最大类别，其注释显著扩展，纳入AlphaMissense与APOGEE 2等评分，后者提供与ACMG/AMP一致的致病性分类。此部分通过数据库内容整合与统计得出全基因组覆盖的结论。
*跨变异类别的注释覆盖*：所有49 704个SNV均具备HGVS命名、功能后果、系统发育保守性、局部约束与群体频率等核心注释；错义变异另有APOGEE 2及通用预测器，tRNA另有nAPOGEE、MitoTIP与PON-mt-tRNA，rRNA另有nAPOGEE。同义、无义及控制区与基因间区3633个SNV目前仅有核心注释。因58个位点重叠计数，分类计数合计49 878而非49 704。该结论由各类别变异数统计得出。
*从蛋白编码向非编码变异扩展致病性注释*：APOGEE 2与nAPOGEE均提供连续评分、后验致病概率及ACMG/AMP对齐分类三重输出。以每位点氨基酸替换的平均APOGEE 2评分构建矩阵，显示替换类型间预测致病性存在显著异质性，涉及甘氨酸与脯氨酸的替换评分偏高。对tRNA与rRNA SNV，tRNA在各替换类型下致病概率持续高于rRNA，转换与颠换间均值相近，提示替换类型本身不足以预测功能影响。在4524个tRNA与7536个rRNA SNV中，ClinVar列为意义未明变异（VUS）或解释冲突者285个，MITOMAP列为“已报告”或“报告冲突”者384个，去重后541个；nAPOGEE对其中80个赋值后验概率>0.99、对45个赋值<0.001，为23.1%的未决变异提供强方向性计算证据，但不构成临床重分类。MLC评分在三类基因间差异显著（Kruskal–Wallis，P = 4.29 × 10
-303），rRNA约束最高，其次tRNA与蛋白编码基因。该部分通过评分矩阵、替换类别比较与约束分布统计得出结论。
*应用：线粒体变异的解读*：以m.5645G > A（MT-TA）为例，该变异在群体数据库缺失、MITOMAP未报告、ClinVar列为与MELAS综合征相关的VUS；其位置保守性可变、MLC为0.58，nAPOGEE预测为可能致病，近期有研究报道其在神经肌肉疾病患者中致病。另一例m.3093C > G（MT-RNR2）在MITOMAP为与MELAS相关的“已报告”变异，群体数据库缺失、位点保守、MLC = 0.74，nAPOGEE评分0.88对应后验概率>0.99；近期一项亚洲心肌病队列研究依ACMG/AMP将其重分类为可能致病。两例说明注释可在实践中组合解读，而非对nAPOGEE性能的独立验证。
*网页界面*：界面支持基因组（m.3307A > G）、蛋白层面（MT-ND1 p.M1L）、基因特异位置查询及基因与区域浏览，注释按可折叠区块组织，并提供字段说明图例；对58个重叠位点按所选基因显示注释；提供VCF注释服务（≤32 MB，比对至NC_012920.1），返回VCF或制表符分隔格式的核心注释；支持批量下载。
**讨论与结论**
线粒体序列变异常为线粒体病病因，此类疾病表型与发病年龄高度可变；高通量测序使新SNV发现速率大增，但数据指数增长使功能解读日益复杂。MitImpact已成为错义变异注释的成熟资源，MitImpact 4将范围扩展至全部线粒体SNV。错义变异因专用预测器而注释最充分，APOGEE 2为当前性能最佳的线粒体专用荟萃预测器；RNA变异工具长期限于tRNA，nAPOGEE填补了rRNA预测空白。作者亦指出局限：当前仅覆盖SNV，不含插入、缺失与大规模重排；注释未考虑异质性水平、组织分布与单倍群背景；同义变异与控制区变异缺乏专用致病性预测器，主要依赖保守性、约束与群体频率。未来版本将扩展至其他变异类型并纳入单倍群信息。研究结论为：通过将上述功能注释整合于单一、持续更新的框架，MitImpact确立了其作为全面资源的地位，可对任一编码或非编码线粒体变异进行潜在功能与临床效应的深入评估，支持临床解读与科研应用；该论文发表于《Nucleic Acids Research》。
