碱基日读 每天读完真实文献,只讲有用的那一段

从单细胞扰动预测到多发性硬化免疫分型,计算方法在三类数据集上的能力边界

2026-09-10·碱基·生信方法

要点 预测模型保守偏差、跨细胞类型任务全线失守;实验锚定越强,计算结论才越可信。

三篇文章都在用计算方法拆解生物系统在受到扰动之后的反应。扰动的形式各不相同——单基因敲除、化合物处理、热旱双重胁迫、体内抗体清除疗法——但它们在同一个问题上汇合:计算分析能走多远,在哪里开始失守。

25个数据集、24个指标构建的评测框架

Science Advances 这篇基准测评选了13种代表性的单细胞扰动响应预测方法,在25个跨物种、跨扰动模态的数据集上测了三类任务:预测未见单基因扰动的表达结果(任务一)、预测组合扰动的相互作用(任务二)、跨细胞类型迁移(任务三)。数据集包含两个原代免疫细胞药物响应资源,覆盖了和临床最相关的场景(PMID: 42715312)。

24个评估指标有意设计成四个层次:表达水平准确性、相对变化量(delta)、差异表达(differential expression, DE)恢复、分布相似度。一个方法在表达水平上看起来很好,在DE恢复上可能毫无信号。这篇文章的核心贡献之一,正是厘清了不同指标探测的是不同层次的生物学问题,而不是对同一个预测问题的冗余总结——这个区分在方法选型上有直接影响。

效应量与所有模型共有的保守偏差

结果里有两条需要分开理解。

小效应扰动下,模型的表达水平一致性最高。原因不是模型好,而是扰动本身接近对照,”几乎没变”本来就是合理预测,任何倾向保守输出的方法都能在这里拿高分。大效应扰动下,delta指标和DE指标才开始有区分力(PMID: 42715312)。

所有13种方法共享一个系统性特征:保守偏差(conservative bias)。其中,微调过的基础模型(fine-tuned foundation model)还进一步压缩方差,对组合扰动中的协同效应系统性低估。这直接影响药物组合筛选——你最不想低估的,恰恰是协同效应最强的那些组合(PMID: 42715312)。

PerturbNet在任务一和任务二的DE特征恢复上表现最好;但任务三——跨细胞类型迁移——没有任何一种方法能稳定泛化,生物一致性在这里压倒了所有技术差异(PMID: 42715312)。这是当前方法集体的硬墙,不是某种方法的特有缺陷。任何声称能跨细胞类型泛化的单细胞扰动预测结果,在独立验证之前,应该先当作待验假说。

双算法流程能走到的边界

水稻热旱双重胁迫的circRNA研究,走的是流程驱动的重分析路线(PMID: 42711459)。作者重新分析了已有的链特异性RNA-seq数据,用CIRI2和CIRCexplorer2双算法联合预测,鉴定出208条分布在全部12条染色体上的高置信circRNA(circular RNA,环形RNA)。Junction reads分析给出了环化强度的梯度谱,从高丰度强回剪位点到被线性转录本背景掩盖的低可信候选,分层清晰。基因组注释显示这些circRNA主要来自外显子区和基因间区,有明显的负链偏向,部分基因还通过可变回剪(alternative back-splicing)产生多个circRNA同工型。

比较分析中,83条仅见于对照,51条仅见于胁迫,74条两种条件共有,差异表达分析最终筛出7条热旱组合胁迫下特异诱导的circRNA(PMID: 42711459)。这7条是目前该方向上可信度最高的实验入口,优先验证时应从junction reads丰度高的候选里挑,而不是从208条里随机取。

但文章随后推断的ceRNA(competing endogenous RNA)调控网络,将osa-miR414、osa-miR1439、osa-miR2919指认为关键拓扑枢纽,靶基因关系和网络拓扑全部基于预测,尚未经实验验证,作者自己也明确说这需要后续实验落地(PMID: 42711459)。

用单细胞扰动评测的结论来校准:连有大量训练数据的单细胞模型,跨细胞类型泛化都全线失守;在数据更稀薄的植物非编码RNA调控网络里,纯计算预测的置信度只会更低。把双算法流程推到数据上限本身是贡献,问题在于过度解读。

ocrelizumab作为体内探针发现的B细胞依赖T细胞状态

