《Forensic Science International: Genetics》：The EDNAP collaborative study on mtDNA heteroplasmy detection using Sanger and Massively Parallel Sequencing technologies
编辑推荐：
这项EDNAP（欧洲DNA分型工作组）协作研究评估了法医实验室使用Sanger、Ion Torrent和Illumina测序平台对线粒体DNA（mtDNA）点异质性（PHP）和长度异质性（LHP）的检测和报告情况。研究人员集中制备了标准化DNA提取物并分发给参
这项EDNAP（欧洲DNA分型工作组）协作研究评估了法医实验室使用Sanger、Ion Torrent和Illumina测序平台对线粒体DNA（mtDNA）点异质性（PHP）和长度异质性（LHP）的检测和报告情况。研究人员集中制备了标准化DNA提取物并分发给参与实验室，以评估实验室间一致性及平台特异性（平台=仪器+软件）效应。原始序列数据显示各平台整体高度一致，观察到的线粒体型（mitotype）差异源于转录错误。PHP检测在不同实验室和技术间总体一致。MPS（大规模并行测序）数据在平台间提供了相对稳定的混合比例。本研究中，若低水平PHP低于实验室特定的次要等位基因频率（MAF）阈值，则通常不予报告，尽管原始数据审核确认了其存在。相比之下，LHP解读显示显著更大的变异。基于Sanger的LHP模式可通过根据当前ISFG（国际法医遗传学会）指南进行审核而协调一致，而基于MPS的LHP结果在测序技术之间和内部均存在差异。对齐工作流程的标准化在一定程度上提高了一致性，但重要差异仍然存在，尤其是在Ion Torrent HVS-I（高变区I）poly(C)-tract数据中。这些差异归因于测序技术、分析软件和人工解读。
**论文解读文章**
**研究背景**
线粒体DNA（mtDNA）异质性指个体内共存多种mtDNA类型，包括序列差异的点异质性（PHP）和大小差异的长度异质性（LHP）。自1983年首次发现以来，mtDNA异质性在法医学中具有重要应用，如罗曼诺夫家族遗骸鉴定和美国田纳西州诉Paul Ware案。早期研究基于Sanger测序，但大规模并行测序（MPS）技术的引入显著降低了PHP检测的次要等位基因频率（MAF）阈值，从10%-20%降至1%-2%。然而，MPS数据可能引入假阳性，如死后损伤、核mtDNA元件（NUMTs）和编码错误，故推荐保守阈值（约10% MAF）。此外，不同测序平台（如Illumina和Ion Torrent）在poly(C)-tract的LHP检测中存在差异，影响法医案例解读。为此，欧洲DNA分型工作组（EDNAP）开展系统性多中心研究，使用标准化对照样本评估不同平台（含仪器和软件）的PHP和LHP检测一致性。
**研究内容与结论**
研究人员制备了5个DNA提取物（来自两个商业细胞系007和9947A、两份颊拭子和一份血液样本），分发给25个法医实验室，使用Sanger、Ion Torrent和Illumina平台进行测序。24个实验室返回结果。主要结论：PHP检测在平台间高度一致，原始数据支持低MAF的PHP存在；LHP解读差异显著，受测序技术、分析软件和人工解读影响，尤其Ion Torrent HVS-I poly(C)-tract数据。标准化对齐流程部分改善一致性，但需更新法医指南。论文发表在《Forensic Science International: Genetics》。
**关键技术方法**
研究人员使用以下主要方法：1）样本来源：两个商业细胞系（007和9947A）的DNA提取物、两份颊拭子（CS-3、CS-5）和一份血液样本（CS-4），均来自志愿者（因斯布鲁克医科大学伦理委员会批准）。2）DNA定量：采用SDquant方法（两种mtDNA靶标69和143 bp，一种核DNA靶标70 bp）在ABI Prism 7500仪器上进行多重PCR定量。3）测序平台：Sanger（3130/3500系统）、Ion Torrent（Ion S5/PGM）和Illumina（MiSeq FGx/MiSeq），各实验室使用自身方法及软件。4）数据分析：组织实验室（Lab 0）收集报告线粒体型和原始数据，由两名经验分析师根据ISFG指南独立审核，并使用基于系统发育对齐原则的SAM2软件进行BAM文件重分析。
**研究结果**
**3.1 DNA定量**
通过SDquant定量，5个DNA提取物在Lab 0中核DNA浓度约0.1 ng/μL，mtDNA拷贝数4298-26155 copies/μL。6个参与实验室的定量结果差异较大（10%-314%），表明运输条件导致可测量DNA量显著损失（26%-77%），但未见额外降解。
**3.2 选定样本的线粒体型**
Lab 0通过Sanger和全线粒体基因组测序确定了5个样本的线粒体型，包括PHP位点（如CS-1的2413Y、CS-2的7861Y、CS-4的12483Y和16093Y、CS-5的234R）和LHP位点（如HVS-I、HVS-II、HVS-III的poly(C)-tract）。
**3.3 报告线粒体型与原始数据评估**
24个实验室返回结果：8个Sanger、13个Ion Torrent、8个Illumina（部分重叠）。非异质性变异的差异极少（共14个转录错误），且原始数据审核确认其存在。
**3.4 点异质性评估**
**3.4.1 Sanger检测**：所有8个实验室报告CS-5的234R，6个报告CS-4的16093Y，其余2个报告主要分子但原始数据确认PHP存在。
**3.4.2 Ion Torrent检测**：所有13个实验室报告CS-4的16093Y（中位AF=58%）和CS-5的234R（中位AF=50.5%），11个报告CS-4的12483Y（中位AF=33.3%）。低MAF的2413Y（~9.8%）和7861Y（~8.0%）仅少数报告。重分析（MAF>5%）后检测一致性提高。
**3.4.3 Illumina检测**：6个实验室报告CS-4的16093Y（中位AF=58.1%），7个报告CS-5的234R（中位AF=48.3%），2个报告CS-4的12483Y（中位AF=28.2%）。重分析中，234R的比值差异因引物覆盖位点234所致。
**3.4.4 MPS间比较**：两种平台在报告和重分析数据中PHP比值相似（IQR较小），但Illumina的234R数据因引物问题出现离群值。
**3.4.5 额外混合位点**：重分析发现低MAF混合，部分为平台特异性伪影（如Ion Torrent在poly(C)-tract附近），其他为背景信号。
**3.5 长度异质性评估**
**3.5.1 HVS-I（约16193位点）**：Sanger数据经审核后一致（CS-3: 16193.1C；CS-5: 16193.2C）。Ion Torrent数据报告与审核结果差异大，包括缺失和插入。Illumina数据审核后以单次C插入为主。
**3.5.2 HVS-II（约309位点）**：Sanger数据一致（CS-2: 309.2C；CS-3: 309.1C）。Ion Torrent报告和审核结果在309.1C和309.2C间波动，Illumina审核后一致（CS-2: 309.2C；CS-3: 309.1C）。Ion Torrent在无LHP样本中误报309del或309c。
**3.5.3 HVS-III（约573位点）**：Sanger数据审核后一致（CS-1: 573.4C；CS-3: 573.3C；CS-5: 573.2C）。Ion Torrent和Illumina审核后结果接近，但仍有差异（如CS-1中Ion Torrent出现573.5C）。
**讨论与结论**
讨论指出，Sanger数据经ISFG指南审核后LHP模式一致，而MPS数据差异源于测序技术、分析软件和人工解读三方面。标准化软件（SAM2）部分改善HVS-II和HVS-III的一致性，但Ion Torrent HVS-I poly(C)-tract仍存在变异。建议更新法医指南，重点关注LHP的验证和报告标准。
结论翻译：这项EDNAP协作研究表明，mtDNA点异质性可在当前法医测序平台（包括Sanger、Ion Torrent和Illumina）上以高一致性检测和报告。相比之下，长度异质性仍显著多变，受测序技术、分析软件和用户解读影响。虽然标准化审核改善了部分MPS数据集的一致性，但某些poly(C)-tract区域仍存在重要差异。研究结果指出，当前主要为Sanger测序开发的mtDNA解读方法已不足以用于常规MPS异质性分析。因此，需要更新法医指南，特别是针对长度异质性的解读和报告，并更强调化验特异性验证、标准化生物信息学工作流程和协调的报告标准。
