IF50.5｜人类4D核组结构功能整合图谱

整合多模态数据，绘制人类胚胎干细胞和成纤维细胞的高分辨率4D核组图谱。

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

研究背景

基因组的三维空间组织与基因功能密切相关。过去二十多年，研究者已注释了人类基因组中的大量功能元件，但基因组如何通过空间折叠调控基因表达仍不完全清楚。染色质在多个尺度上组织：从核小体、染色质环、拓扑关联域（TAD）到核内区室和核体。这些结构动态变化，影响增强子-启动子互作、DNA复制和转录。4D核组计划旨在系统描绘基因组在三维空间和时间中的折叠，并建立结构与功能的因果关系。

科学问题

卡在哪里： 尽管已有多种染色质构象捕获技术，但不同方法在检测基因组折叠特征（如区室、TAD边界、染色质环）时的定量差异和互补性尚不明确。此外，缺乏整合多模态数据（如Hi-C、DamID、TSA-seq、SPRITE）的系统框架，难以全面注释基因组空间状态并关联功能。

本文要回答： 系统比较和基准化多种基因组折叠检测方法，整合多模态数据构建人类4D核组图谱，并解析基因组结构、核体定位与转录、复制等功能的关系。

技术路线

作者主要用了：Hi-C、Micro-C、SPRITE、GAM、ChIA-PET、PLAC-seq、DamID、TSA-seq、Repli-seq、iMARGI、单细胞Hi-C、3D建模。

作者首先在H1胚胎干细胞和HFFc6成纤维细胞中生成并比较了多种染色质互作检测方法（Hi-C、Micro-C、SPRITE、GAM、ChIA-PET、PLAC-seq），评估它们检测区室、TAD边界和染色质环的能力。随后，整合Hi-C、DamID和TSA-seq数据，利用SPIN模型注释基因组空间状态。进一步构建单细胞3D基因组结构模型，分析基因的核内微环境。最后，利用单细胞Hi-C和聚合物模型探索细胞间变异性，并训练深度学习模型预测序列变异对折叠的影响。

图 1　Fig. 1。图1展示了研究设计：左上示意两类互补的基因组折叠检测方法（染色质互作和核体距离），右上比较不同方法在检测区室到环等特征的能力，底部列出整合的多模态数据集。整体框架强调多方法基准化和数据整合。

核心亮点

亮点 1｜多方法基准揭示互补优势

系统比较Hi-C、Micro-C、SPRITE、GAM、ChIA-PET和PLAC-seq发现：SPRITE和Hi-C检测区室化最强，Micro-C动态范围最大，ChIA-PET和PLAC-seq富集转录相关环。各方法检测到的环锚定染色质状态不同，CTCF ChIA-PET主要捕获绝缘子环，Pol II ChIA-PET和H3K4me3 PLAC-seq富集活性启动子环。联合分析得到H1细胞141,365个环，远超单一方法，表明不同方法互补。

图 2　Fig. 2。图2比较了H1和HFFc6细胞中不同方法生成的接触图、区室特征向量和绝缘分数。热图显示各方法捕获相似的区室化模式，但动态范围不同。绝缘分数聚类显示方法间高度相关，SPRITE和GAM相关性较低。

亮点 2｜SPIN状态注释空间区室

整合TSA-seq、DamID和Hi-C数据，利用改进的SPIN模型在H1和HFFc6细胞中鉴定出9种空间状态，从核斑点到核纤层。这些状态与组蛋白修饰、复制时序和染色质相关RNA特征显著关联。例如，活性组蛋白标记（H3K27ac、H3K4me3）在核斑点状态富集，而H3K9me3在核纤层状态富集。SPIN状态为基因组提供了线性空间注释，便于整合其他功能数据。

图 3　Fig. 3。图3展示染色质环的分类。a为UMAP投影的构建流程，b显示H1细胞中141,365个环聚为6类，c显示各类环锚定的染色质状态富集模式，d显示环大小分布，e显示不同方法检测到的环类型比例。

亮点 3｜3D模型揭示基因微环境

利用整合基因组建模（IGM）构建了H1和HFFc6细胞的1000个单细胞3D基因组结构。分析发现，高表达基因分为两类：I类基因靠近核斑点，短且位于基因密集区；II类基因远离核斑点，长且位于基因稀疏区，但具有高CpG密度和活性启动子标记。II类基因更多依赖超长距离增强子互作。细胞类型间差异表达基因的核定位显著不同，如POU3F1在H1中位于核内部，在HFFc6中靠近核边缘。

图 4　Fig. 4。图4展示SPIN状态与功能基因组数据的关联。a热图显示组蛋白修饰、复制时序和caRNA在SPIN状态中的富集，b显示染色体1片段在H1和HFFc6中的平均径向位置，c显示单细胞3D结构示例，d-g展示基因微环境分类。

亮点 4｜单细胞折叠变异性

利用WTC-11多能干细胞单细胞Hi-C数据，结合Higashi插补和SBS聚合物模型，发现单细胞中TAD样结构边界存在变异性。整合分析显示，同一细胞中位于相同区室或TAD样结构域的位点更可能形成强染色质环。例如，RABGAP1L附近的环在环锚点位于相同A区室的单细胞中强度更高。表明环、TAD和区室在单细胞水平上协调变化。

图 5　Fig. 5。图5展示单细胞3D基因组分析。a比较Higashi插补和SBS模型预测的接触图与bulk Hi-C的一致性，b显示三个单细胞的3D模型和接触图，c显示环强度与区室/TAD状态的关系。

亮点 5｜增强子互作与转录广度

分析19,693个蛋白编码基因发现，与更多远端增强子互作的基因转录水平更高。管家基因尤其倾向于与多个增强子互作，且在不同细胞类型中使用的增强子集合高度动态（80.3%的增强子-启动子对为细胞类型特异）。在116个人类组织/细胞系中，增强子互作数量与基因表达广度正相关。此外，位于核纤层相关域内的基因若与增强子互作，则更可能表达，提示局部环出核纤层环境以建立调控互作。

生信可带走

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

• 方法栈： Hi-C、Micro-C、SPRITE、GAM、ChIA-PET、PLAC-seq、DamID、TSA-seq、Repli-seq、iMARGI、单细胞Hi-C、3D建模
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事三维基因组、表观遗传学、单细胞多组学或基因调控研究的科研人员。

带走一句

没有单一方法能全面捕获4D核组特征，应根据研究问题选择互补方法并整合多模态数据；SPIN状态和3D模型为功能注释提供了有力框架。