单细胞 | 转录因子足迹分析

数据加载

在本案例中,将采用之前在轨迹构建案例中已经介绍并处理过的数据集。

复制代码
library(Signac)
library(Seurat)

bone <- readRDS("cd34.rds")
DimPlot(bone, label = TRUE)

要执行足迹分析,必须首先向对象中添加Motif 信息,这包括每个Motif 的精确位置。这一过程可以通过使用"motif"和"packages"这两个包中的函数来实现。

复制代码
library(motifmatchr)
library(JASPAR2020)
library(TFBSTools)
library(BSgenome.Hsapiens.UCSC.hg19)

# extract position frequency matrices for the motifs
pwm <- getMatrixSet(
  x = JASPAR2020,
  opts = list(species = 9606, all_versions = FALSE)
)

# add motif information
bone <- AddMotifs(bone, genome = BSgenome.Hsapiens.UCSC.hg19, pfm = pwm)

Motif 足迹分析

现在可以对任何已知位置信息的Motif 进行足迹分析。通常,这会涵盖基因组中所有该Motif 的实例。也可以设置 in.peaks = TRUE 参数,以便只考虑那些位于分析中峰值区域内的Motif 。Footprint() 函数会收集所有必要的数据,并将其保存在分析结果中。之后,可以使用 PlotFootprint() 函数来绘制这些基序的足迹图。

复制代码
# gather the footprinting information for sets of motifs
bone <- Footprint(
  object = bone,
  motif.name = c("GATA2", "CEBPA", "EBF1"),
  genome = BSgenome.Hsapiens.UCSC.hg19
)

# plot the footprint data for each group of cells
p2 <- PlotFootprint(bone, features = c("GATA2", "CEBPA", "EBF1"))

p2 + patchwork::plot_layout(ncol = 1)

本文由mdnice多平台发布

相关推荐
ChinaRainbowSea14 分钟前
9-2 MySQL 分析查询语句:EXPLAIN(详细说明)
java·数据库·后端·sql·mysql
风象南20 分钟前
SpringBoot基于Java Agent的无侵入式监控实现
java·spring boot·后端
崎岖Qiu26 分钟前
【Spring篇08】:理解自动装配,从spring.factories到.imports剖析
java·spring boot·后端·spring·面试·java-ee
香饽饽~、1 小时前
【第十一篇】SpringBoot缓存技术
java·开发语言·spring boot·后端·缓存·intellij-idea
程序员爱钓鱼2 小时前
Go语言实战指南 —— Go中的反射机制:reflect 包使用
后端·google·go
ℳ₯㎕ddzོꦿ࿐3 小时前
Spring Boot 集成 MinIO 实现分布式文件存储与管理
spring boot·分布式·后端
ai小鬼头7 小时前
百度秒搭发布:无代码编程如何让普通人轻松打造AI应用?
前端·后端·github
考虑考虑7 小时前
@FilterRegistration和@ServletRegistration注解
spring boot·后端·spring
一只叫煤球的猫8 小时前
🔥 同事混用@Transactional和TransactionTemplate被我怼了,三种事务管理到底怎么选?
java·spring boot·后端
你的人类朋友9 天前
(●'◡'●)从Dockerfile快速入门Docker Compose
后端