多组学整合揭示栀子（Gardenia jasminoides J.Ellis）叶片染色质可及性的基因组特征

《Horticulturae》：Integration of Multi-Omics Reveals Genomic Features of Chromatin Accessibility in Gardenia jasminoides J.Ellis Leaves

　　

染色质可及性是顺式调控元件的重要特征，它塑造了植物的基因调控；然而，在传统药用植物如栀子（Gardenia jasminoides J.Ellis，G. jasminoides）中，对它的研究较少。本研究应用转座酶可及染色质测序（Assay for Transposase-Accessible Chromatin using sequencing，ATAC-seq）技术来表征G. jasminoides叶片中可及染色质区域（accessible chromatin regions，ACRs）的全基因组分布和功能特征。本研究共鉴定出26,461个ACRs和13,625个相关基因，研究结果表明ACRs的开放染色质状态与其相关基因的表达水平呈正相关。还发现这些ACRs富集了众多保守的转录因子结合基序（transcription factor binding motifs，TF motifs）。整合两个组蛋白修饰数据集进一步证明，ACRs与激活型组蛋白修饰（包括H₃K4me³和H₃K27ac）密切相关，这些修饰共同调节靶基因的转录。此外，荧光素酶（luciferase，LUC）报告实验验证了两个候选ACRs的转录激活活性。综上所述，这些发现阐明了G. jasminoides中ACRs的基因组特征，并为该物种后续的基因调控分析提供了关键的基因组资源。

染色质可及性是顺式调控元件的重要特征，在植物基因调控中发挥核心作用，然而在传统药用植物如栀子（Gardenia jasminoides J.Ellis）中相关研究仍較为匮乏。栀子作为茜草科（Rubiaceae）常绿灌木，其干燥成熟果实具有镇痛、止血、抗炎、保肝利胆等药用功效，同时兼具观赏和天然色素来源等多重价值，然而其转录调控网络，尤其是转录因子（transcription factor，TF）与非编码调控元件的协同调控机制仍不清楚。为填补该空白，研究人员提出核心问题：可及染色质区域（accessible chromatin regions，ACRs）如何参与栀子基因的转录调控？该研究发表于《Horticulturae》，采用转座酶可及染色质测序（ATAC-seq）技术，结合转录组测序（RNA-seq）、组蛋白染色质免疫沉淀测序（ChIP-seq）以及双荧光素酶（LUC）报告实验，系统表征了栀子叶片中ACRs的基因组分布、表观基因组特征及其转录调控功能。

研究叶片样本采自中国科学院庐山植物园（江西省九江市）的3至4年生栀子（LS016）植株，每个样本混合3片成熟叶，设置两个独立生物学重复。主要技术方法包括：ATAC-seq鉴定开放染色质区域；RNA-seq量化基因表达水平；组蛋白ChIP-seq检测H3K4me3和H3K27ac两种激活型修饰；双荧光素酶报告实验验证候选ACR的转录激活活性。数据处理利用Bowtie2、MACS2、HOMER、FIMO等生物信息学工具完成峰值 calling、注释及基序分析。

3.1 通过ATAC-seq全基因组鉴定G. jasminoides中的ACRs。研究人员对栀子叶片进行ATAC-seq，获得71.4 M高质量reads，其中71.6%以上 uniquely mapping 至参考基因组。共鉴定出26,461个在两个生物学重复中稳定存在的ACRs，关联13,625个蛋白编码基因，分布遍及所有11条染色体。基因组注释显示，ACRs主要位于基因间区（46.62%）和启动子区（25.66%），其余分布于转录终止区（10.51%）、内含子（9.00%）和外显子（8.21%）。ATAC-seq信号在转录起始位点（transcription start site，TSS）附近显著富集，且其强度与基因表达水平（FPKM）呈正相关，表明开放染色质状态与转录活性密切相关。

3.2 ACRs中保守TF motif富集分析。利用PlantTFDB数据库中的已知转录因子结合基序（TF motifs）进行映射，发现ACRs显著富集保守TF motifs，排名前10的家族包括BBR-BPC、MIKC_MADS、AP2、Dof、GATA、ERF、GRAS、LBD、Nin-like和C2H2。进一步构建潜在转录因子调控网络，其中8,899个基因同时含有MYB、MIKC_MADS和C2H2三类motif，提示ACRs可能通过多种TF协同调控下游靶基因。

3.3 侧翼ACRs的组蛋白修饰特征。组蛋白ChIP-seq在叶片中鉴定出18,617个H3K27ac富集峰和28,856个H3K4me3富集峰，两种修饰信号均富集于TSS附近且与基因表达水平正相关。通过k-means聚类将ACRs分为两个主要簇，第一簇以H3K4me3修饰为主，第二簇两种修饰信号均较低。共发现5,652个同时携带H3K4me3和H3K27ac双修饰的ACRs，其相关基因表达水平显著高于仅单修饰或无修饰的ACRs。此外，鉴定出13,063个远端ACRs（distal ACRs，dACRs），其中81.11%缺乏这两种激活型修饰，但双修饰的dACRs与近端ACRs（proximal ACRs，pACRs）共同关联的基因表达水平最高，提示dACRs可能具有增强子样功能。

3.4 候选ACRs的瞬时转染功能验证。基于ACRs、组蛋白修饰及基因表达特征，选取ACR3272和ACR19689两个候选区域，分别克隆至mini35S启动子上游的LUC报告载体中，通过农杆菌介导转化烟草（Nicotiana benthamiana）叶片。双荧光素酶报告实验结果显示，两个候选ACR均显著增强下游荧光素酶活性，与阴性对照（mini35S空载体）相比差异具有统计学意义（p < 0.001），证实了ACRs具有转录激活能力。

讨论部分总结了该研究的重要意义：建立了适用于木本非模式植物的ATAC-seq实验流程，整合多组学数据构建了首个栀子全基因组染色质可及性图谱，系统解析了ACRs的基因组定位、表观遗传特征及其转录调控功能。研究还指出，非编码顺式调控序列在栀子基因调控中发挥关键作用，远端ACRs可能作为增强子样元件协同近端ACRs调控靶基因表达。同时，研究存在局限性：当前数据仅来源于叶片，染色质可及性具有组织特异性；仅检测了两种激活型组蛋白修饰；缺乏跨物种比较分析；仅研究了正常生长条件下的静态染色质状态。未来需扩展至花、果实、根等多组织，检测抑制性修饰（如H₃K27me³）和DNA甲基化（WGBS），并开展与咖啡、金鸡纳等近缘种的比较表观基因组学研究，以及逆境胁迫下的动态染色质重塑分析。

研究结论部分翻译：在本研究中，研究人员使用ATAC-seq和组蛋白ChIP-seq方法系统地表征了G. jasminoides中的全基因组ACRs。此外，研究结果表明ACRs在G. jasminoides基因的转录调控中起着至关重要的作用。总的来说，这项研究为鉴定重要功能基因、阐明转录调控机制以及推进G. jasminoides的分子育种工作提供了必要的基因组资源和理论基础。