单细胞多组学解析线粒体DNA突变嵌合性与动态异质性

《Nature Communications》：Single-cell multi-omic analysis of mitochondrial mutational mosaicism and dynamics

这项研究旨在解决单细胞水平上全面评估线粒体DNA（mtDNA）突变图谱的技术瓶颈。研究人员利用线粒体单细胞ATAC-seq（mtscATAC-seq）和mtDNA高突变的POLGD274A细胞系，引入了单细胞mtDNA每百万碱基对突变数（scmtMPM）和异质性加权线粒体局部约束分数（scwMSS）两项指标，量化了单细胞的突变负荷和体细胞嵌合性。研究表明单个POLGD274A细胞具有复杂的突变图谱，而致病性突变和截断性变异仅以亚阈值水平存在。在人体样本中，该框架揭示了先前未被认识的单细胞水平mtDNA突变景观异质性。该研究为线粒体遗传学、疾病和体细胞嵌合性的单细胞研究提供了新框架。

　　

在细胞这个微小的生命工厂里，线粒体是至关重要的能量供应站，其内部拥有一套独立于细胞核的环状DNA——线粒体DNA（mtDNA）。然而，这套“能量工厂”的蓝图并不总是完美无缺，其突变率比核DNA高出许多，并与衰老、癌症、神经肌肉疾病等多种人类疾病密切相关。与通常认为的不同，mtDNA在每个细胞中拥有成百上千个拷贝，这意味着一个细胞可以同时存在野生型和突变型mtDNA，这种现象被称为“异质性”（heteroplasmy）。突变的比例（变异等位基因频率，VAF）高低，直接影响着细胞的功能和命运。传统上，人们认为细胞需要积累60%-80%的突变mtDNA，其代谢缺陷才会显现出来，但近期遗传学证据显示，即使在更低的VAF水平，表型也可能已经出现。

长期以来，研究人员依靠“群体”测序（bulk sequencing）来研究mtDNA变异，这种方法只能提供一个细胞群体的平均信号。但细胞是异质性的，每个细胞的mtDNA拷贝数、突变种类和比例都可能千差万别。这种“单细胞水平”的突变图谱“马赛克”和动态变化，我们却知之甚少。建立细胞水平的“基因型-表型”关联一直充满挑战，原因就在于这种异质性和细胞状态依赖的mtDNA动态特性。近年来，单细胞测序技术的发展，特别是线粒体单细胞ATAC-seq（mtscATAC-seq），实现了在单次实验中对成千上万个细胞的mtDNA进行基因分型，并能同时获取染色质可及性信息，为我们打开了一扇窥探单细胞线粒体遗传世界的新窗口。但如何系统性地解读单细胞水平的mtDNA遗传变异，量化其整体突变负荷，并评估其功能影响，仍缺乏一个综合性的分析框架。

为了填补这一空白，并深入探索mtDNA突变的基本特性，研究人员在《Nature Communications》上发表了一项研究。他们整合了先进的单细胞多组学技术、具有已知高突变率的细胞模型以及来自线粒体病患者和健康供体的人体样本，引入了一套创新的分析指标，首次在单细胞分辨率上系统描绘了mtDNA的突变景观、动态变化及其生物学意义。

为了开展这项研究，作者主要运用了以下几个关键技术方法：

1. 1.
   线粒体单细胞ATAC-seq（mtscATAC-seq）：该技术可同时对单个细胞的染色质开放区域和mtDNA进行高通量测序，从而在获取细胞表观遗传状态的同时，对每个细胞的整个线粒体基因组进行基因分型。研究人员利用该技术分析了对照HEK293细胞、两种DNA聚合酶γ（POLG）校对缺陷（POLG^D274A）的HEK293细胞系，以及来自2名健康供体和4名线粒体脑肌病伴乳酸酸中毒和卒中样发作（MELAS）综合征患者的外周血单个核细胞（PBMC）。
2. 2.
   细胞哈希标记（Cell Hashing）：用于在混合样本中对不同细胞系或不同处理条件的细胞进行多重标记和后续解卷积，实现在单次测序运行中平行分析多个样本。
3. 3.
   单细胞mtDNA突变负荷量化指标：研究团队引入了两个核心计算指标：单细胞mtDNA每百万碱基对突变数（scmtMPM） 和 异质性加权线粒体局部约束分数（scwMSS）。scmtMPM借鉴了肿瘤基因组学中量化突变负荷的思路，对每个细胞检测到的mtDNA突变数进行测序深度归一化。scwMSS则整合了来自gnomAD数据库的线粒体基因组局部约束（MLC）模型，将每个突变位点的约束分数（反映其功能重要性）乘以其在细胞中的VAF后进行加权求和，从而评估突变集合的潜在功能影响。
4. 4.
   代谢压力模型：使用半乳糖培养细胞，以迫使细胞依赖氧化磷酸化（OXPHOS）供能，从而施加代谢选择压力，观察携带高mtDNA突变负荷的细胞如何适应和响应。

研究结果

1. POLG^D274A细胞中线粒体DNA突变负荷的升高与异质性分布

研究人员首先利用mtscATAC-seq分析了对照（CTRL）和两种POLG^D274A细胞系（KI36和KIA2）。正如预期，突变细胞系的mtDNA测序深度和检测到的mtDNA变异数量均显著高于对照，且突变特征以C>T转换为主，表明复制错误是主要驱动因素。重要的是，单细胞分析揭示了远比传统群体测序所知的更为复杂的突变景观。单个POLG^D274A细胞可携带数百个mtDNA变异，且其VAF分布呈现高度异质性。新产生的“细胞系特异性”变异主要存在于低VAF（<5%），而“亲本共有”的变异则倾向于具有更高的VAF，这与体细胞突变通过遗传漂变和营养分离逐渐积累的预期相符。

