为了研究Type H血管跟间充质干细胞之间的细胞通讯,我需要将Type H血管的细胞提取出来。Type H血管以CD31高表达、Emcn高表达为特点。具体代码如下。pncs是我的seurat数据名。
### 提取某一基因的表达值到metadata中
expr <- pncs@assays$RNA
#Pecam1
gene_expression_Pecam1 <- expr %>% .['Pecam1',] %>% as.data.frame() %>% t()
gene_expression_Pecam1 <- as.data.frame(gene_expression_Pecam1)
colnames(gene_expression_Pecam1) <- 'Pecam1'
gene_expression_Pecam1$cell <- rownames(gene_expression_Pecam1)
#每个细胞表达Pecam1的数值
gene_expression_Pecam1$Pecam1
#Emcn
gene_expression_Emcn <- expr %>% .['Emcn',] %>% as.data.frame() %>% t()
gene_expression_Emcn <- as.data.frame(gene_expression_Emcn)
colnames(gene_expression_Emcn) <- 'Emcn'
gene_expression_Emcn$cell <- rownames(gene_expression_Emcn)
gene_expression_cd31hi <- gene_expression_Pecam1[which(gene_expression_Pecam1$Pecam1 > 2),]
gene_expression_Emcnhi <- gene_expression_Emcn[which(gene_expression_Emcn$Emcn > 2),]
# Assuming gene_expression_Pecam1 and gene_expression_Emcn have 'cell' column
cd31hi_cells <- gene_expression_Pecam1$cell[which(gene_expression_Pecam1$Pecam1 > 2)]
emcnhi_cells <- gene_expression_Emcn$cell[which(gene_expression_Emcn$Emcn > 2)]
# Find the common cells between the two sets
common_cells <- intersect(cd31hi_cells, emcnhi_cells)
# Subset your Seurat object based on common cells
pncs_TypeH <- pncs[, common_cells]
细胞通讯部分待更新。