IF31.7｜瞬时HDAC抑制诱导基因表达与3D基因组折叠的细胞记忆

一次短暂的组蛋白去乙酰化酶抑制，能在胚胎干细胞中留下持久的转录与染色质构象记忆。

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

研究背景

细胞身份由基因调控和表观遗传机制共同塑造，其中染色质状态存储维持转录程序所需的信息。胚胎干细胞具有独特的许可性染色质，激活与抑制性构象常共存，尤其在发育基因周围形成二价或预备状态。这些状态需要精细平衡，既能保持基因低表达，又能为未来激活做好准备。细胞记忆指短暂刺激转化为长期基因表达变化的能力，是发育过程中的重要概念，但不同表观遗传层面间的交互及其对记忆的贡献因功能冗余和实验手段限制而难以研究。三维基因组折叠正逐渐被认为是细胞身份的关键调控者，但其是否参与细胞记忆尚不清楚。

科学问题

卡在哪里： 此前研究难以将基因调控与细胞记忆分离，且缺乏实验手段在急性扰动后同时追踪转录、组蛋白修饰和3D基因组折叠的动态。三维基因组是否以及如何存储短暂刺激的信息，仍是未解之谜。

本文要回答： 本文旨在探究短暂HDAC抑制如何影响转录、组蛋白修饰和基因组折叠，并确定三维基因组组织是否参与细胞记忆的形成。

技术路线

作者主要用了：RNA-seq、ChIP–seq、ATAC–seq、ultradeep Micro-C、生物物理建模。

作者用TSA脉冲处理小鼠胚胎干细胞4小时，诱导全局组蛋白高乙酰化，同时通过优化条件避免细胞周期和乙酰化组的广泛扰动。随后利用多组学手段（RNA-seq、多种组蛋白修饰ChIP–seq、ATAC–seq）刻画染色质状态和转录变化，并用超高深度Micro-C绘制三维基因组互作图谱。为验证记忆效应，作者在TSA洗脱后24小时再次检测，并进一步进行第二轮TSA脉冲以观察重复暴露的后果。此外，通过生物物理聚合物模拟解释染色质刚度变化对区室互作的影响，并在PRC1亚基敲除细胞中验证Polycomb介导的环在记忆中的作用。

图 1　Fig. 1。图1展示TSA处理4小时后多组学变化：染色体18上H3K27ac信号整体升高；典型高乙酰化区域显示H3K27ac和ATAC信号同时增加；差异峰附近基因富集发育相关GO术语；火山图显示大量基因上调（神经发育）和下调（多能性维持）。

核心亮点

亮点 1｜TSA诱导全局H3K27ac与发育基因激活

4小时TSA处理导致H3K27ac在全基因组范围显著升高，活性组蛋白修饰（H3K4me1、H3K4me3）总体增加，抑制性修饰（H3K27me3、H2AK119ub）变化较小。差异ChIP–seq峰的基因注释显示发育过程富集，多能性相关基因下调。RNA-seq证实上调基因与神经谱系发育相关，下调基因与干细胞维持相关。这些变化表明急性HDAC抑制推动细胞退出多能性状态。

图 2　Fig. 2。图2显示Micro-C揭示的3D基因组变化：A–A区室互作减弱，B–B互作增强；H3K27ac在A区室增加更多；差异表达基因TSS主要位于A区室；生物物理模型示意通过增加A区室染色质刚度可重现trans互作增加和区室互作变化；环分析显示CTCF环减弱，非CTCF环增强。

亮点 2｜TSA重塑3D基因组：A区室减弱，非CTCF环增强

超高深度Micro-C显示TSA处理显著增加染色体间（trans）互作，减少顺式（cis）互作，A–A区室互作减弱而B–B互作增强，类似分化细胞特征。生物物理模拟表明，A区室染色质刚度增加大于B区室可解释trans互作比例变化，并预测A区室向核外周移动，免疫荧光证实H3K27ac焦点靠近核膜。环分析发现CTCF介导的环普遍减弱，而非CTCF环（携带活性或抑制性组蛋白标记）增强，差异环位点富集发育基因。

图 3　Fig. 3。图3关联表观基因组与转录变化：上调基因TSS附近H3K27ac和H3K4me1信号增加，下调基因TSS附近活性标记减少；上调基因启动子与最近H3K4me1峰距离更近；环锚点频率显示上调基因TSS附近环数量未显著增加，但下调基因TSS附近环增强。

亮点 3｜基因上调与增强子激活相关，下调与抑制性环相关

上调基因TSS附近H3K27ac和H3K4me1信号增加，且与预存在的启动子-增强子（E–P）接触相关，但E–P接触频率本身未显著改变。下调基因TSS附近活性标记减少，但未获得抑制性组蛋白修饰；相反，下调基因启动子周围环强度增加，部分形成新的H3K9me3相关环。这表明上调可能由增强子过度激活驱动，而下调与抑制性染色质环形成有关。

图 4　Fig. 4。图4展示恢复实验：TSA洗脱24小时后H3K27ac、H3K4me1和ATAC信号基本恢复；H3K27ac峰重新富集于TSS；火山图显示仅少数基因仍失调；胃泌素样体实验表明TSA脉冲后转录组大部分恢复，但Sox2阳性区域扩大。

亮点 4｜恢复后部分基因和构象保持记忆

TSA洗脱24小时后，组蛋白乙酰化和大多数转录变化恢复，但仍有164个基因维持差异表达。更长时间追踪显示差异基因数量随时间增加。3D基因组折叠未完全恢复：cis–trans比例部分恢复，A–A互作在cis中仍减少，B–B互作保持增强。局部区室切换（如F11–Klkb1–Cyp4v3位点）在基因表达恢复后仍维持A区室身份。增强的环在恢复后仍保持强度。这些结果表明基因组构象携带了过去状态的记忆。

图 5　Fig. 5。图5显示恢复后3D基因组仍保留记忆：鞍图显示cis中A–A互作未完全恢复，B–B互作保持增强；F11–Klkb1–Cyp4v3位点区室切换在恢复后仍维持；Pitx2位点显示不完全的构象恢复；差异环的pile-up显示TSA中增强的环在恢复后仍保持强度。

亮点 5｜重复TSA暴露加剧转录记忆，Polycomb环介导持续下调

第二轮TSA处理后，恢复更不完全，767个基因持续失调，且与发育过程相关。分层分析显示部分基因在两次处理中逐渐加重失调，提示细胞记忆。持续下调的基因启动子区域富集H3K27me3，且Polycomb环在恢复后仍增强。在PCGF2/PCGF4敲除细胞中破坏Polycomb介导的空间聚类，使持续下调基因数量减少超过两倍，证明PRC1介导的染色质环是TSA诱导的持续基因下调所必需的。

生信可带走

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

• 方法栈： RNA-seq、ChIP–seq、ATAC–seq、ultradeep Micro-C、生物物理建模
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 关注表观遗传记忆、3D基因组和干细胞状态转换的研究者，尤其是做多组学整合分析的生信人员。

带走一句

3D基因组折叠不仅是基因调控的读出，还能存储短暂扰动的记忆；分析表观遗传扰动时，应同时关注染色质构象的持久变化。