知识积累---空间数据结合图像分割实现单细胞级解卷积分析

作者,Evil Genius

今日参考文献

中国人发表的,但是,

这里多说一句,大家有条件多支持一下正版,盗版如果横行,说实话没有动力做下去了,这是大多数原创者的心声,当然了,我做学生的时候,也没钱,资料很多也是盗版,工作了以后也买卖过盗版, 现在不过不再跟盗版打交道了,知道了其中的辛酸。

当然了,量力而行,别跟我读书似的一个月就500的补贴,还非要买几千上万块钱的资料,家里条件很好,这可以,否则就叫败家。

花了钱的,好好学习,变成自己的,才是目的。

知识积累

  • 对空间转录组而言,在固定位置获取的相同大小的spot不能精确捕获随机位置的单细胞,因此本质上不能在单细胞水平上描述转录组(无论是10X Visium 、HD、华大、百迈客都有这样的问题)。
最佳的解决方案:将空间转录组与其匹配的基于组织学图像的形态对齐,恢复缺失的细胞,从而实现真正的单细胞水平和全片尺度的反卷积、卷积和聚类低分辨率和高分辨率spot。
  • 目前可用的空间分辨转录组学技术有两大类:(i)原位杂交或原位测序技术,如seqFISH、MERFISH、STARmap和FISSEQ;(ii)原位捕获技术,如空间转录组学(ST)、SLIDE-seq、ZipSeq和HDST。
  • 迄今为止,最主要的商用测序技术是10X Visium,这是一种基于spot的空间转录组学,在2018年从ST技术进一步发展而来。
基于spot的空间转录组学的主要技术限制之一在于spot的低分辨率,它通常覆盖多个细胞。目前的10X Visium spot尺寸为直径55 μm,spot覆盖细胞的数量根据生物组织的不同在1到30之间。为了提高分辨率,BayesSpace提出利用贝叶斯统计绘制空间转录组数据中的邻域结构,将分辨率提高到subspot水平。xfuse的方法使用深度生成模型将空间转录组与组织学图像数据相结合,可以在微米分辨率下表征解剖特征的转录组。还有一些方法,如stLearn和SpaGCN,将组织学图像整合到空间上平滑基因表达,并改善了连续空间模式的聚类。除了计算加强之外,空间转录组学技术的分辨率也在不断提高,例如,新版本的10X Visium HD,达到了2um 的精度。
然而,细胞大小的spot分辨率不等于单细胞水平。与单细胞测序技术不同的是,相同大小和固定位置的spot无法精确捕获实际大小不同、位置随机的细胞,仍然可以部分覆盖同时接近的多个细胞,这使得基于spot的空间基因表达无论spot大小如何,本质上都无法达到单细胞水平。这种单细胞水平的缺乏不能通过目前的方法来解决,仅通过计算或技术上提高spot分辨率,阻碍了空间转录组学中单个细胞空间组织和细胞类型特异性基因表达变异的表征。此外,在Visium平台上,在面积为6500 × 6500μm2的载玻片上分布了4992个直径为55 μm的spot,两个spot中心之间的距离为100 μm,在~70%的面积内未测量基因表达。
与新鲜冷冻组织相比,FFPE组织的染色图像更好地保留了细胞形态,因此可以更准确地分割细胞核。
将空间转录组数据与匹配的基于组织学图像的细胞核形态信息相结合,以弥合spot分辨率和单细胞水平之间的差距。

结果1、单细胞分辨率不等于单细胞水平

在10X Visium平台上,将组织切片放置在载玻片上,在载玻片上用条形码逆转录物(RT)引物在固定坐标上分析相同大小的spot。在组织渗透过程中,mRNA分子垂直向下扩散到固体表面,并在原位spot内与RT引物局部杂交。进一步提取cDNA mRNA复合物用于文库制备和测序读出。与单细胞RNA测序(scRNAseq)不同的是,单细胞技术可以精确地单独捕获每个细胞并描绘整个细胞的基因表达,而预定位置的spot极有可能同时部分覆盖近距离的多个细胞。然而,这些空间数据不能仅仅通过计算或技术上将spot大小缩小到单细胞大小来解决,即单细胞分辨率不是单细胞水平。

10X HD数据仍然有很大的机会部分覆盖多个细胞(这个在课程上已经讲过了)。一个spot,即使在非常小的尺寸上,仍然很可能同时覆盖多个细胞,每个细胞覆盖的很小一部分,使得spot仍然捕获细胞类型的混合基因表达。
因此,spot分辨率和单细胞水平之间的上述差距提出了新的计算挑战:分别针对低分辨率和高分辨率spot的“单细胞反卷积”和“单细胞卷积”。其中的区别在于单细胞反卷积是指将大的低分辨率spot分裂成单个细胞,而单细胞卷积是指将相邻的高分辨率spot覆盖的小块组装成一个完整的单个细胞。

结果2、STIE工作流程概述

  • spot水平的基因表达和基于匹配组织学图像的细胞核分割。


STIE的基本原理是,给定细胞类型转录组特征,可以从每个spot的基因表达估计细胞类型比例,同时也可以从组织学图像中通过分割细胞/细胞核,提取细胞/细胞核形态特征,将单个细胞聚类成groups,定位每个单个细胞并计算每个spot覆盖的细胞比例。STIE有两种模式:(1)单细胞反STIE有两种模式:(1)单细胞反卷积/卷积:根据细胞类型转录组特征,STIE将低/高分辨率spot反卷积/卷积成单细胞。转录组特征可以从现有的不同组织的单细胞图谱中定义。(2)单细胞聚类(Signature-free 单细胞反卷积/卷积):在没有细胞类型转录组特征的情况下,STIE可以在单细胞水平上进行聚类。聚类算法建立在反卷积/卷积的基础上:给定聚类数量和聚类转录组特征的初始值,STIE迭代改进基于单细胞反卷积/卷积的聚类特征,并重新估计特征,直到迭代收敛。
最后,根据spot基因表达和细胞形态,被spot覆盖的细胞被分配细胞类型,其中也包括那些位于捕获区域的spot大小之外但位于扩大的spot区域范围内的细胞。因此,STIE填补了spot间隙区域中缺失的细胞。

结果3、基因表达和核形态的整合使单细胞水平的反卷积/卷积在空间转录组学中得以实现(小鼠脑、乳腺)

结果4、STIE将空间转录组聚类提高到单细胞水平

结果5、STIE揭示了单细胞水平的空间转录组学见解

与单细胞基因表达谱不同,spot不能总是捕获完整的细胞,导致spot边界细胞的信息丢失或外部细胞的污染。因此,提出了第一个问题:被spot捕获的基因表达是来自spot内的所有细胞,还是来自spot外的更多细胞?
其次,与传统的细胞型反褶积方法不同,STIE将细胞核形态信息与基因表达结合起来。
基于细胞类型反卷积的细胞类型共定位。在给定基因表达的情况下,可以估计每个spot内的细胞类型比例。

结果6、STIE能够研究空间分辨的细胞-细胞相互作用

最后,感受一下效果

图像分割前

图像分割后


必须要结合图像分割

示例代码在GitHub - zhushijia/STIE: Spatial Transcriptome Image and Expression integration enables single cell level spatial data analysis

生活很好,有你更好

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

推荐阅读更多精彩内容