《International Journal of Molecular Sciences》：Trichostatin A-Induced H3K9ac Hyperacetylation Selectively Enriches Promoters in Adipogenic Human Periodontal Ligament Stem Cells
编辑推荐：
脂肪形成（Adipogenesis）受到与染色质可及性相关的组蛋白修饰调控，这些修饰包括组蛋白H3第9位赖氨酸乙酰化（H3K9ac），这是一种与转录激活相关的标记。其动态变化由组蛋白乙酰转移酶和组蛋白去乙酰化酶（HDACs）调控。尽管Trichostatin
脂肪形成（Adipogenesis）受到与染色质可及性相关的组蛋白修饰调控，这些修饰包括组蛋白H3第9位赖氨酸乙酰化（H3K9ac），这是一种与转录激活相关的标记。其动态变化由组蛋白乙酰转移酶和组蛋白去乙酰化酶（HDACs）调控。尽管Trichostatin A（TSA）能够诱导全局性的高乙酰化，但其在人细胞成脂分化过程中与调控元件处H3K9ac差异性富集的关联尚不清楚。因此，本研究描绘了TSA诱导的人牙周膜干细胞（hPDLSCs）中H3K9ac富集的基因组分布特征。为此，研究人员对hPDLSCs进行TSA预处理（100 nM，72 h）后，再诱导其成脂分化14天。通过免疫荧光评估细胞形态、脂滴半径以及H3K9ac和H3K9me2的核内分布，并通过ChIP-seq评估H3K9ac的基因组分布。采用非配对双尾Student's t检验（p < 0.05）和MACS3峰值检测（q = 0.001）评估统计学显著性。TSA使脂滴半径减小了1.69倍，并使核内H3K9ac信号增强了10.17倍，同时使H3K9me2信号降低了1.38倍。此外，与对照组（364个）相比，TSA处理组有752个启动子表现出更强的H3K9ac富集，这些启动子关联的基因包括成脂相关基因如PPARGC1B和AGPAT2，以及染色质重塑基因SETD5和KDM2A。总体而言，这些发现表明TSA改变了成脂表型，并在与代谢和染色质重塑程序相关的启动子处增加了H3K9ac的富集。
研究背景与目的
脂肪形成（Adipogenesis）是一个受精密调控的细胞分化过程，对于脂肪组织的稳态和可塑性至关重要，其失调与肥胖和2型糖尿病密切相关。这一过程涉及基因表达的协同改变，其中表观遗传机制，特别是组蛋白翻译后修饰中的乙酰化和去乙酰化，发挥着核心作用。组蛋白H3第9位赖氨酸乙酰化（H3K9ac）作为关键的转录激活标记，与开放的染色质状态直接相关，而其拮抗标记H3K9me2/3则与转录抑制相关。两者之间的平衡深刻影响着成脂相关基因的转录状态。尽管已知H3K9ac在成脂分化中扮演重要角色，但现有证据多来自小鼠模型，其在人类成体干细胞中的基因组分布特征仍缺乏深入的表征。人牙周膜干细胞（hPDLSCs）因其较高的成脂分化潜能，成为从表观遗传学角度研究脂肪形成的理想模型。前期研究已证实，组蛋白去乙酰化酶（HDAC）抑制剂Trichostatin A（TSA）预处理能增强hPDLSCs的成脂分化，并伴随H3K9ac蛋白水平的持续升高，在第14天达到峰值。然而，这种全局性的核内H3K9ac增加是否意味着该修饰在特定基因组区域（尤其是成脂基因启动子）的选择性富集，还是仅仅表现为无区域特异性的均匀增加，这一问题尚未明确。因此，本研究旨在利用TSA作为H3K9ac高乙酰化模型，通过ChIP-seq技术系统表征TSA预处理后hPDLSCs在成脂分化第14天时H3K9ac的全基因组分布特征，验证TSA诱导的高乙酰化是否能导致H3K9ac在成脂相关基因启动子处的选择性富集，从而鉴定出与成脂表型变化相关的候选调控区域和基因。
主要关键技术方法
研究人员采用了以下关键技术方法：首先，通过免疫荧光（Immunofluorescence）结合共聚焦显微镜定量分析了H3K9ac和H3K9me2的核内分布及荧光强度；其次，运用染色质免疫沉淀测序（ChIP-seq）技术，在单一样本（来源于单个供体的hPDLSCs细胞系）中对H3K9ac修饰的全基因组分布进行初始图谱绘制；最后，通过生物信息学分析流程，包括使用Bowtie2比对至人类参考基因组（hg38）、MACS3进行峰值检测（Peak calling）、ChIPseeker进行基因组注释、GREAT进行基因邻近关联分析，以及MACS3 bdgdiff进行差异富集分析，并结合ShinyGO和STRING进行功能富集与蛋白质互作网络分析。
研究结果
2.1. TSA对hPDLSC成脂分化的影响
通过明场显微镜观察和形态计量学分析，研究人员发现TSA预处理后的hPDLSCs在成脂诱导第14天，其脂滴形态与对照组存在显著差异。TSA处理组的脂滴平均半径（2.04 ± 0.49 μm）较对照组（3.44 ± 1.98 μm）更小，但每个细胞内的脂滴数量更多（范围41-241 vs. 35-126），呈现出聚集的小脂滴表型。
2.2. TSA改变hPDLSC成脂分化过程中H3K9ac和H3K9me2的核内分布
免疫荧光结果显示，TSA处理72小时后，H3K9ac的核内荧光信号即变得更强且更均一。在成脂分化第14天，TSA处理组的H3K9ac核内平均荧光强度（14.24 AFU）远高于对照组（1.40 AFU）。相反，H3K9me2的信号在TSA处理后减弱，其平均荧光强度从对照组的0.44 AFU降至0.32 AFU。这表明TSA诱导了持久的表观遗传重塑，伴随着H3K9ac的全局性富集和H3K9me2的减少。
2.3. 质量指标与比对验证
通过读取去重、比对率、文库复杂度等指标评估了ChIP-seq文库质量。经下采样均衡测序深度后，对照组和TSA组分别获得约2590万和2600万的比对读数，并分别识别出24,749和27,283个峰值，FRiP分数分别为19.62%和23.90%，表明数据质量可靠。
2.4. H3K9ac富集区域的基因组分布与基因关联
基因组注释显示，在两个实验组中，H3K9ac富集区域主要定位于近端启动子（距转录起始位点[TSS]≤1 kb），TSA处理组中启动子占比更高（55.10% vs. 49.74%）。H3K9ac信号在TSS周围高度集中，且TSA组信号更强。GREAT分析表明，绝大多数富集区域都与至少一个基因相关联，其中大部分关联距离TSS在5-500 kb范围内。
2.5. 差异富集的H3K9ac启动子及其关联基因的功能特征
差异富集分析鉴定出，与对照组（364个）相比，TSA处理组有更多的启动子（752个）显示出更高的H3K9ac富集。这752个启动子中有311个与蛋白编码基因相关。基于hPDLSCs特异性背景的GO功能富集分析显示，这些基因显著富集于氧化磷酸化、细胞呼吸、质子跨膜转运和ATP合成等代谢过程，这与脂肪细胞成熟所需的能量需求一致。STRING蛋白互作网络分析揭示了包含MAPK1、PPARGC1B、WNT10B、SLC2A4以及多个线粒体组分在内的模块化拓扑结构。
2.6. 成脂、能量代谢和染色质重塑相关基因的局部H3K9ac富集构型
研究人员选取了14个代表性基因座进行局部H3K9ac分布的可视化分析，并根据信号拓扑结构将其归纳为三种描述性架构：（i）主导性局灶峰架构，如AAMDC、SETD5、ELOVL7等，其富集信号尖锐地限制在基因5'端；（ii）分段富集架构，如PPARGC1B、WNT10B、SLC2A4、VSTM2A，表现为多个离散的、连续的富集区域沿基因座分布；（iii）广泛富集结构域，如LMF1、CTTN、MBD3、AGPAT2，其特征是H3K9ac信号连续覆盖了启动子、CpG岛和广泛的基因内区。所有选定基因座的启动子区域均检测到正的log2倍数变化值，确认了其定量上的差异富集。
讨论与结论
讨论部分总结：本研究表明，TSA预处理不仅改变了hPDLSCs成脂分化过程中的H3K9ac核水平，还导致了其全基因组分布的靶向性重组，而非简单的全局性增加。这种区域特异性的H3K9ac富集，特别是在752个启动子的选择性增加，支持了TSA诱导的高乙酰化与成脂定向相关的假说。第14天的基因组图谱捕捉到了活跃的染色质状态，而非终末成熟期的稳定状态。研究发现H3K9ac主要富集于近端启动子，与转录活性染色质标志一致。不同基因座呈现出的局灶峰、分段富集和广泛结构域等多种局部构型，反映了基因调控的复杂性，可能与CpG岛密度、顺式调控元件多样性等因素有关。TSA处理导致的更小、更多数量的脂滴表型，并非发育停滞的标志，而是与快速从头脂质合成和高代谢需求相一致，这与ChIP-seq数据中发现的氧化磷酸化、线粒体生物合成等过程的富集相互印证。研究还重点讨论了VSTM2A、PPARGC1B、AGPAT2等候选基因在早期成脂定向和脂质代谢中的潜在作用。同时，研究人员也坦诚指出了研究的局限性，包括ChIP-seq仅使用单一生学样本、未进行spike-in归一化、细胞来源于单一供体以及TSA的非特异性作用等，强调这些结果为后续的功能验证提供了初步的表观遗传框架。
研究结论翻译：TSA预处理与成脂表型的改变以及H3K9ac和H3K9me2的丰度变化相关联，促进了hPDLSCs中与成脂程序相连的启动子处H3K9ac的富集。总之，这些发现支持H3K9ac是与成脂启动相关的一个关键表观遗传特征，并为识别调控区域和候选基因提供了一个初步的基因组框架，可为再生医学和代谢研究领域的后续功能验证及靶向表观遗传学研究提供信息。
