黑树莓果实成熟过程中DNA甲基化与转录调控动态互作的多组学解析

《Communications Biology》：Multi-omics analyses reveal dynamic interactions between DNA methylation and transcriptional regulation during black raspberry ripening

果实成熟是一个复杂的发育过程，涉及基因表达、代谢和表观遗传调控的协调变化。DNA甲基化动态如何影响黑树莓（Rubus occidentalis）果实成熟过程中的基因转录和代谢物积累仍知之甚少。本研究进行了发育序列的转录组、甲基化组和代谢组分析；分析结果揭示了黑树莓果实成熟过程中广泛的表观遗传重编程。果实组织整体DNA甲基化水平高于叶片，尤其在CHG和CHH序列背景中更为显著。成熟过程中，所有胞嘧啶背景下的局部甲基化在启动子区域逐步降低，而CHG和CHH甲基化在富含转座子的区域增加。两种主要的甲基化转变——启动子低甲基化和CHH高甲基化——调控参与成熟过程基因的转录变化。具有激活潜力的甲基结合转录因子可能促进与CHH高甲基化相关的转录激活。多组学整合揭示花青素积累与一致网络内生物合成基因和调控基因的表达平行协调。异染色质中非CG甲基化升高与介导染色质压缩的组蛋白变体转录增加同时发生，提示染色质重塑因子精细调节甲基转移酶的可及性。本研究强调了全基因组和位点特异性表观遗传重编程，并证明了黑树莓果实成熟过程中DNA甲基化与转录调控之间的协调互作。

研究背景方面，肉质果实成熟伴随颜色、质地、糖含量、香气和风味等生理变化，是受精密调控的发育过程。果实传统上分为呼吸跃变型和非呼吸跃变型，其中呼吸跃变型果实（如番茄）的成熟机制研究较为深入，而草莓、葡萄和柑橘等非呼吸跃变型果实的成熟机制尚不清楚。DNA甲基化作为重要的表观遗传修饰，参与调控分化、基因转录、转座子沉默和基因组印记等过程。在番茄中，转录因子CNR基因因启动子高甲基化而沉默，RIN作为MADS-box转录因子激活众多成熟相关基因，而DNA去甲基化酶SlDML2的缺失会延迟成熟；苹果和梨中MYB位点DNA甲基化改变与花青素生物合成变化相关。然而，黑树莓作为一种富含花青素、具有优异抗氧化能力的非呼吸跃变型果实，其DNA甲基化如何调控果实成熟过程中的基因表达尚缺乏系统研究。此前Rubus属成熟研究主要集中于转录组和代谢谱分析，鉴定了一些参与花青素生物合成的转录因子，但DNA甲基化在果实成熟中的调控作用基本未被探索。因此，研究人员开展了黑树莓果实发育系列的多组学分析，旨在揭示DNA甲基化动态与转录调控在果实成熟中的协调互作机制，这对于理解非呼吸跃变型果实成熟的表观遗传调控具有重要意义。

在技术方法方面，研究人员以黑树莓NC 349品系为材料，采集五个成熟阶段（绿、红、软、转色、黑）的果实及叶片，每个阶段设置生物学重复。关键技术方法包括：全基因组重亚硫酸盐测序（BS-Seq）进行单碱基分辨率DNA甲基化谱分析，覆盖CG、CHG和CHH三种序列背景；RNA-Seq转录组测序分析基因表达动态；高效液相色谱-二极管阵列检测（HPLC-DAD）进行靶向代谢组分析检测花青素类化合物；采用CGmapTools鉴定差异甲基化区域（DMRs）；运用DESeq2鉴定差异表达基因（DEGs）；通过加权基因共表达网络分析（WGCNA）构建基因调控网络；利用Chop-PCR和qRT-PCR进行候选基因验证。

