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

从UMI碰撞到植物细胞壁难题:三项组学研究指向的同一个方法学教训

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

要点 组学流程假设失效时仍在产出数字,使用前先把前提假设对照自己的数据检查一遍。

工具超出说明书范围之后,通常还是会给你一个数。RNA-seq的UMI去重流程、蛋白质组学的内标追踪算法、植物空间转录组的分析管线,失效模式各不相同,共同点只有一个:流程不会告诉你假设已经失效。这正是以下三篇工作各自面对的问题。

UMI碰撞与矩估计矫正

UMI(unique molecular identifier,唯一分子标识符)是RNA-seq建库时给每条原始转录本打上的随机序列标签,PCR扩增后靠它区分来自不同原始分子的读数。这套逻辑的前提是每条转录本拿到的标签在当前批次里是唯一的。

序列越短,可能的组合数越少。文库量较大时,不同转录本碰到同一个标签(碰撞)的概率就不可忽视。标准去重流程把碰撞的两条当一条,系统性地低估基因表达量,而且流程不报错,输出的计数矩阵看起来和没有碰撞时完全一样(PMID: 42741932)。

Agyemang等的解法是矩估计量(method-of-moments estimator)。UMI的经验分布是非均匀的,有些序列因碱基组成偏差而更常出现,利用这一特性可以反推碰撞造成的低估幅度并矫正。作者证明,配合矫正估计量,较短的UMI定量结果与长UMI相当,同时降低测序和合成成本(PMID: 42741932)。

这项工作以数据集分析和模拟为主,不同平台、建库方案、细胞类型下碰撞率的系统性差异还需要更多实验验证。在采用更短UMI之前,先在自己的文库规模下估算一下碰撞率,是合理的谨慎。

BrScan-MS1与MS1层溴素指纹追踪

蛋白质组学(MS-based proteomics)实验里,向样品中加入内标肽段是校准保留时间、监测仪器漂移的标准做法。传统方法靠MS/MS鉴定来追踪内标,前提是内标信号强到足以触发二级碎裂。

样品量极少时(单细胞蛋白质组)或基质复杂时(血浆),内标在MS/MS层常常消失或不可靠,校准失败,仪器漂移就进了结果,而流程通常不会明确报告追踪失败(PMID: 42742076)。

Li等引入的BrIS(brominated internal standards,溴化内标)绕过了这一依赖。溴具有特征性同位素指纹,含溴肽段在MS1谱图上呈现出可辨认的特征信号。BrScan-MS1算法利用这个指纹,加上BrIS肽段保守的洗脱顺序,只在MS1层就能定位内标,不需要触发MS/MS。

在DDA稀释系列里,BrScan-MS1在每个含BrIS的重复样品中都识别出了全部内标肽段,常规MS/MS方法在低浓度下会丢失。血浆基质的长时间DIA分析里,BrIS保留时间漂移更小、定量更稳,均优于iRT肽段。涉及三个细胞系的单细胞蛋白质组实验里,低投入量下BrIS仍然有效,细胞系特异性的蛋白质组模式大体得以保留,对照组差异表达模式的保留也比iRT稍好(PMID: 42742076)。

从设计逻辑看,这不只是换了一组参考肽段,而是把MS1层的化学特征,也就是溴的同位素模式,做成了可用于追踪的信号,让内标追踪不再依赖需要触发鉴定才能完成的步骤。

植物细胞壁与空间转录组的系统性挑战

空间转录组(spatial transcriptomics)技术在保留组织结构的同时做基因表达测量,主流平台几乎全是为动物组织开发的。

植物打破了这些平台的前提假设。细胞壁刚性、化学组成多样,直接影响样品制备方式、细胞与亚细胞的分割精度,以及信号检测效率和多模态数据整合(PMID: 42740657)。结果是许多协议只针对特定物种或组织类型,换个材料就容易复现不了。

这篇综述(Jhu等,The Plant Cell,PMID: 42740657)把分散在各实验室的问题整理成了结构化的瓶颈图谱,涵盖样品制备、细胞与亚细胞分割、信号检测、多模态数据整合四个环节,并提出了社区驱动的协议共享与基准测试路线图。综述本身不提供新数据,这是它的局限,也是当前阶段的现实。方法尚未成熟,能说清楚问题在哪、哪些策略值得尝试,已经有用。

对打算做植物空间组学的读者,务实的预期是从样品制备开始,每个环节都可能需要重新优化。不报错不等于对,只等于还没发现错在哪。

矩估计与溴素内标共享的设计起点

这三件事共享一个结构。方法在假设成立时工作良好,假设失效后仍在产出数字,却没有任何提示。

UMI矩估计量和BrIS的解法,起点都不是换一个更好的算法,而是先把依赖的假设可见化,然后设计出不依赖它的路径。矩估计量把碰撞率量化出来;BrScan-MS1用MS1层的物理特征绕开了MS/MS鉴定这个前提。

植物空间组学还处于更早的阶段,连问题的全貌都在整理中。但揭示的东西一样。方法从设计对象迁移到新的系统时,失效不一定发出声音。

文献来源
  1. Unlocking the Full Potential of Spatial Omics in Plants: Practical Challenges, Solutions, and a Path Forward.,The Plant cell(PMID: 42740657)
  2. Unique Molecular Identifiers Don't Need to be Unique: A Collision-Aware Estimator for RNA-Seq Quantification.,Journal of computational biology : a journal of computational molecular cell biology(PMID: 42741932)
  3. BrIS: Visualizing Internal Standards at the MS1 Level via Bromine Fingerprint for Sensitive and Robust Proteomics Quality Control.,Angewandte Chemie (International ed. in English)(PMID: 42742076)

RNA-seq定量UMI去重蛋白质组质控空间转录组植物细胞壁

同主题文章

相关工具

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

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

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

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

凝胶条带大小估算(DNA bp / 蛋白 kDa,半对数标准曲线)

输入 marker 各带的大小与迁移距离,拟合半对数标准曲线估算未知条带;并指出两个常见错误——在两条带之间按大小线性内插(永远高估),以及对跑出 marker 范围的条带给估计值(那是外推)。

动物给药量计算器(mg/kg 剂量 → 注射体积 / 批量备药)

输入剂量 (mg/kg)、体重与药液浓度,算出注射体积;支持批量备药并给出各途径最大注射量参考。

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