《Signal Transduction and Targeted Therapy》：RMzyme: regulations of RNA-modifying enzymes in humans
编辑推荐：
本研究为解决人类RNA修饰蛋白(RMPs)的调控机制与功能尚不明确的问题，系统地分析了378个涵盖基因组、转录组、表观转录组、蛋白质组及翻译后修饰(PTM)的多组学数据集。研究人员构建了名为RMzyme的平台，在63种人体组织/细胞类型中，揭示了RNA修饰的动态景观、细胞类型特异性模式及其在急性髓系白血病(AML)等疾病中的关键作用，例如确认了ALKBH5在m6A动态调控和白血病干细胞自我更新中的核心功能。该研究整合了多组学数据，提供了关于RMPs及其下游效应的深刻见解，为理解RNA修饰在人类疾病分子机制中的作用提供了宝贵资源。
在生命这部由DNA编写的宏伟剧本之外，还存在着一个精细而活跃的“后记”层——RNA修饰。如同在写完的信件上进行批注，这些化学修饰在不改变RNA碱基序列的前提下，动态调控着RNA的命运，深刻影响着基因表达、细胞功能乃至人类健康。迄今为止，科学家们已发现了超过170种RNA化学修饰，构成了复杂的“表观转录组”景观。研究表明，RNA修饰及其调控蛋白(RNA-modifying proteins, RMPs)与多种人类疾病，尤其是癌症的发生发展密切相关。然而，尽管相关数据库不断涌现，但对于RMPs如何在各种生理和病理情境中行使其功能，其自身又受到怎样的调控，我们仍知之甚少。为了拨开这层迷雾，一项发表在《Signal Transduction and Targeted Therapy》杂志上的研究，通过整合海量多组学数据，绘制了迄今为止最全面的人类RMPs调控图谱，并构建了一个名为RMzyme的知识平台，为我们深入理解RNA修饰的奥秘打开了一扇新的大门。
为了系统地阐明RMPs的调控网络与功能，研究团队汇集并深入分析了来自63种人类组织的378个多组学数据集，涵盖了基因组学、转录组学（包括单细胞水平）、表观转录组学、蛋白质组学以及蛋白质的翻译后修饰(Post-translational modifications, PTMs)数据。研究者对来自42个细胞系的表观转录组数据进行重新处理，识别了超过200万个RNA修饰位点，并分析了RMPs扰动实验（如敲低、过表达）引起的差异修饰峰和差异表达基因。在单细胞层面，研究者利用非负矩阵分解(Nonnegative matrix factorization, NMF)算法，基于179个RMPs的表达谱，在264个单细胞RNA测序数据集中鉴定了RNA修饰相关的细胞亚型。此外，研究者还从文献和公共数据库中收集、整理了RMPs的PTM信息，并利用临床蛋白质组肿瘤分析联盟(CPTAC)的蛋白质组学数据，分析了RMPs在多种癌症中的蛋白质表达、PTM水平及蛋白质相互作用网络。
RNA修饰在细胞系中的分布图谱
研究首先在42个细胞系中系统描绘了8种主要RNA修饰（包括m6A、m5C、m7G等）的转录组-wide分布。分析发现了超过200万个修饰位点，其中m6A最为丰富。不同的修饰类型显示出独特的区域偏好性，例如m1A主要富集在5‘非翻译区(UTR)和起始密码子附近，而m6A则广泛分布于mRNA的5’UTR、编码区(CDS)、3‘UTR以及非编码RNA中。这些分布模式揭示了不同RNA修饰可能具有各自独特的生物合成机制和功能。
不同条件下RNA修饰峰变化的检测与注释
通过分析RMPs扰动（如METTL3敲低、ALKBH5过表达）后的转录组和表观转录组数据，研究者鉴定了大量差异表达基因和差异修饰峰。以急性髓系白血病细胞系为例，在ALKBH5过表达的NOMO1细胞中，ANXA1基因的m6A甲基化和表达均显著降低；而在METTL3敲低的MOLM13细胞中，FPR1基因的m6A水平下降，表达上升。功能富集分析显示，这些受影响的基因参与了mRNA加工、翻译起始、免疫调控等重要通路。研究还发现，这些差异修饰位点常与微小RNA(miRNA)靶位点、RNA结合蛋白(RBP)结合区域或体细胞单核苷酸变异(SNV)位点重叠，暗示了RNA修饰可能与这些转录后调控层面存在复杂的相互作用。
单细胞水平RNA修饰的调控模式
研究整合了264个单细胞转录组数据集，利用NMF算法鉴定出7538个由RMP表达定义的细胞亚簇。分析发现，m6A相关调节酶在大多数细胞类型中广泛表达，而m7G和m5C相关酶则分别在胚胎细胞和髓系细胞中富集，表现出细胞类型特异性。在急性髓系白血病中，与健康个体相比，患者体内METTL14+B细胞、ALKBH5+干细胞等RNA修饰相关细胞亚群的比例显著升高。对ALKBH5+干细胞的功能分析显示，其差异表达基因富集在内质网蛋白质加工、泛素介导的蛋白酶体降解、自噬等通路。细胞通讯分析进一步揭示了ALKBH5+干细胞通过特定的配体-受体对（如MIF--CD74+CXCR4）与肿瘤微环境细胞进行活跃交流。实验验证表明，ALKBH5在白血病干细胞中高表达，其过表达能促进细胞增殖、周期进程和ZNF503表达，而敲低则损害干细胞自我更新，证实了其关键作用。
与RNA修饰和人类疾病相关的细胞类型特异性体细胞变异
研究在120个单细胞转录组数据集的疾病样本中检测了超过1000万个单细胞体细胞变异，并将其中位于RNA修饰位点附近的变异与RMP介导的NMF细胞簇相关联。在急性髓系白血病中，NPM1是最常发生突变（如p.E105K错义突变）的基因之一，且突变主要富集在造血干细胞、多能祖细胞等早期细胞中。分析还发现，ALKBH5是与体细胞变异关联最频繁的RMP，其相关的变异在粒细胞-单核细胞祖细胞(GMP)和干细胞中尤为常见，提示m6A去甲基化过程的破坏可能与特定细胞类型的恶性转化有关。此外，研究观察到某些组织对（如乳腺-甲状腺、肺-胰腺）共享较高比例的RNA修饰相关变异位点，这可能反映了它们在发育起源或生理调控上的相似性。
RNA修饰酶的蛋白质表达谱、翻译后修饰与突变图谱
研究者从文献和数据库中收集了179个人类RMPs上总计28种类型、11378个PTM位点的信息，其中磷酸化是最常见的修饰类型，其次是泛素化和乙酰化。绝大多数RMPs具有多种PTM类型，表明了其调控的复杂性。利用CPTAC癌症蛋白质组学数据，研究鉴定了与RMPs在RNA、蛋白质及PTM位点水平上显著相关的潜在互作蛋白。例如，m5C甲基转移酶NSUN2在肿瘤中的乙酰化、磷酸化和泛素化位点水平均上调。研究还发现，RMPs的突变频率与癌症的代谢活性相关，例如ALKBH5和EIF3A在代谢活跃的肿瘤中突变更频繁。值得注意的是，在子宫内膜癌中，研究者发现YTHDC1基因的一个移码突变(p.N97Mfs*13)位于一个已知的乙酰化位点(K90)附近，这种PTM与突变的共现可能通过稳定截短蛋白并损害其识别m6A的能力，从而促进肿瘤发生。
综合以上结果，本研究通过大规模多组学整合分析，极大地深化了我们对人类RNA修饰酶调控网络的理解。研究揭示了RNA修饰，特别是m6A，在转录组中的广泛分布和动态特性，并阐明了RMPs扰动如何通过改变特定靶基因（如ANXA1, ZNF503）的修饰状态和表达水平，影响诸如翻译起始、免疫应答等关键细胞通路。在单细胞分辨率下，研究首次系统描绘了RNA修饰模式的细胞类型特异性，并鉴定出在疾病（如急性髓系白血病）中异常扩增的、由特定RMP（如ALKBH5）定义的干细胞样亚群，这些细胞具有独特的代谢特征和细胞间通讯网络，可能是疾病进展的驱动因素。此外，研究将体细胞遗传变异、蛋白质翻译后修饰与RNA修饰景观联系起来，发现了它们之间潜在的功能交叉，例如某些突变可能通过影响RMPs的PTM状态或直接位于修饰位点附近，从而干扰正常的RNA代谢过程。
这项研究的核心贡献在于构建了RMzyme这一综合性知识库和平台。它不仅是首个利用数百个多组学数据集对多种RNA修饰类型及其调节酶进行全面表征的资源，更重要的是，它从分子、细胞到组织水平，整合了遗传、表观转录、蛋白质互作等多维度信息，为揭示RMPs在生理和病理状态下的复杂调控机制提供了前所未有的视角。该资源有望加速针对RNA修饰通路的生物标志物发现和靶向治疗策略的开发，例如针对ALKBH5等关键RMP的小分子抑制剂，从而为癌症等多种人类疾病的精准医学研究开辟新的道路。
