单细胞测序专题集合单细胞转录组生信工具

velocyto||sc-RNA速率:一种细胞轨迹推断方法

2019-10-25  本文已影响0人  周运来就是我
18年nature发了一篇单细胞方法类文章,讲得就是如何利用RNA velocity来做细胞发育路径的推断。

RNA速度详细显示了神经元和其他的细胞如何在大脑发育和成熟时获得它们的特定功能。我们特别兴奋的是,这种新方法有望协助揭示出大脑通常如何发育,同时也为破解人类大脑发育障碍(比如精神分裂症和自闭症)中发生的问题提供线索。
---卡罗林斯卡医学院分子系统生物学教授Sten Linnarsson

RNA丰度是单个细胞状态的有力指标。单细胞RNA测序可以定量、准确、灵敏、高通量地揭示RNA丰度。然而,这种方法只能在某个时间点捕获静态快照,这对分析胚胎发生或组织再生等随时间变化的现象提出了挑战。在这里,我们证明了RNA的速率——基因表达状态的时间导数——可以通过区分普通单细胞RNA测序中未剪接和剪接的mrna来直接估计。RNA速率是一种高维的载体,可以在数小时的时间尺度上预测单个细胞的未来状态。我们在神经嵴谱系中验证了它的准确性,展示了它在多个已发表的数据集和技术平台上的应用,揭示了发育中的小鼠海马的分支谱系树,并检测了人类胚胎大脑中的转录动力学。我们期待RNA速度将极大地帮助分析发育谱系和细胞动力学,特别是在人类。

在发育过程中,分化发生在数小时至数天的时间尺度上,这与典型的mRNA半衰期相当。初生(未剪接的)和成熟(剪接的)mRNA的相对丰度可以用来估计基因剪接和降解的速率,而不需要进行代谢标记,如 bulk显示的那样。我们推断在单细胞RNA测序(RNA-seq)数据中可以检测到类似的信号,并可以揭示动态过程中整个转录组的变化速率和方向。

所有常见的单细胞RNA-seq协议都依赖于 oligo-dT 引物来富集聚腺苷酸化的mRNA分子,检测基于SMART-seq2的单细胞RNA-seq数据集C1, inDrop和1 10x Genomics Chromium protocols,我们发现15-25%的reads包含未剪接的内含子序列。单细胞水平(~20%) 与bulk4中先前的观测结果一致(14.6%)。大多数这样的读数来自于内含子区域内的二次启动位置基因组铬文库,我们还发现了大量不一致的启动从更常见的内孔- polyt序列,这可能是在PCR扩增过程中,通过对第一链cDNA引物产生的。

为了量化前体的abun- dance与成熟mRNA之间的时间依赖关系,我们假设了一个简单的转录动力学模型,其中剪接mRNA丰度(RNA速度)的一阶导数由未剪接mRNA产生的剪接mRNA与mRNA降解的平衡决定。在动态过程中,转录速率的增加导致unspliced mRNA的快速增加,其次是增加subse——quent拼接mRNA,直到达到一个新的稳定状态。

原理优点复杂,用法很简单,安装完之后(velocyto.py||install
按照官网示例教程一步一步走下去就好了。

提供了R语言版的和Python版的,可以根据自己的喜好来安装使用。

Usage: velocyto [OPTIONS] COMMAND [ARGS]...

        Options:
          --version  Show the version and exit.
          --help     Show this message and exit.

        Commands:
          run            Runs the velocity analysis outputting a loom file
          run10x         Runs the velocity analysis for a Chromium Sample
          run-dropest    Runs the velocity analysis on DropEst preprocessed data
          run-smartseq2  Runs the velocity analysis on SmartSeq2 data (independent bam file per cell)
          tools          helper tools for velocyto


velocyto
velocyto-notebooks
RNA velocity of single cells
velocyto.py||install
RNA velocity | RNA速率
Nature:测量RNA速度可预测单个细胞的未来状态和最终命运
What happens to the mRNA after it has completed its job at the ribosome and the protein has been released?
Database for mRNA Half-Life of 19 977 Genes Obtained by DNA Microarray Analysis of Pluripotent and Differentiating Mouse Embryonic Stem Cells
Messenger RNA Half-Life Measurements in Mammalian Cells

上一篇下一篇

猜你喜欢

热点阅读