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

细胞解卷积十种算法评测:参考集错配让最好的方法r值跌到0.17

2026-11-02·碱基·生信方法

要点 解卷积算法选型不能只看基准分数,参考集与样本的生物学匹配程度才是真正的性能上限。

参考集的选择比算法本身更重要,但在多数论文的方法部分,这句话只占一行。

批量RNA测序(bulk RNA-seq)样本量大、成本低,代价是把组织里所有细胞的信号混在一起。细胞解卷积(cellular deconvolution)是一条出路——用单细胞RNA测序(scRNA-seq)的参考数据反推各细胞类型在 bulk 样本里的比例,绕开物理分选。这个方向近年方法涌现,系统横向比较一直不够。

MuSiC、DWLS 和 BayesPrism 在基准条件下领跑

Malik 等人对十种主流的基于参考集的解卷积算法做了系统评测,包括 MuSiC、DWLS、BayesPrism、CIBERSORTx、SCDC、BisqueRNA、DISSECT、TAPE、Scaden 和 scpDeconv,用两个独立数据集跑了六组实验(PMID: 42762696)。选这十个的理由是架构多样且使用普遍,传统统计模型和深度学习方法都有代表。

基准条件下,MuSiC、DWLS 和 BayesPrism 表现最强,每样本均值 Pearson 相关系数 r 大于 0.95、Lin 一致性相关系数(CCC)大于 0.95(PMID: 42762696)。深度学习方法整体变异更大。

这三组数字看起来很漂亮。问题在“基准条件”这四个字里。

参考集限定为单一发育阶段后 r 最低跌至 0.17

真实分析中,参考 scRNA-seq 数据和待分析的 bulk 样本往往来自不同实验室、不同发育阶段或不同疾病状态。Malik 等人专门模拟了这种错配:把参考集限定为单一发育阶段,看各算法性能怎么变。

各方法均值 Pearson r 落在 0.17 到 0.58 之间(PMID: 42762696)。

基准 r > 0.95,参考集错配后最低跌到 0.17。不是渐变,是断崖。

BayesPrism 在参考集错配实验里平均鲁棒性最高;测序深度那组实验里,SCDC 和 DWLS 在四个深度水平上最稳定(PMID: 42762696)。没有哪个算法在全部六组实验里都排第一。

这意味着方法部分写“我们使用 MuSiC 进行解卷积”本身信息量不足。参考集来自哪个组织、哪个发育阶段,与样本的生物学距离有多远,才是决定结果可不可信的那几行字。

这份评测的代码已公开在 GitHub 上,做到了基本可复现。但两个数据集覆盖的生物学范围有限,结论能否推到更复杂的组织类型或疾病状态,还需要更多验证。

TopClust 用拓扑数据分析绕开参数主观性

从 bulk 解卷积往上走一层,scRNA-seq 数据本身的聚类分析面临一个结构性问题:大多数方法需要用户手动指定参数,不同参数导致不一致的细胞类型划分,可重复性差(PMID: 41981916)。

Hu 等人提出的 TopClust 用高维密度图替代主观的参数选择。具体做法:先用共享最近邻模型(shared nearest neighbor model)构建高维密度图,再用拓扑数据分析(topological data analysis)在图上找局部极大值。每个极大值是一个簇的顶峰细胞(summit cell),顶峰附近密度较高的核心细胞群是峰值细胞(peak cells)(PMID: 41981916)。聚类的数量和边界由拓扑特征决定,不需要调参。

在肺癌细胞系和荧光激活细胞分选(fluorescence-activated cell sorting,FACS)分选的免疫细胞两个基准数据集上,TopClust 的结果与已知细胞类型吻合良好(PMID: 41981916)。单独分析峰值细胞,可以减少刻画每个簇的转录变异,改善细胞类型注释的准确度。

方法还提供了一个差异分析工具:计算顶峰细胞对之间的方向向量,把基因投影到这个向量上,识别最能区分两类细胞的基因,作为传统差异表达分析(differential expression analysis)的补充(PMID: 41981916)。

验证目前只在两个相对简单的细胞体系上做过。在细胞类型高度复杂的组织或罕见细胞比例极低的情况下是否同样可靠,还需要更多数据支持。

胞内凝胶化在单细胞层面的均一性挑战

同期一篇综述从湿实验方向切入细胞状态调控。Chen 等人梳理了胞内凝胶化(intracellular gelation)技术——在细胞质内原位构建三维聚合物网络,形成“人工细胞骨架(artificial cytoskeleton)”,直接调控细胞内的物理化学微环境(PMID: 42746902)。相比从细胞外部作用的传统生物材料,这种“由内向外”的策略可以在不改变外部环境的前提下干预细胞状态。综述覆盖的范围包括触发机制、材料体系、生物学效应和医学应用。

然而单细胞层面的异质性同样是这条路线的障碍:各细胞之间凝胶化效率不一致。加上目前缺乏体内靶向技术,绝大多数工作停留在体外;人工细胞骨架的长期安全性也尚无定论(PMID: 42746902)。均一性、靶向性、安全性,这三条限制都是可靠应用这项技术的前提,目前都没有解决。

scRNA-seq细胞解卷积算法评测拓扑数据分析bulk RNA-seq

同主题文章

相关工具

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

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

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

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

DNA 序列翻译成蛋白质(六框翻译)在线工具

标准遗传密码表,支持单框与六框翻译、ORF 查找。

质谱质量误差计算器:ppm 误差与质量窗口在线换算

输入实测与理论 m/z,立刻得到 ppm 误差和 mDa 偏差;或由理论质量和允许 ppm 范围反算 XIC 质量提取窗口,适用于 Orbitrap、Q-TOF 等高分辨质谱日常数据分析。

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