2019年5月 week2 文献学习: Identifying cell populations with scRNASeq

2019-05-19 week2文献学习

纠结着要看什么文献,因为最近要开始准备托福考试,然后很想放弃一周一篇的文献学习,想想还是要学习的呢。然后刚好单细胞天地推送的一篇文献 Identifying cell populations with scRNASeq,就刚好拿起来读一读,加上最近自己用seurat 分析数据,也算是总结一下。

摘要

scRNASeq已成为量化单个细胞转录组的有效方法。 然而,来自scRNASeq 实验的数据通常既有噪声又高维度,使得计算分析变得非常重要。 在这里,我们提供了不同实验方案的概述和计算分析的最流行的方法,主要专注于识别生物学重要基因的方法,将数据投射到较低维度并将数据聚类到假定的细胞群体中。 最后,我们讨论了鉴定细胞类型或细胞状态的验证和生物解释的方法。

文章主要综述了四方面的内容
1)实验方面:scRNAseq实验设计以及实验过程中产生的噪音;
2)如何解决数据的高纬度问题:一是降维,二是进行特征选择
3) 如何将数据聚类到假定的细胞类群
4) 如何鉴定细胞类型或者细胞状态

文章主要内容

一、实验方面

1. scRNAseq主要是分成三个方面

1)单细胞分离

1) 细胞分离的过程需要将样品进行解离,然后分选到PCR的单独孔中,或者利用droplets或者microwells或者微流控捕获单个细胞
2)文库制备
文库制备需要反转录和扩增
3) 测序

2)分选细胞的方法

对于一次性要获得很多细胞的研究,可以使用基于液滴的方法,比如In-drop, drop-seq以及10X Chromium,可以一次性捕获成千上万细胞。确定细胞类型所需的最低测序深度可以为25000-50000reads/cells。或者使用微孔的方法替代droplet的方法。但是这些方法需要在细胞裂解前加上barcode,因此只支持3/5端测序。
如果细胞量不大的话,可以使用流式分选,这些方法捕获细胞的性价比比较高,检出率也比较高,不仅可以支持3/5端测序,也可以支持全长转录本测序

===收获:明白两种不同的分选细胞的策略,然后可以支持的测序方法也是不同的,才明白了实验室师姐所用的分选细胞的方法。

3)单细胞测序中存在的问题

1)Doublet的问题
Doublet是指一个液滴或者微孔中存在2个或者多个细胞,会对细胞类型的鉴定过程产生影响。
2)batch effect
批次效应是不同时间或者不同人员制备的实验重复之间的实验效率或者不同细胞状态而产生的。可以使用一些算法比如quantile,SVA的Combat, RUVs去除批次效应。

想要消除批次效应需要通过仔细的实验设计,将每个实验条件分散到各个实验批次中,做到一视同仁。【不大理解这样真的能消除批次效应吗??】

2.scRNAseq中的技术噪音

解决办法:UMIs和已知浓度的外源RNA分子(spike-ins)

1)UMI

UMI是反向添加到每个cDNA的5或者3端,长度是4-10bp的分子,作用是区分哪些reads来自于同一个cDNA分子,然后用来估算原始的分子数量。但是存在一个问题UMI和转录本的5/3端结合后进行测序,会存在丢失isoform信息的情况,捕获的遗传变异较少等,评价等位基因较难等。

2)spike-ins

标准的spike-ins 是ERCC组织指定的一段细菌序列。存在的问题是:ERCC的捕获效率低于内源性mRNA,具有较高的技术变异性,有时比内源性基因的含量还多,spike-ins的计数受到生物条件的影响,因此有时候会失去作为control的优势。如果使用spike-ins,比对应该将spike-ins序列和参考基因组序列合并作为共同参考基因组。

3)多重测序

多重测序也是产生技术噪音的一个原因,因为它会导致不同细胞之间的reads不在一个层次上,使用标准化可以纠正不同细胞之间的的测序深度的影响。
SCnorm可以解释不同测序深度对基因不同表达水平的影响。

3)收获

== 还是不明白spike-ins是啥,作用是啥,不是很明白,是为了作为对照吗?但是知道了如果使用spike-ins,比对应该将spike-ins序列和参考基因组序列合并作为共同参考基因组。
== SCnorm可以解释不同测序深度对基因不同表达水平的影响。

二、数据维度的处理

1.降维

方法主要有PCA(线性降维),tSNE(非线性降维,主要用于数据的可视化),DM方法(非线性降维,主要用于分析细胞的连续发展)

2.特征选择

移除掉信息量较少的基因,减少数据噪音,提高数据运行速度
M3Drop;HVG(High variabel genes);spike-in based methods; correlated expression.

一般的话是先进行特征选择,然进行降维,然后再进行细胞聚类

三、非监督聚类鉴定细胞类型

方法主要有k均值聚类,层次聚类,基于密度的聚类方法,图聚类方法。
聚类中最关键的一个点在于要聚多少类。
许多的聚类工具可以在ASAP中找到的~,是个web工具。

四、细胞类型的鉴定

1)计算角度
对数据运用多种聚类方法然后比较得到一个比较好的聚类结果。
计算方法主要是为了提高结果的可靠性,然后细胞类型的自动化鉴定未来我觉得可行的,并且是流行的。

2)实验角度
主要是基于Marker的方法~

五、总结

这篇文章算是帮我把分析思路再次理清了下,这次主要的收获在于对于UMI有了比较深的一个理解,对于特征选择和降维两者之间有了比较清楚地认识,这两种方法主要都是用于对于数据维度的处理,只是处理方法不一样,特征选择的话主要是移除掉信息量较低的基因。

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