《NAR Genomics and Bioinformatics》：Rp-Bp and the Ribotools box for translatome profiling
编辑推荐：
研究人员提出Ribotools“box”，用于分析配对的核糖体图谱（Ribo-seq）与RNA测序（RNA-seq）数据。Ribotools可独立使用，也可与Rp-Bp联合用于Ribo-seq质量控制（QC）、开放阅读框（ORF）发现、翻译效率（TE）及差异表
研究人员提出Ribotools“box”，用于分析配对的核糖体图谱（Ribo-seq）与RNA测序（RNA-seq）数据。Ribotools可独立使用，也可与Rp-Bp联合用于Ribo-seq质量控制（QC）、开放阅读框（ORF）发现、翻译效率（TE）及差异表达分析。在本研究中，研究人员使用Ribotools刻画人诱导多能干细胞（hiPSCs）向心肌细胞（CMs）的分化，揭示了一个高度连贯的转变：从可塑的、信号驱动的状态向稳定的、结构定型的心肌细胞状态过渡，并经由多层翻译控制分阶段推进。结合互补的蛋白质组学，所得数据集为研究hiPSC心脏谱系定向模型提供了丰富的多组学资源。此外，研究人员提出了一个可重复且可扩展的流程，用于基准测试翻译效率分析工具。Ribotools与Rp-Bp可在Bioconda上获得，有助于使研究更易获取且可重复。详细文档、测试数据集和教程见https://ribotools.readthedocs.io与https://rp-bp.readthedocs.io。
研究背景：核糖体图谱（Ribo-seq；基于RNA测序的RNA翻译读出）通过分离并深度测序核糖体保护RNA片段（RPFs；核糖体足迹），在亚密码子分辨率提供全基因组翻译组快照。自2010年以来，Ribo-seq已广泛用于多种生物、组织或细胞，揭示由数千个未注释开放阅读框（ORFs；可翻译阅读框）构成的“暗翻译组”。ORF鉴定、表征和注释，以及组织、阶段或疾病特异性翻译调控研究持续受到关注；相关公共数据库、服务器、标准化注释和指南亦在发展。已有工具多利用翻译核糖体的3-nt周期性特征识别Ribo-seq ORFs，其中包括Rp-Bp；Rp-Bp已用于需要无监督选择读长以稳健估计周期性的研究。当前问题在于：翻译效率分析（TEA）领域缺乏稳健、模块化、端到端（E2E）、积极维护、文档完善，并可将特征按基因调控层级进行原则性分类的方案。为此，研究人员提出Ribotools，与Rp-Bp组合形成Ribo-seq质量控制（QC）、ORF发现和TEA的E2E解决方案；提供TEA工具基准流程，并与现有ORF发现工具比较；同时生成hiPSCs向CMs分化的时间过程数据集，以解析心脏谱系定向中的翻译调控。
研究人员开展了工具基准比较、Ribo-seq/RNA-seq配对分析及蛋白质组学验证，得出Ribotools为最佳全方位工具，并可揭示hiPSC向CM分化中从早期谱系程序缓冲、定向时翻译推进到后期信号与重塑通路抑制的动态轨迹；所生成ORF图谱和多组学数据为研究hiPSC心脏发育模型提供资源与蓝图。
关键技术方法（≤250字）：研究人员采用Lukas Cyganek提供的hiPSC系UMGi014-C clone 14（isWT1.14），将其分化为心室CMs，并在分化第0、1、5、7、17、31天采集配对样本，开展RNA-seq与Ribo-seq，同时补充质谱蛋白质组学。分析端以GRCh38和Ensembl release 102指导从头转录组注释，使用Rp-Bp预测Ribo-seq ORFs，使用Ribotools完成QC、丰度估计、TE/差异表达与调控分类，并提供可重复的TEA工具基准流程。
研究结果：
Ribotools stands out as the strongest all-around tool：通过第一项基准，研究人员比较Ribotools、Riborex、anota2seq和deltaTE的检测性能，发现Ribotools与deltaTE相当，并在多数测试中略优于anota2seq；Riborex在翻译特异性比较中的表面优势反映其未区分TE交互评分而非真正优越，在考虑所有调控类时表现最差，原因在于不能检测丰度调控基因。分类准确性方面，各分类工具对abundance和background类高且一致，对translation和buffering类中等，anota2seq对translation基因回收较低。第二项基准使用新数据集子集进行标签置换，在零假设下FDR < 0.05时假阳性为零，支持真实数据中显著基因的特异性，但受n = 2重复数限制。定性比较ORF发现工具显示，Ribotools/Rp-Bp是少数积极维护且有广泛在线文档和教程的选项，此点仅与ORFik可比；仅Ribotools/Rp-Bp提供交互式可视化仪表板；ORFik的ORF翻译预测依赖启发式监督学习，而Rp-Bp基于概率模型进行无偏de novo发现，提供完整后验分布而非点分类。结论是Ribotools为唯一E2E、模块化且积极维护的TEA工具，整体表现最佳，并与Rp-Bp构成稳健的Ribo-seq分析、QC和ORF发现工作流。
Transcriptome and translatome of hiPSC differentiation into cardiomyocytes：通过Ribotools分析hiPSCs向CMs分化的配对转录组与翻译组，样本覆盖hiPSCs（第0天）、早期中胚层（第1天）、心脏中胚层（第5天）、心脏祖细胞（第7天）、早期CMs（第17天）和CMs（第31天），获得心脏谱系定向的全局RNA翻译快照，并表明hiPSC-CMs是研究心脏谱系定向的合理模型。
Quality control：通过QC仪表板评估各Ribo-seq库，每个文库平均产生1050万条比对 reads，其中约670万条周期性 reads用于后续分析；reads主要长度为28或29 nt，平均84%比对到注释编码序列（CDSs）的主阅读框；长度分布和metagene图谱可用于检查读长周期性，并可观察第五密码子位置翻译起始后暂停等翻译特征，说明数据质量高。
ORF discovery：通过Rp-Bp和ORF发现仪表板，研究人员发现Ribo-seq ORFs主要位于Ensembl注释CDSs内，但每个发育阶段平均还翻译1861个上游开放阅读框（uORFs）、1425个注释为非编码转录本中的ORF（ncORFs）和约1218个注释CDS变体；从头转录组使研究人员预测到392个完全新颖、未表征的ORFs；1941个Ribo-seq ORFs匹配GENCODE Ribo-seq ORFs共识集，其中约95%为未注释小开放阅读框（sORFs），多为uORFs，例如GATA4 uORF及SLC2A3上游的c12riboseqorf16，后者被Ribotools鉴定为第17天后翻译下调。
Downstream analyses using ORF counts：通过Ribotools丰度估计流程获得Ribo-seq ORF计数并进行主成分分析（PCA；降维可视化），发现hiPSCs向CMs进展在转录组和翻译组水平均按可比发育顺序区分；第0天和第1天细胞轮廓高度相似并表达多能性和早期中胚层标志物，第5天和第7天由心脏定向中胚层进展为心脏祖细胞，第17天后主要表达CM标志物；心脏定向hiPSCs的丰度谱显示MYH7:MYH6、TNII3:TNNI1、MYL2:MYL7和IRX4:NR2F2比值由第5、7天较低转为第17、31天较高，确认由胎儿样向成熟样CMs转变。
Translational control during hiPSC differentiation into cardiomyocytes：通过Ribotools，研究人员显示CM分化中的翻译调控遵循动态轨迹，并具有定向、执行、控制和稳定等多层控制。第1天出现协调变化，介导细胞分化、细胞命运和发育相关广谱细胞反应。第5天细胞退出多能性并发生谱系定向时经历重大重编程；早期CM分化中翻译相关基因主要发生翻译下调，与从增殖性向谱系特异性程序转变时翻译机器协调抑制一致。基因集富集分析（GSEA；基因集合富集检验）显示rRNA加工（NES ?2.8）、翻译延伸（NES ?2.6）、翻译起始（NES ?2.6）和细胞质翻译（NES ?2.2）负富集，其中大量翻译下调核糖体蛋白基因属于核心富集集；这些变化可能反映细胞获得谱系特异性翻译程序时翻译能力和选择性的转变，是否涉及核糖体组成或特化改变仍待确定。同时，关键发育心脏通路在转录上激活但在翻译上被缓冲，提示“时序和剂量”机制精细调节心脏蛋白表达；第5天被缓冲通路富集于心脏隔膜形态发生（NES 2.3）、心肌组织（NES 2.2）和心内膜垫（NES 2.1）。差异蛋白质组学部分确认上述观察，上调蛋白仍主要富集于替代谱系程序。第7天，心脏和肌源性基因由缓冲转为翻译推进，表明翻译和蛋白水平的谱系限制及分化程序主动执行；心脏、脂质代谢和核转运过程的持续缓冲提示功能与代谢网络仍在精细调节和分阶段激活。第17天向成熟转变时，细胞由响应/适应状态转向更稳定的功能状态，在翻译水平主动减少重塑和信号可塑性并关闭替代命运决定程序；尽管翻译水平发生细化和重组，蛋白丰度变化仍主要与心脏谱系和功能程序相关。第31天系统达到稳定状态，蛋白产出似被主动维持而非动态重塑；在无变化条件下，翻译调控主要服务蛋白稳态，多数通路通过缓冲机制保持稳定；RNA丰度变化由翻译效率补偿（显著交互项），尽管转录水平显著变化，RPF水平仍稳定。多数发育阶段差异蛋白质组学与Ribotools分类预期广泛一致。
讨论部分总结：研究人员指出，RNA-seq与Ribo-seq联合分析具有挑战性且需要可靠工具；Ribotools box是开源、基于配置、模块化的软件包，并由Rp-Bp支撑。与朴素比率方法不同，TE变化通过在RNA-seq和Ribo-seq计数拟合的负二项广义线性模型（GLM；此处为负二项模型）中检验交互项推断。该框架保证方法可靠性，同时弥补TEA领域实践缺口：专用QC工具侧重QC和可视化而非完整分析，仅TEA工具面向现成计数数据且缺乏Ribo-seq特异性QC，且多数已弃用或过时。Ribotools是少数文档完善且积极维护的E2E工具之一，通过易用而稳健的统计框架“开箱”提供QC、ORF发现和TEA；其价值不在引入全新算法，而在提供连贯实用的E2E解决方案。未来工作可包括改进QC指标和可视化、交互式TEA仪表板及特定协议支持。尽管有这些限制，Ribotools在检测、分类及TEA与ORF发现工具整体比较中表现最佳；基准并非另一次穷尽比较，而是精简、务实、可演进的框架，完整工作流开放且完全可重复。研究人员还指出，翻译在人心脏发生细胞命运决定中的积极作用早已被认识，Ribo-seq近年越来越多用于心脏重编程和大规模心脏翻译组表征；本研究展示Rp-Bp和Ribotools如何简化复杂数据集分析并提供心脏谱系定向机制见解。研究人员承认Ribo-seq ORFs的蛋白质组学支持仍有限，尤其是sORFs，这是公认挑战；研究目标是证明Ribotools不仅简化翻译组图谱分析，也促进ORF发现早期阶段。
研究结论部分翻译：研究人员的目的是证明Ribotools不仅简化翻译组图谱分析，而且促进ORF发现的早期阶段。这些广泛数据为将hiPSC-CMs作为模型系统进行研究提供了强大资源，并提供了潜在新型调控性uORFs、替代翻译事件以及更多内容的来源。
