《Horticulture Research》：Distal cis-regulatory elements in 3′ UTRs orchestrate gene expression via single-gene domain architecture in Brassica rapa
编辑推荐：
顺式调控元件（cis-regulatory elements, CREs）是精确基因调控的基础。尽管位于基因间区的远端CREs已被广泛研究，但驻留在3′非翻译区（3′ untranslated regions, 3′ UTRs）中的那些CREs的功能景观仍很大
顺式调控元件（cis-regulatory elements, CREs）是精确基因调控的基础。尽管位于基因间区的远端CREs已被广泛研究，但驻留在3′非翻译区（3′ untranslated regions, 3′ UTRs）中的那些CREs的功能景观仍很大程度上未被探索。在本研究中，研究人员生成了白菜（Brassica rapa）叶片在两个发育阶段（莲座期和生殖期）的综合表观基因组图谱，整合了ATAC-seq、针对五种组蛋白修饰（H3K4me3、H3K27ac、H3K27me3、H3K36me3、H3K9me2）的CUT&Tag数据以及BL-Hi-C数据。研究人员鉴定了位于3′ UTRs内的6362个CREs，定义为转录终止位点（transcription termination sites, TTS）上下游0.5 Kb范围内的可及染色质区域。这些3′ UTR内的CREs根据组蛋白修饰被功能分类为：H3K27ac标记的、H3K27me3标记的和二价（bivalent）CREs。利用分辨率高达100 bp的BL-Hi-C，研究人员发现这些3′ UTR内的CREs促进了单基因域（single-gene domains）的形成。发育阶段间的比较分析揭示，这些远端CREs通过调节染色质相互作用来协调阶段特异性的基因表达，这一机制以开花调控因子BrFLC2为例得到阐明。双荧光素酶报告实验验证了BrFLC2远端CRE的转录激活活性。总之，研究人员的发现确立了3′ UTR嵌入的远端CREs作为关键调控元件，通过白菜中的单基因域架构塑造基因表达。
研究背景与科学问题
顺式调控元件（cis-regulatory elements, CREs）是实现基因精确时空表达调控的基因组基础。以往对CREs的研究多集中于启动子近端区域或基因间区的远端增强子/沉默子，而对位于3′非翻译区（3′ untranslated regions, 3′ UTRs）内部的远端CREs的功能认知严重不足。白菜（Brassica rapa）是全球重要的蔬菜作物，也是进化表观基因组学的优良模式物种。然而，由于此前常规Hi-C技术分辨率有限（通常为2-5 Kb），无法解析单个基因尺度的精细三维染色质结构如单基因域（single-gene domains），因此白菜中是否存在3′ UTR嵌入的远端CREs、它们如何参与更高阶染色质架构的形成以及这种架构如何在发育过程中动态调控基因表达，均是未解之谜。为此，研究人员开展了本项多组学整合研究。
主要研究方法
研究人员采用了以下关键技术方法：（1）ATAC-seq用于全基因组染色质可及性图谱绘制；（2）CUT&Tag技术检测五种组蛋白修饰（H3K4me3、H3K27ac、H3K27me3、H3K36me3、H3K9me2）；（3）BL-Hi-C技术（桥连连接子Hi-C，使用HpyCH4V和AluI双酶切）构建超高分辨率染色质构象捕获图谱，分辨率可达100 bp；（4）双荧光素酶报告实验进行候选CRE的功能验证。所有实验材料均取自白菜品种Chiifu的幼叶，分为莲座期（40天苗龄）和生殖期（低温处理30天后恢复室温10天）两个发育阶段。
研究结果
1. 白菜叶片表观基因组图谱构建
研究人员对莲座期和生殖期白菜叶片生成了12套全基因组表观遗传数据集，包括五种组蛋白修饰、染色质可及性和转录组数据，各两个生物学重复。所有重复间Pearson相关系数超过0.97。通过多元隐马尔可夫模型（Hidden Markov Model, HMM）将基因组划分为12种功能不同的染色质状态，其中State 8由H3K27ac主导定义为远端转录激活元件，State 10主要由H3K27me3标记代表远端转录抑制元件，State 9同时富集H3K27ac和H3K27me3为二价染色质状态。
2. 远端CREs的鉴定与特征分析
基于莲座期ATAC-seq数据，共鉴定出14498个近端CREs和6362个远端CREs（定义为TTS上下游0.5 Kb内的可及染色质区域）。从头基序分析显示近端和远端CREs均富集相似的转录因子结合基序。表观修饰谱对比表明，近端CREs主要发挥启动子功能，而远端CREs兼具增强子或沉默子潜能。
3. 远端CREs的功能分类
整合H3K27ac和H3K27me3数据将远端CREs分为三类：1929个H3K27ac标记的CREs、609个H3K27me3标记的CREs和2256个二价CREs。基因表达分析显示，H3K27ac标记CREs关联的基因表达水平显著高于二价状态，而H3K27me3标记CREs关联的基因表达则显著更低。GO富集分析揭示H3K27ac标记基因富集于刺激响应通路，H3K27me3标记基因富集于生物合成调控通路。
4. 精细三维结构揭示远端CREs参与单基因域形成
BL-Hi-C数据产生约290百万有效接触，染色体内部与染色体间接触比为4.4±1，质量优异。与常规Hi-C相比，BL-Hi-C在500 bp、200 bp甚至100 bp分辨率下均能清晰解析单基因域结构。通过按基因边界可及性将基因分为四类，堆积接触矩阵可视化显示两端均有可及CREs的基因其基因内互作强度最高（平均653），而两端均无的基因最低（平均321）。染色质可及性与基因内互作强度呈正相关（Spearman R=0.77），基因表达水平与基因内互作强度同样正相关（Spearman R=0.77）。叶片极性发育基因BrKAN2.3作为典型案例展示了清晰的单基因域结构。
5. 远端CREs通过单基因域调控基因表达
发育阶段比较分析识别出1447个增强和2195个减弱的远端可及染色质区域（distal accessible chromatin regions, dACRs）。增强型dACRs关联的基因表达变化倍数显著更高，且活跃组蛋白修饰（H3K4me3、H3K27ac、H3K36me3）显著富集，而抑制性修饰H3K27me3在减弱型dACRs中更丰富。以开花调控基因BrFLC2为例，莲座期该基因形成强健的单基因域，伴随活跃组蛋白修饰和高表达；进入生殖期后域结构减弱，涉及远端CRE的互作从109次降至50次，伴随表达下降和H3K27me3富集增加。双荧光素酶报告实验证实BrFLC2的远端CRE具有转录激活活性。全基因组层面，上调基因的远端CRE互作频率变化显著高于下调基因。
讨论与结论
讨论部分指出，本研究首次在白菜中系统观察到单基因域结构并将其强度与远端CRE可及性和基因表达水平定量关联。BL-Hi-C技术以约260倍覆盖度实现了与Micro-C相当的百碱基分辨率，成为紧凑植物基因组精细染色质构象分析的高效工具。植物中增强子-启动子互作距离通常较短（几Kb到几十Kb），这与动物中超过100 Kb的长程互作形成对比，反映了植物高基因密度和紧凑基因组的特点。BrFLC2的案例揭示了发育过程中染色质架构、组蛋白修饰与基因表达的协同重编程。研究也指出了局限性：材料均为幼叶，无法研究茎尖分生组织特异性表观变化；异源体系验证需进一步通过CRISPR-Cas9等内源编辑手段确认；单基因域形成的分子机制尚待阐明。
研究结论翻译如下：总之，本研究呈现了白菜的综合表观基因组景观，鉴定了数千个先前未表征的位于3′ UTRs内的远端CREs，并揭示了它们通过单基因域架构协调基因表达的作用。这些发现确立了3′ UTR嵌入的远端CREs作为关键的调控元件，并为未来旨在精细调控芸薹属作物基因表达的功能研究提供了宝贵资源。
