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

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,直到达到一个新的稳定状态。

  • a 剪接计数和未剪接计数是通过包含内含子序列的单独计数读数来估计的。对于基于唯一分子标识符(UMI)的协议,与给定分子相关的多次读取被分组(带有星号的方框)。饼状图显示了未拼接分子的典型部分。

  • b 转录模式的动态,捕捉转录(α),拼接(β)和降解率(γ)参与生产unspliced (u)和拼接(s) mRNA的产物。

  • c 模型的解决方案b作为时间的函数,显示unspliced和拼接动力学反应步骤中α的变化。

  • d ,相图显示与c相同的解(固体曲线)。稳定状态不同的转录率α值落在坡的对角线γ(虚线)。未剪接的mRNA水平高于或低于这个比例表示增加(红色底纹)或分别减少(蓝色底纹)基因的表达。24小时内小鼠肝脏中环状rna相关基因的剪接(s)和非剪接(u) mrna的丰度。未剪接的mRNA可预测下一次剪接的时间点。

  • f,g 观察到一对环状驱动基因的相位图:Fgf1 (f)和Cbs(g)。每个点的昼夜节律时间用时钟符号表示(与e中的对应)。冲对角线显示了稳定的关系,为预测的γ

  • h,表达状态在未来时间t的变化,如模型预测的那样,显示在前两个主要成分(PCs)的空间中,概括了沿着昼夜节律周期的进展。每个圆圈显示观察到的表达式状态,箭头指向未来状态的位置,从速度估计中推断出来。

原理优点复杂,用法很简单,安装完之后(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

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 212,029评论 6 492
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,395评论 3 385
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 157,570评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,535评论 1 284
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 65,650评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 49,850评论 1 290
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,006评论 3 408
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,747评论 0 268
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,207评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,536评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,683评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,342评论 4 330
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,964评论 3 315
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,772评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,004评论 1 266
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,401评论 2 360
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,566评论 2 349

推荐阅读更多精彩内容

  • 知识的学习没有一蹴而就,没有捷近,扎实的学习是唯一的捷近。 一篇RNA-seq分析流程的综述,全面而详细!深度好文...
    dandanwu90阅读 48,271评论 5 152
  • 欢迎关注oddxix RNA-seq是高通量测序中最常见的一种应用,本期视频介绍其:1.方法原理2.生物信息分析 ...
    oddxix阅读 4,682评论 1 44
  • 记录: 昨天和教育负责人交流后,收获真的很多,其中最大一点是术有专攻! 要很清楚你要教哪个科目,甚至是哪个版块,如...
    新灵清阅读 172评论 0 2
  • 我好像蛮喜欢简书新版的界面蛮好的的。 正文分割线 曾经有一个人在网上围棋室下赢了所有世界高手,让所有人都吃惊和着迷...
    止行阅读 485评论 0 0
  • 实现效果:新闻列表标题,在标题最后一个字符加已读未读标签。 用富文本标签实现的,在富文本标签里面添加文本附件图片插...
    骑马纵天下阅读 1,710评论 2 1