PacBio HiFi长读长测序实现帕金森病线粒体DNA异质性高精度检测

《BMC Genomics》：Long-read sequencing enables high-accuracy mitochondrial heteroplasmy detection in Parkinson’s disease

背景：低频异质性线粒体DNA（mitochondrial DNA, mtDNA）变异与衰老及神经系统疾病（包括帕金森病（Parkinson's disease, PD））相关。采用PacBio HiFi长读长进行靶向深度mtDNA测序，有望在全线粒体基因组范围内以高精度解析异质性。 方法：为验证Vega PacBio测序在检测mtDNA异质性中的效能，研究人员分析了四种由两种mtDNA单倍型组成的预定义混合物。研究人员生成了覆盖整个线粒体基因组的单一长片段PCR扩增子，并以预定义比例（次要混合单倍型组分：5%、2%、1%和0.1%）混合这些扩增子。使用Mutserve2进行变异检测，并通过计算预期与检测变异之间比较所得的F1分数评估准确性。将全长mtDNA PacBio测序应用于探索来自五名LRRK2 p.Gly2019Ser变异携带者（n=3名PD患病者，n=2名未患病携带者）的成纤维细胞传代中的异质性。使用线性混合模型纵向评估mtDNA异质性水平和变异负荷的变化，作为原理验证。 结果：单扩增子方法实现了全长单倍型解析，避免了与重叠PCR策略相关的扩增偏倚。预定义混合物在5%至1%异质性水平间的F1分数为1.0，在0.1%时仍保持较高水平（0.91）。在针对0.1%混合物特别降低的检测阈值下，观察到n=10/62个与Illumina参考不一致的额外变异，但该混合物中的敏感性仍保持非常高的1.00。检测到的次要变异与预期异质性水平高度吻合，平均变异水平分别为0.057（5%）、0.022（2%）、0.011（1%）和0.001（0.1%）。在十二次成纤维细胞传代中，研究人员观察到较少的mtDNA异质性变异。在这项小型原理验证实验中，患病者（n=3）的异质性变异负荷高于未患病者（n=2）LRRK2变异携带者，且与年龄增长相关。值得注意的是，研究人员观察到异质性变异呈现不同模式，部分变异在传代过程中异质性水平升高，部分则降低。 结论：PacBio HiFi测序结合单扩增子策略，可实现准确的全长mtDNA异质性检测和纵向分析，为研究疾病中线粒体变异和动态变化提供了有价值的工具。

论文解读：PacBio HiFi长读长测序在帕金森病线粒体DNA异质性检测中的应用

一、研究背景与问题

线粒体DNA（mitochondrial DNA, mtDNA）在生殖系和体细胞中持续发生突变与选择，形成mtDNA异质性，即同一细胞、组织或生物体内共存多种不同的线粒体基因型。mtDNA的损伤修复效率低于核DNA，其突变率约为核DNA的10倍，因此并非所有mtDNA拷贝共享相同变异，仅存在于部分拷贝中的变异被称为异质性变异。已有研究表明，mtDNA异质性不仅与衰老密切相关，还涉及多种疾病，包括帕金森病（Parkinson's disease, PD）等神经退行性疾病。因此，准确表征和定量mtDNA异质性变异对于理解线粒体生物学及mtDNA异质性在疾病中的作用至关重要。

目前，下一代短读长测序技术已实现高深度mtDNA分析，被视为异质性检测的金标准。然而，这些方法依赖约50-300bp的读长，可能导致覆盖不均、比对偏倚、易受核线粒体插入片段（NUMTs）污染影响，并限制全长mtDNA单倍型的直接解析。近年研究已展示长读长测序在mtDNA分析中的价值，包括结构变异表征、从头组装、单倍型重建及靶向多扩增子mtDNA测序。但既往研究多聚焦于基因组组装、缺失检测或较高水平的异质性。研究人员此前利用牛津纳米孔测序（Oxford Nanopore Technology, ONT）的基准研究表明，变异水平高于5%时方可实现准确检测。尽管ONT的化学和碱基识别模型已有显著改进，但针对低频mtDNA异质性检测的系统性基准评估仍然有限。PacBio HiFi测序提供高精度长读长，是检测低频变异的理想工具，然而其在低变异水平（<5%）下的系统基准评估仍较为缺乏。此外，全长mtDNA分子的深度测序常因使用多个重叠PCR扩增子而受到影响。

二、研究内容与意义

本研究旨在探索PacBio HiFi测序结合单一全长长片段PCR扩增方法检测mtDNA异质性的性能。通过预定义单倍型混合物，研究人员评估了低变异水平下的检测准确性，并将该工作流程应用于LRRK2 p.Gly2019Ser变异携带者来源的成纤维细胞传代中异质性的纵向动态分析。该研究为长读长测序在mtDNA异质性检测中的应用提供了系统验证，并首次展示了PacBio HiFi测序结合单扩增子策略在全长mtDNA异质性检测和纵向分析中的高精度性能，为研究线粒体变异及其在疾病中的作用提供了有力工具。

三、主要技术方法

