Argonaute 1对拟南芥信使RNA 5-甲基胞嘧啶修饰的新型调控作用

《The Plant Genome》：The novel regulatory influence of Argonaute 1 on messenger RNA 5-methylcytosine modification in Arabidopsis thaliana

5-甲基胞苷（m5C）修饰在RNA代谢中发挥关键作用。TRM4B（tRNA特异性甲基转移酶4B）是目前已知的唯一介导mRNA（信使RNA）甲基化的植物m5C机器组分。迄今为止，除了一些简并序列基序外，尚未发现调控位点特异性m5C修饰的调控因子。借鉴动物中的相关发现，研究人员利用拟南芥ago1突变体，研究了AGO1（Argonaute 1）在m5C调控中的作用。m5C-RIP-seq（m5C RNA免疫沉淀结合深度测序）分别在野生型（WT）和ago1中鉴定出高置信度的5981个和4559个m5C峰。这些峰中有相当一部分在WT和ago1之间存在重叠，表明其具有高度的基因特异性。与WT相比，ago1表现出366个超甲基化峰和168个低甲基化峰。虽然超甲基化和低甲基化都可能与基因表达的上调或下调相关，但具有相对复杂结构的活跃转录基因优先受到m5C修饰的影响。相当一部分m5C相关基因在功能上参与mRNA监控和核糖体生物发生。有趣的是，在WT中鉴定的超过60%的m5C峰与miRNA（微小RNA）结合位点相关联。此外，基于AGO1-RIP-seq实验，在超过40%的m5C峰对应的转录本上可以检测到AGO1结合信号。免疫沉淀联合质谱分析验证了TRM4B与TRANSPORTIN1（TRN1）之间的相互作用，TRN1是一种已知的AGO1互作蛋白，也参与miRNA通路。这些发现支持了一个模型，即miRNA引导的AGO1复合物在介导位点特异性m5C修饰中发挥关键作用。

一、研究背景

RNA修饰是真核生物基因表达调控的重要层面。5-甲基胞苷（m⁵C）是继N⁶-甲基腺苷（m⁶A）之后最普遍的RNA甲基化标记之一，广泛参与RNA加工、结构重塑、稳定性控制、翻译调控及亚细胞定位等多种RNA调控过程。在模式植物拟南芥中，TRM4B（tRNA特异性甲基转移酶4B）是目前唯一已知的介导mRNA（信使RNA）和长链非编码RNA（lncRNA）m⁵C修饰的“写入器”。然而，关于植物m⁵C调控机制的基础认知仍然十分有限。例如，除了几个简并的序列基序（如HACCR和CUYCUYC）被认为可能指导位点特异性m⁵C修饰外，尚无明确的调控因子被发现。鉴于m⁵C修饰仅发生在0.02%-0.09%的胞嘧啶上，这些简并的顺式元件显然不足以赋予精确RNA甲基化所需的高序列特异性。此前在人类mRNA中的一项开创性研究发现，m⁵C位点与Argonaute（AGO）蛋白的结合区域在全转录组水平上存在显著关联。AGO蛋白是RNA干扰（RNAi）通路的核心效应器，由微小RNA（miRNA）或小干扰RNA引导进行序列依赖性靶向。鉴于RNAi和m⁵C修饰在真核生物中高度保守，探索这两个进化保守通路之间潜在的功能串扰具有重要意义。拟南芥基因组编码10个AGO家族成员，其中AGO1是最具特征的旁系同源物，主要与miRNA结合以介导序列特异性的转录后基因调控。因此，本研究旨在系统探究AGO1在拟南芥RNA m⁵C沉积中的潜在调控作用。

二、主要技术方法

研究人员利用拟南芥弱功能突变体ago1-25及其野生型（WT）Col-0的16日龄幼苗，进行了以下关键技术分析：（1）m⁵C-RIP-seq（m⁵C RNA免疫沉淀结合深度测序），每个基因型设置三个生物学重复，用于全转录组m⁵C峰鉴定和定量比较；（2）AGO1-RIP-seq（AGO1 RNA免疫沉淀结合深度测序），使用AGO1特异性抗体对2周龄WT幼苗进行两个生物学重复的实验，以鉴定AGO1直接结合的转录本；（3）小RNA测序（sRNA-seq），每个基因型三个生物学重复，用于分析miRNA的表达变化；（4）免疫沉淀联合质谱（IP-MS），利用过表达AGO1和TRN1-3×FLAG的转基因植株，鉴定AGO1复合物的潜在互作蛋白。

三、研究结果

3.1 ago1突变对全转录组m⁵C修饰的影响

