空间转录组数据可视化介绍--loupe browser
,
就像10X genomics单细胞转录组测序一样,做完10X Visium 空间转录组测序拿到分析结果后,我们想对感兴趣的基因进行空间位置表达可视化、或者想特殊标注感兴趣的某个cluster,又或者我们想根据已知的某种细胞类型特异性高表达marker基因的表达标注这些细胞类型大致部位,抑或想自主分析2个不同cluster间差异表达基因绘制热图,等等。
嗯,像小编这种生信代码小白只能通过loupe browser来实现。loupe browser除了可以对10X genomics单细胞转录组数据进行可视化外,对于10X Visium空间转录分析结果中的cloupe文件也可进行个性化可视化操作。
数据及软件下载
使用10X官网的10X Visium空间转录组数据进行演示,数据下载链接:
https://support.10xgenomics.com/spatial-gene-expression/datasets。以Human Breast Cancer (Block A Section 1)数据为例。
上下滑动查看图片
loupe browser界面介绍
首先打开loupe browser软件后,选择上面下载的cloupe文件打开,之后的界面如下图所示。下面多图框中的第三张图对loupe browser界面做了详细介绍。有一点补充说明:如果要输出图片,对于组织切片显示不同cluster时,需要将切片用鼠标拖到软件正中,这样输出的图片中主图在中间比较美观,位置大家可以自主拖动后输出尝试。
上下滑动查看图片
右侧的界面介绍,我们下面作详细介绍。官方操作指南见Tutorial - Spatial Gene Expression,https://support.10xgenomics.com/spatial-gene-expression/software/visualization/latest/analysis。
我们可以选择以空间转录组最常见的这种Spatial方式(即组织切片)展示,还是将这些spot降维聚类的t-SNE和UMAP散点图形式(10X genomics单细胞转录组常见的可视化图)
上下滑动查看图片
如何标记感兴趣的基因或细胞的位置
本次使用的是乳腺癌切片空间转录组数据,那么乳腺癌中乳腺上皮细胞和乳腺腺体、以及NK T细胞的分布,我们可以尝试标记一下。
首先,我们从cellmarker数据库(http://biocc.hrbmu.edu.cn/CellMarker/)中查询一下相关细胞的marker基因。Mammary gland和Mammary epithelium的marker基因查找,大家按照下方截图所示类比查找即可。
上下滑动查看图片
接着,我们进行loupe browser可视化这些marker基因的表达。比Mammary gland Cancer stem cell,marker为CD49f、CD61、ESA,symbol分别为ITGA6、ITGB3、FLOT2。我们还可以针对选择某个cluster再根据某种规则细分类,通过Filters功能实现,大致方法类型,选中某个cluster后,根据某个或某些基因表达量高低划分。不再截图详述,大家可自行尝试。
上下滑动查看图片
基因表达热图及差异基因分析
上下滑动查看图片
上下滑动查看图片
推荐阅读
空间转录组ST & 单细胞转录组测序联合揭示胰腺导管腺癌的组织结构