《Genes》：Characterization of Large-Scale Mitochondrial DNA Rearrangements Using Complementary Short- and Long-Read Sequencing Approaches
编辑推荐：
背景：大规模线粒体DNA（mitochondrial DNA, mtDNA）重排（包括缺失或重复）是原发性线粒体疾病的重要病因之一。在临床实践中，基于短读长测序（short-read sequencing, SRS）的检测方法通常通过识别区域性的mtDNA拷贝
背景：大规模线粒体DNA（mitochondrial DNA, mtDNA）重排（包括缺失或重复）是原发性线粒体疾病的重要病因之一。在临床实践中，基于短读长测序（short-read sequencing, SRS）的检测方法通常通过识别区域性的mtDNA拷贝数变化来发现这些重排，但一般无法直接区分缺失与重复。因此，mtDNA重复可能被错误归类为缺失，从而导致不准确的分子诊断和基因型-表型关联判断。本研究旨在评估长读长测序（long-read sequencing, LRS）对于明确表征大规模mtDNA重排的效用。方法：两名此前经SRS临床检测报告携带单一大规模mtDNA缺失的个体，因其临床表现与典型mtDNA缺失综合征不符而接受重新评估。研究人员对已有的SRS数据进行了再分析，并分别采用Pacific Biosciences（PacBio）和Oxford Nanopore Technologies两种长读长测序平台独立表征了可疑的重排。结果：SRS数据的再分析表明，所报告的缺失更符合mtDNA重复的特征。这一发现在其中一例病例中得到了PacBio和Oxford Nanopore两种LRS方法的证实。将该病例中的“缺失”重新分类为“重复”，显著改善了基因型-表型的解释，并更好地解释了该患者表现出的非典型临床特征。然而，在第二例病例中，由于mtDNA结构变异的低异质性水平，LRS未能捕获足够多的长读长以提供重复的直接证据，因此该病例仍无定论。结论：SRS与LRS在大规模mtDNA重排的评估中提供了互补的信息。SRS能够灵敏地识别mtDNA大片段缺失，但无法区分mtDNA缺失与重复。另一方面，LRS可以提供直接的结构信息并确认重复的存在，但对低水平的异质性结构变异可能存在灵敏度不足。一种整合两种方法并结合临床相关性及对潜在不一致发现进行再分析的综合策略，有望改进mtDNA重排的表征与解读。
研究背景与科学问题
线粒体DNA（mitochondrial DNA, mtDNA）是细胞能量代谢的核心，编码氧化磷酸化（oxidative phosphorylation, OXPHOS）系统的关键组分。大规模mtDNA重排，包括千碱基级别的缺失和重复，已被广泛报道。单一大规模mtDNA缺失（single large-scale mtDNA deletion, SLSMD）通常导致一系列重叠的临床特征，如Pearson综合征（PS）、Kearns–Sayre综合征（KSS）及慢性进行性眼外肌麻痹（chronic progressive external ophthalmoplegia, CPEO）。相比之下，大规模mtDNA重复虽在KSS和肌病表型患者中有所报道，但其临床特征尚未被充分阐明。
目前，临床上普遍采用短读长测序（short-read sequencing, SRS）来检测mtDNA结构变异，其原理是通过分析不同区域的测序深度差异来推断拷贝数的增减。然而，标准分析流程存在一个根本性局限：当重复分子由一个全长mtDNA分子和一个相应的缺失分子首尾相接形成串联重复时，其断裂点序列与单纯缺失产生的序列完全相同。因此，SRS虽然能敏感地识别异常mtDNA群体的存在，却无法区分缺失与重复。这种误判可能导致错误的分子诊断和基因型-表型关联，进而影响预后评估和治疗决策。本研究正是为了解决这一临床诊断困境，评估长读长测序（long-read sequencing, LRS）在明确表征这类复杂mtDNA重排中的价值。
研究方法概述
研究对象来自费城儿童医院线粒体医学前沿项目（Mitochondrial Medicine Frontier Program, MMFP）的两名患者，均因临床表现不符合典型SLSMD综合征而接受进一步遗传学分析。研究人员采用了互补的测序策略：首先对已有SRS数据进行再分析，随后独立应用Pacific Biosciences（PacBio）HiFi测序和Oxford Nanopore Technologies（ONT）长读长测序平台对可疑重排进行表征。对于ONT测序，研究人员采用了靶向PCR扩增和PvuII线性化两种文库制备策略。此外，研究人员开发了一套基于断裂点感知参考重建的计算分析方法，用于从长读长数据中鉴定串联重复结构。
研究结果
病例1的mtDNA分析
该病例此前通过BluePrint Genetics的捕获法SRS检测报告了一个m.3577_15546del11970的大片段缺失，血液样本中异质性为60.5%。然而，患者的临床表现并不符合典型的KSS。研究人员的后续分析在尿液样本中通过长范围PCR联合SRS确认了相同的断裂点，并通过数字微滴PCR（digital droplet PCR, ddPCR）测得尿液样本中缺失异质性为43%。值得注意的是，研究人员还在此病例中检测到一个新的新生同质性错义变异m.5983G>A（p.Gly27Glu），位于MT-COI基因，该变异高度保守且功能预测工具提示可能影响蛋白质功能，但被归类为意义未明的变异。
PacBio LRS的结果彻底改变了对此结构变异的认知：LRS揭示，未被删除的mtDNA片段（标记为S区和E区）实际上是以串联重复的形式存在的。携带此重复的全长mtDNA分子跨度超过20kb，且所有捕获到断裂点的读段均显示为重复结构。因此，原本被判定为m.3577_15546缺失的结构，实际上是全长mtDNA分子与一个缺失分子融合后形成的m.15546_3577dup4600串联重复。断裂点处存在一段11bp的微同源序列（ccccctcccca）。
病例2的mtDNA分析
该病例此前通过GeneDx的检测报告了一个位于口腔拭子样本中的m.79_13637del13559缺失，异质性低于15%。研究人员在肌肉样本中通过长范围PCR联合SRS确认了相同断裂点。然而，对单核苷酸多态性（SNP）的分析揭示了关键线索：位于m.13637_79区域内的SNP呈现约50%的异质性，且这些SNP属于J1c7单倍群，而患者本身的单倍群为N1b1b。这表明在该区域存在额外的mtDNA片段，而非单纯的缺失。ddPCR定量显示此重复的异质性约为7%，且在母亲的口腔拭子样本中未检出。
然而，PacBio和ONT的LRS结果均未能提供确凿证据。PacBio仅获得一条覆盖重复区段的读段，ONT也仅有一条跨越全长的读段，但包含一段150bp的非预期插入。为了验证该读段是否为真实重复或人工产物，研究人员使用PvuII限制性内切酶对mtDNA进行线性化处理，结果显示仅有野生型单倍群，未发现任何支持串联重复的证据。因此，尽管SRS数据强烈提示存在m.13637_79dup3011的串联重复，但LRS因结构变异异质性水平过低而无法最终确认，该病例的诊断仍悬而未决。
讨论与结论
本研究的讨论部分指出，大规模mtDNA缺失与重复的形成机制尚不完全清楚，可能涉及复制滑脱、拷贝选择重组或双链断裂修复障碍等。MITOMAP和MitoBreak数据库收录了大量单一缺失，但关于mtDNA重复的报道极为有限，这很大程度上是由于传统检测方法（如Southern印迹、长范围PCR、限制性片段长度多态性分析）通量低、敏感性差，而SRS又无法区分二者。
病例1的发现极具启示性：LRS显示几乎所有异常mtDNA分子均为重复而非缺失。鉴于既往研究表明，携带mtDNA重复的细胞系无明显呼吸链缺陷，而携带缺失的细胞系则表现出严重功能障碍，这很好地解释了为何病例1虽有高达60.5%的报告异质性，却未表现出典型的SLSMD症状。研究人员据此提出，该病例严重的神经系统表现可能主要归因于新发现的MT-COI基因错义变异，而mtDNA重复本身可能是一个对表型影响较小的独立结构发现。
病例2则凸显了当前技术的另一项挑战。虽然SRS数据强烈指向重复，但由于结构变异异质性极低（约7%），即使采用高深度的LRS也无法捕获足够数量的跨重复全长读段。该额外mtDNA片段的起源尚不明确——它不属于患者自身的单倍群，在人核基因组中也未找到高度同源的序列，且其母系样本中不存在。尽管不能完全排除父系mtDNA传递的可能性，但与已知的双亲遗传案例不同，本病例的嵌合仅影响一个3011bp的片段，提示其形成机制可能与众不同。
研究结论
综合两个病例，本研究清晰地展示了SRS与LRS在表征大规模mtDNA重排中的互补价值。SRS能灵敏地识别mtDNA拷贝数异常，但无法区分缺失与重复；LRS则能提供直接的结构信息，但在低异质性结构变异面前敏感性受限。病例1证明了被SRS误判为缺失的mtDNA重排，经LRS评估后实为占主导地位的重复结构；病例2则警示，即使是LRS也可能不足以解析低水平且潜在复杂的mtDNA重排。这些发现强调了一种整合两种测序方法、结合正交实验手段和临床相关性的综合策略的重要性，这对于改进mtDNA重排的表征和解读至关重要。
