scATAC-seq 细胞(based scRNA-seq)注释

1. 首先根据scRNA-seq找到细胞类群的marker 基因;比如下面三大类群细胞(虚线圈起来的);

比如:Myeloid cell: 20 个marker 基因。

2. 根据scATAC-seq 计算细胞之间的相似度;按照peak 开放性;

tSNE图,先根据开放性peak,对所有细胞进行聚类分析;然后随机挑取一个细胞(Myeloid,Maker 基因promoter开放),根据开放性,计算其他所有细胞和这个细胞的相似性,并在图中用 "point" 大小表明。紧接着保留相似性最高的n个细胞,用下图的线衡量他们之间关联;假设cellA 和 cellB 、cellC 关联.

3. 对于scATAC-seq里面的cellA, 假设计算与Myeloid cell 的influence score;

influence score= M(cell A)/N; 当score 大于一定值得时候,ATAC-seq里面的cellA ,就可以被推断是不是Myeloid cell.

M 指的是:在cellA里面, 落在promoter of Myeloid cell marker genes的ATAC-seq peaks

N 指的是: 在cellA里面 total ATAC-seq peaks。

4. 由于cellB 和cellC里面ATAC-seq peak由于数据稀疏问题,没有办法检测在Myeloid cell marker genes的ATAC-seq信息。但是根据 scATAC-seq里面:cellA 和 cellB 、cellC 关联. 因此也可以推断cellB cellC 也是Myeloid cell;这个可以被认为是diffuse influence.

ref:

https://genomebiology.biomedcentral.com/articles/10.1186/s13059-021-02279-1

https://static-content.springer.com/esm/art%3A10.1038%2Fs41586-021-03674-1/MediaObjects/41586_2021_3674_MOESM1_ESM.pdf

  • 2
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值