CellChat：细胞间相互作用分析利器
前言
我们知道，细胞间信息传递方式一个是细胞表面配受体的相互作用，另一个通过细胞产生的可溶性小分子，即细胞因子。在单细胞数据分析中下游，有时候我们想看某几种细胞类型之间的相互作用，就有人推荐我们做一个配受体分析。那什么是配受体？我们在文章Cell-Cell Interaction Database|| 单细胞配受体库你还在文章的附录里找吗？中提到配受体其实是细胞的特定蛋白，蛋白追溯到基因表达上就是基因对。
Inference and analysis of cell-cell communication using CellChat
Suoqin Jin, Christian F. Guerrero-Juarez, Lihua Zhang, Ivan Chang, Peggy Myung, Maksim V. Plikus, Qing Nie
bioRxiv 2020.07.21.214387; doi: https://doi.org/10.1101/2020.07.21.214387
今天，我们就用CellChat来分析一下我们的PBMC数据，看看配受体分析的一般流程。
除了从任何给定的scRNA-seq数据推断细胞间通信外，CellChat还提供了进一步的数据探索、分析和可视化功能。
它能够分析细胞与细胞间的通讯，以获得细胞发展轨迹上的连续状态。
该方法结合社会网络分析、模式识别和多种学习方法，可以定量地描述和比较推断出的细胞间通信网络。
它提供了一个易于使用的工具来提取和可视化推断网络信息。例如，它可以随时预测所有细胞群的主要信号输入和输出，以及这些细胞群和信号如何协调在一起实现功能。
它提供了几个可视化输出，以方便用户引导的直观数据解释。
1devtools::install_github("sqjin/CellChat")
CellChat需要两个输入:
一个是细胞的基因表达数据，
另一个是细胞标签(即细胞标签)。
对于基因表达数据矩阵，基因应该在带有行名的行中，cell应该在带有名称的列中。CellChat分析的输入是均一化的数据(Seurat@assay$RNA@data)。如果用户提供counts数据，可以用normalizeData函数来均一化。对于细胞的信息，需要一个带有rownames的数据格式作为CellChat的输入。
这两个文件在我们熟悉的Seurat对象中是很容易找到的，一个是均一化之后的数据，一个是细胞类型在metadata中。那么就让我们开始chat之旅吧。
数据配置
首先，我们加载包和引入实例数据。
1library(CellChat)
2library(ggplot2)
3library(ggalluvial)
4library(svglite)
5library(Seurat)
6library(SeuratData)
7options(stringsAsFactors = FALSE)
我们用Seurat给出的pbmc3k.final数据集，大部分的计算已经存在其对象中了：
1pbmc3k.final
2
3An object of class Seurat
413714 features across 2638 samples within 1 assay
5Active assay: RNA (13714 features, 2000 variable features)
6 2 dimensional reductions calculated: pca, umap
7
8 pbmc3k.final@commands$FindClusters # 你也看一看作者的其他命令，Seurat是记录其分析过程的。
9Command: FindClusters(pbmc3k.final, resolution = 0.5)
10Time: 2020-04-30 12:54:53
11graph.name : RNA_snn
12modularity.fxn : 1
13resolution : 0.5
14method : matrix
15algorithm : 1
16n.start : 10
17n.iter : 10
18random.seed : 0
19group.singletons : TRUE
20verbose : TRUE
按照我们刚才说的，我们在Seurat对象中提出CellChat需要的数据：
1data.input <- pbmc3k.final@assays$RNA@data
2identity = data.frame(group =pbmc3k.final$seurat_annotations , row.names = names(pbmc3k.final$seurat_annotations)) # create a dataframe consisting of the cell labels
3unique(identity$group) # check the cell labels
4
5[1] Memory CD4 T B CD14+ Mono NK CD8 T Naive CD4 T FCGR3A+ Mono DC Platelet
6Levels: Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
创建一个Cell Chat对象。
1cellchat <- createCellChat(data = data.input)
2cellchat
3
4An object of class CellChat
5 13714 genes.
6 2638 cells.
7
8summary(cellchat)
9
10Length Class Mode
11 1 CellChat S4
S4 类学会了吗？
在学习单细胞数据分析工具的时候，在知道了要干嘛之后，第二步就是看数据格式，俗称：单细胞数据格式。我们在听说你的单细胞对象需要一个思维导图？，曾给出一个简单的可视化数据结构的方法：导图。
1library(mindr)
2(out <- capture.output(str(cellchat)))
3out2 <- paste(out, collapse="\n")
4mm(gsub("\\.\\.@","# ",gsub("\\.\\. ","#",out2)),type ="text")
当然，我们可以用str来看，就是有点冗长：
1> str(cellchat)
2Formal class 'CellChat' [package "CellChat"] with 14 slots
3 ..@ data.raw : num[0 , 0 ]
4 ..@ data :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
5 .. .. ..@ i : int [1:2238732] 29 73 80 148 163 184 186 227 229 230 ...
6 .. .. ..@ p : int [1:2639] 0 779 2131 3260 4220 4741 5522 6304 7094 7626 ...
7 .. .. ..@ Dim : int [1:2] 13714 2638
8 .. .. ..@ Dimnames:List of 2
9 .. .. .. ..$ : chr [1:13714] "AL627309.1" "AP006222.2" "RP11-206L10.2" "RP11-206L10.9" ...
10 .. .. .. ..$ : chr [1:2638] "AAACATACAACCAC" "AAACATTGAGCTAC" "AAACATTGATCAGC" "AAACCGTGCTTCCG" ...
11 .. .. ..@ x : num [1:2238732] 1.64 1.64 2.23 1.64 1.64 ...
12 .. .. ..@ factors : list()
13 ..@ data.signaling: num[0 , 0 ]
14 ..@ data.scale : num[0 , 0 ]
15 ..@ data.project : num[0 , 0 ]
16 ..@ net : list()
17 ..@ netP : list()
18 ..@ meta :'data.frame': 0 obs. of 0 variables
19Formal class 'data.frame' [package "methods"] with 4 slots
20 .. .. ..@ .Data : list()
21 .. .. ..@ names : chr(0)
22 .. .. ..@ row.names: int(0)
23 .. .. ..@ .S3Class : chr "data.frame"
24 ..@ idents :Formal class 'factor' [package "methods"] with 3 slots
25 .. .. ..@ .Data : int(0)
26 .. .. ..@ levels : chr(0)
27 .. .. ..@ .S3Class: chr "factor"
28 ..@ DB : list()
29 ..@ LR : list()
30 ..@ var.features : logi(0)
31 ..@ dr : list()
32 ..@ options : list()
我们把metadata信息加到CellChat对象中，这个写法跟Seurat很像啊。
1cellchat <- addMeta(cellchat, meta = identity, meta.name = "labels")
2cellchat <- setIdent(cellchat, ident.use = "labels") # set "labels" as default cell identity
3levels(cellchat@idents) # show factor levels of the cell labels
4
5[1] "Naive CD4 T" "Memory CD4 T" "CD14+ Mono" "B" "CD8 T" "FCGR3A+ Mono" "NK"
6
7groupSize <- as.numeric(table(cellchat@idents)) # number of cells in each cell group
8
9[1] 697 483 480 344 271 162 155 32 14
导入配受体数据库
CellChat提供了人和小鼠的配受体数据库，分别可以用CellChatDB.human,CellChatDB.mouse来导入。来看一下这个数据库的结构吧。
1CellChatDB <- CellChatDB.human
2(out3 <- capture.output(str(CellChatDB)))
3out4 <- paste(out3, collapse="\n")
4mm(gsub("\\$","# ",gsub("\\.\\. ","#",out4)),type ="text")
这个数据库的信息是很全面的：
1> colnames(CellChatDB$interaction)
2 [1] "interaction_name" "pathway_name" "ligand" "receptor" "agonist" "antagonist" "co_A_receptor"
3 [8] "co_I_receptor" "evidence" "annotation" "interaction_name_2"
4> CellChatDB$interaction[1:4,1:4]
5 interaction_name pathway_name ligand receptor
6TGFB1_TGFBR1_TGFBR2 TGFB1_TGFBR1_TGFBR2 TGFb TGFB1 TGFbR1_R2
7TGFB2_TGFBR1_TGFBR2 TGFB2_TGFBR1_TGFBR2 TGFb TGFB2 TGFbR1_R2
8TGFB3_TGFBR1_TGFBR2 TGFB3_TGFBR1_TGFBR2 TGFb TGFB3 TGFbR1_R2
9TGFB1_ACVR1B_TGFBR2 TGFB1_ACVR1B_TGFBR2 TGFb TGFB1 ACVR1B_TGFbR2
10> head(CellChatDB$cofactor)
11 cofactor1 cofactor2 cofactor3 cofactor4 cofactor5 cofactor6 cofactor7 cofactor8 cofactor9 cofactor10 cofactor11 cofactor12
12ACTIVIN antagonist FST
13ACTIVIN inhibition receptor BAMBI
14ANGPT inhibition receptor 1 TIE1
15ANGPT inhibition receptor 2 PTPRB
16BMP antagonist NBL1 GREM1 GREM2 CHRD NOG BMP3 LEFTY1 LEFTY2
17BMP inhibition receptor BAMBI
18 cofactor13 cofactor14 cofactor15 cofactor16
19ACTIVIN antagonist
20ACTIVIN inhibition receptor
21ANGPT inhibition receptor 1
22ANGPT inhibition receptor 2
23BMP antagonist
24BMP inhibition receptor
25> head(CellChatDB$complex)
26 subunit_1 subunit_2 subunit_3 subunit_4
27Activin AB INHBA INHBB
28Inhibin A INHA INHBA
29Inhibin B INHA INHBB
30IL12AB IL12A IL12B
31IL23 complex IL12B IL23A
32IL27 complex IL27 EBI3
33> head(CellChatDB$geneInfo)
34 Symbol Name EntrezGene.ID Ensembl.Gene.ID MGI.ID Gene.group.name
35HGNC:5 A1BG alpha-1-B glycoprotein 1 ENSG00000121410 MGI:2152878 Immunoglobulin like domain containing
36HGNC:37133 A1BG-AS1 A1BG antisense RNA 1 503538 ENSG00000268895 Antisense RNAs
37HGNC:24086 A1CF APOBEC1 complementation factor 29974 ENSG00000148584 MGI:1917115 RNA binding motif containing
38HGNC:7 A2M alpha-2-macroglobulin 2 ENSG00000175899 MGI:2449119 C3 and PZP like, alpha-2-macroglobulin domain containing
39HGNC:27057 A2M-AS1 A2M antisense RNA 1 144571 ENSG00000245105 Antisense RNAs
40HGNC:23336 A2ML1 alpha-2-macroglobulin like 1 144568 ENSG00000166535 C3 and PZP like, alpha-2-macroglobulin domain containing
其实是记录了许多许多受配体相关的通路信息，不像有的配受体库只有一个基因对。这样，我们就可以更加扎实地把脚落到pathway上面了。在CellChat中，我们还可以先择特定的信息描述细胞间的相互作者，这个可以理解为从特定的侧面来刻画细胞间相互作用,比用一个大的配体库又精细了许多呢。
1CellChatDB.use <- subsetDB(CellChatDB, search = "Secreted Signaling") # use Secreted Signaling for cell-cell communication analysis
2cellchat@DB <- CellChatDB.use # set the used database in the object
有哪些可以选择的侧面呢？
1> unique(CellChatDB$interaction$annotation)
2[1] "Secreted Signaling" "ECM-Receptor" "Cell-Cell Contact"
预处理
对表达数据进行预处理，用于细胞间的通信分析。首先在一个细胞组中识别过表达的配体或受体，然后将基因表达数据投射到蛋白-蛋白相互作用(PPI)网络上。如果配体或受体过表达，则识别过表达配体和受体之间的相互作用。
1cellchat <- subsetData(cellchat) # subset the expression data of signaling genes for saving computation cost
2future::plan("multiprocess", workers = 4) # do parallel 这里似乎有一些bug，在Linux上居然不行。de了它。
3cellchat <- identifyOverExpressedGenes(cellchat)
4 cellchat <- identifyOverExpressedInteractions(cellchat)
5cellchat <- projectData(cellchat, PPI.human)
相互作用推断
然后，我们通过为每个相互作用分配一个概率值并进行置换检验来推断生物意义上的细胞-细胞通信。
1# cellchat <- computeCommunProb(cellchat) 注意这个函数如果你可以用就用，这个是作者的。
2mycomputeCommunProb <-edit(computeCommunProb) # computeCommunProb内部似乎有一些bug，同一套数据在window10上没事，到了Linux上有报错。发现是computeExpr_antagonist这个函数有问题，(matrix(1, nrow = 1, ncol = length((group))))，中应为(matrix(1, nrow = 1, ncol = length(unique(group))))？ 不然矩阵返回的不对。de了它。
3environment(mycomputeCommunProb) <- environment(computeCommunProb)
4cellchat <- mycomputeCommunProb(cellchat) # 这儿是我de过的。
推测细胞间在信号通路水平上的通讯。我们还通过计算与每个信号通路相关的所有配体-受体相互作用的通信概率来推断信号通路水平上的通信概率。
注:推测的每个配体-受体对的细胞间通信网络和每个信号通路分别存储在“net”和“netP”槽中。
我们可以通过计算链路的数量或汇总通信概率来计算细胞间的聚合通信网络。
1cellchat <- computeCommunProbPathway(cellchat)
2cellchat <- aggregateNet(cellchat)
让我们看看这结果。
1> cellchat@netP$pathways
2 [1] "TGFb" "NRG" "PDGF" "CCL" "CXCL" "MIF" "IL2" "IL6" "IL10" "IL1" "CSF"
3[12] "IL16" "IFN-II" "LT" "LIGHT" "FASLG" "TRAIL" "BAFF" "CD40" "VISFATIN" "COMPLEMENT" "PARs"
4[23] "FLT3" "ANNEXIN" "GAS" "GRN" "GALECTIN" "BTLA" "BAG"
5> head(cellchat@LR$LRsig)
6 interaction_name pathway_name ligand receptor agonist antagonist co_A_receptor co_I_receptor
7TGFB1_TGFBR1_TGFBR2 TGFB1_TGFBR1_TGFBR2 TGFb TGFB1 TGFbR1_R2 TGFb agonist TGFb antagonist TGFb inhibition receptor
8TGFB1_ACVR1B_TGFBR2 TGFB1_ACVR1B_TGFBR2 TGFb TGFB1 ACVR1B_TGFbR2 TGFb agonist TGFb antagonist TGFb inhibition receptor
9TGFB1_ACVR1C_TGFBR2 TGFB1_ACVR1C_TGFBR2 TGFb TGFB1 ACVR1C_TGFbR2 TGFb agonist TGFb antagonist TGFb inhibition receptor
10TGFB1_ACVR1_TGFBR1 TGFB1_ACVR1_TGFBR1 TGFb TGFB1 ACVR1_TGFbR
11WNT10A_FZD1_LRP5 WNT10A_FZD1_LRP5 WNT WNT10A FZD1_LRP5 WNT agonist WNT antagonist WNT activation receptor WNT inhibition receptor
12WNT10A_FZD2_LRP5 WNT10A_FZD2_LRP5 WNT WNT10A FZD2_LRP5 WNT agonist WNT antagonist WNT activation receptor WNT inhibition receptor
13 evidence annotation interaction_name_2
14TGFB1_TGFBR1_TGFBR2 KEGG: hsa04350 Secreted Signaling TGFB1 - (TGFBR1+TGFBR2)
15TGFB1_ACVR1B_TGFBR2 PMID: 27449815 Secreted Signaling TGFB1 - (ACVR1B+TGFBR2)
16TGFB1_ACVR1C_TGFBR2 PMID: 27449815 Secreted Signaling TGFB1 - (ACVR1C+TGFBR2)
17TGFB1_ACVR1_TGFBR1 PMID: 29376829 Secreted Signaling TGFB1 - (ACVR1+TGFBR1)
18WNT10A_FZD1_LRP5 KEGG: hsa04310; PMID: 23209157 Secreted Signaling WNT10A - (FZD1+LRP5)
19WNT10A_FZD2_LRP5 KEGG: hsa04310; PMID: 23209159 Secreted Signaling WNT10A - (FZD2+LRP5)
可视化
在推断细胞-细胞通信网络的基础上，CellChat为进一步的探索、分析和可视化提供了各种功能。
通过结合社会网络分析、模式识别和多种学习方法的综合方法，t可以定量地描述和比较推断出的细胞-细胞通信网络。
它提供了一个易于使用的工具来提取和可视化推断网络的高阶信息。例如，它可以随时预测所有细胞群的主要信号输入和输出，以及这些细胞群和信号如何协调在一起实现功能。
你可以使用层次图或圈图可视化每个信号通路。如果使用层次图可视化通信网络，请定义vertex.receiver，它是一个数字向量，给出作为第一个层次结构图中的目标的细胞组的索引。我们可以使用netVisual_aggregate来可视化信号路径的推断通信网络，并使用netVisual_individual来可视化与该信号路径相关的单个L-R对的通信网络。
在层次图中，实体圆和空心圆分别表示源和目标。圆的大小与每个细胞组的细胞数成比例。边缘颜色与信源一致。线越粗，信号越强。这里我们展示了一个MIF信号网络的例子。所有显示重要通信的信令路径都可以通过cellchat@netP$pathways访问。
1>cellchat@netP$pathways
2 [1] "TGFb" "NRG" "PDGF" "CCL" "CXCL" "MIF" "IL2" "IL6" "IL10" "IL1"
3[11] "CSF" "IL16" "IFN-II" "LT" "LIGHT" "FASLG" "TRAIL" "BAFF" "CD40" "VISFATIN"
4[21] "COMPLEMENT" "PARs" "FLT3" "ANNEXIN" "GAS" "GRN" "GALECTIN" "BTLA" "BAG" 1levels(cellchat@idents)
2vertex.receiver = seq(1,4) # a numeric vector
3# check the order of cell identity to set suitable vertex.receiver
4#cellchat@LR$LRsig$pathway_name
5#cellchat@LR$LRsig$antagonist
6pathways.show <- "MIF"
7# netVisual_aggregate(cellchat, signaling = pathways.show, vertex.receiver = vertex.receiver, vertex.size = groupSize) # 原函数
8mynetVisual_aggregate(cellchat, signaling = pathways.show, vertex.receiver = vertex.receiver, vertex.size = groupSize) 原函数这里似乎有一个和igraph相关的小问题在不同igraph可能会表现bug，不巧我遇到了，de了它。
经典的配受体圈图：
1mynetVisual_aggregate(cellchat, signaling = c("MIF"), layout = "circle", vertex.size = groupSize,pt.title=20,vertex.label.cex = 1.7)
计算和可视化每个配体-受体对整个信号通路的贡献度。
1netAnalysis_contribution(cellchat, signaling = pathways.show)
识别细胞群的信号转导作用，通过计算每个细胞群的网络中心性指标，CellChat允许随时识别细胞间通信网络中的主要发送者、接收者、调解者和影响者。
1cellchat <- netAnalysis_signalingRole(cellchat, slot.name = "netP") # the slot 'netP' means the inferred intercellular communication network of signaling pathways
1netVisual_signalingRole(cellchat, signaling = pathways.show, width = 12, height = 2.5, font.size = 10)
识别特定细胞群的全局通信模式和主要信号。除了探索单个通路的详细通讯外，一个重要的问题是多个细胞群和信号通路如何协调运作。CellChat采用模式识别方法来识别全局通信模式以及每个小群的关键信号。
识别分泌细胞外向交流模式。随着模式数量的增加，可能会出现冗余的模式，使得解释通信模式变得困难。我们选择了5种模式作为默认模式。一般来说，当模式的数量大于2时就可以认为具有生物学意义。
1nPatterns = 5
2# 同样在这里遇到了bug，难道说是我没有安装好吗，de了它。
3# cellchat <- myidentifyCommunicationPatterns(cellchat, pattern = "outgoing", k = nPatterns)
4myidentifyCommunicationPatterns <- edit(identifyCommunicationPatterns)
5environment(myidentifyCommunicationPatterns) <- environment(identifyCommunicationPatterns)
6cellchat <- myidentifyCommunicationPatterns(cellchat, pattern = "outgoing", k = nPatterns)
1# Visualize the communication pattern using river plot
2netAnalysis_river(cellchat, pattern = "outgoing")
1# Visualize the communication pattern using dot plot
2netAnalysis_dot(cellchat, pattern = "outgoing")
识别目标细胞的传入(incoming)通信模式。
1netAnalysis_river(cellchat, pattern = "incoming")
1netAnalysis_dot(cellchat, pattern = "incoming")
作为结尾有大量的空间，我们得以先看看cellchat配受体推断的结构是如何的。
1> head(cellchat@LR$LRsig)
2 interaction_name pathway_name ligand receptor agonist antagonist co_A_receptor co_I_receptor
3TGFB1_TGFBR1_TGFBR2 TGFB1_TGFBR1_TGFBR2 TGFb TGFB1 TGFbR1_R2 TGFb agonist TGFb antagonist TGFb inhibition receptor
4TGFB1_ACVR1B_TGFBR2 TGFB1_ACVR1B_TGFBR2 TGFb TGFB1 ACVR1B_TGFbR2 TGFb agonist TGFb antagonist TGFb inhibition receptor
5TGFB1_ACVR1C_TGFBR2 TGFB1_ACVR1C_TGFBR2 TGFb TGFB1 ACVR1C_TGFbR2 TGFb agonist TGFb antagonist TGFb inhibition receptor
6TGFB1_ACVR1_TGFBR1 TGFB1_ACVR1_TGFBR1 TGFb TGFB1 ACVR1_TGFbR
7WNT10A_FZD1_LRP5 WNT10A_FZD1_LRP5 WNT WNT10A FZD1_LRP5 WNT agonist WNT antagonist WNT activation receptor WNT inhibition receptor
8WNT10A_FZD2_LRP5 WNT10A_FZD2_LRP5 WNT WNT10A FZD2_LRP5 WNT agonist WNT antagonist WNT activation receptor WNT inhibition receptor
9 evidence annotation interaction_name_2
10TGFB1_TGFBR1_TGFBR2 KEGG: hsa04350 Secreted Signaling TGFB1 - (TGFBR1+TGFBR2)
11TGFB1_ACVR1B_TGFBR2 PMID: 27449815 Secreted Signaling TGFB1 - (ACVR1B+TGFBR2)
12TGFB1_ACVR1C_TGFBR2 PMID: 27449815 Secreted Signaling TGFB1 - (ACVR1C+TGFBR2)
13TGFB1_ACVR1_TGFBR1 PMID: 29376829 Secreted Signaling TGFB1 - (ACVR1+TGFBR1)
14WNT10A_FZD1_LRP5 KEGG: hsa04310; PMID: 23209157 Secreted Signaling WNT10A - (FZD1+LRP5)
15WNT10A_FZD2_LRP5 KEGG: hsa04310; PMID: 23209159 Secreted Signaling WNT10A - (FZD2+LRP5)
16> head(cellchat@dr)
17list()
18> head(cellchat@data)
196 x 2638 sparse Matrix of class "dgCMatrix"
20 [[ suppressing 70 column names 'AAACATACAACCAC', 'AAACATTGAGCTAC', 'AAACATTGATCAGC' ... ]]
21
22AL627309.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
23AP006222.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
24RP11-206L10.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
25RP11-206L10.9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
26LINC00115 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
27NOC2L . . . . . . . . . . . 1.646272 . . . . . . . . 1.398186 . . . . . . . . . . . . 1.89939 . . . . . . . 1.36907 1.721224 . . . . . . . . .
28
29AL627309.1 . . . . . . . . . . . . . . . . . . ......
30AP006222.2 . . . . . . . . . . . . . . . . . . ......
31RP11-206L10.2 . . . . . . . . . . . . . . . . . . ......
32RP11-206L10.9 . . . . . . . . . . . . . . . . . . ......
33LINC00115 . . . . . . . . . . . . . . . . . . ......
34NOC2L . . . 1.568489 1.678814 . 1.253835 . . 3.791113 . . . . . . . . ......
35
36 .....suppressing 2568 columns in show(); maybe adjust 'options(max.print= *, width = *)'
37 ..............................
38> head(cellchat@idents)
39[1] Memory CD4 T B Memory CD4 T CD14+ Mono NK Memory CD4 T
40Levels: Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
41> head(cellchat@meta)
42 labels
43AAACATACAACCAC Memory CD4 T
44AAACATTGAGCTAC B
45AAACATTGATCAGC Memory CD4 T
46AAACCGTGCTTCCG CD14+ Mono
47AAACCGTGTATGCG NK
48AAACGCACTGGTAC Memory CD4 T
49> head(cellchat@netP$pathways)
50[1] "TGFb" "NRG" "PDGF" "CCL" "CXCL" "MIF"
51> head(cellchat@netP$prob)
52[1] 0 0 0 0 0 0
53> head(cellchat@netP$centr)
54$TGFb
55$TGFb$outdeg
56 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
570.000000e+00 5.798502e-07 2.634094e-05 0.000000e+00 1.108822e-06 9.977646e-06 9.953461e-06 2.840617e-07 3.475282e-06
58
59$TGFb$indeg
60 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
610.000000e+00 1.002762e-05 1.384499e-05 0.000000e+00 7.596075e-06 1.270618e-05 5.256794e-06 5.744824e-07 1.713913e-06
62
63$TGFb$hub
64 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
65 0.00000000 0.02278982 1.00000000 0.00000000 0.04484954 0.37878876 0.37787064 0.01116456 0.13193619
66
67$TGFb$authority
68 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
69 0.00000000 0.74712407 1.00000000 0.00000000 0.56314554 0.86435263 0.37969073 0.04280264 0.11659336
70
71$TGFb$eigen
72 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
73 0.01217244 0.31304003 1.00000000 0.01217244 0.25802457 0.58202001 0.37843282 0.02320534 0.12622971
74
75$TGFb$page_rank
76 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
77 0.02054795 0.13742492 0.21555291 0.02054795 0.11208641 0.31212523 0.09458943 0.02724384 0.05988138
78
79$TGFb$betweenness
80 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
81 0 0 24 0 0 10 0 0 0
82
83$TGFb$flowbet
84[1] 0.000000e+00 4.342669e-06 2.862661e-05 0.000000e+00 6.752863e-06 2.460332e-05 1.254051e-05 1.032200e-06 6.967716e-06
85
86$TGFb$info
87[1] 0.00000000 0.16628670 0.19401551 0.00000000 0.12870372 0.18191312 0.16895822 0.03556505 0.12455769
88
89$NRG
90$NRG$outdeg
91 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
921.116774e-10 1.024289e-10 2.194763e-10 5.436629e-11 5.792191e-11 1.166520e-10 4.634672e-11 1.511780e-11 1.629172e-12
93
94$NRG$indeg
95 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
960.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 7.256165e-10
97
98$NRG$hub
99 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
100 0.508835533 0.466696996 1.000000000 0.247709130 0.263909627 0.531501345 0.211169583 0.068881216 0.007422998
101
102$NRG$authority
103 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
1044.163336e-17 4.163336e-17 4.163336e-17 4.163336e-17 4.163336e-17 4.163336e-17 4.163336e-17 4.163336e-17 1.000000e+00
105
106$NRG$eigen
107 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
108 0.36342198 0.33332567 0.71422288 0.17691953 0.18849029 0.37961042 0.15082215 0.04919654 1.00000000
109
110$NRG$page_rank
111 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
112 0.01666667 0.01666667 0.01666667 0.01666667 0.01666667 0.01666667 0.01666667 0.01666667 0.86666667
113
114$NRG$betweenness
115 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
116 0 0 0 0 0 0 0 0 0
117
118$NRG$flowbet
119[1] 0 0 0 0 0 0 0 0 0
120
121$NRG$info
122[1] 0 0 0 0 0 0 0 0 0
123
124$PDGF
125$PDGF$outdeg
126 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
1272.117157e-10 5.254122e-10 1.830680e-09 0.000000e+00 3.046756e-10 1.195279e-09 6.457814e-10 1.492427e-10 0.000000e+00
128
129$PDGF$indeg
130 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
1319.596760e-10 7.355168e-10 1.375790e-09 0.000000e+00 4.145239e-10 1.028332e-09 2.501300e-10 9.881712e-11 0.000000e+00
132
133$PDGF$hub
134 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
135 0.09759699 0.32222056 1.00000000 0.00000000 0.18684898 0.65291566 0.35275497 0.08152314 0.00000000
136
137$PDGF$authority
138 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
1399.058608e-01 6.942716e-01 1.000000e+00 2.363558e-17 3.912788e-01 6.197010e-01 2.361036e-01 7.182571e-02 2.363558e-17
140
141$PDGF$eigen
142 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
143 0.44237332 0.51181753 1.00000000 0.07396075 0.29250188 0.67517921 0.29135234 0.07823533 0.07396075
144
145$PDGF$page_rank
146 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
147 0.15091590 0.12046482 0.24044555 0.02054795 0.07685927 0.27934926 0.05452706 0.03634225 0.02054795
148
149$PDGF$betweenness
150 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
151 1 0 18 0 0 5 0 0 0
152
153$PDGF$flowbet
154[1] 8.857166e-10 1.204604e-09 4.049689e-09 0.000000e+00 8.517939e-10 3.745196e-09 1.048193e-09 4.458839e-10 0.000000e+00
155
156$PDGF$info
157[1] 0.16144948 0.14611532 0.20300365 0.00000000 0.10956327 0.17885050 0.14080069 0.06021709 0.00000000
158
159$CCL
160$CCL$outdeg
161 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
1621.682814e-04 6.442088e-04 9.328993e-04 9.764691e-05 4.601953e-03 1.067399e-05 2.613615e-03 5.048297e-05 2.374245e-04
163
164$CCL$indeg
165 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
1661.013085e-03 1.208426e-03 4.952297e-04 5.869028e-04 3.900117e-03 1.125963e-04 1.773075e-03 7.483047e-05 1.929230e-04
167
168$CCL$hub
169 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
170 0.050735945 0.193934101 0.210819077 0.029282445 1.000000000 0.003249727 0.551908511 0.013914236 0.052892139
171
172$CCL$authority
173 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
174 0.30095289 0.35990610 0.14750945 0.17431275 1.00000000 0.03323390 0.45558530 0.02222082 0.04989215
175
176$CCL$eigen
177 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
178 0.17214374 0.27750548 0.17861545 0.09964869 1.00000000 0.01772801 0.50285164 0.01802822 0.05152917
179
180$CCL$page_rank
181 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
182 0.08937815 0.10366984 0.05186354 0.05878616 0.41583926 0.02465234 0.19773523 0.02202754 0.03604793
183
184$CCL$betweenness
185 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
186 0 0 0 0 56 0 0 0 0
187
188$CCL$flowbet
189[1] 6.253950e-04 1.206020e-03 1.184412e-03 4.216339e-04 7.464863e-03 7.286026e-05 3.851205e-03 1.024123e-04 5.393918e-04
190
191$CCL$info
192[1] 0.13488584 0.13862093 0.12659975 0.11726949 0.15963716 0.03961851 0.15306688 0.04024833 0.09005310
193
194$CXCL
195$CXCL$outdeg
196 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
1970.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 2.948861e-08
198
199$CXCL$indeg
200 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
2016.251119e-09 5.660697e-09 4.984283e-09 2.735102e-09 2.997064e-09 3.851281e-09 2.461799e-09 4.823805e-10 6.488065e-11
202
203$CXCL$hub
204 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
205 0 0 0 0 0 0 0 0 1
206
207$CXCL$authority
208 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
209 1.00000000 0.90554935 0.79734257 0.43753795 0.47944431 0.61609465 0.39381731 0.07716707 0.01037905
210
211$CXCL$eigen
212 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
213 0.5394037 0.4884566 0.4300895 0.2360096 0.2586140 0.3323237 0.2124265 0.0416242 1.0000000
214
215$CXCL$page_rank
216 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
217 0.1198308 0.1181000 0.1161172 0.1095240 0.1102919 0.1127960 0.1087229 0.1029205 0.1016966
218
219$CXCL$betweenness
220 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
221 0 0 0 0 0 0 0 0 0
222
223$CXCL$flowbet
224[1] 0 0 0 0 0 0 0 0 0
225
226$CXCL$info
227[1] 0.12994155 0.12702636 0.12305974 0.10129559 0.10488823 0.11427509 0.09707279 0.03583427 0.16660638
228
229$MIF
230$MIF$outdeg
231 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
2320.0012989751 0.0039272021 0.0006234461 0.0006401726 0.0005135156 0.0002049902 0.0003848437 0.0001321595 0.0000000000
233
234$MIF$indeg
235 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
2360.0005188736 0.0008184262 0.0007859180 0.0035144980 0.0009227472 0.0008137752 0.0001170739 0.0002339928 0.0000000000
237
238$MIF$hub
239 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
2404.252550e-01 1.000000e+00 2.238501e-01 2.095786e-01 1.680262e-01 7.360549e-02 1.160678e-01 4.315756e-02 2.774719e-18
241
242$MIF$authority
243 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
2442.769140e-01 4.020539e-01 2.249636e-01 1.000000e+00 3.209851e-01 2.590427e-01 6.228011e-02 7.151140e-02 4.690529e-18
245
246$MIF$eigen
247 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
248 0.40785736 1.00000000 0.28217435 0.81714092 0.31062247 0.21639268 0.11882053 0.07323643 0.01492405
249
250$MIF$page_rank
251 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
252 0.02128513 0.02564654 0.12732754 0.50874503 0.11392566 0.11715499 0.01911107 0.04839913 0.01840491
253
254$MIF$betweenness
255 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
256 0 10 0 17 14 0 11 0 0
257
258$MIF$flowbet
259[1] 0.0010772004 0.0004430504 0.0013253722 0.0018828736 0.0013830050 0.0007361476 0.0002572374 0.0005064761 0.0000000000
260
261$MIF$info
262[1] 0.10896205 0.16504074 0.11589344 0.17947163 0.13549734 0.10051455 0.12311142 0.07150883 0.00000000
每个pattern有outgoing和ingoing两种。
1> head(cellchat@netP$pattern$outgoing$pattern$cell)
2 CellGroup Pattern Contribution
31 Naive CD4 T Pattern 1 9.182571e-01
42 Memory CD4 T Pattern 1 8.643879e-01
53 CD14+ Mono Pattern 1 6.958107e-04
64 B Pattern 1 8.943340e-01
75 CD8 T Pattern 1 8.497941e-02
86 FCGR3A+ Mono Pattern 1 2.351798e-05
9> head(cellchat@netP$pattern$outgoing$pattern$signaling)
10 Pattern Signaling Contribution
111 Pattern 1 TGFb 1.509635e-08
122 Pattern 2 TGFb 5.851347e-01
133 Pattern 3 TGFb 2.021400e-01
144 Pattern 4 TGFb 4.466321e-08
155 Pattern 5 TGFb 2.127253e-01
166 Pattern 1 NRG 3.333424e-01
17> head(cellchat@netP$pattern$outgoing$data)
18 TGFb NRG PDGF CCL CXCL MIF IL2 IL6 IL10 IL1 CSF IL16 IFN-II
19Naive CD4 T 0.00000000 0.5088355 0.1156487 0.036567375 0 0.33076349 1.000000000 0.21361180 0.017388599 1.043256e-04 0.0006363636 0 0.004454402
20Memory CD4 T 0.02201327 0.4666970 0.2870039 0.139985939 0 1.00000000 0.948036204 0.22211580 1.000000000 1.150654e-04 0.0006048585 0 0.004707477
21CD14+ Mono 1.00000000 1.0000000 1.0000000 0.202718122 0 0.15875069 0.000000000 0.09461735 0.005818249 1.000000e+00 0.0010788329 0 0.005461241
22B 0.00000000 0.2477091 0.0000000 0.021218579 0 0.16300984 0.009150461 0.02181469 0.003863723 2.876928e-05 0.0002110580 0 0.001720322
23CD8 T 0.04209499 0.2639096 0.1664276 1.000000000 0 0.13075865 0.475620565 0.12534217 0.527133566 4.519162e-05 0.0003131413 0 0.003303116
24FCGR3A+ Mono 0.37878860 0.5315013 0.6529157 0.002319449 0 0.05219751 0.000000000 0.03752352 0.253673778 7.630358e-05 1.0000000000 0 0.004745991
25 LT LIGHT FASLG TRAIL BAFF CD40 VISFATIN COMPLEMENT PARs FLT3 ANNEXIN GAS
26Naive CD4 T 1.0000000 0.0000000 0.12801302 0.00000000 1.987539e-04 0.0052298348 0 1.0000000 0 1.0000000000 0.3515932720 0.02399186
27Memory CD4 T 0.8516886 1.0000000 0.85744830 0.09989685 2.286423e-04 1.0000000000 0 0.9403386 0 0.6925428133 1.0000000000 0.03584303
28CD14+ Mono 0.0512085 0.0000000 1.00000000 1.00000000 1.000000e+00 0.0080996253 0 0.8803694 0 0.0006179983 0.7171291990 0.02706222
29B 0.5629699 0.0000000 0.06312626 0.00000000 8.393504e-05 0.0003093270 0 0.3587101 0 0.0003490343 0.0003780528 0.01054186
30CD8 T 0.1842115 0.0000000 0.08407400 0.00000000 6.513411e-05 0.0008636328 0 0.5033253 1 0.0004055095 0.4595993742 0.01898338
31FCGR3A+ Mono 0.0832080 0.2745868 0.63644930 0.93360412 3.279022e-01 0.0044454725 1 0.3187685 0 0.0002367928 0.2119665274 0.01193921
32 GRN GALECTIN BTLA BAG
33Naive CD4 T 0.0000000 0.0000000 0.0000000 1.0000000
34Memory CD4 T 0.0000000 0.0000000 1.0000000 0.9388102
35CD14+ Mono 1.0000000 0.8983294 0.0000000 0.7920962
36B 0.0000000 0.0000000 0.5998942 0.4454517
37CD8 T 0.0000000 0.0000000 0.0000000 0.4831780
38FCGR3A+ Mono 0.1277283 1.0000000 0.2785847 0.3247730
39> cellchat@net
40$prob
41, , TGFB1_TGFBR1_TGFBR2
42
43 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
44Naive CD4 T 1.222691e-11 1.692462e-09 2.264589e-09 4.620186e-12 1.291360e-09 1.960243e-09 8.655394e-10 9.634429e-11 2.629338e-10
45Memory CD4 T 2.270338e-09 3.142597e-07 4.204920e-07 8.578932e-10 2.397814e-07 3.639734e-07 1.607142e-07 1.788942e-08 4.882008e-08
46CD14+ Mono 2.719456e-08 3.763876e-06 5.036034e-06 1.027602e-08 2.871745e-06 4.358185e-06 1.924748e-06 2.142640e-07 5.844517e-07
47B 3.582287e-12 4.958639e-10 6.634879e-10 1.353639e-12 3.783474e-10 5.743193e-10 2.535890e-10 2.822731e-11 7.703534e-11
48CD8 T 1.736672e-09 2.403890e-07 3.216497e-07 6.562368e-10 1.834175e-07 2.784145e-07 1.229360e-07 1.368429e-08 3.734375e-08
49FCGR3A+ Mono 1.030133e-08 1.425741e-06 1.907620e-06 3.892565e-09 1.087800e-06 1.650808e-06 7.290809e-07 8.116246e-08 2.213748e-07
50NK 1.027623e-08 1.422259e-06 1.902958e-06 3.883081e-09 1.085141e-06 1.646755e-06 7.272983e-07 8.096435e-08 2.208291e-07
51DC 1.112404e-09 1.539695e-07 2.060130e-07 4.203442e-10 1.174768e-07 1.783002e-07 7.873805e-08 8.764877e-09 2.391283e-08
52Platelet 3.590036e-09 4.966840e-07 6.644667e-07 1.356569e-09 3.789035e-07 5.745492e-07 2.539314e-07 2.827603e-08 7.699129e-08
53
54, , TGFB1_ACVR1B_TGFBR2
55
56 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
57Naive CD4 T 8.440868e-12 3.075855e-10 4.274550e-10 3.315884e-12 2.229500e-10 3.681800e-10 1.570298e-10 1.834846e-11 4.756356e-11
58Memory CD4 T 1.567331e-09 5.711352e-08 7.937122e-08 6.157055e-10 4.139808e-08 6.836461e-08 2.915780e-08 3.407004e-09 8.831697e-09
59CD14+ Mono 1.877381e-08 6.841044e-07 9.506996e-07 7.375044e-09 4.958626e-07 8.188301e-07 3.492476e-07 4.080909e-08 1.057772e-07
60B 2.473037e-12 9.011756e-11 1.252374e-10 9.715001e-13 6.532072e-11 1.078708e-10 4.600719e-11 5.375799e-12 1.393535e-11
61CD8 T 1.198914e-09 4.368838e-08 6.071417e-08 4.709778e-10 3.166702e-08 5.229470e-08 2.230394e-08 2.606152e-09 6.755695e-09
62FCGR3A+ Mono 7.111536e-09 2.591388e-07 3.601247e-07 2.793674e-09 1.878326e-07 3.101709e-07 1.322948e-07 1.545849e-08 4.006799e-08
63NK 7.094210e-09 2.585072e-07 3.592468e-07 2.786868e-09 1.873748e-07 3.094142e-07 1.319723e-07 1.542082e-08 3.997016e-08
64DC 7.679496e-10 2.798377e-08 3.888915e-08 3.016789e-10 2.028365e-08 3.349550e-08 1.428628e-08 1.669323e-09 4.327040e-09
65Platelet 2.478389e-09 9.030435e-08 1.254923e-07 9.736029e-10 6.545434e-08 1.080686e-07 4.610002e-08 5.386992e-09 1.395861e-08
66
67, , TGFB1_ACVR1C_TGFBR2
68
69 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
70Naive CD4 T 0 0 0 0 0 0 0 0 0
71Memory CD4 T 0 0 0 0 0 0 0 0 0
72CD14+ Mono 0 0 0 0 0 0 0 0 0
73B 0 0 0 0 0 0 0 0 0
74CD8 T 0 0 0 0 0 0 0 0 0
75FCGR3A+ Mono 0 0 0 0 0 0 0 0 0
76NK 0 0 0 0 0 0 0 0 0
77DC 0 0 0 0 0 0 0 0 0
78Platelet 0 0 0 0 0 0 0 0 0
79
80, , TGFB1_ACVR1_TGFBR1
81
82 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
83Naive CD4 T 1.189109e-11 3.873823e-10 4.965448e-10 4.544316e-12 2.846539e-10 4.232163e-10 1.881368e-10 2.208094e-11 6.074030e-11
84Memory CD4 T 2.207939e-09 7.192901e-08 9.219821e-08 8.437887e-10 5.285441e-08 7.858227e-08 3.493314e-08 4.099982e-09 1.127813e-08
85CD14+ Mono 2.644406e-08 8.614599e-07 1.104207e-06 1.010590e-08 6.330077e-07 9.410925e-07 4.183726e-07 4.910372e-08 1.350587e-07
86B 3.484106e-12 1.135035e-10 1.454883e-10 1.331490e-12 8.340397e-11 1.240029e-10 5.512432e-11 6.469743e-12 1.779699e-11
87CD8 T 1.688886e-09 5.501955e-08 7.052374e-08 6.454268e-10 4.042910e-08 6.010862e-08 2.672085e-08 3.136136e-09 8.626777e-09
88FCGR3A+ Mono 1.001622e-08 3.262944e-07 4.182392e-07 3.827814e-09 2.397636e-07 3.564543e-07 1.584664e-07 1.859898e-08 5.115538e-08
89NK 9.993118e-09 3.255413e-07 4.172736e-07 3.818983e-09 2.392101e-07 3.556306e-07 1.581005e-07 1.855606e-08 5.103702e-08
90DC 1.081809e-09 3.524210e-08 4.517296e-08 4.134256e-10 2.589626e-08 3.850073e-08 1.711561e-08 2.008818e-09 5.525462e-09
91Platelet 3.490750e-09 1.137069e-07 1.457441e-07 1.334030e-09 8.355043e-08 1.241924e-07 5.521981e-08 6.481446e-09 1.781969e-08
92
93, , WNT10A_FZD1_LRP5
94
95 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
96Naive CD4 T 0 0 0 0 0 0 0 0 0
97Memory CD4 T 0 0 0 0 0 0 0 0 0
98CD14+ Mono 0 0 0 0 0 0 0 0 0
99B 0 0 0 0 0 0 0 0 0
100CD8 T 0 0 0 0 0 0 0 0 0
101FCGR3A+ Mono 0 0 0 0 0 0 0 0 0
102NK 0 0 0 0 0 0 0 0 0
103DC 0 0 0 0 0 0 0 0 0
104Platelet 0 0 0 0 0 0 0 0 0
105
106, , WNT10A_FZD2_LRP5
107
108 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
109Naive CD4 T 0 0 0 0 0 0 0 0 0
110Memory CD4 T 0 0 0 0 0 0 0 0 0
111CD14+ Mono 0 0 0 0 0 0 0 0 0
112B 0 0 0 0 0 0 0 0 0
113CD8 T 0 0 0 0 0 0 0 0 0
114FCGR3A+ Mono 0 0 0 0 0 0 0 0 0
115NK 0 0 0 0 0 0 0 0 0
116DC 0 0 0 0 0 0 0 0 0
117Platelet 0 0 0 0 0 0 0 0 0
118
119, , WNT10A_FZD3_LRP5
120
121 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
122Naive CD4 T 0 0 0 0 0 0 0 0 0
123Memory CD4 T 0 0 0 0 0 0 0 0 0
124CD14+ Mono 0 0 0 0 0 0 0 0 0
125B 0 0 0 0 0 0 0 0 0
126CD8 T 0 0 0 0 0 0 0 0 0
127FCGR3A+ Mono 0 0 0 0 0 0 0 0 0
128NK 0 0 0 0 0 0 0 0 0
129DC 0 0 0 0 0 0 0 0 0
130Platelet 0 0 0 0 0 0 0 0 0
131
132, , WNT10A_FZD6_LRP5
133
134 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
135Naive CD4 T 0 0 0 0 0 0 0 0 0
136Memory CD4 T 0 0 0 0 0 0 0 0 0
137CD14+ Mono 0 0 0 0 0 0 0 0 0
138B 0 0 0 0 0 0 0 0 0
139CD8 T 0 0 0 0 0 0 0 0 0
140FCGR3A+ Mono 0 0 0 0 0 0 0 0 0
141NK 0 0 0 0 0 0 0 0 0
142DC 0 0 0 0 0 0 0 0 0
143Platelet 0 0 0 0 0 0 0 0 0
144
145, , WNT10B_FZD1_LRP5
146
147 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
148Naive CD4 T 0 0 0 0 0 0 0 0 0
149Memory CD4 T 0 0 0 0 0 0 0 0 0
150CD14+ Mono 0 0 0 0 0 0 0 0 0
151B 0 0 0 0 0 0 0 0 0
152CD8 T 0 0 0 0 0 0 0 0 0
153FCGR3A+ Mono 0 0 0 0 0 0 0 0 0
154NK 0 0 0 0 0 0 0 0 0
155DC 0 0 0 0 0 0 0 0 0
156Platelet 0 0 0 0 0 0 0 0 0
157
158, , WNT10B_FZD2_LRP5
159
160 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
161Naive CD4 T 0 0 0 0 0 0 0 0 0
162Memory CD4 T 0 0 0 0 0 0 0 0 0
163CD14+ Mono 0 0 0 0 0 0 0 0 0
164B 0 0 0 0 0 0 0 0 0
165CD8 T 0 0 0 0 0 0 0 0 0
166FCGR3A+ Mono 0 0 0 0 0 0 0 0 0
167NK 0 0 0 0 0 0 0 0 0
168DC 0 0 0 0 0 0 0 0 0
169Platelet 0 0 0 0 0 0 0 0 0
170
171, , WNT10B_FZD3_LRP5
172
173 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
174Naive CD4 T 0 0 0 0 0 0 0 0 0
175Memory CD4 T 0 0 0 0 0 0 0 0 0
176CD14+ Mono 0 0 0 0 0 0 0 0 0
177B 0 0 0 0 0 0 0 0 0
178CD8 T 0 0 0 0 0 0 0 0 0
179FCGR3A+ Mono 0 0 0 0 0 0 0 0 0
180NK 0 0 0 0 0 0 0 0 0
181DC 0 0 0 0 0 0 0 0 0
182Platelet 0 0 0 0 0 0 0 0 0
183
184, , WNT10B_FZD6_LRP5
185
186 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
187Naive CD4 T 0 0 0 0 0 0 0 0 0
188Memory CD4 T 0 0 0 0 0 0 0 0 0
189CD14+ Mono 0 0 0 0 0 0 0 0 0
190B 0 0 0 0 0 0 0 0 0
191CD8 T 0 0 0 0 0 0 0 0 0
192FCGR3A+ Mono 0 0 0 0 0 0 0 0 0
193NK 0 0 0 0 0 0 0 0 0
194DC 0 0 0 0 0 0 0 0 0
195Platelet 0 0 0 0 0 0 0 0 0
196
197, , WNT16_FZD1_LRP5
198
199 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
200Naive CD4 T 0 0 0 0 0 0 0 0 0
201Memory CD4 T 0 0 0 0 0 0 0 0 0
202CD14+ Mono 0 0 0 0 0 0 0 0 0
203
204 [ reached getOption("max.print") -- omitted 6 row(s) and 114 matrix slice(s) ]
205
206$pval
207, , TGFB1_TGFBR1_TGFBR2
208
209 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
210Naive CD4 T 1.00 1.00 1.00 1.00 1.00 1.00 1.00 0.99 0.56
211Memory CD4 T 1.00 0.67 0.39 1.00 0.33 0.00 0.15 0.44 0.01
212CD14+ Mono 0.87 0.00 0.00 0.80 0.00 0.00 0.00 0.00 0.00
213B 1.00 1.00 0.98 1.00 0.99 0.95 0.95 0.99 0.69
214CD8 T 1.00 0.36 0.04 0.99 0.07 0.00 0.00 0.44 0.00
215FCGR3A+ Mono 0.73 0.00 0.00 0.67 0.00 0.00 0.00 0.00 0.00
216NK 0.74 0.00 0.00 0.70 0.00 0.00 0.00 0.01 0.00
217DC 0.70 0.20 0.21 0.68 0.22 0.00 0.10 0.26 0.01
218Platelet 0.52 0.00 0.00 0.48 0.00 0.00 0.00 0.00 0.00
219
220, , TGFB1_ACVR1B_TGFBR2
221
222 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
223Naive CD4 T 1.00 1.00 1.00 1.00 1.00 1.00 1.00 0.99 0.94
224Memory CD4 T 1.00 0.73 0.39 1.00 0.48 0.00 0.24 0.46 0.02
225CD14+ Mono 0.87 0.00 0.00 0.78 0.00 0.00 0.00 0.00 0.00
226B 1.00 1.00 0.99 1.00 0.99 0.97 0.96 0.99 0.92
227CD8 T 0.92 0.39 0.04 0.93 0.16 0.00 0.00 0.45 0.00
228FCGR3A+ Mono 0.71 0.00 0.00 0.67 0.00 0.00 0.00 0.00 0.00
229NK 0.75 0.00 0.00 0.70 0.00 0.00 0.00 0.01 0.00
230DC 0.66 0.21 0.21 0.64 0.23 0.00 0.10 0.26 0.01
231Platelet 0.42 0.00 0.00 0.42 0.00 0.00 0.00 0.00 0.00
232
233, , TGFB1_ACVR1C_TGFBR2
234
235 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
236Naive CD4 T 1 1 1 1 1 1 1 1 1
237Memory CD4 T 1 1 1 1 1 1 1 1 1
238CD14+ Mono 1 1 1 1 1 1 1 1 1
239B 1 1 1 1 1 1 1 1 1
240CD8 T 1 1 1 1 1 1 1 1 1
241FCGR3A+ Mono 1 1 1 1 1 1 1 1 1
242NK 1 1 1 1 1 1 1 1 1
243DC 1 1 1 1 1 1 1 1 1
244Platelet 1 1 1 1 1 1 1 1 1
245
246, , TGFB1_ACVR1_TGFBR1
247
248 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
249Naive CD4 T 1.00 1.00 1.00 1.00 1.00 1.00 1.00 0.99 0.95
250Memory CD4 T 1.00 0.75 0.46 1.00 0.38 0.00 0.22 0.47 0.02
251CD14+ Mono 0.88 0.00 0.00 0.80 0.00 0.00 0.00 0.00 0.00
252B 1.00 1.00 1.00 1.00 0.99 0.97 0.98 0.99 0.91
253CD8 T 0.92 0.38 0.05 0.92 0.05 0.00 0.00 0.46 0.00
254FCGR3A+ Mono 0.71 0.00 0.00 0.67 0.00 0.00 0.00 0.00 0.00
255NK 0.76 0.00 0.00 0.71 0.00 0.00 0.00 0.02 0.00
256DC 0.66 0.21 0.23 0.63 0.23 0.00 0.12 0.25 0.01
257Platelet 0.40 0.00 0.00 0.40 0.00 0.00 0.00 0.00 0.00
258
259, , WNT10A_FZD1_LRP5
260
261 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
262Naive CD4 T 1 1 1 1 1 1 1 1 1
263Memory CD4 T 1 1 1 1 1 1 1 1 1
264CD14+ Mono 1 1 1 1 1 1 1 1 1
265B 1 1 1 1 1 1 1 1 1
266CD8 T 1 1 1 1 1 1 1 1 1
267FCGR3A+ Mono 1 1 1 1 1 1 1 1 1
268NK 1 1 1 1 1 1 1 1 1
269DC 1 1 1 1 1 1 1 1 1
270Platelet 1 1 1 1 1 1 1 1 1
271
272, , WNT10A_FZD2_LRP5
273
274 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
275Naive CD4 T 1 1 1 1 1 1 1 1 1
276Memory CD4 T 1 1 1 1 1 1 1 1 1
277CD14+ Mono 1 1 1 1 1 1 1 1 1
278B 1 1 1 1 1 1 1 1 1
279CD8 T 1 1 1 1 1 1 1 1 1
280FCGR3A+ Mono 1 1 1 1 1 1 1 1 1
281NK 1 1 1 1 1 1 1 1 1
282DC 1 1 1 1 1 1 1 1 1
283Platelet 1 1 1 1 1 1 1 1 1
284
285, , WNT10A_FZD3_LRP5
286
287 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
288Naive CD4 T 1 1 1 1 1 1 1 1 1
289Memory CD4 T 1 1 1 1 1 1 1 1 1
290CD14+ Mono 1 1 1 1 1 1 1 1 1
291B 1 1 1 1 1 1 1 1 1
292CD8 T 1 1 1 1 1 1 1 1 1
293FCGR3A+ Mono 1 1 1 1 1 1 1 1 1
294NK 1 1 1 1 1 1 1 1 1
295DC 1 1 1 1 1 1 1 1 1
296Platelet 1 1 1 1 1 1 1 1 1
297
298, , WNT10A_FZD6_LRP5
299
300 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
301Naive CD4 T 1 1 1 1 1 1 1 1 1
302Memory CD4 T 1 1 1 1 1 1 1 1 1
303CD14+ Mono 1 1 1 1 1 1 1 1 1
304B 1 1 1 1 1 1 1 1 1
305CD8 T 1 1 1 1 1 1 1 1 1
306FCGR3A+ Mono 1 1 1 1 1 1 1 1 1
307NK 1 1 1 1 1 1 1 1 1
308DC 1 1 1 1 1 1 1 1 1
309Platelet 1 1 1 1 1 1 1 1 1
310
311, , WNT10B_FZD1_LRP5
312
313 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
314Naive CD4 T 1 1 1 1 1 1 1 1 1
315Memory CD4 T 1 1 1 1 1 1 1 1 1
316CD14+ Mono 1 1 1 1 1 1 1 1 1
317B 1 1 1 1 1 1 1 1 1
318CD8 T 1 1 1 1 1 1 1 1 1
319FCGR3A+ Mono 1 1 1 1 1 1 1 1 1
320NK 1 1 1 1 1 1 1 1 1
321DC 1 1 1 1 1 1 1 1 1
322Platelet 1 1 1 1 1 1 1 1 1
323
324, , WNT10B_FZD2_LRP5
325
326 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
327Naive CD4 T 1 1 1 1 1 1 1 1 1
328Memory CD4 T 1 1 1 1 1 1 1 1 1
329CD14+ Mono 1 1 1 1 1 1 1 1 1
330B 1 1 1 1 1 1 1 1 1
331CD8 T 1 1 1 1 1 1 1 1 1
332FCGR3A+ Mono 1 1 1 1 1 1 1 1 1
333NK 1 1 1 1 1 1 1 1 1
334DC 1 1 1 1 1 1 1 1 1
335Platelet 1 1 1 1 1 1 1 1 1
336
337, , WNT10B_FZD3_LRP5
338
339 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
340Naive CD4 T 1 1 1 1 1 1 1 1 1
341Memory CD4 T 1 1 1 1 1 1 1 1 1
342CD14+ Mono 1 1 1 1 1 1 1 1 1
343B 1 1 1 1 1 1 1 1 1
344CD8 T 1 1 1 1 1 1 1 1 1
345FCGR3A+ Mono 1 1 1 1 1 1 1 1 1
346NK 1 1 1 1 1 1 1 1 1
347DC 1 1 1 1 1 1 1 1 1
348Platelet 1 1 1 1 1 1 1 1 1
349
350, , WNT10B_FZD6_LRP5
351
352 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
353Naive CD4 T 1 1 1 1 1 1 1 1 1
354Memory CD4 T 1 1 1 1 1 1 1 1 1
355CD14+ Mono 1 1 1 1 1 1 1 1 1
356B 1 1 1 1 1 1 1 1 1
357CD8 T 1 1 1 1 1 1 1 1 1
358FCGR3A+ Mono 1 1 1 1 1 1 1 1 1
359NK 1 1 1 1 1 1 1 1 1
360DC 1 1 1 1 1 1 1 1 1
361Platelet 1 1 1 1 1 1 1 1 1
362
363, , WNT16_FZD1_LRP5
364
365 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
366Naive CD4 T 1 1 1 1 1 1 1 1 1
367Memory CD4 T 1 1 1 1 1 1 1 1 1
368CD14+ Mono 1 1 1 1 1 1 1 1 1
369
370 [ reached getOption("max.print") -- omitted 6 row(s) and 114 matrix slice(s) ]
371
372$count
373 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
374Naive CD4 T 4 9 15 5 11 21 12 14 6
375Memory CD4 T 13 21 22 9 22 31 21 21 13
376CD14+ Mono 12 20 25 12 23 28 26 28 14
377B 3 6 11 4 6 17 9 11 6
378CD8 T 7 13 22 7 15 27 20 19 12
379FCGR3A+ Mono 12 25 28 12 22 33 26 30 15
380NK 10 19 24 9 20 26 21 23 12
381DC 13 24 25 13 21 32 22 26 18
382Platelet 2 6 10 2 10 11 10 11 9
383
384$sum
385 Naive CD4 T Memory CD4 T CD14+ Mono B CD8 T FCGR3A+ Mono NK DC Platelet
386Naive CD4 T 5.235731e-04 6.742952e-04 3.909235e-04 7.501420e-04 5.434838e-04 3.885489e-04 1.610436e-04 4.846562e-05 4.413932e-06
387Memory CD4 T 1.007867e-03 1.385925e-03 6.727733e-04 1.319087e-03 1.129907e-03 6.201049e-04 4.244407e-04 1.029006e-04 2.323768e-05
388CD14+ Mono 2.212146e-04 3.583798e-04 1.213175e-03 5.313253e-04 5.061446e-04 5.027468e-04 2.294104e-04 8.682125e-05 2.022770e-05
389B 1.301160e-05 9.973032e-05 1.565374e-04 3.703069e-04 1.646528e-04 2.057724e-04 4.275688e-05 2.459992e-05 3.097154e-06
390CD8 T 7.640382e-04 9.283023e-04 4.849123e-04 6.086610e-04 1.986549e-03 1.788599e-04 8.787072e-04 5.912427e-05 9.023021e-05
391FCGR3A+ Mono 1.374292e-04 2.766033e-04 4.453398e-04 1.984605e-04 1.309001e-04 2.772841e-04 6.165247e-05 3.351834e-05 9.078602e-07
392NK 4.436511e-04 4.983154e-04 3.013077e-04 3.858570e-04 1.078647e-03 9.820542e-05 4.720637e-04 3.638077e-05 4.795777e-05
393DC 3.642583e-05 8.053200e-05 1.016134e-04 9.111682e-05 6.074735e-05 6.164358e-05 2.886705e-05 1.000832e-05 1.323708e-06
394Platelet 2.580361e-05 3.406017e-05 1.414725e-05 1.492857e-05 9.745813e-05 3.867913e-06 4.425967e-05 2.105407e-06 4.930773e-06
395
396head(cellchat@netP$similarity)
397head(cellchat@net$count)
398head(cellchat@net$prob)
399head(cellchat@net$sum)
400head(cellchat@DB)
401head(cellchat@var.features)
github 仓库在：https://github.com/sqjin/CellCha
References
[1] Cell-Cell Interaction Database|| 单细胞配受体库你还在文章的附录里找吗？: https://www.jianshu.com/p/49613adce465
[2] 听说你的单细胞对象需要一个思维导图？: https://www.jianshu.com/p/7560f4fd0d77