Spatial mitochondrial lineage tracing uncovers a premetastatic niche and microenvironment programmed fate switching in osteosarcoma

技术摘要：空间线粒体谱系追踪揭示骨肉瘤的前转移生态位与微环境程序化命运转换

问题陈述
骨肉瘤的进展以复杂的时空动态为特征，然而控制肿瘤细胞命运决策的谱系关系仍不为人所知。虽然单细胞 RNA 测序（scRNA-seq）揭示了瘤内异质性，但它将细胞从其原生组织架构中分离出来，抹去了对理解细胞邻域至关重要的空间坐标。相反，空间转录组学（ST）保留了架构，但仅捕捉静态快照，无法区分相邻细胞是共享共同的克隆起源，还是通过趋同演化汇聚到了相似的状态。此外，现有的谱系追踪方法通常依赖于转基因操作（如 CRISPR 条形码），这使得它们无法应用于原发性人类临床标本。迫切需要利用内源性示踪剂整合空间、时间和谱系历史，以确定促转移状态是由瞬时微环境可塑性产生的，还是由决定性的克лоnal 锁定（clonal locking）驱动的。

方法论
作者在一个单一的、未经治疗的、晚期非转移性人类骨肉瘤标本上，整合了三个分析维度——空间、时间和谱系：

1. 单细胞与空间转录组学： 对 11,203 个细胞进行 scRNA-seq，以定义细胞群体并推断拷贝数变异（CNVs），从而区分恶性与非恶性细胞。在同一样本上进行空间转录组学（10x Visium）分析，将基因表达映射到组织架构上。
2. 伪时间重建： 作者应用了 cell2fate（一种贝叶斯 RNA 速率模型）来推断转录动力学和基因模块，重建了肿瘤状态和基质重构的时间进程。
3. 内源性谱系追踪： 为了在无需遗传操作的情况下解析克隆祖先，作者利用了自然发生的体细胞线粒体 DNA（mtDNA）变异。他们通过 MASTER 方案从 Visium cDNA 文库中富集了线粒体转录本。
4. 空间与系统发育分析：
   • 空间相干性： 应用基于 Moran's I 的空间自相关过滤器对线粒体变异等位基因频率（VAFs）进行处理，以确认克隆空间相干性并过滤掉深度驱动的伪影。
   • 谱系重建： 使用结合了 Wright-Fisher 漂变模型和隐马尔可夫树（WF-HMT）的 MitoDrift 框架来重建系统发育树，并模拟随机线粒体分离。
   • 命运转换建模： 拟合离散状态连续时间马尔可夫模型（ARD），以估计细胞状态（COL3A1⁺、ALPL⁺、THY1⁺ 等）之间的定向转换速率。
   • 相互作用分析： 使用 CellChat 和 SpaCET 在空间背景下分析细胞间通信，以识别信号枢纽。

关键结果

• 肿瘤异质性与 CNV 分析： scRNA-seq 识别出 12 个截然不同的群体，其中肿瘤细胞（66.5%）占据主导地位。恶性亚群包括 ALPL⁺（成骨细胞型）、THY1⁺（间充质/干细胞样）和 COL3A1⁺（祖细胞/基质重构型）细胞。CNV 分析证实了这些亚群的恶性起源，其中 THY1⁺ 细胞表现出独特的 16 号染色体增益。
• 逐步命运分叉： 伪时间重建揭示了一个逐步级联过程：COL3A1⁺ 祖细胞分叉为两个截然不同的谱系：ALPL⁺ 成骨谱系和 THY1⁺ 间充质谱系。这种分叉在时间上与髓系驱动的重构模块的峰值活性相吻合。
• 差异化转录程序：
  • 向 ALPL⁺ 谱系的转变是由 EMT 程序和向氧化磷酸化转变的代谢切换驱动的。
  • 向 THY1⁺ 谱系的转变以 NF-κB 介导的炎症重编程、抗原处理下调和代谢重构为特征，这与由 EMT 驱动的 ALPL⁺ 路径截然不同。
• THY1_Endo 生态位的发现： 空间转录组学识别出一个特定的促转移生态位，其中 THY1⁺ 肿瘤细胞与 PLVAP⁺ 功能障碍内皮细胞共定位。该 THY1_Endo 生态位富含 PTN（多效蛋白）和 NOTCH 信号，形成了一个空间受限的信号枢纽。细胞间相互作用分析证实，THY1⁺ 细胞与内皮细胞及特定巨噬细胞亚群之间的信号相互作用最为强烈，这种关系在其他肿瘤亚型中并未观察到。
• 谱系-状态耦合与定向命运锁定：
  • 空间相干性： 线粒体变异显示出显著的空间自相关性（Z = 10.07, p ≈ 4×10⁻²⁴），证实了 mtDNA VAFs 在 Visium 分辨率下保留了足够的克隆信息。
  • 不可逆性： 使用 WF-HMT 模型进行的系统发育重建揭示了显著的转换速率不对称。从 COL3A1⁺ 到 THY1⁺ 的速率为 724.08，而反向转换速率仅为 4.09（约 177 倍的不对称性）。相比之下，COL3A1⁺ 到 ALPL⁺ 的转换速率较低（0.47）。
  • 克隆独立性： 低于 0.1 的 Blomberg's K 值表明，THY1⁺ 和 ALPL⁺ 状态多次独立地起源于不同的 COL3A1⁺ 祖细胞，然而一旦建立，THY1⁺ 命运基本上是不可逆的。
  • 突变特征： 研究识别出了分支特异性的线粒体突变，存在于高纯度的 THY1⁺ 和 ALPL⁺ 克隆中，这表明表型锁定之前存在突变启动过程。

意义与主张
本文提出了一个骨肉瘤进展的“时间-空间-谱系”接力模型。作者声称，转移性的 THY1⁺ 间充质命运不仅受细胞自主程序的控制，还需要通过与功能障碍内皮细胞直接共定位来获得严格的空间许可，从而形成 THY1_Endo 生态位。

主要主张包括：

1. 方法学验证： 本研究证明了体细胞线粒体变异可以作为强大的内源性谱系示踪剂，用于研究原发性人类肿瘤，能够在无需转基因操作的情况下解析克隆历史和命运决策。
2. 机制见解： 研究确立了 THY1⁺ 状态是一种不可逆的、克隆锁定的终末分化状态，是由从 COL3A1⁺ 祖细胞进行的定向转换驱动的，而非一种瞬时的可塑状态。
3. 生态位依赖性： THY1⁺ 细胞的转移能力取决于由 THY1_Endo 生态位提供的空间受限信号环境（PTN/NOTCH）。
4. 临床意义： THY1_Endo 生态位被提议为一个切实的空间靶点。作者建议，治疗策略应旨在同时破坏谱系承诺（不可逆的命运）和生态位支持（空间信号），以拦截转移级联反应。

研究结论指出，整合 RNA 速率、空间转录组学和线粒体谱系追踪，提供了一个理解并可能拦截癌症演化的多维框架，弥合了单细胞状态、空间组织、时间动态与克隆祖先之间的鸿沟。