《Frontiers in Plant Science》：Plant-specific function of H3K9me3 as a permissive chromatin mark during Arabidopsis seed germination
编辑推荐：
本研究通过整合表观基因组与转录组学分析，揭示了在拟南芥种子萌发过程中，传统抑制性组蛋白标记H3K9me3扮演了关键的植物特异性角色。研究发现，H3K4me3和H3K9me3在常染色质区域共定位，与基因活性正相关，且H3K9me3在种子从后熟状态向萌发刺激状态转换期间显著增加，与翻译、能量代谢和细胞分裂等关键基因的激活密切相关。休眠释放时，休眠相关基因(如DOG1)伴随H3K4me3和H3K9me3水平的靶向降低而下调。这挑战了动物中H3K9me3作为抑制性异染色质标记的传统观点，表明其在植物萌发中可能作为促进转录重编程的活性或许可性表观遗传标记。
种子萌发是植物生命周期中的一个关键发育转变，其精准调控对确保在不利环境中的生存至关重要。该过程与内部休眠状态紧密相连，后者在种子成熟过程中形成，并通过干燥贮藏（后熟）和层积处理逐渐解除。长期以来，赤霉素(GAs)和脱落酸(ABA)等激素及其下游转录网络的互作被认为是调控休眠与萌发的核心机制。然而，表观遗传机制，特别是组蛋白甲基化在这些状态转变中的作用尚不十分清楚。
为了探究与萌发潜力相关的表观遗传景观，研究团队聚焦于三种组蛋白甲基化标记：通常代表激活标记的H3K4me3、代表抑制标记的H3K27me3，以及功能尚不明确的H3K9me3。研究在两种拟南芥生态型——代表浅休眠的哥伦比亚型(Col)和代表深休眠的海角佛得角型(Cvi)中，对新鲜采收(FH)、后熟(AR)和萌发刺激(GS)三种种子状态进行了全基因组染色质免疫沉淀测序(ChIP-seq)分析，并整合了已发表的转录组(RNA-seq)数据。
全基因组水平H3K4me3、H3K9me3和H3K27me3在休眠、后熟及萌发刺激种子中的分布特征
ChIP-seq分析揭示了三种组蛋白甲基化标记的不同分布模式。H3K4me3主要富集在常染色质的基因区，H3K27me3则广泛分布于基因区和基因间区，符合其已知的转录抑制功能。与动物中H3K9me3作为组成型异染色质标记的经典角色不同，本研究发现拟南芥种子中的H3K9me3主要富集于常染色质区域，常与H3K4me3和转录活跃的基因共定位。全基因组分布图进一步证实，H3K4me3和H3K9me3在常染色质富集，而异染色质相关标记H3K9me2则局限于着丝粒周围异染色质。组蛋白甲基化荟萃分析显示，H3K4me3和H3K9me3显著富集在转录起始位点(TSS)附近，而H3K27me3则在基因体上广泛沉积。对非转座元件(TE)基因、TE基因和TE片段的分析表明，H3K9me3主要与非TE基因关联，在TE相关序列中水平显著较低，这支持了H3K9me3在植物中主要与活跃转录的编码基因相关的观点。
组蛋白甲基化标记的组合发生
对三种组蛋白标记在三种种子状态下的全基因组相关性分析显示，H3K4me3与H3K9me3在所有状态下均呈现一致的正相关，而H3K4me3与H3K27me3无显著相关。在基因水平上，最常见的类别是同时被H3K4me3和H3K9me3标记的基因，其次是仅被H3K4me3标记的基因。大多数H3K9me3峰与H3K4me3共定位。值得注意的是，在Col和Cvi中，同时被H3K4me3和H3K9me3标记的基因比例从AR到GS状态有所增加。这些发现暗示H3K4me3和H3K9me3在调节染色质结构和基因表达方面存在功能上的相互作用。
组蛋白甲基化对FH、AR和GS种子转录输出的影响
整合ChIP-seq和RNA-seq数据的分析表明，H3K4me3和H3K9me3水平与基因表达水平呈正相关，而H3K27me3水平则呈负相关。这表明植物中的H3K9me3可能以一种类似于H3K4me3的、转录许可或支持性的方式发挥作用。进一步分析不同组蛋白标记组合对基因表达的影响发现，同时被H3K4me3和H3K9me3标记的基因，其转录水平显著高于仅被其中一种标记的基因。相反，同时被H3K4me3/H3K27me3或H3K9me3/H3K27me3标记的基因，其转录水平则低于仅被H3K4me3或H3K9me3标记的基因。这说明H3K4me3和H3K9me3与转录激活相关并可能起协同作用，而H3K27me3的存在则抵消了它们的激活效应，从而有助于微调种子中的基因表达。
种子状态生理和发育转变过程中组蛋白甲基化的动态变化
对组蛋白甲基化水平在FH、AR和GS阶段间相关性的分析显示，在Col中，所有标记在FH和AR之间都表现出强相关性，而GS与更早种子状态（FH或AR）的相关性对于H3K4me3和H3K9me3有所下降，表明萌发期间活性染色质特征发生了重组。相比之下，H3K27me3在所有种子状态间保持了相对较强的相关性，表明这种抑制性标记从休眠状态到萌发具有更高的稳定性。Cvi中的总体相关性模式与Col相似。荟萃分析显示，H3K4me3和H3K27me3水平在两种生态型的所有种子状态下保持相对恒定，而H3K9me3则表现出更动态的行为：其在FH到AR转变期间水平维持，但在AR到GS转变期间显著增加。这表明H3K9me3的增加可能是拟南芥萌发启动相关的保守表观遗传特征。
差异峰分析进一步支持了上述发现。在FH到AR转变期间，H3K4me3和H3K9me3的差异峰数量相对有限，但在AR到GS转变期间显著增加。值得注意的是，在AR到GS转变期间，显示H3K9me3增加的区域数量远超过减少的区域。这表明H3K9me3在萌发相关的转录激活中可能扮演关键角色。与相关性分析一致，H3K27me3在AR和GS之间的差异峰数量较少，提示H3K27me3介导的转录抑制在萌发中可能作用较小。总之，全基因组组蛋白甲基化景观在FH到AR转变期间基本保持稳定，但在AR到GS转变期间，大量基因组区域的H3K9me3显著增加。这表明虽然染色质状态在休眠释放期间得以维持，但H3K9me3可能在促进萌发期间的转录重编程中扮演突出角色。
组蛋白甲基化动态与种子状态转变期间转录输出的关系
差异表达基因(DEG)分析显示，AR到GS转变期间的DEG数量远多于FH到AR转变，这与差异H3K4me3和H3K9me3峰的模式一致。将所有DEG汇集并通过k均值聚类分为五个簇(C1至C5)后，基因本体(GO)分析揭示，在Col和Cvi的C1和C2簇（在AR到GS转变期间下调的基因）中，富集了与水分剥夺、低氧和ABA信号等胁迫相关过程的GO术语。相反，在C4和C5簇（在AR到GS转变期间上调的基因）中，则富集了与翻译、脂肪酸生物合成、细胞分裂、分化和发育生长等相关的GO术语。
为了确定组蛋白甲基化在这些簇中的作用，研究量化了所有种子状态下这些簇内基因的H3K4me3、H3K9me3和H3K27me3水平。总体而言，表达动态模式与H3K4me3和H3K9me3水平的变化相关。值得注意的是，H3K9me3动态与Col和Cvi的C4和C5簇（基因在FH和AR状态表达相对较低，但在AR到GS转变期间被激活）的转录组动态存在显著相关性。相比之下，H3K4me3在此转变期间与基因表达变化的相关性较H3K9me3更为温和。对C4和C5簇基因在GS状态下的H3K4me3和H3K9me3沉积分析发现，最大比例的基因（Col中68%，Cvi中76%）同时被H3K4me3和H3K9me3标记。其中，相当一部分含有AR < GS差异H3K9me3峰。DEG与差异甲基化基因(DMG)的比较进一步支持了这些发现。值得注意的是，约20%在AR到GS转变期间上调的基因含有AR < GS差异H3K9me3峰。此外，在AR到GS转变期间，被差异H3K4me3或H3K9me3峰标记的基因的转录水平与H3K4me3和H3K9me3水平呈正相关。相比之下，在FH到AR转变期间识别的DMG则显示相对恒定或仅有轻微变化的转录水平。这些结果进一步支持了H3K9me3可能作为关键表观遗传调节因子，调节AR到GS转变期间的转录重编程，并可能促进萌发必需基因表达的观点。
萌发期间富含H3K9me3基因的生物学功能
对C4和C5簇基因根据H3K9me3状态进行分类后的GO分析显示，不同组别富集的GO术语各异，提示可能的功能分化。在Col中，具有AR < GS差异H3K9me3峰的基因主要与翻译和核糖体生物合成等GO术语相关。相反，具有非差异（共有）H3K9me3峰的基因则富集了与氧化磷酸化、有氧呼吸和ATP生物合成等能量代谢相关的GO术语。缺乏H3K9me3峰的基因组显示出与细胞周期相关的强烈富集。在Cvi中也观察到类似的趋势，携带AR < GS差异H3K9me3峰的基因同样主要富集于翻译相关GO术语。这些结果表明，功能多样的基因在萌发期间受到不同的表观遗传调控。
后熟期间组蛋白甲基化和转录输出的变化
尽管在干燥种子中活跃转录可能不会发生，但研究仍鉴定出大量FH到AR转变期间的DEG。其中，共同下调的基因包括编码种子贮藏蛋白、脂质生物合成酶的基因，以及与休眠或萌发相关的转录因子编码基因，如DOG1的反式激活因子bZIP67和萌发负调控因子PEI1/TZF6。组蛋白甲基化分析显示，大多数AR抑制基因被H3K4me3标记，一小部分被H3K9me3标记，仅有相对较小部分被H3K27me3标记。在AR状态，观察到H3K9me3峰数量显著减少。组蛋白甲基化水平的比较揭示了在Col和Cvi中，相当一部分AR抑制基因存在差异富集。总体而言，组蛋白甲基化的变化是适度的。这些结果表明，FH到AR转变期间的转录抑制与H3K4me3和H3K9me3标记的靶向减少有关。
共同上调的基因则与对非生物和生物胁迫的响应相关。与AR下调基因类似，大多数AR上调基因被H3K4me3或H3K9me3标记，而非H3K27me3。AR上调基因中的组蛋白甲基化变化总体上不如AR下调基因中明显。对随机选择的非差异表达基因的分析表明，绝大多数非DEG表现出稳定的组蛋白甲基化水平，没有显著偏向于增加或减少。这进一步支持了FH到AR转变期间的转录抑制与H3K4me3和H3K9me3标记的靶向减少相关，而非全局性组蛋白甲基化模式转变的结果。
结论
本研究揭示，传统上被视为抑制性组蛋白标记的H3K9me3，在拟南芥种子中与转录正相关，并在萌发期间显著增加。组蛋白甲基化谱在FH到AR期间基本保持稳定，但休眠释放伴随着一部分下调基因的H3K4me3和H3K9me3减少。这些发现支持了H3K9me3在植物中作为一种活性或许可性染色质标记，支持萌发期间转录重编程的特异性角色。更广泛地说，它们凸显了组蛋白修饰功能在界级水平上的潜在差异，拓展了我们对真核生物表观遗传多样性的理解。
