Uridylation regulates mRNA decay directionality in fission yeast
尿苷化调控裂殖酵母mRNA降解方向
期刊《Nature Communications》，2024年第15卷，文章编号8359，DOI:10.1038/s41467-024-50824，2024年9月27日在线发表
1 摘要总结
真核细胞质mRNA存在5'→3'、3'→5'两条外切降解通路，poly(A)尾尿苷化修饰广泛存在，但全局mRNA周转中的功能尚不清晰。本研究开发全基因组3'-RACE（gw3'-RACE）测序技术，以裂殖酵母为模式生物解析尿苷化调控机制：细胞存在两类尿苷化修饰，分别由两种末端尿苷转移酶催化；Cid1负责缩短poly(A)尾加尿苷，Cid16对无尾RNA进行寡尿苷修饰。poly(A)尾尿苷化可招募Lsm1-7复合物，推动5'→3'降解；敲除Cid1导致过度去腺苷，mRNA转向3'→5'降解途径。该研究阐明尿苷化通过双向调控，决定mRNA降解通路选择，完善真核mRNA周转精细调控网络。
2 中文关键词
裂殖酵母，mRNA降解，poly(A)尾，尿苷化，末端尿苷转移酶，Cid1，Cid16，5'-3'降解，3'-5'降解，全基因组3'-RACE
3 研究目的
1. 建立可同时捕获有尾、无尾RNA的gw3'-RACE高通量测序方法，完整解析全转录组3'末端修饰图谱；
2. 区分裂殖酵母两类RNA尿苷化修饰（poly(A)U修饰、寡U修饰），鉴定对应的催化酶Cid1、Cid16；
3. 阐明poly(A)尾尿苷化对mRNA降解通路（5'→3' / 3'→5'）的定向调控作用；
4. 解析尿苷化缺失后过度去腺苷、3'→5'降解代偿的分子机制；
5. 结合Lsm1、Ski、Dis3l2突变株，构建尿苷化介导mRNA降解通路切换完整调控模型；
6. 明确高尿苷化靶基因功能特征，关联糖、氨基酸代谢通路基因表达稳态调控。
4 研究思路
1. 方法学构建：开发gw3'-RACE建库测序流程，可同时捕获poly(A)、poly(A)U、寡U、无尾四类3'末端，与PAL-seq交叉验证尾长数据（图1）；
2. 单/双敲菌株构建：构建Δcid1、Δcid16、Δcid1Δcid16突变株，gw3'-RACE测序对比各类3'末端丰度变化，区分两种酶底物偏好（图3）；
3. 尾长、修饰定量分析：统计不同长度poly(A)尾尿苷化比例，关联基因拷贝数、半衰期、功能分类（图2）；
4. 5'降解通路验证：构建Δlsm1、Δedc1、Δxrn1突变株，检测尿苷化尾积累，验证Lsm1-7是尿苷下游效应因子（图4）；
5. 3'降解通路解析：敲除Ski2、Dis3l2，联合Δcid1双突变，分析无尾、寡U RNA积累，明确3'→5'代偿通路（图5）；
6. 生长动力学检测：采用Bioscreen高通量仪测定所有单/双/三突变株生长曲线，关联降解通路失衡的生理表型（图4c、图5c）；
7. 多组学验证：RNA-seq差异转录、SWATH定量蛋白组、Northern印迹测定mRNA半衰期，关联尿苷化频率与基因稳态（图4h-l）；
8. 整合全部测序、酶活、生长、组学数据，绘制尿苷化调控mRNA降解通路模型（图6）。
5 研究亮点
1. 创新gw3'-RACE技术，一次性捕获所有类型RNA 3'末端修饰，突破传统测序只能检测poly(A)尾局限；
2. 首次在裂殖酵母区分两种功能分化尿苷转移酶：Cid1专一修饰短poly(A)尾，Cid16催化无RNA寡尿苷，底物完全分离；
3. 揭示尿苷化双重调控功能：既招募Lsm1-7促进5'→3'降解，又阻止mRNA过度去腺苷、抑制3'→5'降解；
4. 证明尿苷缺失时Ski介导外体、Dis3l2两条3'→5'通路代偿激活，形成降解通路平衡切换机制；
5. 发现糖代谢、氨基酸合成类mRNA高频率发生尿苷修饰，修饰调控其表达动态平衡；
6. 多突变株遗传上位分析完整搭建RNA 3'修饰—两条降解通路的调控层级模型。
6 可延伸的研究方向
1. 体外重组Cid1/Cid16蛋白，解析两种酶识别不同RNA底物的结构分子机制；
2. 拓展人细胞样本，验证同源TUT4/TUT7是否存在与酵母保守的降解调控逻辑；
3. 筛选小分子抑制剂特异性抑制Cid1，用于调控细胞mRNA稳态、肿瘤基因表达；
4. 动态胁迫（营养饥饿、氧化应激）下测序，探究胁迫如何重编程全转录组尿苷化图谱；
5. 单分子测序观测活细胞mRNA尾动态尿苷修饰、降解时序过程；
6. 结合核糖体图谱，解析尿苷化如何联动翻译与mRNA降解的交叉调控。
7 测量数据、对应图表与研究意义
1. gw3'-RACE方法原理图、标准品尾长检测密度图（图1a-c）：验证测序技术精准定量poly(A)尾长度，PAL-seq相关性r=0.79，证明方法可靠；
2. 全转录组3'末端分类、尾长-尿苷化关联散点图（图1d、图2）：明确短poly(A)尾更易发生尿苷修饰，高表达代谢基因修饰频率最高；
3. cid1/cid16突变株3'末端丰度、尾长分布（图3）：直接证明Cid1主导poly(A)U修饰，Cid16负责寡U修饰；
4. Δlsm1等5'通路突变株尾长、转录组、蛋白组数据（图4）：证实尿苷修饰通过Lsm1复合物启动脱帽5'降解；
5. Ski2/Dis3l2单/多突变末端分布与生长曲线（图5）：阐明尿苷缺失时无尾RNA增多，3'→5'通路代偿激活；
6. 通路调控模型示意图（图6）：可视化尿苷化双向平衡两条mRNA降解通路的完整分子流程；
7. Bioscreen生长曲线（图4c、图5c）：量化各降解通路突变株生长缺陷，验证通路失衡带来细胞生理损伤。
8 整体研究结论
1. 自主研发gw3'-RACE测序，可全面解析含poly(A)U、寡U、无尾在内所有RNA 3'末端修饰，是尾修饰研究核心工具；
2. 裂殖酵母存在功能分工的两种尿苷转移酶：Cid1优先修饰短poly(A)尾生成poly(A)U；Cid16针对完全去腺苷无尾RNA添加寡U；
3. poly(A)尾尿苷化具备双重调控功能：一是招募Lsm1-7复合物，激活脱帽并导向5'→3'降解；二是阻断持续去腺苷，避免mRNA被3'→5'外切酶降解；
4. 敲除Cid1后poly(A)尿苷化大幅下降，mRNA发生过度去腺苷，无尾RNA积累，Ski/Dis3l2介导的3'→5'降解通路代偿增强；
5. 糖酵解、氨基酸代谢相关mRNA尿苷修饰频率显著高于其他基因，尿苷化维持其转录产物动态平衡；
6. 两条降解通路存在互补代偿关系：正常条件5'→3'降解占主导，尿苷修饰缺失时3'→5'通路成为主要降解途径，二者由3'端尿苷化修饰决定平衡。
9 芬兰Bioscreen仪器微生物生长曲线数据详细研究意义
本研究使用Bioscreen C全自动微生物检测仪，每15 min采集OD600，全程恒温振荡，批量测定野生型、cid/cid16/lsm1/ski2/dis3l2单、双、三重突变株液体生长曲线，研究意义如下：
1. 高通量同步上百株突变平行对照，统一接种浓度与培养基，消除分批摇瓶人工误差，直观对比不同RNA降解通路缺陷的生长差异；
2. 连续动态记录迟滞期、对数期、稳定期，区分轻微生长迟缓与重度增殖缺陷，定量评估尿苷化缺失、降解通路紊乱对细胞生理影响程度；
3. 实现遗传上位表型判定：Δcid1单株生长轻微受损，Δcid1Δski2双株生长严重停滞，直观证明3'→5'通路代偿会加重细胞生存压力；
4. 区分完全培养基（YES）与限定合成培养基（EMM）条件下表型差异，明确营养环境会放大RNA降解失衡缺陷；
5. 生长表型与gw3'-RACE、RNA-seq测序数据联动，建立“酶缺失→3'末端修饰改变→降解通路切换→细胞生长缺陷”完整证据链；
6. 标准化时序OD数据可拟合生长速率、倍增时间等动力学参数，为不同突变株缺陷程度提供定量统计学依据。