查看原文
其他

系统学习单细胞转录组测序scRNA-Seq(一)

豆豆花花 生信星球 2022-06-07

 今天是生信星球陪你的第308天


   大神一句话,菜鸟跑半年。我不是大神,但我可以缩短你走弯路的半年~

   就像歌儿唱的那样,如果你不知道该往哪儿走,就留在这学点生信好不好~

   这里有豆豆和花花的学习历程,从新手到进阶,生信路上有你有我!

刘小泽写于19.3.15

好了,终于可以学新知识了,单细胞测序是目前研究的热点,那么我重新回来探究一下,首先问自己几个问题

都是转录组测序,有什么不一样?

看到转录组我们都不陌生,有参无参也都做了许多了,但是对一个细胞测序,为什么呢?10X给出了一个说明

https://www.10xgenomics.com/solutions/single-cell/

image.png

传统的是测一群细胞(Bulk RNAseq),然后我们定量是基于所有细胞内的基因平均表达量进行统计,这样得到的是组织之间的差异性,也就是说只能看到大体的趋势,但是对于继续深入研究还是存在比较大的背景噪音。对于癌症或者系统发育相关的研究,往往是想看到某些特殊的细胞发挥的作用,也是相当于一步步筛选。

单细胞更加精准,最终得到的是每个细胞型的表达量,从而可以帮助理解组织异质性,鉴定罕见的细胞类型,检测细胞组分的变化,例如在癌症异质性、胚胎发育、细胞对药物的反应、神经细胞分型、差异表达通路等方面都发挥重要作用

怎么去做单细胞研究?

我们常规的转录组测序,需要提取样本RNA、建库测序,但是对于细胞来讲,材料提取是个大问题,但是目前主要有两个获取细胞的方法:Smart-seq和Drop-seq

Smart-seq的特点就是测序细胞量少(因为一次通过一个细胞),但是测序的reads数多;Drop-seq(10X 公司的方法)一次可以测许多细胞,但是测得reads数都不是很多

image.png
image.png

获得了样本,一般公司会直接给你表达矩阵,也就是一般的转录组测序利用featureCount获得的结果,自己需要先进行归一化(例如CPM)

image.png

当然,这仅仅是许多资料中的一个流程

  • For differential expression:

  • –   SCDE: http://pklab.med.harvard.edu/scde/index.html

  • –   SCA: https://github.com/RGLab/SingleCellAssay

  • –   MAST: https://github.com/RGLab/MAST

  • –   SAMseq: http://cran.r-project.org/web/packages/samr

  • For clustering etc.:

  • –   Monocle: https://github.com/cole-trapnell-lab/monocle-release

  • –   Rtsne: http://cran.r-project.org/web/packages/Rtsne

  • –   Sincell: http://master.bioconductor.org/packages/devel/bioc/html/sincell.html

  • –   scLVM: https://github.com/PMBio/scLVM

  • –   BASiCS: https://github.com/catavallejos/BASiCS

  • –   Pagoda: http://pklab.med.harvard.edu/scde

  • –   Seurat toolkit: http://www.satijalab.org/seurat.html

  • –   Sincera pipeline: https://research.cchmc.org/pbge/sincera.html 


点击底部的“阅读原文”,获得更好的阅读体验哦😻

初学生信,很荣幸带你迈出第一步。

我们是生信星球,一个不拽术语、通俗易懂的生信知识平台。由于是2018年新号,竟然没有留言功能。需要帮助或提出意见请后台留言、联系微信或发送邮件到jieandze1314@gmail.com,每一条都会看到的哦~



您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存