研究结果方面，第一，黑树莓果实发育和成熟过程中的甲基化组特征：研究人员对五个发育阶段果实和叶片进行了单碱基分辨率DNA甲基化谱分析，发现果实组织整体DNA甲基化水平显著高于叶片，尤其在CHG和CHH背景中差异超过20%。随着成熟推进，CG和CHG对称位点甲基化逐渐下降，CHH甲基化位点数在ST2达到峰值后下降。第二，果实成熟过程中局部DNA甲基化模式改变：CG背景下低甲基化DMR多于高甲基化DMR，而非CG背景下高甲基化DMR占优势。大多数DMR发生在相同位点，低甲基化DMR长度随果实发育增加，表明去甲基化主要发生在相同位点。第三，果实成熟过程中特定位点的渐进性主动DNA去甲基化：ST5-vs-leaf低甲基化DMR在所有胞嘧啶背景中随成熟进展持续下降，RoDML3表达量最高且为番茄SlDML2的直接同源基因，提示其可能是成熟过程中局部去甲基化的主要执行者。第四，低甲基化DMR主要位于启动子和基因体区域：三个胞嘧啶背景的低甲基化DMR多数位于基因区，其中超过半数位于启动子；高甲基化DMR主要位于基因间区，非CG高甲基化DMR与LTR反转录转座子富集区域共定位。第五，软果阶段发生广泛的转录组重编程：转录组分析将五个阶段分为早期组（ST1-ST2）和晚期组（ST3-ST5），ST2到ST3转变伴随最显著的转录变化（2191个基因上调、3914个基因下调），上调基因富集于脂肪酸生物合成、糖酵解、类黄酮生物合成和ABA信号传导等功能。第六，部分基因活性与DNA甲基化变化相关：CG/CHG/CHH低甲基化和CHH高甲基化相关DEG数量远超CG或CHG高甲基化相关DEG，启动子或基因体低甲基化DMR与基因上调和下调均相关。第七，促进甲基化基因表达和DNA甲基化的表观遗传调控因子转录动态：甲基结合转录因子和蛋白呈现早期表达组（如RoLEC2、RoSUVH1/3）和晚期表达组（如RoMBD9和RoDNAJ2）两种模式。第八，花青素化合物主要在红色至黑色阶段合成：代谢组分析检测到89种代谢物包括7种花青素，分为红色阶段（ST2）大幅增加组和软果阶段（ST3）上升组。第九，花青素积累与生物合成基因及其调控转录因子表达相关：12个酶编码基因转录本与花青素水平强相关（r≥0.8），RobHLH3与花青素含量强正相关，RoSPL9强负相关。第十，花青素生物合成基因及其转录因子在基因网络中的互连：WGCNA将14457个基因组织为14个调控模块，模块6和13与七种花青素高度关联，包含RobHLH3、RoBBX22、RoMYB10-1和RoTTG1四个正相关转录因子。第十一，果实基因组异染色质中非CG甲基化水平高于叶片：CHG和CHH甲基化在异染色质中果实显著高于叶片，而在真染色质区域未观察到此模式。第十二，DNA甲基转移酶和压缩组蛋白基因的同时升高表达提示成熟过程中独特的染色质动态：RoCMT2在果实晚期活跃，RoDRM2和大多数RdDM途径基因主要在早期表达，RoCMT3从早期到晚期均有活性；H2A.W变体转录本在成熟晚期增加。第十三，多种染色质重塑因子在果实成熟过程中转录活跃：44个重塑因子分为早期高活性组（包括RoDDM1）和晚期上调组，多个基因家族在两组中均有代表，提示亚功能化以适应阶段特异性功能。

讨论部分总结了以下要点：动态甲基化变化在生殖过程中对增强转座子沉默和维持基因组稳定性具有关键作用，黑树莓成熟可能涉及类似机制。研究人员采用较宽松的p值阈值（p<0.05）可能增加假阳性DMR数量，这是本研究的局限性之一，但通过整合DMRs与DEGs增强了生物学支持。内含子DMR的普遍存在提示转录后调控的可能性，内含子甲基化影响可变剪接和转录本异构体产生。启动子区域的局部CG去甲基化伴随CHG和CHH同时去甲基化，与转录改变相关并富集果实成熟功能。RoDML3和RoNOR在整个果实发育过程中高表达，RoNOR启动子CG去甲基化随果实成熟发生，与番茄中NOR直接靶向激活SlDML2形成正反馈循环一致。与番茄不同，RoCNR、RoTDR4和RoRIN随成熟进展均显示强转录抑制，这可能归因于呼吸跃变型和非呼吸跃变型果实之间的差异。CHH高甲基化与基因激活和抑制均相关，与草莓中的研究一致。黑树莓果实异染色质中非CG甲基化显著升高，表明果实特异性异染色质重塑程序，与拟南芥胚乳、花粉和胚胎中的发现相似。许多染色质重塑因子基因在不同阶段受转录动态调控，可能 contribu
te to 升高的非CG甲基化。研究结论翻译如下：本研究揭示了 orchestrating 黑树莓果实成熟和花青素生物合成的遗传和表观遗传程序。研究强调了全基因组和位点特异性表观遗传重编程，并证明了黑树莓果实成熟过程中DNA甲基化与转录调控之间的协调互作。