m⁵C-RIP-seq在WT和ago1中分别鉴定出5981个和4559个高置信度m⁵C峰。差异分析发现ago1中存在366个超甲基化峰和168个低甲基化峰。所有m⁵C峰在五条染色体上大致均匀分布，但在着丝粒和端粒区域明显缺失。序列基序富集分析恢复了富含CU和AC的基序，其中AC富集基序的富集依赖于功能性AGO1。基因组特征注释显示，约48%-61%的m⁵C峰位于编码区（CDS）。值得注意的是，超甲基化峰主要富集于3′非翻译区（UTR），而低甲基化峰则显著富集于5′UTR。

3.2 m⁵C峰相关基因的结构与表达

m⁵C修饰优先靶向含有更多外显子、结构更复杂的基因，并且这些基因的表达水平显著高于非甲基化基因。在ago1中，超甲基化和低甲基化均可与基因表达的上调或下调相关。通过整合差异甲基化峰与差异表达基因，研究人员鉴定出一批与发育和胁迫响应相关的候选基因，例如调控细胞骨架组织的AT4G29340、细胞壁组织相关的AT3G54580、冷驯化和盐耐受相关的AT5G15970等，这些基因的m⁵C修饰变化与其表达水平的改变呈现相关性。

3.3 m⁵C峰相关基因的功能分析

功能富集分析表明，m⁵C相关基因显著富集于RNA结合、加工、剪接等核心RNA相关功能类别。京都基因与基因组百科全书（KEGG）通路分析揭示，这些基因在mRNA监控通路和核糖体生物发生通路中显著过度代表。具体而言，在ago1中超甲基化的15个核糖体相关基因编码不同的大亚基和小亚基蛋白，而6个mRNA监控通路相关基因包括外显子连接复合体、蛋白磷酸酶2A、多聚腺苷酸结合蛋白及切割和多聚腺苷酸化特异性因子的组分。尽管这些基因的m⁵C水平发生变化，但其稳态mRNA水平并未检测到统计学显著变化，提示这些超甲基化标记可能被未知的m⁵C“读取器”识别，从而在不影响mRNA丰度的情况下发挥下游转录后调控效应。

3.4 miRNA参与AGO1介导的位点特异性m⁵C调控

通过TargetFinder算法预测miRNA结合位点（MBS），发现在ago1中，59.29%的超甲基化峰和55.36%的低甲基化峰与至少一个预测的MBS相关联。在WT中，超过60%的m⁵C峰与MBS相关联，这一比例在独立发表的数据集中得到验证（61.71%）。MBS与m⁵C峰均显著富集于3′UTR区域。AGO1-RIP-seq实验提供了直接的生化证据，显示在WT中，约47.58%的m⁵C峰位于至少一个重复支持的AGO1结合转录本上；对于超甲基化和低甲基化峰，该比例分别为43.17%和42.26%。超几何分布检验证实，差异甲基化基因与AGO1结合基因之间的重叠具有高度统计显著性。此外，小RNA测序数据显示，多个与差异m⁵C峰相关的miRNA在ago1突变体中显著下调。IP-MS实验未检测到AGO1与TRM4B的直接相互作用，但鉴定出TRN1（一种已知的AGO1互作蛋白）与TRM4B之间存在可靠的互作信号。同时，ago1突变不影响TRM4B的转录水平。

四、讨论与结论

本研究通过全转录组m⁵C-RIP-seq在WT和ago1突变体中鉴定出数千个高置信度m⁵C峰，证实mRNA m⁵C修饰在拟南芥转录组中广泛存在。比较甲基化组分析发现数百个统计显著的差异m⁵C峰，直接证明AGO1是全转录组m⁵C沉积的关键调控因子。超甲基化峰和低甲基化峰在UTR区域呈现截然不同的分布模式，提示其可能受不同机制调控。高度转录且结构复杂的基因优先受到m⁵C修饰，这可能源于更高的底物可用性或m⁵C本身对转录本稳定性的促进作用。功能分析揭示，大量超甲基化基因参与mRNA监控和核糖体生物发生，表明m⁵C对转录本稳定性和翻译的调节效应可通过这些核心调控通路被功能放大。

综合多组学证据——包括m⁵C甲基化组图谱、MBS共定位分析和AGO1 RIP-seq——研究人员提出一个工作模型：miRNA引导的AGO1效应复合物作为拟南芥位点特异性m⁵C沉积的关键上游调控因子。该模型提出两种非互斥的机制：其一，AGO1复合物在其靶标MBS上的占据产生局部空间位阻，阻止m⁵C机器接近和修饰邻近的m⁵C位点；其二，AGO1作为分子支架，通过与m⁵C酶促机器的相互作用（经中间蛋白TRN1桥接），主动招募m⁵C“写入器”或“擦除器”在其结合区域附近沉积或移除m⁵C。在该模型中，miRNA通过序列互补性将AGO1复合物引导至特定区域，而这些区域优先富集于3′UTR。本研究为植物转录组中位点特异性m⁵C沉积的分子机制提供了重要的新见解。