注释用的参考数据库
![](https://img.haomeiwen.com/i25071788/238185eec6664c91.png)
细胞类型鉴定
![](https://img.haomeiwen.com/i25071788/274bef4abb3cd8b4.png)
思路
代码在cell_marker_annotation.R
1、查看所有marker基因
![](https://img.haomeiwen.com/i25071788/503fd53d22c5b602.png)
0亚群很多是T细胞的marker,可以初步注释为T细胞
1亚群也像T或者NK细胞
0和1很相似,可以计算这个亚群的差异基因,代码在sub_analysis_scRNA
2、计算特定两组细胞之间的差异基因,决定要不要把两个亚群合并
![](https://img.haomeiwen.com/i25071788/8ffa94fded9410b3.png)
这里看的是cluster0和3的差异,下表按照差异大小排序,
![](https://img.haomeiwen.com/i25071788/68971e9e2dfb4a2e.png)
发现cluster0和3差异大的都是线粒体基因,说明本身基因差异不大,那就可以合并;接着看看cluster0,3和cluster1的差异
![](https://img.haomeiwen.com/i25071788/877a6fef7dc8051a.png)
![](https://img.haomeiwen.com/i25071788/47dc661d6a6f6f50.png)
这些差异基因是有意义的,那么1也许可以不合并进去
免疫细胞注释可以用singleR
![](https://img.haomeiwen.com/i25071788/651dc2e7cc83717b.png)
注释小结
![](https://img.haomeiwen.com/i25071788/a1a2e7c2f243ee0f.png)
cellmarker优于panglaoDB优于其他数据库,但只能解决一部分问题,还要靠文献查找
可以把每个亚群找到的marker用在线工具做功能富集,初步看一下功能
![](https://img.haomeiwen.com/i25071788/ad6ad2c1df4e47ac.png)
![](https://img.haomeiwen.com/i25071788/8f2f32e3c547f651.png)
![](https://img.haomeiwen.com/i25071788/7b781efdddc2f75d.png)
![](https://img.haomeiwen.com/i25071788/51bddec42a246a3f.png)
网友评论