研究人员采用两阶段研究设计。第一阶段为混合物验证实验，使用来自欧洲血统个体的血液DNA，经长片段PCR生成覆盖完整线粒体基因组的单一扩增子，将两种已知单倍型（主要组分D4e1'3与次要组分J1c2）的PCR产物按5%、2%、1%和0.1%的预设比例混合。变异检测使用Mutserve2工具，以Illumina NextSeq测序数据作为金标准，通过F₁分数评估准确性，并利用Haplocheck检测样本污染、Himito评估NUMTs污染。第二阶段为生物学应用，纳入五名LRRK2 p.Gly2019Ser变异携带者（3名PD患病者，2名未患病携带者）来源的成纤维细胞系，在传代1、4、8和12时提取DNA并进行PacBio Vega平台测序。测序数据经Filtlong过滤、Minimap2比对至rCRS参考基因组，使用Mutserve2进行变异检测，并采用线性混合效应模型评估异质性变异负荷随传代的变化。

四、研究结果

PacBio HiFi测序生成高质量全长mtDNA读长，超过95%的读长达到Q30碱基质量。过滤后平均覆盖深度为83,180X（范围64,035X-107,279X），覆盖高度均匀。单扩增子方法实现了全长单倍型解析，避免了重叠PCR策略相关的扩增偏倚。

PacBio工作流程在混合物模型中实现高精度mtDNA变异检测：在5%至1%的混合物水平，预期与检测变异完全一致，F₁分数为1.0，且未检测到与Illumina参考不一致的额外变异。在0.1%混合物中应用降低的检测阈值后，所有预期次要单倍型变异均被恢复，但检测到10个额外变异，F₁分数为0.91。检测到的次要变异水平与预期高度吻合，平均变异水平分别为0.057（5%）、0.022（2%）、0.011（1%）和0.001（0.1%）。

通过逐步降采样评估不同测序深度下的检测性能，结果显示在大于1500X的测序深度下，混合物比例≥1%的F₁分数范围为0.98-1.00，敏感性为1.00。5%混合物在所有测试深度（包括约300X）均保持F₁分数和敏感性为1.00。

未检测到系统性测序错误且测序噪声极低：额外变异仅在0.1%混合物中出现，且非特定位点的系统性错误，部分位于低复杂度区域如短同聚物。常见变异在约64,000X深度下变异水平保持为1，降采样至约300X后仍保持高度一致（≥0.997）。

读长水平单倍型分配验证了全长mtDNA分相能力：在所有预定义混合物中，次要单倍型读长的观察比例与预期混合比例及Mutserve2报告的变异水平高度吻合，信息性变异在单个测序分子上持续共现，证实了全长线粒体单倍型的重建。

原理验证实验中，研究人员纵向评估了五名LRRK2 p.Gly2019Ser变异携带者来源成纤维细胞中的线粒体异质性。在12次传代中观察到异质性变异数量减少，线性混合模型确认随时间推移异质性变异减少（传代12对传代1：β=-3.2，SE=1.3，p=0.026）。探索性分析显示，患病者来源成纤维细胞的异质性变异负荷高于未患病携带者（β=5.0，SE=1.0，p=1.0×10??），且与采样年龄较大相关（β=1.5，SE=0.3，p=0.001）。变异特异性分析揭示了异质性动态的异质性模式，部分变异水平逐渐降低，部分逐渐升高。例如，在患病者来源细胞系L-15,695中，chrM:7642A（MT-CO2，同义变异）和chrM:14118G（MT-ND5，同义变异）水平随传代下降，而chrM:4818A（MT-ND2，p.E117K）水平逐渐升高。利用Mitorsaw和Sniffles2未检测到等位基因频率>1%的结构变异。

五、讨论与结论

讨论部分指出，PacBio HiFi工作流程在1%变异水平下实现100%的一致性，且在高深度测序下可支持大规模多重分析。0.1%混合物中检测到的额外变异可能反映残留测序和比对噪声，而非特定于0.1%混合物的变异，不应被解读为建立0.1%的一般分析检测限。该研究首次将PacBio HiFi测序与单倍型混合物模型整合，证明了其在1%-5%异质性水平下的高精度检测能力，优于此前ONT测序在相同基准框架下>5%的检测限。研究局限性包括长片段PCR可能引入引物脱落风险、无法评估PCR采样效应和聚合酶错误、PCR富集去除天然DNA甲基化、以及结构变异检测性能未系统评估。成纤维细胞传代中异质性变异的动态变化提示选择压力或克隆扩增过程塑造了线粒体异质性，与既往单细胞研究结论一致。

结论：PacBio HiFi测序结合单扩增子策略可实现准确的全长mtDNA异质性检测和纵向分析。基准评估证明异质性变异可高精度检测至1%，0.1%检测保持高敏感性但受参考不一致变异影响。将该工作流程应用于LRRK2 p.Gly2019Ser变异携带者成纤维细胞模型，揭示了随时间变化的动态、变异特异性异质性变化，以及较高异质性变异负荷与LRRK2相关PD之间的关联。该研究为研究线粒体异质性及其在衰老和疾病中的作用提供了有力工具。