详细
markers=list(
epithelial=c(#'KRT14', 'KRT17', 'KRT6A', 'KRT5',
'EPCAM', 'KRT8', 'KRT18','KRT19', #'KRT16', 'KRT6B', 'KRT15', 'KRT6C', 'KRTCAP3',
#'SFN', 'PROM1', 'ALDH1A1',
'CD24'),
fibro=c("MME",'FAP', 'PDPN', "VCAN", #'COL1A2', 'DCN', 'COL3A1',
'COL6A1'),
endo =c('PECAM1', 'VWF', 'ENG'),
#myocytes=c('ACTA1', 'ACTN2', 'MYL2', 'MYH2'),
"T" =c("PTPRC", #'CD2', 'TRAC','CD7',
'CD3D', 'CD3E', #'CD3G',
'CD4', 'CD8A', 'CD8B', 'SELL',
'FOXP3', "CTLA4", "CD274", "PRF1", "IFNG", "RUNX3", "ZBTB7B"),
NK=c('NKG7', 'GNLY', 'GZMA', 'GZMB', "GZMK", "GZMH","GZMM",
"NCAM1", 'FGFBP2', 'CX3CR1'),
B = c( 'CD79A', 'CD79B', 'CD19', 'MS4A1', 'VPREB3', 'SLAMF7','BLNK', 'FCRL5'),
Plasma=c('IGHG1', 'MZB1', 'SDC1'),
mono=c('CD14', "FCGR3A"),
macro=c('CD163', 'CD68', 'FCGR2A', 'CSF1R', "SPP1"),
M1=c('NOS2', 'IRF5', 'PTGS2'),
M2=c( 'VSIG4', 'MS4A4A'), #'CD163',
DC=c('ITGAX','CD40', 'CD80', 'CD83', 'CCR7', "IL7R", "LAMP3", "CD1C"),
pDC=c("FSCN1","CCL19","IDO1", "IDO2", 'LILRA4', 'IRF7', 'IL3RA'), #MZB1
MHCII=c("HLA-DQA1", "HLA-DQB1", "HLA-DQB2", "CD74"),
#mast = c('CMA1', 'MS4A2', 'TPSAB1', 'TPSB2', 'CPA3'),
#eryth=c("HBA1", "HBA2", "HBB", "HBD"),
neu=c("MPO","ELANE", "LTF","CAMP", "CMTM2","PROK2","FCGR3B"),
#stem=c('CD34', 'KIT', 'SOCS2', 'SOCS1'),
#cycle=c('MKI67', 'CCND1', 'CCND3', 'CCNE1', 'CCNE2', 'CCNB1', 'CCNB2'),
#APA =c('NUDT21', "CPSF6", "CSTF2"),
meta=c("GAPDH", 'ACTB', "SCCPDH")
)
pdf( paste0(outputRoot, keyword, "_05_2.Manual.marker.DotPlot.pdf"), width=26, height=pdf_height)
DotPlot(scObj_2, features = markers, cluster.idents = T, cols=c('lightgrey', "red")) + RotatedAxis()+
ggtitle(keyword)
dev.off()
简短
cell_marker.Simple=list(
"T Cells"= c('PTPRC', 'CD3D', 'CD3E', "CD3G", 'CD8A', "CD8B", "CD4", "IL7R", "FOXP3"),
"B cells" = c('CD19', 'CD79A', "CD79B", 'MS4A1'), #[CD20]
"Plasma cells" = c('IGHG1', 'MZB1', 'SDC1'), #'CD79A'
"Monocytes and macrophages" = c('CD68', 'CD163', 'CD14'),
"NK Cells" = c('FGFBP2', 'FCGR3A', 'CX3CR1'),
"Photoreceptor cells" = c('RCVRN'),
"Fibroblasts" = c('FGF7', 'MME'),
"Endothelial cells" = c('PECAM1', 'VWF'),
"epi or tumor" = c('EPCAM', 'KRT19', 'PROM1', 'ALDH1A1', 'CD24')
#"immune" =c('PTPRC')
#"epithelial/cancer" = c( 'EPCAM'),
#"stromal" =c('MME', 'PECAM1')
)
DotPlot(tmp, features = cell_marker.Simple, cols = c("lightgrey", "red"), cluster.idents = T)+RotatedAxis()
# More gene list:
https://gitee.com/dawnEve/bioToolKit/blob/dev/R_scripts/single_script/do_Seurat.script.R