《SCIENCE ADVANCES》：Integrative analysis of mRNA stability regulation uncovers a metastasis-suppressive program in breast cancer
编辑推荐：
本研究针对乳腺癌基因表达异质性中转录后调控机制尚不清楚的问题，开发了整合RNA序列特征与RNA结合蛋白表达的深度学习框架GreyHound，以解析mRNA稳定性调控网络。研究发现RBP RBMS3通过稳定包含AUA基序的靶向mRNA（特别是TXNIP）发挥抑癌作用，其缺失导致mRNA不稳定并与不良临床预后及增强的转移潜力相关。体内外实验证实，RBMS3通过调节TXNIP mRNA稳定性驱动了关键的抑制转移程序。这项工作揭示了一个新的乳腺癌转录后调控机制，并展示了可解释的RNA动态模型在疾病调控程序发现中的能力。
在癌症的复杂图景中，基因表达的异常是核心特征。长久以来，科学家们主要从基因本身的突变和表现遗传修饰层面探寻癌症异质性的根源。然而，越来越多的证据表明，在基因完成转录、信使RNA（mRNA）诞生之后，其命运的调控——即转录后调控——同样扮演着至关重要的角色。特别是在乳腺癌中，癌细胞转移是导致患者死亡的首要原因，而驱动转移的基因调控网络变化仍未完全明晰。mRNA的稳定性，作为决定其最终丰度和功能的关键一环，如何在不同的乳腺癌细胞中变化，又由哪些“幕后推手”精确调控，进而影响疾病的进程？这些问题的答案，如同一片尚待探索的“暗物质”领域，对理解乳腺癌的进展和开发新的治疗策略具有迫切的科学意义。以往的研究通常将影响mRNA稳定的顺式（cis，指RNA序列本身）因素和反式（trans，指RNA结合蛋白等细胞因子）因素分开探讨，但生命系统的精妙在于其复杂的相互作用。为了系统性地照亮这片“暗物质”，一项发表于《SCIENCE ADVANCES》的研究应运而生。研究人员开发了一个名为“GreyHound”的革命性计算框架，它像一位训练有素的“猎犬”，能够整合RNA序列信息与细胞中RNA结合蛋白的表达谱，精准预测并解析不同细胞状态下mRNA稳定性的调控逻辑。通过这条“猎犬”的引领，研究团队成功锁定了一个以前未被充分认识的抑制乳腺癌转移的关键调控轴，为理解癌症的转录后调控机制打开了新的大门。
为开展此项研究，作者综合运用了多项前沿技术。首先，他们使用SLAM-seq（硫醇链接烷基化代谢RNA测序）技术在六个代表不同亚型（管腔型、HER2阳性、三阴性）的乳腺癌细胞系中系统测量了mRNA的动态（合成与降解）。其次，他们开发了GreyHound这个可解释的多模态深度学习框架，其架构整合了卷积神经网络处理RNA序列和一个预训练的变分自编码器处理RBP表达数据，用以预测mRNA稳定性并识别候选的转录后调控网络。关键的实验验证技术包括：利用短发卡RNA（shRNA）在细胞中进行基因敲低；通过改进的CLIP-seq（紫外交联免疫沉淀测序）在全转录组范围绘制RBMS3蛋白的RNA结合位点；使用双荧光报告系统验证特定序列元件对mRNA稳定性的影响；在免疫缺陷的NSG小鼠中进行尾静脉注射肺定植实验以评估体内转移能力；以及采用双向导CRISPR干扰（CRISPRi）文库进行体内功能性筛选，以寻找RBMS3下游的关键效应因子。临床关联分析则利用了来自METABRIC、TCGA等多个大型乳腺癌患者队列的基因表达与生存数据。
研究结果
不同乳腺癌细胞系模型中的mRNA动态存在异质性
为了系统表征不同乳腺癌模型和亚型中的RNA动态，研究者在六个乳腺癌细胞系中进行了SLAM-seq分析。结果显示，虽然转录速率在很大程度上与基因表达相关，但仍有相当一部分基因表达的变异无法仅用转录变化来解释。通过估算mRNA降解速率，研究者发现尽管全局降解模式在不同细胞系间大体保守，但仍有数百个转录本的降解速率存在显著异质性。这种异质性具有亚型特异性，例如三阴性乳腺癌（TNBC）模型倾向于聚在一起。这些发现表明，mRNA稳定性是导致乳腺癌异质性的一个重要因素，并激励研究者进一步探索其背后的决定因素。
mRNA稳定性的整合建模揭示乳腺癌中的转录后异质性
为了理解导致乳腺癌异质性的转录后调控基础，研究者开发了GreyHound框架。该模型能同时评估转录本稳定性的顺式和反式决定因素。模型在测试集上表现出预测值与测量值之间良好的相关性。通过模型的解构分析，研究者识别出在不同细胞状态下重要性可变的序列区域和RBP。对高可变序列区域的基序分析，结合RBP表达与靶转录本稳定性之间的相关性验证，构建了一个RBP-基序相互作用网络。其中，RNA结合蛋白RBMS3与一个富AUA序列元件的相互作用被确定为最显著的调控轴。分析显示，携带该AUA元件的转录本（即RBMS3的候选调节子）的降解速率与RBMS3的表达呈负相关，这提示RBMS3可能是乳腺癌中mRNA稳定性的潜在调节因子。
RBMS3在乳腺癌细胞中作为mRNA稳定性的转录后调节因子发挥作用
在GreyHound模型预测的基础上，研究者通过实验探究了RBMS3的功能。在MDA-MB-231细胞中敲低RBMS3后，RNA测序分析发现，下调转录本的3‘非翻译区（3’UTR）中AUA基序显著富集。通过CLIP-seq实验，研究者在全基因组范围绘制了RBMS3的结合位点，证实其广泛结合于数千个位点，并且这些位点同样显著富集AUA基序。SLAM-seq实验进一步证明，在RBMS3敲低后，含有AUA基序的转录本降解速率增加，但其转录输出未受显著影响，表明RBMS3主要通过转录后机制发挥作用。使用双荧光报告系统的功能验证实验表明，RBMS3结合位点足以在RBMS3存在时增加报告mRNA的稳定性，而这种效应在RBMS3敲低的细胞中消失。这些结果证实了RBMS3通过直接结合AUA基序来稳定其靶转录本。
RBMS3及其下游靶向调节子与乳腺癌进展相关
基因集富集分析显示，RBMS3结合的转录本显著富集于血管内皮生长因子和转化生长因子-β等与癌症进展和转移相关的通路。对多个独立乳腺癌患者数据集（如METABRIC、TCGA）的分析表明，RBMS3表达降低与不良临床预后（如更短的无病生存期和无复发生存期）显著相关。此外，RBMS3表达在肿瘤组织中普遍低于配对正常组织，并且在更高分级、更晚分期的肿瘤中表达更低。这些关联在基底样/Claudin-low亚型中尤为明显，提示RBMS3在这些侵袭性更强的乳腺癌亚型中可能扮演重要角色。
调控RBMS3影响乳腺癌的转移和侵袭
为了在功能上建立RBMS3与转移潜力的因果关系，研究者在免疫缺陷小鼠中进行了肺定植实验。敲低RBMS3的MDA-MB-231细胞和HCC1806细胞在尾静脉注射后，表现出显著增强的肺转移定植能力。相反，在高度转移的MDA-LM2细胞中过表达RBMS3，则能显著抑制肺转移。细胞增殖实验排除了增殖差异的影响。进一步的Transwell侵袭实验表明，敲低RBMS3能显著增强MDA-MB-231细胞的侵袭能力。这些结果确立了RBMS3在三阴性/基底样乳腺癌中作为转移和侵袭抑制因子的角色。
RBMS3-TXNIP调控轴抑制乳腺癌转移
为了阐明RBMS3抑制转移的下游分子机制，研究者从RBMS3直接结合的靶基因中，筛选出可能介导其表型的关键效应因子。通过一个包含13个候选基因的双向导CRISPRi体内筛选，研究者发现，敲低硫氧还蛋白相互作用蛋白（TXNIP）能重现RBMS3敲低所导致的增强肺定植表型，而不影响细胞增殖。临床数据分析显示，TXNIP表达与RBMS3表达呈正相关，并且低TXNIP表达同样与患者的不良预后相关。在一组90例乳腺癌临床样本中也证实，TXNIP表达随着疾病分期进展而降低。机制上，RBMS3敲低导致TXNIP的mRNA和蛋白水平下降，但并未影响其前体mRNA水平，表明RBMS3在转录后水平调节TXNIP。体内 epistasis（上位性）实验表明，同时敲低RBMS3和TXNIP并未产生比单独敲低TXNIP更强的转移表型，支持了RBMS3位于TXNIP上游、通过稳定TXNIP mRNA来发挥抑制转移作用的线性调控模型。
研究结论与意义
本研究通过整合mRNA稳定性测量、可解释的深度学习建模和多层次的实验验证，系统性地揭示了一个此前未知的乳腺癌转移抑制调控程序。核心结论是：RNA结合蛋白RBMS3通过直接结合靶转录本3‘UTR中的AUA富集序列元件，发挥稳定mRNA的作用。其关键下游靶点TXNIP的mRNA稳定性受到RBMS3的正向调控。RBMS3-TXNIP轴的破坏（表现为RBMS3或TXNIP的低表达）会导致转录本不稳定、TXNIP水平下降，进而显著增强乳腺癌细胞的侵袭能力和体内转移定植潜力，且与乳腺癌患者的不良临床预后密切相关。
这项研究的重大意义在于多个层面。首先，在生物学机制上，它明确地将一个具体的转录后调控机制（由特定RBP介导的mRNA稳定性调节）与癌症转移这一致命表型直接联系起来，为理解乳腺癌转移的分子基础增添了新的维度。其次，在方法论上，成功开发的GreyHound框架提供了一种强大的、可推广的研究范式，能够从复杂的多组学数据中整合顺式与反式信息，发现疾病背景下具有功能意义的RNA-蛋白质互作网络，这对于研究其他依赖转录后调控的复杂疾病具有重要参考价值。最后，在临床转化潜力上，RBMS3和TXNIP的表达水平显示出显著的预后价值，有望作为识别高转移风险乳腺癌的生物标志物。虽然目前的研究尚未将其确立为直接的治疗靶点，但阐明的RBMS3-TXNIP轴为未来开发旨在恢复该通路功能、抑制转移的新型治疗策略提供了坚实的理论基础和精准的分子靶点。
