拟时序分析

简要内容
--拟时序分析的定义
--结果解读
--------做出结果(cluster、state)
--------确定起始点
--------做start····end的图(并且映射到tsne中又一个图)
--------热图和感兴趣的-单个基因的动力学变化趋势图
--做拟时序分析的适合的情况
--------细胞图谱
--------case control的特征转变
--------重编程或转变轨迹
--文章实例
--------little tips(遇到三叉口的拟时序分析两种方法)


首先做细胞注释(在project中有了,不再赘述)

一、拟时序分析是什么

通过构建细胞间变化轨迹预测细胞随时间的变化
模拟时间的细胞序列的分析
时序中细胞发生的变化
也不一定非要设置不同的时间段去做实验
细胞本身存在拟时序变化,细胞是有变化的,才可以做拟时序分析

二、结果解读

image.png

主成分一和二,根据高可变基因进行排序,
每个点是一个细胞
左图按cluster排序,线就是进化轨迹
交点1、2,是有节点的(说明基因有表达量的差异)
右上按state排序,每个分叉一个state,一个表达模式
cluster:得到拟时序分析结果,将cluster映射到结果中
state:根据结果把每个分叉点做命名(表达模式不同),如下

image.png

计算机没法判断起始点,要通过生物学意义(基因表达和功能)去确定起始点,
法一:比如一个case,一个control(对照组),一部分枝杈上case的细胞数多(control少),一部分的枝杈上control的细胞数多(case少),那么起点就是control,细胞由control向case转化


image.png

法二:判断每个state或cluster中什么功能的细胞多,比如乃鹦鹉的到激活型的,,,,某个state中乃鹦鹉的细胞数目多,而另一个state中激活型的细胞数目多,则乃鹦鹉的那个state为起始点

判断起始点后,做下面这张图

image.png

这图与上面的那个cluster映射的图一样的形状,但是下面这个图时根据外面认为的确定一个起始点,颜色越黑的是起始点,越黄,随时间越靠后,然后将这个时间顺序映射到tsne中去,黑色的细胞不断往浅色的细胞转换,

image.png

时间顺序得到之后,接下来可以做热图和动力学变化趋势图

热图


image.png

image.png

热图最左边到最右边就是从起始点到后面的这个动力变化过程中基因的表达量(基因变化趋势)变化(最右边那一列是基因名),左边那一列是将所有的基因分为了三块(三种表达模式),绿色、澄色、蓝色,比如绿色的表达模式,那些高可变基因的表达量逐渐减少(表达量多少的图是右边的那个标签,红色最多,蓝色最少),澄色的部分的表达量也是随着时间的流逝,细胞转化过程中,逐渐减少。蓝色的部分是逐渐增多。
还可以针对感兴趣的某个基因做单个基因在里面的变化趋势


image.png

这个图中纵坐标代表表达量,图中每个点代表每个细胞,基因在这些细胞中的表达量
image.png

左上的基因趋势不变后来升高,

右上的基因先降低后升高(代表橙色部分)
左下基因一直下降(代表绿色部分)
右下基因一直升高(有可能对应热图中的蓝色部分)


image.png

三、文章中如何应用

1做细胞图谱(无组别)

比如细胞有两个亚型A 和 B,两者之间存在一个转变关系,那我们就去做一个拟时序分析,


image.png

2、control组的细胞往case组细胞特征转变

两个组别,研究药物处理或其他处理,导致表达谱发生了变化,一种细胞类型去做拟时序分析,否则无意义

3、某一个细胞类型下面的重编程或转变轨迹

比如T细胞的CD4的乃鹦鹉的细胞往TH的转变
比如CD8的乃鹦鹉T细胞往CTL毒性方面的转变
或者干细胞的重编程,人为的特定指定干细胞往哪个方向转变

这些方面都比较适合做拟时序分析

四、例子文章

1、文章一

肿瘤的

首先看以下这个图


image.png

这个图中都是导管细胞,做完拟时序分析后,左边的为正常的导管细胞,恶化的导管都是右边的,所以人为的指定左边的正常导管细胞为起始点,判断从左往后的肿瘤变化轨迹,然后得到了热图


image.png

左边正常,右边恶化,分为P5和P6两种模式,P5的基因表达量都升高,P6的基因先升高后面有点降低
image.png

这两个图是一个图

得到了逐渐升高的基因,就可以去做这些基因的功能分析了(可以挑一些基因去做,也可以差异基因去做)

2、文章二

类似于重编程


image.png

image.png

分为早期、中边、晚期
最右下角的那个图是单个基因的动态变化趋势,看逐渐升高还是降低,

3、文章三

细胞图谱的


image.png

image.png

A中颜色深的是起始点,浅的是终点
为什么下面不是终点,上面是终点?因为上面的线段长(时间长)


image.png

image.png

c中表示从start到end的基因的变化趋势

遇到这种三角的图(三岔口的拟时序分析)有两种处理方式


image.png

1、直接描绘start---end的变化趋势,像本文中的c

2、因为分支比较明显,可以研究这两者之间的差异


image.png

image.png

红黄这两种变化之间的差异

比如这篇重编程的文章


image.png

image.png

下面这条是sucessful 代表干细胞的重编程功能是成功的,得到成功的成纤维,上面那条是重编程失败的细胞的变化趋势,然后绘制热图


image.png

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