IF31.7｜单细胞Micro-C揭示多增强子枢纽调控基因表达

scMicro-C以5kb分辨率解析单细胞三维基因组，发现启动子-增强子条纹结构及多增强子枢纽。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

哺乳动物基因组中，增强子通过形成染色质环与启动子相互作用，调控基因的时空表达。然而，增强子与启动子之间的三维空间组织方式，尤其是在单细胞水平上的动态变化，仍不清楚。此前的研究主要依赖群体平均的Hi-C或Micro-C数据，无法解析单个细胞内的染色质折叠异质性。单细胞Hi-C技术虽然能够重建单细胞三维基因组结构，但受限于酶切位点，分辨率通常只有20kb，不足以研究精细的增强子-启动子相互作用。Micro-C利用微球菌核酸酶切割核小体间DNA，可获得核小体分辨率，但尚未应用于单细胞。

科学问题

卡在哪里： 单细胞三维基因组研究受限于分辨率，无法解析精细的增强子-启动子相互作用。群体平均的Micro-C虽然分辨率高，但掩盖了细胞间的异质性。此外，多增强子如何协同调控基因表达，其三维组织方式在单细胞水平上缺乏直接证据。

本文要回答： 开发高分辨率的单细胞三维基因组技术，解析增强子-启动子相互作用的结构基础，并揭示多增强子协同调控的机制。

技术路线

作者主要用了：scMicro-C、单细胞三维基因组重建、聚合物模拟、bulk Micro-C、ChIP-seq、ABC模型。

作者改进了bulk Micro-C流程，通过引入SDS提高连接效率、优化MNase消化程度、省略生物素富集并采用META全基因组扩增，建立了scMicro-C方法。随后对GM12878细胞进行单细胞测序，利用Dip-C算法重建单细胞三维基因组结构，并系统分析了启动子-增强子相互作用。结合bulk Micro-C、ChIP-seq和聚合物模拟，揭示了启动子-增强子条纹（PES）的形成机制及多增强子枢纽的组装。

图 1　Fig. 1。图1展示了scMicro-C的开发流程和优化结果。a为实验流程示意图。b显示SDS处理提高了连接效率，染色质片段长度分布更优。c显示SDS处理使单细胞接触数增加8.1倍。d显示不同MNase浓度下的染色质消化程度。e和f显示随着MNase浓度增加，接触比例和独特接触数增加，且scMicro-C优于Dip-C。g显示不同MNase浓度下染色质结构一致。h显示scMicro-C保留核小体占位和TF足迹，而scHi-C不能。i和j显示ensemble scMicro-C与bulk Micro-C在A/B区室、TAD和染色质环上高度一致。

核心亮点

亮点 1｜scMicro-C实现5kb分辨率

作者通过引入SDS处理、优化MNase消化和采用META扩增，开发了scMicro-C。与原始Micro-C相比，SDS处理使单细胞检测到的染色质接触数提高了8.1倍。在GM12878细胞中，724个高质量细胞平均获得835k个独特接触。聚合所有单细胞数据得到的ensemble scMicro-C可达到5kb分辨率，与bulk Micro-C在A/B区室、TAD和染色质环等特征上高度一致。单细胞三维基因组重建显示，约38.8%的细胞可解析5kb分辨率结构，且与bulk Micro-C接触频率高度反相关（Pearson r=-0.907）。

图 2　Fig. 2。图2展示了scMicro-C重建的单细胞三维基因组结构。a显示一个GM12878细胞的单倍型填补接触图和对应的5kb分辨率三维结构。b显示单细胞三维距离矩阵与bulk Micro-C接触频率矩阵在20kb、10kb和5kb分辨率下高度反相关。c显示散点图定量了这种反相关关系（Pearson r=-0.931、-0.922、-0.907）。d和e显示嵌套CTCF环在单细胞结构中倾向于同时形成。f展示了动态的环形成事件。

亮点 2｜发现启动子-增强子条纹结构

在bulk Micro-C数据中，作者发现许多基因的启动子处存在一种线状条纹结构，命名为启动子-增强子条纹（PES）。PES从转录起始位点（TSS）延伸至基因下游，与多个H3K27ac标记的增强子重叠。在GM12878中，共鉴定出819个PES基因，占4077个活跃长基因的20%。PES基因表达水平显著高于非PES基因（P=9.88×10-9），且启动子处cohesin结合更强。通过分析cohesin敲除的mESC Micro-C数据，发现PES形成依赖于cohesin介导的loop extrusion。

图 3　Fig. 3。图3展示了PES基因的特征和形成机制。a显示PES基因的pile-up分析结果，显示从TSS延伸的条纹信号。b显示PES基因表达水平显著高于非PES基因。c显示PES基因体上H3K27ac信号更强。d显示PES基因启动子处RAD21结合更强。e显示cohesin敲除后PES减弱，而WAPL敲除后PES延伸。f为PES形成的模型示意图。

亮点 3｜单细胞结构揭示多增强子枢纽

利用scMicro-C重建的单细胞三维结构，作者发现PES基因的启动子与多个增强子之间存在动态的相互作用。以EBF1基因为例，约30%的单细胞结构中启动子与下游增强子形成环状接触。更重要的是，在单个细胞中，多个增强子可同时与启动子相互作用，形成多增强子枢纽。随着枢纽中增强子数量增加，增强子间网络更紧密，染色质结构更紧凑（回转半径减小）。这种多增强子枢纽在IRF2等其他PES基因中也存在，表明其普遍性。

图 4　Fig. 4。图4展示了EBF1基因座的单细胞E-P相互作用和多增强子枢纽。a显示bulk Micro-C和ensemble scMicro-C接触图，显示PES和多个增强子。b显示按TSS相互作用排序的单细胞接触图，显示动态的条纹形成。c为EBF1基因座示意图。d显示启动子与各增强子形成环的比例。e展示单细胞三维结构中动态的E-P环。f显示多增强子同时与启动子相互作用的比例。g展示多增强子枢纽的三维可视化。h显示增强子数量增加导致更紧密的增强子网络。i显示增强子数量增加导致染色质更紧凑。

亮点 4｜聚合物模拟揭示形成机制

通过聚合物模拟，作者发现仅靠cohesin介导的loop extrusion和CTCF边界不足以形成PES和E-P/E-E环。需要引入增强子对cohesin的弱捕获以及启动子-增强子间的分子亲和力，才能重现实验观察到的条纹和环结构。模拟还显示，只有在强亲和力和cohesin捕获同时存在时，才能形成多增强子枢纽。这些结果表明，cohesin介导的loop extrusion促进启动子扫描下游增强子，而分子亲和力则稳定多增强子与启动子的空间聚集。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： scMicro-C、单细胞三维基因组重建、聚合物模拟、bulk Micro-C、ChIP-seq、ABC模型
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事三维基因组、基因调控和单细胞多组学研究的科研人员。

带走一句

scMicro-C将单细胞三维基因组分辨率提升至5kb，使精细的E-P相互作用研究成为可能，为解析多增强子协同调控提供了新工具。