IF31.7｜活细胞成像揭示cohesin介导的环挤出动力学在人类细胞中高度均一

TAD锚点每小时接触一次，持续6–19分钟，cohesin动力学参数跨区域高度一致。

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

研究背景

哺乳动物基因组被折叠成数千个拓扑关联结构域（TAD），这些结构由cohesin介导的环挤出形成，并被收敛取向的CTCF位点所界定。Hi-C等群体测序方法揭示了TAD的静态特征，但无法直接观测环挤出的动态过程。体外实验测得cohesin的挤出速度为0.5–1 kb/s，然而在活细胞中，由于存在染色质环境、转录等多种因素，实际速度可能不同。此前在mESCs中的研究显示TAD锚点接触短暂且罕见，但仅局限于单个位点，且mESCs与分化细胞在染色质状态和细胞周期上存在差异，因此亟需在人类分化细胞中系统研究环挤出动力学。

科学问题

卡在哪里： 此前对环挤出的研究主要依赖体外实验或mESCs中的单个位点，缺乏在人类分化细胞中跨多个内源基因组区域的系统定量。同时，cohesin动力学参数（密度、驻留时间、速度）是否在不同基因组区域存在局部调控，以及CTCF位点分布如何影响TAD动力学，仍不清楚。

本文要回答： 在活的人类细胞中直接观测多个内源TAD锚点的运动，定量环挤出动力学参数，并探究这些参数在不同基因组区域间是否一致，以及CTCF位点如何决定TAD动力学。

技术路线

作者主要用了：CRISPR基因组编辑、活细胞三维成像、Capture Micro-C、聚合物模拟。

作者利用CRISPR在HCT116细胞中将TetOx96和CuOx150重复阵列分别插入到三个不同TAD（L1、L2、T1）的锚点，并通过TetR-splitGFP和CymR-Halo标记实现双色荧光成像。同时构建了RAD21-mAID降解系统，通过添加auxin条件性去除cohesin，以区分cohesin依赖与非依赖的染色质运动。通过长时间三维追踪锚点间距离，结合聚合物模拟，系统比较不同区域的环挤出动力学参数。

图 1　Fig. 1。图1展示了实验系统的设计：在HCT116细胞中通过CRISPR将TetOx96和CuOx150阵列插入到三个不同TAD的锚点，并用TetR-splitGFP和CymR-Halo标记。图中对比了有无auxin处理下的锚点距离分布，显示去除cohesin后距离增大，验证了系统的有效性。

核心亮点

亮点 1｜TAD锚点接触短暂且频繁

通过定义‘近端状态’（锚点距离持续小于阈值），作者发现TAD锚点平均每小时接触1.1–1.2次，每次持续6–19分钟，占G1期总时间的7–26%。相比之下，auxin去除cohesin后，近端状态比例下降3.1–3.9倍，频率降至约0.4次/小时，寿命缩短约1.3倍。这些结果表明TAD是高度动态的结构，锚点接触由cohesin介导的环挤出驱动，但接触本身是瞬时的。

图 2　Fig. 2。图2量化了TAD锚点的近端状态。a显示原始距离时间序列及分割出的近端状态（蓝色条）。b–d分别展示近端状态的比例、频率和寿命，auxin处理显著降低比例和频率，但寿命缩短幅度较小。e展示通过模型估计的三种状态比例，显示TAD几乎无开放状态。

亮点 2｜TAD几乎持续处于挤出状态

利用基于锚点距离分布的模型，作者将染色质状态分为开放、挤出和闭合三种。在未处理细胞中，L1和L2的开放状态比例中位数为0%，T1为11%，而挤出状态比例分别为76%、79%和73%。在RAD21降解后，挤出状态完全消失，开放状态升至89–94%。这表明TAD内部几乎始终有cohesin在挤出环，而非处于完全松弛状态，支持多cohesin同时作用的模型。

图 3　Fig. 3。图3解释了闭合速率（CR）与cohesin分子速度的区别，并展示了如何从距离时间序列中提取CR。通过对齐近端状态开始前的数据并拟合线性下降，作者估计了各区域的CR。随机打乱数据作为对照，证实了下降源于过程性挤出。

亮点 3｜锚点闭合速率约为0.1 kb/s

通过对齐近端状态开始前的距离时间序列并拟合线性下降段，作者估计了闭合速率（CR）：L1和L2为0.07 kb/s，T1为0.16 kb/s，half TAD为0.09 kb/s。随机打乱时间序列后，约90%的bootstrap样本不再显示线性下降，证实了该下降源于过程性的环挤出。结合模拟，作者推断cohesin的分子速度低于0.5 kb/s，显著慢于体外测得的0.5–1 kb/s，提示体内存在减速因素。

图 4　Fig. 4。图4展示了聚合物模拟的参数设置和结果。a列出变化的cohesin参数（密度、驻留时间、速度）。b显示四个区域CTCF位点的驻留时间分布。c和d展示不同参数组合下的接触图谱和p(s)曲线，说明cohesin参数对TAD结构的影响。

亮点 4｜cohesin动力学参数跨区域一致

通过聚合物模拟，作者系统扫描cohesin密度（1–40 Mb⁻¹）和驻留时间（2–33 min），并比较模拟与实验的Micro-C接触图谱和成像数据。结果显示，四个基因组区域的最佳参数组合高度重叠，支持统一的cohesin密度（约12–18 Mb⁻¹）和驻留时间（约15–25 min）。区域间的动力学差异主要源于CTCF位点的位置和结合强度，而非cohesin参数的局部调节。

图 5　Fig. 5。图5比较了模拟与实验的偏差。a和b分别展示在两种假设速度下，不同参数组合与Micro-C和成像数据的吻合程度。c显示四个区域的最佳参数区域重叠，支持统一的cohesin动力学。d和e通过模拟不同速度下的闭合速率，推断体内cohesin速度低于0.5 kb/s。

亮点 5｜TAD内部存在多个同时挤出的环

模拟显示，在最佳参数下，TAD锚点之间平均由2.5–3.7个内部环连接，未挤出的DNA长度在112–161 kb之间，尽管TAD大小相差2.7倍，但有效基因组距离仅相差1.4倍。此外，TAD外部的DNA序列也可能参与锚点间的最短路径，表明邻近区域的CTCF和cohesin对局部染色质相互作用有重要影响。

生信可带走

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

• 方法栈： CRISPR基因组编辑、活细胞三维成像、Capture Micro-C、聚合物模拟
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 研究三维基因组、cohesin/CTCF功能、或希望将环挤出参数纳入定量模型的学者。

带走一句

环挤出动力学参数可能在全基因组范围内均一，TAD差异主要由CTCF位点决定，这为基于CTCF分布的预测模型提供了基础。