多发性硬化(multiple sclerosis, MS)研究给出了一个不同层次的范例(PMID: 42714077)。作者把抗CD20疗法ocrelizumab明确定位为”体内扰动探针”,在开始治疗的患者中纵向采集高维多模态免疫数据,结合机器学习流程解析连续免疫细胞状态的变化。连续状态,而非传统离散分群——这个方法论选择让分析得以捕捉跨越分群边界的共享生物学过程,也是这套流程区别于常规scRNA-seq分析的地方。

他们找到一个慢性激活、促炎-细胞毒性、具有CNS归巢特性的T细胞状态,B细胞清除后被选择性削减。这个状态在多个独立数据集中得到印证:在临床侵袭性复发患者外周血中富集,在CSF来源的克隆性扩增T细胞里可见(其中一个亚群具有已证实的EB病毒特异性),也出现在终末期MS慢性炎性病灶边缘(PMID: 42714077)。这条证据链同时把B细胞疗法奏效但机制不清楚的局面推进了一步——B细胞的作用可能不是直接的,而是通过维持这个病理性T细胞状态间接实现的。

natalizumab实验提供了关键的交叉验证。这个靶向整合素α4(integrin α4)、阻断淋巴细胞跨血脑屏障迁移的疗法,反而导致这个T细胞状态在外周血中富集。逻辑自洽:ocrelizumab切断B细胞对这个T细胞状态的激活信号,状态减少;natalizumab封住CNS入口,状态无法入脑而积压在血里。两种疗法作用于同一轴线的不同节点,效果方向相反,机制完全一致(PMID: 42714077)。

纵向数据来自接受治疗的患者,不是随机对照,建立的是强关联而非完整因果链。这条局限需要记在心里。但多数据集独立验证加双疗法交叉印证,已经把证据强度推到了同类研究的高位。

三道关口

把三篇放在一起,能整理出一条实用的判断线索。

效应量是第一道。扰动效应小时,表达水平指标没有区分力,任何保守方法都能拿高分;大效应下delta和DE才开始有意义(PMID: 42715312)。做扰动分析之前先评估效应量,再决定信哪套指标。

跨细胞类型是第二道,目前全线失守。训练集与测试集细胞类型不同,当前所有方法都说不准(PMID: 42715312)。水稻研究里从已知胁迫条件外推到完整调控网络,也在同一道关口上(PMID: 42711459)。

实验锚定是第三道,也是最能拉开差距的一道。MS研究展示了当前能达到的最高可信度:纵向采集、多数据集独立验证、两种疗法交叉印证,让分析从产生假说升级到支持干预设计。作者指出,直接操控这个T细胞状态可能是靶向慢性神经炎症的关键机会,也可能为减少大范围免疫细胞清除策略提供依据(PMID: 42714077)。

水稻那篇的7条差异表达circRNA是当前最高置信的实验入口,三个miRNA枢纽是待验假说;MS那篇发现的T细胞状态,验证链条已足够长,直接干预实验是下一步。

单细胞扰动预测基准评测circRNA多发性硬化T细胞状态

同主题文章

相关工具

多重检验校正在线计算:BH (FDR)、BY 与 Bonferroni

粘贴一列 p 值,一次得到 BH 的 q 值、BY 与 Bonferroni 校正值,三者并列对照。

log2FC 与倍数变化 (Fold Change) 在线互转

log2FC ↔ 倍数、上调下调百分比,附常用阈值对照表。

米氏方程 Km 与 Vmax 计算(非线性拟合 + 三种线性化对照)

输入若干组底物浓度与反应速率,用非线性拟合求 Km 与 Vmax,并把 Lineweaver-Burk、Hanes-Woolf、Eadie-Hofstee 三种线性化的结果并排列出——它们的分歧本身就是数据质量的读数。

DNA 连接反应计算器(载体与插入片段用量 / Ligation)

输入线性化载体的用量与大小、插入片段大小及目标摩尔比,自动计算连接所需插入片段质量与各组分 fmol 数,适用于 T4 连接酶黏端和平端克隆。

没有广告,不收集数据,每个公式都验算过。维护它需要持续投入——如果帮你省了时间,请杯咖啡。
请杯咖啡 了解更多 →