IF50.5｜核小体纤维拓扑引导转录因子结合增强子

转录因子组合不是随机搜索增强子，而是沿着核小体纤维上的“路标”定向移动。

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

研究背景

细胞身份由多种转录因子（TF）共同结合到细胞类型特异性基因的增强子上决定。尽管TF能识别可及染色质内的特定DNA基序，但仅凭基序信息无法解释TF如何选择增强子。在重编程过程中，先锋TF如OCT4、SOX2、KLF4（OSK）可结合关闭染色质并启动后续激活，但组合结合如何定位到远端增强子仍不清楚。已有研究表明，单个TF可在单核小体上识别基序，但组合结合是否依赖更高阶的染色质结构尚不明确。

科学问题

卡在哪里： 此前研究多关注单个TF在单核小体上的基序识别，但无法解释多个TF如何协同定位到特定增强子。单核小体上的基序排列不能预测组合结合事件，提示可能存在更高阶的染色质纤维拓扑信号，但该信号的性质和作用机制未知。

本文要回答： 本文旨在揭示TF组合如何利用核小体纤维的3D拓扑和基序排列（“路标元素”）引导其结合到细胞类型特异性增强子。

技术路线

作者主要用了：ChIP-seq、ATAC-seq、MNase-seq、Micro-C、RNA-seq、EMSA、重编程系统。

作者构建了四种不同的TF组合（OSKM、GETM、GETMR、BS9G4M）诱导小鼠胚胎成纤维细胞重编程，系统比较了早期和终末阶段的TF结合、染色质可及性、核小体定位和3D基因组结构。通过整合ChIP-seq、ATAC-seq、MNase-seq和Micro-C数据，在单核小体分辨率下分析TF结合模式，并利用基序分析、报告基因实验和H1水平扰动验证了核小体纤维拓扑在TF组合结合中的作用。

图 1　Fig. 1。图1展示了四种重编程组合（OSKM、GETM、GETMR、BS9G4M）的实验设计和RNA-seq主成分分析，显示GETMR重编程轨迹分叉为iPS和iTS。热图显示OSK基序在开放和关闭染色质核小体上的密度分布，揭示先锋TF在关闭染色质核小体上具有独特的基序读取模式。

核心亮点

亮点 1｜单核小体基序无法解释组合结合

作者比较了OSK、GET等TF在单核小体上的基序分布，发现组合结合（combo-nucs）与单独结合（solo-nucs）的基序排列模式显著不同，但组合结合事件中不同TF基序的共发生频率与随机期望无显著差异（P(OSK)≈P(O)P(S)P(K)）。这表明单核小体上的基序识别只能解释TF的单独结合，组合结合需要更高阶的染色质结构信息。

图 2　Fig. 2。图2显示OSK在早期重编程中结合到大的核小体阵列（平均6个核小体，约7 kb），且基序在阵列内外呈方向特异性分布。排序后的热图显示SOX2基序方向性将阵列分为两类，KLF4基序在边界形成方向相反的峰。重编程完成后OSK结合位点侧向移位至邻近增强子，伴随H3K27ac沉积和H1排除。

亮点 2｜核小体阵列上的基序语法具有方向性

OSK在早期重编程中结合到平均含6个核小体、跨度约7 kb的核小体阵列上。这些阵列内外的OSK基序呈现方向特异性聚集，尤其SOX2基序的方向性可将阵列分为两类。KLF4基序在阵列中心和远端边界形成方向相反的峰。这种基序排列模式在0.8-2.8 kb大小的阵列中最为明显，且伴随核小体相位化，提示核小体纤维上的基序语法可能引导OSK向阵列边界聚集。

图 3　Fig. 3。图3展示Nanog位点的OSK结合从早期（路标元素）到终末（增强子）的移位。双荧光报告基因实验显示，路标元素方向翻转后，eGFP在重编程早期被激活，而tdTomato直到重编程完成才表达，证明路标方向性对增强子激活至关重要。

亮点 3｜路标元素引导OSK侧向移位至增强子

在重编程完成时，OSK结合位点从初始核小体阵列的边界侧向移位至邻近的增强子区域（约500 bp），伴随染色质可及性和H3K27ac增加。报告基因实验显示，将Nanog增强子-启动子之间的路标元素方向翻转后，eGFP报告基因在重编程早期被激活，而tdTomato报告基因直到重编程完成才表达，证明路标元素的方向性对于增强子激活至关重要。

图 4　Fig. 4。图4显示Micro-C分析揭示OSK核小体阵列边界之间存在染色质环相互作用，且环锚点接触强度在重编程后降低。H1.4敲低或过表达改变核小体纤维折叠构象，并抑制OSK结合和重编程效率，表明染色质纤维拓扑调控TF组合结合。

亮点 4｜OSK路标元素位于染色质环内

Micro-C分析显示，早期重编程中OSK核小体阵列的边界之间存在强烈的染色质环相互作用，且环锚点处的接触强度在重编程完成后显著降低。H1.4敲低或过表达改变了核小体纤维的折叠构象（从zig-zag螺旋变为untwisted ladder或更扭曲的螺旋），并显著抑制OSK结合和重编程效率，表明染色质纤维拓扑而非整体可及性调控TF组合结合。

图 5　Fig. 5。图5显示GET结合到较小且富含H1和核小体间接触的核小体阵列，TFAP2C基序位于阵列边界。GET阵列在TAD边界形成条纹，并与TS细胞增强子通过环连接。H1.4扰动阻断GET介导的重编程，提示GET通过染色质引导的易位到达增强子。

亮点 5｜GET结合高度互联的路标元素

与OSK不同，GET（GATA3、EOMES、TFAP2C）结合到较小（3-5个核小体）且富含H1和核小体间接触的核小体阵列上。TFAP2C基序（回文序列）位于阵列边界，使GET阵列分为两类。GET阵列在拓扑相关域边界形成条纹，并与TS细胞增强子通过环显著连接。H1.4扰动同样阻断GET介导的重编程，提示GET通过染色质引导的易位机制到达增强子。

亮点 6｜MYC与不同TF组合遵循不同基序语法

MYC在OSK阵列中显示E-box基序的中心缺失和边界富集，但在GET阵列中其结合完全不依赖E-box。AlphaFold-Multimer预测和EMSA/免疫共沉淀实验证实MYC与TFAP2C存在直接蛋白-蛋白相互作用，表明MYC与GET的结合由蛋白互作驱动，而非基序识别。

生信可带走

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

• 方法栈： ChIP-seq、ATAC-seq、MNase-seq、Micro-C、RNA-seq、EMSA、重编程系统
• 公开数据： GSE201852、GSE98124、PRJNA20689
• 代码： 公开仓库 https://github.com/Rdatatable/data.table、https://github.com/YaqiangCao/cLoops2、https://github.com/broadinstitute/picard；本地已克隆：data.table、FastQC、cLoops2、pairtools、distiller-nf
• 谁该点开原文： 研究转录因子组合、染色质结构和细胞重编程的生信与实验人员。

带走一句

分析TF组合结合时，不要只看单核小体基序，要关注核小体阵列的3D拓扑和基序方向性，它们可能才是决定增强子选择的关键。