单细胞挖掘思路单细胞流程

单细胞笔记15-轨迹推断之RNA velocity

2023-03-20  本文已影响0人  江湾青年

之前大多数 single cell 技术都是静态的刻画 cell state(即都是snapshot),没有能够引入细胞随时间的变化,living cell imaging tracking 的技术虽然能提供时间信息,但是能够覆盖到gene 数量是非常少的(通常也就不到10)。而发育分化过程中 cell state 是动态变化的,因此怎么能够理解 cell state transition,怎么将这些 snapshot 拼成 video 就是一个很有意思且很重要的问题。

18年 “RNA velocity of single cells” 首先引入 RNA velocity 这个概念,充分利用现有的数据来得到全新的概念。核心的思想就是,利用 splicing 这个生物过程来引入时间这个要素,从而推断出 cell state 的 velocity(速率)

RNA velocity

具体是什么意思呢,首先从一个例子说起,当我们看到下图坤坤跳舞的中间状态的时候,我们ikun脑海中是可以串联起这几个动作的。我们可以从其中一个 snapshot 知道坤坤的运动状态,并且预测运动下一秒的状态。本质原因是因为我们对跳舞这个过程有一个动力学的认识。

image.png

那么我们回到细胞内的转录,我们只需要通过中心法则中依赖于时间的生物过程,并对这个过程进行动力学刻画即可以得到这些 snapshot 的未来方向,那么应该使用什么生物过程呢?首先需要满足两个条件,1. 依赖于时间,2. 现有的数据中应该可以直接得到

答案就是 splicing! 首先 splicing 肯定依赖于时间,其次现有的测序 reads 天然的可以 mapping 到 intro 和 exon, 从而得到 unspliced 和 spliced 的RNA。

Gene 在转录成 mRNA 的过程中会经过 splicing,而 splicing 的动力学我们是通过 ODE 来 model 的,因此我们便可以通过 unspliced mRNA (u)和 spliced RNA(s) ,以及动力模型( ODE model) 来计算 cell state 的 velocity 。

v2-58e3f688b1eaf1e8cd3d20ec36112672_1440w.webp.png

总结:RNA velocity通过计算剪接和未剪接的mRNA的比率,来恢复定向的动态信息。这种方法不仅量化基因活性,而且量化它们在单个细胞中的变化(RNA 速率),开辟了研究细胞分化的新方法。


参考

https://zhuanlan.zhihu.com/p/351104418
https://zhuanlan.zhihu.com/p/359100974

上一篇下一篇

猜你喜欢

热点阅读