027-62435310 | 
service@benagen.com | 
“长polyA尾更稳定？” 科学家打脸传统观点：RNA调控没这么简单！
文章标题：Nanopore direct RNA sequencing of human transcriptomes reveals the complexity of mRNA modifications and crosstalk between regulatory features
发表期刊：Cell Genomics
发表时间：2025.5
影响因子：11.1
研究介绍
近年来，RNA修饰在基因表达调控中的作用日益受到关注，尤其是N6-甲基腺苷（m⁶A）作为mRNA中最常见的修饰，参与调控mRNA的剪接、稳定性、翻译和降解等过程。然而，m⁶A与其他RNA调控特征（如多聚腺苷酸尾长度、mRNA丰度和选择性剪接）之间的相互作用机制尚未完全阐明。传统研究方法通常依赖短读长测序技术，难以同时捕获RNA修饰和全长转录本信息。本研究利用牛津纳米孔直接RNA测序（dRNA-seq）技术，首次在人类白血病细胞系中同步分析了转录组和表观转录组，揭示了m⁶A修饰与多聚腺苷酸尾长度、mRNA丰度及剪接事件之间的复杂互作网络，为RNA调控的机制研究提供了新的视角和数据资源。
研究结果
1. dRNA-seq全面解析白血病细胞的转录组与表观转录组特征
研究通过牛津纳米孔dRNA-seq技术，对白血病细胞系MOLM13 的4个对照样本（scramble shRNA）和2个METTL3敲低（KD）样本进行测序（图1A）。使用m6Anet软件分析，实现了单分子水平的m⁶A化学计量分析，在四个生物学重复中鉴定出高置信度的m⁶A修饰位点（图1B-C），其中75%的位点在至少两个重复中被检测到，表明数据的高度可重复性。进一步分析显示，m⁶A修饰在转录本上的分布呈现典型的5'UTR富集和剪接位点附近缺失的特征（图1D-E），且外显子长度与m⁶A修饰频率呈正相关（图1F）。
此外，METTL3（m⁶A甲基转移酶）敲除后，m⁶A修饰水平显著降低，尤其是DRACH基序（图1G-H），且修饰比例高的位点受影响更明显（图1I-J）。基因功能分析表明，m⁶A修饰减少的基因主要富集于核糖体生物发生和RNA剪接等通路（图1K），并在CDKN2C和RARA等基因中验证了修饰信号的下降（图1L-M）。
图1 dRNA-seq全面解析白血病细胞的转录组与表观转录组特征
2. m⁶A修饰与mRNA丰度、稳定性及翻译活性的关联分析
研究发现，m⁶A修饰与mRNA丰度之间存在微弱的负相关性（图2A-B），而修饰位点数量和总修饰比例较高的转录本表现出更低的mRNA稳定性（图2E-H）。值得注意的是，m⁶A修饰的转录本与核糖体结合频率呈负相关（图2I-L），提示m⁶A可能通过抑制翻译调控基因表达。METTL3敲低后，77%表达下调基因和85%表达上调基因均与m⁶A修饰减少相关（图2M），且这些基因分别富集于细胞周期调控和代谢通路（图2N-O），表明m⁶A的调控作用具有显著的基因特异性。
图2 m⁶A修饰与mRNA丰度、稳定性及翻译活性的关联分析
3. polyA尾长度的全基因组特征及其与mRNA丰度的关系
通过dRNA-seq，研究首次在单分子水平上精确测量了polyA尾的长度（图3A-C），发现其分布范围从几十到近400 nt，中位长度为92 nt。根据polyA尾长度，转录本可分为短、中、长三组（图3D），其中短尾转录本更倾向于高表达（图3J-K），且与长尾转录本相比，具有更长的半衰期（图3L-M），挑战了“长polyA尾的转录本更稳定且表达水平更高”的传统观点。此外，长尾转录本的3'UTR和CDS区域更长，且AU含量更高（图3E-F）。功能分析显示，短尾基因富集于线粒体呼吸链相关通路，而长尾基因则与GTP酶活性调控相关（图3G-I）。
图3 polyA尾长度的全基因组特征及其与mRNA丰度的关系
4. m⁶A修饰与polyA尾长度的动态关联
研究揭示了m⁶A修饰与polyA尾长度之间的复杂关系：虽然修饰的转录本整体polyA尾更长（图4A-D），但修饰比例高的转录本polyA尾长度轻微增加（图4B）。METTL3敲除后，全局m⁶A水平下降5%-10%，其中186个转录本的polyA尾长度显著增加，而100个转录本的polyA尾长度减少（图4E-G）。功能分析表明，polyA尾长增加的基因与病毒感染和蛋白质折叠通路相关，而polyA尾长减少的基因则参与细胞周期调控（图4H-I），暗示m⁶A可能通过调控polyA尾长度影响特定通路的mRNA命运。
图4 m⁶A修饰与polyA尾长度的动态关联
5. METTL3介导的m⁶A修饰调控选择性剪接
METTL3敲低导致选择性剪接事件显著增加，尤其是外显子跳跃和内含子保留（图5A-C）。这些事件涉及的基因富集于细胞质翻译和核苷酸代谢通路（图5E）。进一步分析显示，m⁶A修饰比例在20%-60%的转录本中剪接事件最多（图5F-G），且修饰位点与剪接位点的空间距离呈现特定分布模式（图5H），表明m⁶A可能通过局部修饰影响剪接因子的结合效率。
图5 METTL3介导的m⁶A修饰调控选择性剪接
研究总结 本研究利用纳米孔dRNA-seq技术，首次在单分子水平上同步解析了mRNA修饰、polyA尾长度和剪接事件之间的动态互作，揭示了m⁶A修饰在转录本稳定性、翻译调控和剪接中的多层次作用。研究不仅挑战了“长polyA尾增强mRNA稳定性”的传统观点，还提出了m⁶A修饰的调控效果高度依赖其修饰比例和基因组背景。这些发现为理解RNA修饰的生物学功能提供了全新视角，并为后续研究提供了宝贵的数据资源。未来，随着测序技术的进步和计算工具的优化，dRNA-seq有望在RNA表观遗传学领域发挥更重要的作用。 参考文献： Kim Y, Saville L, O’Neill K, et al. Nanopore direct RNA sequencing of human transcriptomes reveals the complexity of mRNA modifications and crosstalk between regulatory features[J]. Cell Genomics, 2025.
