3D-DNA 结构改变如何塑造 IDH1 与 DNMT3A 突变型 AML？

HiC 比甲基化数据多测了什么

这篇论文要回答的是：两种表观遗传调控因子 IDH1 和 DNMT3A 的突变，是否会重塑白血病细胞中 DNA 的物理折叠，以及这种重塑是否会产生单靠甲基化数据看不到的靶点。已知这两种突变都会全局改变 DNA 甲基化，而 IDH1 突变此前已被与 DNA 环形成改变和癌基因上调联系起来。本研究把这一问题扩展到全基因组，并并排比较这两种基因型。

核心测量手段是染色质构象捕获，即 HiC：它读出细胞核内哪些相距较远的 DNA 片段在物理上相互接触，从而生成全基因组互作图谱。作者从这些图谱中提取两个层次的结构。区室是大尺度区域，标记为活跃（A）或不活跃（B）；DNA 环是较短距离的接触，常锚定在 CTCF 结合位点，把增强子带到启动子附近。环分数（loop score）量化某个环存在的强度，分数升高意味着接触更频繁，分数下降意味着接触减弱。

细胞模型与患者样本中的共享环与基因型特异环

作者把 IDH1 p.R132H 或 DNMT3A p.R882H 分别导入 KG1a 髓系白血病细胞，构建细胞模型；这些细胞缺乏其他已知驱动突变。他们使用了 11 个 CRISPR 修饰克隆：3 个 IDH1 野生型对照、3 个 IDH1 突变克隆、2 个 DNMT3A 野生型对照和 3 个 DNMT3A 突变克隆。同时分析了 8 例初诊 AML 患者样本——3 例两基因均野生型、3 例携带 IDH1 p.R132H/C、1 例携带 DNMT3A p.R882H/C、1 例同时携带两种突变——这些样本作为探索性和验证性数据集。

在患者样本中共鉴定出 9,839 个环，其中 3,617 个为各基因型共有。环的获得部分共享（888 个环在两种突变型中均获得），部分特异：934 个仅在 IDH1 突变样本中获得，1,378 个仅在 DNMT3A 突变样本中获得。丢失也呈同样模式：IDH1 突变样本丢失 724 个，DNMT3A 突变样本丢失 872 个，两组共同丢失 1,426 个。唯一一例双突变患者同时携带与 IDH1 和 DNMT3A 相关的获得环和丢失环，作者将此视为支持突变特异性结构特征、而非单一通用白血病特征。

论文还把环变化与甲基化联系起来：高甲基化的 CTCF 锚点在丢失的环中富集，低甲基化的环则在获得的环中富集。这是甲基化状态与环状态之间的相关性，并未证明甲基化导致了环变化。

改变后的环锚点携带差异表达基因

为了把结构与转录联系起来，作者考察了差异表达基因的位置。在 IDH1 或 DNMT3A 突变型 KG1a 细胞中表达发生变化的基因，显著更多地位于丢失或获得环的锚点上，而不是两种基因型共有环的锚点上。锚点至少携带一个差异表达基因的环，其环分数在基因型间的变化幅度也更大。

在 KG1a 中，两种基因型的环分数变化都与基因表达变化正相关，且变化方向在细胞系与患者之间大体一致：IDH1 突变型 KG1a 与患者样本中环分数变化重叠的基因有 53%，DNMT3A 突变样本中有 60%，同时表现出相似的表达变化。重叠只是部分的，这决定了细胞系结果能外推多远。

在 KG1a 与患者共享的改变环中，作者列出的上调基因包括 IDH1 突变样本中的 EXT1、CHD7 和 IGF1R，以及 DNMT3A 突变样本中的 IKZF3、ETS1、SEPTIN9 和 IRS2。对于 IGF1R，基因位点周围相互作用强度增加在 KG1a 和 AML 样本中均得到确认，并与下游 CTCF 锚点处 CpG 甲基化降低同时出现。作者把环改变呈现为转录变化的潜在驱动因素；现有证据是跨基因型和样本类型的一致关联，而非此阶段直接的因果检验。

药物相关性筛选把环基因列为靶点

为了把环相关基因转化为候选脆弱点，作者使用了 BeatAML 2.0 数据集，该数据集将 300 多种化合物的离体药物反应谱与 805 例 AML 患者的基因表达配对；其中 543 例样本同时具备两类数据。由于突变特异性药物反应数据有限，他们在全部样本中计算药物 AUC 值与环基因表达的相关性，而不是在各基因型内部计算，然后把 Spearman 系数与来自 KG1a 和患者样本的基因表达 log2 倍数变化合并为秩积（rank product）。