2. 有害mtDNA变异的异质性动态受限

进一步将变异按功能分类（同义、错义、截断、致病性）分析发现，同义和错义变异的VAF分布范围很广，部分可达到较高水平。然而，致病性变异和截断变异（如无义突变、移码突变）的VAF分布则高度偏斜，几乎全部维持在很低的水平（通常远低于60-80%的经典生化阈值）。这些有害变异在细胞群中也表现出较低的香农熵值，表明其异质性动态受到限制，提示在长期培养过程中受到了负向选择。这一发现表明，POLG^D274A相关的线粒体缺陷可能并非由单个或少数几个高VAF的有害变异驱动，而是由遍布整个基因组、成百上千个低VAF变异的累积负担所导致。

3. 引入量化单细胞mtDNA突变负荷的新框架：scmtMPM与scwMSS

鉴于单个细胞可能携带大量分布各异的mtDNA变异，研究团队提出了scmtMPM和scwMSS两项指标，用于量化单细胞的整体mtDNA突变负荷。应用该框架分析POLG^D274A细胞，发现其scmtMPM和scwMSS值均显著高于对照，且两者在突变细胞中呈现良好的相关性。该框架还能解析不同功能基因群（如编码复合物I、III、IV、V，tRNA，rRNA的基因）的突变负荷，揭示了突变在基因组范围内的分布并非完全随机，不同基因群间存在一定的相关性。

4. 半乳糖压力诱导全局性适应而非克隆选择

当用半乳糖处理POLG^D274A细胞以施加OXPHOS依赖的代谢压力时，细胞出现了明显的适应反应，包括mtDNA拷贝数的显著增加，以及染色质可及性在GDF15、CCL20等与线粒体/内质网应激和炎症信号通路相关基因上的变化。然而，无论是群体还是单细胞水平，主要mtDNA变异的VAF分布并未发生显著偏移，也未观察到明显的基于特定mtDNA变异谱的克隆选择。这表明，在这种普遍高突变负荷但缺乏单个主导性致病突变的背景下，细胞通过全局性的拷贝数扩增和转录重编程来应对代谢压力，而非清除特定有害变异。

5. 人体免疫细胞中线粒体DNA体细胞嵌合性的异质性

将该分析框架应用于健康供体和MELAS患者的PBMC，成功量化了人体内单细胞水平的mtDNA突变负荷。结果显示，scmtMPM和scwMSS评分存在年龄相关性和供体/细胞类型特异性异质性。在MELAS患者中，携带致病性mt.3243A>G变异（位于MT-TL1基因）的细胞显示出更高的突变负荷。有趣的是，当从评分中剔除该致病变异后，年轻患者的双峰分布消失，但老年患者仍保持较高的突变负荷，提示随着年龄增长，新的体细胞突变不断累积。此外，通过scMMB评分，研究还成功识别出在特定个体和细胞类型中富集、可能具有功能影响的复杂I基因变异（如mt.10599G>A， mt.10270T>C， mt.10398A>G），这些变异在免疫细胞亚群中显示出克隆性扩张的迹象，表明了细胞谱系对特定mtDNA变异的耐受性或偏好性。

研究结论与意义

这项研究通过整合先进的单细胞多组学技术和创新的计算框架，首次在单细胞分辨率上系统性地描绘和量化了线粒体DNA的突变景观、动态异质性和累积负荷。研究不仅证实了POLG突变细胞中存在远超预期的体细胞mtDNA变异复杂性，揭示了有害变异在负向选择下被限制在低水平的动态规律，更重要的是，它建立了一套可广泛应用于任何单细胞mtDNA测序数据的分析标准——scmtMPM和scwMSS指标。这套框架使得研究人员能够超越对单个变异的关注，从“基因组突变负荷”和“潜在功能影响”两个维度，评估单个细胞的线粒体遗传状态。

该研究的核心意义在于：

1. 1.
   提供了新工具与新视角：scmtMPM和scwMSS为线粒体遗传学、衰老研究、线粒体病和癌症生物学等领域提供了强大的单细胞分析工具，使得比较不同细胞、不同个体、不同疾病状态的线粒体遗传“健康度”成为可能。
2. 2.
   深化了对疾病机制的理解：研究表明POLG相关缺陷可能与广泛的低水平突变累积负担有关，而非少数高负荷致病突变。在人体中，该框架揭示了致病性mtDNA变异如何影响体细胞突变景观的演变，并识别出具有潜在病理意义的细胞类型特异性克隆扩张。
3. 3.
   连接基因型与表型：通过同时获取染色质可及性数据，研究展示了如何将线粒体基因型与核基因组调控响应关联起来，为理解“核-线粒体”通讯在应对线粒体应激中的作用提供了范例。
4. 4.
   具有广泛的临床应用潜力：该框架可用于识别疾病相关的mtDNA变异模式，评估线粒体突变负荷在衰老和疾病进展中的风险，并可能在未来助力于线粒体相关疾病的诊断、预后判断和治疗策略开发。

总而言之，这项研究标志着我们对线粒体遗传的认识从“群体平均”迈入了“单细胞解析”的新时代，为在单细胞水平上探索线粒体在健康与疾病中的核心作用奠定了方法论基础，并开辟了新的研究方向。