查看原文
其他

scHLAcount || 单细胞转录组HLA等位基因分析

周运来 单细胞天地 2022-06-07

分享是一种态度

作者 |  周运来

男,

一个长大了才会遇到的帅哥,

稳健,潇洒,大方,靠谱。

一段生信缘,一棵技能树,

一枚大型测序工厂的螺丝钉,

一个随机森林中提灯觅食的津门旅客。




scHLAcount允许我们使用个性化的参考基因组计算HLA I类基因HLA-A、B和C的单细胞转录组序列数据中的分子数;和HLA II类基因DPA1, DPB1, DRA1, DRB1, DQA1, DQB1。可以使用由替代方法确定的提供的HLA类型,也可以使用此工具分析HLA类型,然后根据这些调用进行量化。

scHLAcount可用于观察HLA基因等位基因的特异性表达。通过将细胞特异性HLA计数mapping在t-SNE图谱上,寻找单倍型完全缺失的簇,它还可以用来评估杂合性的缺失。一般的HLA表达损失也可以用scHLAcount来评估,并且在这个任务上比默认的Cell Ranger表现得更好,特别是当样本具有与参考不符的HLA单倍型时。

scHLAcount 是10X 官方出品,也是单细胞转录组数据挖掘的另一个例子,其实我们仅看HLA表达量的话,我们写过如何快速地从单细胞数据中观察HLA基因表达模式https://www.jianshu.com/p/cc9eca31dffd

但是用scHLAcount可以分析HLA等位基因的情况,当然是从bam文件中calling出来的。

这个也给我们释放了一个信号:单细胞转录组数据挖掘越来越多要从比对结果BAM文件中来挖了。

github 地址:https://github.com/10XGenomics/scHLAcount


往期回顾

提高单细胞分析准确度的工具之一:Self-assembling manifolds

单细胞转录组高级分析八:整合V(D)J数据

空白云服务器的一些生物信息学设置






如果你对单细胞转录组研究感兴趣,但又不知道如何入门,也许你可以关注一下下面的课程




看完记得顺手点个“在看”哦!


生物 | 单细胞 | 转录组丨资料每天都精彩

长按扫码可关注

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存