在 IDH1 突变型 AML 中，48 个重叠环基因被排序；在 DNMT3A 突变型 AML 中为 57 个。对 IDH1 突变样本，ST8SIA6 和 ADAMTS3 与 LSD1 抑制剂 GSK-2879552 的相关性排名最高（Spearman 系数 −0.432，p = 4.09 × 10⁻⁵ 和 −0.289，p = 0.0076），IGF1R 则出现在与 IGF1R 抑制剂 BMS-754807 相关性最高的候选基因中（系数 −0.424，p = 0.0057）。对 DNMT3A 突变样本，TBC1D16 和 FMN1 分别与 ralimetinib（系数 0.381，p = 0.0133）和 HDAC 抑制剂 panobinostat（系数 −0.306，p = 1.28 × 10⁻⁷）的相关性排名最高。

按环基因靶点数量排序，IDH1 突变型 AML 中 venetoclax 和 BMS-754807 居前，DNMT3A 突变型 AML 中则是 elesclomol 和 panobinostat。BeatAML 数据还显示，IDH1 突变患者对 BMS-754807 比 IDH1 野生型患者更敏感，ralimetinib 也位列 DNMT3A 突变患者反应更好的药物之中。这些是在混合患者集合中的相关性排序，因此只是提名候选药物，并未确立基因型特异性疗效。

细胞系测试支持两个基因型特异候选药物

随后，作者用 WST-1 活力实验在 CRISPR 模型中直接测试了被提名的药物。IDH1 突变克隆对 BMS-754807 单药治疗敏感性增加，而 BMS-754807 与 IDH1 特异性抑制剂 ivosidenib 联用在 Loewe 协同模型下产生协同效应。对 GSK-2879552 的敏感性可见，但幅度不大。

DNMT3A 突变细胞对 ralimetinib 的敏感性在单药和与 venetoclax 联用时均增加。而 DNMT3A 野生型与突变型细胞对 BMS-754807 和 elesclomol 的敏感性差异在体外仅属轻微。

一项特异性检查强化了环与药物之间的联系：用 CRISPR 敲除 IGF1R 和 TBC1D16 的 CTCF 锚点，不仅逆转了基因表达，也逆转了相应的药物反应。该实验把候选脆弱点与特定锚点而非仅与突变本身联系起来，不过它是在细胞系模型中完成的。

结论成立的范围与变薄之处

患者数据集规模小且不均衡：仅 8 例样本，其中只有 1 例 DNMT3A 突变和 1 例双突变，作者将其描述为探索性和验证性数据。患者样本在基因表达和 3D-DNA 改变上也差异很大，这与 AML 的复杂图景一致。作者把依赖 KG1a 作为主要细胞模型列为一项局限，同时指出 KG1a 与患者之间存在共享改变和重叠的转录模式，以此说明该机制并非模型特有。

药物相关性这一步也是间接的。由于突变特异性药物反应数据有限，相关性是在全部 543 例样本中计算的，而不是在基因型内部计算，因此某个环基因与药物 AUC 的相关性本身并不能表明该药物在突变细胞中通过该环起作用。真正把 BMS-754807 和 ralimetinib 从排序候选推进为可检验假设的，是体外协同实验和锚点敲除结果。

为何 IGF1R 环对 IDH1 突变型 AML 最重要

这篇论文最独特的论断不是 3D 结构在 AML 亚型之间存在差异，而是某一特定结构改变指向某一特定药物组合。在 IDH1 突变型 AML 中，IGF1R 位点相互作用强度增加，其下游 CTCF 锚点 CpG 甲基化降低，表达上调；同样的细胞对 BMS-754807 敏感，且当 IGF1R 抑制剂与 ivosidenib 联用时敏感性进一步增强。作者指出，IGF1R 上调此前未在 IDH1 突变型 AML 中报道过，并将其与 MAPK 信号激活这一可能的生存机制联系起来。

对 DNMT3A 突变型 AML，对应的候选药物是 ralimetinib，依据是 TBC1D16 的排序以及突变细胞系中单药和 venetoclax 联用的敏感性。作者认为这尤其值得测试，因为 DNMT3A 突变患者对 venetoclax 的反应存在差异。因此，数据支持的实践优先级是一个基因型分野：IDH1 突变疾病用 IGF1R 抑制联合 ivosidenib，DNMT3A 突变疾病用 p38/MAPK 抑制；两者都仍需临床评估，尚不能作为既定疗法。

本文引用的文献