Inference and analysis of cell-cell communication using CellChat-Melange
Nasir Mahmood Abbasi
2023-09-12
1. load libraries
2. Load Seurat Object
#Load Seurat Object Integrated Object
load("Melange_v3.Robj")
MelangeAn object of class Seurat
59684 features across 81220 samples within 3 assays
Active assay: SCT (23076 features, 3000 variable features)
3 layers present: counts, data, scale.data
2 other assays present: RNA, ADT
2 dimensional reductions calculated: pca, umapDefaultAssay(Melange) <- "SCT"
3. Data input & processing and initialization of CellChat object
#------------------------------------------------------------------------------
# **Part I: Data input & processing and initialization of CellChat object**
#------------------------------------------------------------------------------
# Loading and Preparing Data
Melange <- AddMetaData(object = Melange,
metadata = paste0("Cluster_", Melange$seurat_clusters),
col.name = 'Clusters')
#Melange[["RNA"]] <- JoinLayers(Melange[["RNA"]])
# Extract the normalized data
#data_input <- GetAssayData(Melange, assay = "RNA", slot = "data")
data.input <- Melange@assays$SCT$data
#1. Create CellChat object
cellchat <- createCellChat(object = data.input, meta = Melange@meta.data, group.by = "Clusters")[1] "Create a CellChat object from a data matrix"Warning: The 'meta' data does not have a column named `samples`. We now add this column and all cells are assumed to belong to `sample1`! Set cell identities for the new CellChat object
The cell groups used for CellChat analysis are Cluster_0, Cluster_1, Cluster_10, Cluster_11, Cluster_12, Cluster_13, Cluster_14, Cluster_15, Cluster_16, Cluster_17, Cluster_18, Cluster_19, Cluster_2, Cluster_20, Cluster_21, Cluster_22, Cluster_23, Cluster_24, Cluster_25, Cluster_26, Cluster_27, Cluster_28, Cluster_29, Cluster_3, Cluster_4, Cluster_5, Cluster_6, Cluster_7, Cluster_8, Cluster_9 #2. Set the ligand-receptor interaction database
cellchatDB <- CellChatDB.human # or CellChatDB.mouse
showDatabaseCategory(cellchatDB)
# Show the structure of the database
dplyr::glimpse(cellchatDB$interaction)Rows: 3,233
Columns: 28
$ interaction_name <chr> "TGFB1_TGFBR1_TGFBR2", "TGFB2_TGFBR1_TGFBR2", "TGFB3_TGFBR1_TGFBR2", "TGFB1_ACVR1B_TGFBR2", "TGFB…
$ pathway_name <chr> "TGFb", "TGFb", "TGFb", "TGFb", "TGFb", "TGFb", "TGFb", "TGFb", "TGFb", "TGFb", "TGFb", "TGFb", "…
$ ligand <chr> "TGFB1", "TGFB2", "TGFB3", "TGFB1", "TGFB1", "TGFB2", "TGFB2", "TGFB3", "TGFB3", "TGFB1", "TGFB2"…
$ receptor <chr> "TGFbR1_R2", "TGFbR1_R2", "TGFbR1_R2", "ACVR1B_TGFbR2", "ACVR1C_TGFbR2", "ACVR1B_TGFbR2", "ACVR1C…
$ agonist <chr> "TGFb agonist", "TGFb agonist", "TGFb agonist", "TGFb agonist", "TGFb agonist", "TGFb agonist", "…
$ antagonist <chr> "TGFb antagonist", "TGFb antagonist", "TGFb antagonist", "TGFb antagonist", "TGFb antagonist", "T…
$ co_A_receptor <chr> "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "…
$ co_I_receptor <chr> "TGFb inhibition receptor", "TGFb inhibition receptor", "TGFb inhibition receptor", "TGFb inhibit…
$ evidence <chr> "KEGG: hsa04350", "KEGG: hsa04350", "KEGG: hsa04350", "PMID: 27449815", "PMID: 27449815", "PMID: …
$ annotation <chr> "Secreted Signaling", "Secreted Signaling", "Secreted Signaling", "Secreted Signaling", "Secreted…
$ interaction_name_2 <chr> "TGFB1 - (TGFBR1+TGFBR2)", "TGFB2 - (TGFBR1+TGFBR2)", "TGFB3 - (TGFBR1+TGFBR2)", "TGFB1 - (ACVR1B…
$ is_neurotransmitter <lgl> FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE,…
$ ligand.symbol <chr> "TGFB1", "TGFB2", "TGFB3", "TGFB1", "TGFB1", "TGFB2", "TGFB2", "TGFB3", "TGFB3", "TGFB1", "TGFB2"…
$ ligand.family <chr> "TGF-beta", "TGF-beta", "TGF-beta", "TGF-beta", "TGF-beta", "TGF-beta", "TGF-beta", "TGF-beta", "…
$ ligand.location <chr> "Extracellular matrix, Secreted, Extracellular space", "Extracellular matrix, Secreted, Extracell…
$ ligand.keyword <chr> "Disease variant, Signal, Reference proteome, Extracellular matrix, Disulfide bond, Secreted, Mit…
$ ligand.secreted_type <chr> "growth factor", "growth factor", "cytokine;growth factor", "growth factor", "growth factor", "gr…
$ ligand.transmembrane <lgl> FALSE, FALSE, TRUE, FALSE, FALSE, FALSE, FALSE, TRUE, TRUE, FALSE, FALSE, TRUE, FALSE, FALSE, FAL…
$ receptor.symbol <chr> "TGFBR2, TGFBR1", "TGFBR2, TGFBR1", "TGFBR2, TGFBR1", "TGFBR2, ACVR1B", "TGFBR2, ACVR1C", "TGFBR2…
$ receptor.family <chr> "Protein kinase superfamily, TKL Ser/Thr protein kinase", "Protein kinase superfamily, TKL Ser/Th…
$ receptor.location <chr> "Cell membrane, Secreted, Membrane raft, Cell surface, Cell junction, Tight junction", "Cell memb…
$ receptor.keyword <chr> "Membrane, Secreted, Disulfide bond, Kinase, Transmembrane, Differentiation, ATP-binding, Recepto…
$ receptor.surfaceome_main <chr> "Receptors", "Receptors", "Receptors", "Receptors", "Receptors", "Receptors", "Receptors", "Recep…
$ receptor.surfaceome_sub <chr> "Act.TGFB;Kinase", "Act.TGFB;Kinase", "Act.TGFB;Kinase", "Act.TGFB;Kinase", "Act.TGFB;Kinase", "A…
$ receptor.adhesome <chr> "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "…
$ receptor.secreted_type <chr> "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "…
$ receptor.transmembrane <lgl> TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, T…
$ version <chr> "CellChatDB v1", "CellChatDB v1", "CellChatDB v1", "CellChatDB v1", "CellChatDB v1", "CellChatDB …# use a subset of CellChatDB for cell-cell communication analysis
cellchatDB.use <- subsetDB(cellchatDB, search = "Secreted Signaling", key = "annotation") # use Secreted Signaling
# set the used database in the object
cellchat@DB <- cellchatDB.use
#------------------------------------------------------------------------------
# **Preprocessing the expression data for cell-cell communication analysis**
#------------------------------------------------------------------------------
#3. Preprocess the expression data
cellchat <- subsetData(cellchat)
cellchat <- identifyOverExpressedGenes(cellchat)
cellchat <- identifyOverExpressedInteractions(cellchat)
| | 0 % ~calculating
|+ | 1 % ~00s
|++ | 2 % ~00s
|++ | 4 % ~00s
|+++ | 5 % ~00s
|+++ | 6 % ~00s
|++++ | 7 % ~00s
|+++++ | 8 % ~00s
|+++++ | 9 % ~00s
|++++++ | 11% ~00s
|++++++ | 12% ~00s
|+++++++ | 13% ~00s
|++++++++ | 14% ~00s
|++++++++ | 15% ~00s
|+++++++++ | 16% ~00s
|+++++++++ | 18% ~00s
|++++++++++ | 19% ~00s
|++++++++++ | 20% ~00s
|+++++++++++ | 21% ~00s
|++++++++++++ | 22% ~00s
|++++++++++++ | 24% ~00s
|+++++++++++++ | 25% ~00s
|+++++++++++++ | 26% ~00s
|++++++++++++++ | 27% ~00s
|+++++++++++++++ | 28% ~00s
|+++++++++++++++ | 29% ~00s
|++++++++++++++++ | 31% ~00s
|++++++++++++++++ | 32% ~00s
|+++++++++++++++++ | 33% ~00s
|++++++++++++++++++ | 34% ~00s
|++++++++++++++++++ | 35% ~00s
|+++++++++++++++++++ | 36% ~00s
|+++++++++++++++++++ | 38% ~00s
|++++++++++++++++++++ | 39% ~00s
|++++++++++++++++++++ | 40% ~00s
|+++++++++++++++++++++ | 41% ~00s
|++++++++++++++++++++++ | 42% ~00s
|++++++++++++++++++++++ | 44% ~00s
|+++++++++++++++++++++++ | 45% ~00s
|+++++++++++++++++++++++ | 46% ~00s
|++++++++++++++++++++++++ | 47% ~00s
|+++++++++++++++++++++++++ | 48% ~00s
|+++++++++++++++++++++++++ | 49% ~00s
|++++++++++++++++++++++++++ | 51% ~00s
|++++++++++++++++++++++++++ | 52% ~00s
|+++++++++++++++++++++++++++ | 53% ~00s
|++++++++++++++++++++++++++++ | 54% ~00s
|++++++++++++++++++++++++++++ | 55% ~00s
|+++++++++++++++++++++++++++++ | 56% ~00s
|+++++++++++++++++++++++++++++ | 58% ~00s
|++++++++++++++++++++++++++++++ | 59% ~00s
|++++++++++++++++++++++++++++++ | 60% ~00s
|+++++++++++++++++++++++++++++++ | 61% ~00s
|++++++++++++++++++++++++++++++++ | 62% ~00s
|++++++++++++++++++++++++++++++++ | 64% ~00s
|+++++++++++++++++++++++++++++++++ | 65% ~00s
|+++++++++++++++++++++++++++++++++ | 66% ~00s
|++++++++++++++++++++++++++++++++++ | 67% ~00s
|+++++++++++++++++++++++++++++++++++ | 68% ~00s
|+++++++++++++++++++++++++++++++++++ | 69% ~00s
|++++++++++++++++++++++++++++++++++++ | 71% ~00s
|++++++++++++++++++++++++++++++++++++ | 72% ~00s
|+++++++++++++++++++++++++++++++++++++ | 73% ~00s
|++++++++++++++++++++++++++++++++++++++ | 74% ~00s
|++++++++++++++++++++++++++++++++++++++ | 75% ~00s
|+++++++++++++++++++++++++++++++++++++++ | 76% ~00s
|+++++++++++++++++++++++++++++++++++++++ | 78% ~00s
|++++++++++++++++++++++++++++++++++++++++ | 79% ~00s
|++++++++++++++++++++++++++++++++++++++++ | 80% ~00s
|+++++++++++++++++++++++++++++++++++++++++ | 81% ~00s
|++++++++++++++++++++++++++++++++++++++++++ | 82% ~00s
|++++++++++++++++++++++++++++++++++++++++++ | 84% ~00s
|+++++++++++++++++++++++++++++++++++++++++++ | 85% ~00s
|+++++++++++++++++++++++++++++++++++++++++++ | 86% ~00s
|++++++++++++++++++++++++++++++++++++++++++++ | 87% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++ | 88% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++ | 89% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 91% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 92% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++ | 93% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 94% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 95% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++++ | 96% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=00s
| | 0 % ~calculating
|+ | 1 % ~00s
|++ | 2 % ~00s
|++ | 4 % ~00s
|+++ | 5 % ~00s
|+++ | 6 % ~00s
|++++ | 7 % ~00s
|+++++ | 8 % ~00s
|+++++ | 9 % ~00s
|++++++ | 11% ~00s
|++++++ | 12% ~00s
|+++++++ | 13% ~00s
|++++++++ | 14% ~00s
|++++++++ | 15% ~00s
|+++++++++ | 16% ~00s
|+++++++++ | 18% ~00s
|++++++++++ | 19% ~00s
|++++++++++ | 20% ~00s
|+++++++++++ | 21% ~00s
|++++++++++++ | 22% ~00s
|++++++++++++ | 24% ~00s
|+++++++++++++ | 25% ~00s
|+++++++++++++ | 26% ~00s
|++++++++++++++ | 27% ~00s
|+++++++++++++++ | 28% ~00s
|+++++++++++++++ | 29% ~00s
|++++++++++++++++ | 31% ~00s
|++++++++++++++++ | 32% ~00s
|+++++++++++++++++ | 33% ~00s
|++++++++++++++++++ | 34% ~00s
|++++++++++++++++++ | 35% ~00s
|+++++++++++++++++++ | 36% ~00s
|+++++++++++++++++++ | 38% ~00s
|++++++++++++++++++++ | 39% ~00s
|++++++++++++++++++++ | 40% ~00s
|+++++++++++++++++++++ | 41% ~00s
|++++++++++++++++++++++ | 42% ~00s
|++++++++++++++++++++++ | 44% ~00s
|+++++++++++++++++++++++ | 45% ~00s
|+++++++++++++++++++++++ | 46% ~00s
|++++++++++++++++++++++++ | 47% ~00s
|+++++++++++++++++++++++++ | 48% ~00s
|+++++++++++++++++++++++++ | 49% ~00s
|++++++++++++++++++++++++++ | 51% ~00s
|++++++++++++++++++++++++++ | 52% ~00s
|+++++++++++++++++++++++++++ | 53% ~00s
|++++++++++++++++++++++++++++ | 54% ~00s
|++++++++++++++++++++++++++++ | 55% ~00s
|+++++++++++++++++++++++++++++ | 56% ~00s
|+++++++++++++++++++++++++++++ | 58% ~00s
|++++++++++++++++++++++++++++++ | 59% ~00s
|++++++++++++++++++++++++++++++ | 60% ~00s
|+++++++++++++++++++++++++++++++ | 61% ~00s
|++++++++++++++++++++++++++++++++ | 62% ~00s
|++++++++++++++++++++++++++++++++ | 64% ~00s
|+++++++++++++++++++++++++++++++++ | 65% ~00s
|+++++++++++++++++++++++++++++++++ | 66% ~00s
|++++++++++++++++++++++++++++++++++ | 67% ~00s
|+++++++++++++++++++++++++++++++++++ | 68% ~00s
|+++++++++++++++++++++++++++++++++++ | 69% ~00s
|++++++++++++++++++++++++++++++++++++ | 71% ~00s
|++++++++++++++++++++++++++++++++++++ | 72% ~00s
|+++++++++++++++++++++++++++++++++++++ | 73% ~00s
|++++++++++++++++++++++++++++++++++++++ | 74% ~00s
|++++++++++++++++++++++++++++++++++++++ | 75% ~00s
|+++++++++++++++++++++++++++++++++++++++ | 76% ~00s
|+++++++++++++++++++++++++++++++++++++++ | 78% ~00s
|++++++++++++++++++++++++++++++++++++++++ | 79% ~00s
|++++++++++++++++++++++++++++++++++++++++ | 80% ~00s
|+++++++++++++++++++++++++++++++++++++++++ | 81% ~00s
|++++++++++++++++++++++++++++++++++++++++++ | 82% ~00s
|++++++++++++++++++++++++++++++++++++++++++ | 84% ~00s
|+++++++++++++++++++++++++++++++++++++++++++ | 85% ~00s
|+++++++++++++++++++++++++++++++++++++++++++ | 86% ~00s
|++++++++++++++++++++++++++++++++++++++++++++ | 87% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++ | 88% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++ | 89% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 91% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 92% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++ | 93% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 94% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 95% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++++ | 96% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=00s
| | 0 % ~calculating
|+ | 1 % ~00s
|++ | 2 % ~00s
|++ | 3 % ~00s
|+++ | 4 % ~00s
|+++ | 5 % ~00s
|++++ | 6 % ~00s
|++++ | 7 % ~00s
|+++++ | 8 % ~00s
|+++++ | 9 % ~00s
|++++++ | 10% ~00s
|++++++ | 11% ~00s
|+++++++ | 12% ~00s
|+++++++ | 13% ~00s
|++++++++ | 14% ~00s
|++++++++ | 15% ~00s
|+++++++++ | 16% ~00s
|+++++++++ | 17% ~00s
|++++++++++ | 18% ~00s
|++++++++++ | 19% ~00s
|+++++++++++ | 20% ~00s
|+++++++++++ | 21% ~00s
|++++++++++++ | 22% ~00s
|++++++++++++ | 23% ~00s
|+++++++++++++ | 24% ~00s
|+++++++++++++ | 25% ~00s
|++++++++++++++ | 26% ~00s
|++++++++++++++ | 27% ~00s
|+++++++++++++++ | 28% ~00s
|+++++++++++++++ | 29% ~00s
|++++++++++++++++ | 30% ~00s
|++++++++++++++++ | 31% ~00s
|+++++++++++++++++ | 32% ~00s
|+++++++++++++++++ | 33% ~00s
|++++++++++++++++++ | 34% ~00s
|++++++++++++++++++ | 35% ~00s
|+++++++++++++++++++ | 36% ~00s
|+++++++++++++++++++ | 37% ~00s
|++++++++++++++++++++ | 38% ~00s
|++++++++++++++++++++ | 39% ~00s
|+++++++++++++++++++++ | 40% ~00s
|+++++++++++++++++++++ | 41% ~00s
|++++++++++++++++++++++ | 42% ~00s
|++++++++++++++++++++++ | 43% ~00s
|+++++++++++++++++++++++ | 44% ~00s
|+++++++++++++++++++++++ | 45% ~00s
|++++++++++++++++++++++++ | 46% ~00s
|++++++++++++++++++++++++ | 47% ~00s
|+++++++++++++++++++++++++ | 48% ~00s
|+++++++++++++++++++++++++ | 49% ~00s
|++++++++++++++++++++++++++ | 51% ~00s
|++++++++++++++++++++++++++ | 52% ~00s
|+++++++++++++++++++++++++++ | 53% ~00s
|+++++++++++++++++++++++++++ | 54% ~00s
|++++++++++++++++++++++++++++ | 55% ~00s
|++++++++++++++++++++++++++++ | 56% ~00s
|+++++++++++++++++++++++++++++ | 57% ~00s
|+++++++++++++++++++++++++++++ | 58% ~00s
|++++++++++++++++++++++++++++++ | 59% ~00s
|++++++++++++++++++++++++++++++ | 60% ~00s
|+++++++++++++++++++++++++++++++ | 61% ~00s
|+++++++++++++++++++++++++++++++ | 62% ~00s
|++++++++++++++++++++++++++++++++ | 63% ~00s
|++++++++++++++++++++++++++++++++ | 64% ~00s
|+++++++++++++++++++++++++++++++++ | 65% ~00s
|+++++++++++++++++++++++++++++++++ | 66% ~00s
|++++++++++++++++++++++++++++++++++ | 67% ~00s
|++++++++++++++++++++++++++++++++++ | 68% ~00s
|+++++++++++++++++++++++++++++++++++ | 69% ~00s
|+++++++++++++++++++++++++++++++++++ | 70% ~00s
|++++++++++++++++++++++++++++++++++++ | 71% ~00s
|++++++++++++++++++++++++++++++++++++ | 72% ~00s
|+++++++++++++++++++++++++++++++++++++ | 73% ~00s
|+++++++++++++++++++++++++++++++++++++ | 74% ~00s
|++++++++++++++++++++++++++++++++++++++ | 75% ~00s
|++++++++++++++++++++++++++++++++++++++ | 76% ~00s
|+++++++++++++++++++++++++++++++++++++++ | 77% ~00s
|+++++++++++++++++++++++++++++++++++++++ | 78% ~00s
|++++++++++++++++++++++++++++++++++++++++ | 79% ~00s
|++++++++++++++++++++++++++++++++++++++++ | 80% ~00s
|+++++++++++++++++++++++++++++++++++++++++ | 81% ~00s
|+++++++++++++++++++++++++++++++++++++++++ | 82% ~00s
|++++++++++++++++++++++++++++++++++++++++++ | 83% ~00s
|++++++++++++++++++++++++++++++++++++++++++ | 84% ~00s
|+++++++++++++++++++++++++++++++++++++++++++ | 85% ~00s
|+++++++++++++++++++++++++++++++++++++++++++ | 86% ~00s
|++++++++++++++++++++++++++++++++++++++++++++ | 87% ~00s
|++++++++++++++++++++++++++++++++++++++++++++ | 88% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++ | 89% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++ | 90% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 91% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 92% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++ | 93% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++ | 94% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 95% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 96% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=00s
The number of highly variable ligand-receptor pairs used for signaling inference is 675 4. Inference of cell-cell communication network
#------------------------------------------------------------------------------
# **Part II: Inference of cell-cell communication network**
#------------------------------------------------------------------------------
#------------------------------------------------------------------------------------
# **Compute the communication probability and infer cellular communication network**
#------------------------------------------------------------------------------------
#4. Compute communication probability and infer cellular communication network
cellchat <- computeCommunProb(cellchat)triMean is used for calculating the average gene expression per cell group.
[1] ">>> Run CellChat on sc/snRNA-seq data <<< [2024-09-12 15:29:38.877459]"
|
| | 0%
|
|= | 0%
|
|= | 1%
|
|== | 1%
|
|== | 2%
|
|=== | 2%
|
|=== | 3%
|
|==== | 3%
|
|==== | 4%
|
|===== | 4%
|
|====== | 5%
|
|======= | 5%
|
|======= | 6%
|
|======== | 6%
|
|======== | 7%
|
|========= | 7%
|
|========= | 8%
|
|========== | 8%
|
|========== | 9%
|
|=========== | 9%
|
|============ | 9%
|
|============ | 10%
|
|============= | 10%
|
|============= | 11%
|
|============== | 11%
|
|============== | 12%
|
|=============== | 12%
|
|=============== | 13%
|
|================ | 13%
|
|================= | 14%
|
|================== | 14%
|
|================== | 15%
|
|=================== | 15%
|
|=================== | 16%
|
|==================== | 16%
|
|==================== | 17%
|
|===================== | 17%
|
|====================== | 18%
|
|======================= | 19%
|
|======================== | 19%
|
|======================== | 20%
|
|========================= | 20%
|
|========================= | 21%
|
|========================== | 21%
|
|========================== | 22%
|
|=========================== | 22%
|
|=========================== | 23%
|
|============================ | 23%
|
|============================= | 23%
|
|============================= | 24%
|
|============================== | 24%
|
|============================== | 25%
|
|=============================== | 25%
|
|=============================== | 26%
|
|================================ | 26%
|
|================================ | 27%
|
|================================= | 27%
|
|================================== | 28%
|
|=================================== | 28%
|
|=================================== | 29%
|
|==================================== | 29%
|
|==================================== | 30%
|
|===================================== | 30%
|
|===================================== | 31%
|
|====================================== | 31%
|
|====================================== | 32%
|
|======================================= | 32%
|
|======================================== | 32%
|
|======================================== | 33%
|
|========================================= | 33%
|
|========================================= | 34%
|
|========================================== | 34%
|
|========================================== | 35%
|
|=========================================== | 35%
|
|=========================================== | 36%
|
|============================================ | 36%
|
|============================================= | 37%
|
|============================================== | 37%
|
|============================================== | 38%
|
|=============================================== | 38%
|
|=============================================== | 39%
|
|================================================ | 39%
|
|================================================ | 40%
|
|================================================= | 40%
|
|================================================== | 41%
|
|=================================================== | 41%
|
|=================================================== | 42%
|
|==================================================== | 42%
|
|==================================================== | 43%
|
|===================================================== | 43%
|
|===================================================== | 44%
|
|====================================================== | 44%
|
|====================================================== | 45%
|
|======================================================= | 45%
|
|======================================================== | 46%
|
|========================================================= | 46%
|
|========================================================= | 47%
|
|========================================================== | 47%
|
|========================================================== | 48%
|
|=========================================================== | 48%
|
|=========================================================== | 49%
|
|============================================================ | 49%
|
|============================================================= | 50%
|
|============================================================== | 51%
|
|=============================================================== | 51%
|
|=============================================================== | 52%
|
|================================================================ | 52%
|
|================================================================ | 53%
|
|================================================================= | 53%
|
|================================================================= | 54%
|
|================================================================== | 54%
|
|=================================================================== | 55%
|
|==================================================================== | 55%
|
|==================================================================== | 56%
|
|===================================================================== | 56%
|
|===================================================================== | 57%
|
|====================================================================== | 57%
|
|====================================================================== | 58%
|
|======================================================================= | 58%
|
|======================================================================= | 59%
|
|======================================================================== | 59%
|
|========================================================================= | 60%
|
|========================================================================== | 60%
|
|========================================================================== | 61%
|
|=========================================================================== | 61%
|
|=========================================================================== | 62%
|
|============================================================================ | 62%
|
|============================================================================ | 63%
|
|============================================================================= | 63%
|
|============================================================================== | 64%
|
|=============================================================================== | 64%
|
|=============================================================================== | 65%
|
|================================================================================ | 65%
|
|================================================================================ | 66%
|
|================================================================================= | 66%
|
|================================================================================= | 67%
|
|================================================================================== | 67%
|
|================================================================================== | 68%
|
|=================================================================================== | 68%
|
|==================================================================================== | 68%
|
|==================================================================================== | 69%
|
|===================================================================================== | 69%
|
|===================================================================================== | 70%
|
|====================================================================================== | 70%
|
|====================================================================================== | 71%
|
|======================================================================================= | 71%
|
|======================================================================================= | 72%
|
|======================================================================================== | 72%
|
|========================================================================================= | 73%
|
|========================================================================================== | 73%
|
|========================================================================================== | 74%
|
|=========================================================================================== | 74%
|
|=========================================================================================== | 75%
|
|============================================================================================ | 75%
|
|============================================================================================ | 76%
|
|============================================================================================= | 76%
|
|============================================================================================= | 77%
|
|============================================================================================== | 77%
|
|=============================================================================================== | 77%
|
|=============================================================================================== | 78%
|
|================================================================================================ | 78%
|
|================================================================================================ | 79%
|
|================================================================================================= | 79%
|
|================================================================================================= | 80%
|
|================================================================================================== | 80%
|
|================================================================================================== | 81%
|
|=================================================================================================== | 81%
|
|==================================================================================================== | 82%
|
|===================================================================================================== | 83%
|
|====================================================================================================== | 83%
|
|====================================================================================================== | 84%
|
|======================================================================================================= | 84%
|
|======================================================================================================= | 85%
|
|======================================================================================================== | 85%
|
|======================================================================================================== | 86%
|
|========================================================================================================= | 86%
|
|========================================================================================================== | 87%
|
|=========================================================================================================== | 87%
|
|=========================================================================================================== | 88%
|
|============================================================================================================ | 88%
|
|============================================================================================================ | 89%
|
|============================================================================================================= | 89%
|
|============================================================================================================= | 90%
|
|============================================================================================================== | 90%
|
|============================================================================================================== | 91%
|
|=============================================================================================================== | 91%
|
|================================================================================================================ | 91%
|
|================================================================================================================ | 92%
|
|================================================================================================================= | 92%
|
|================================================================================================================= | 93%
|
|================================================================================================================== | 93%
|
|================================================================================================================== | 94%
|
|=================================================================================================================== | 94%
|
|=================================================================================================================== | 95%
|
|==================================================================================================================== | 95%
|
|===================================================================================================================== | 96%
|
|====================================================================================================================== | 96%
|
|====================================================================================================================== | 97%
|
|======================================================================================================================= | 97%
|
|======================================================================================================================= | 98%
|
|======================================================================================================================== | 98%
|
|======================================================================================================================== | 99%
|
|========================================================================================================================= | 99%
|
|========================================================================================================================= | 100%
|
|==========================================================================================================================| 100%[1] ">>> CellChat inference is done. Parameter values are stored in `object@options$parameter` <<< [2024-09-12 15:37:26.883988]"cellchat <- filterCommunication(cellchat, min.cells = 10)
#------------------------------------------------------------------------------
# **Infer the cell-cell communication at a signaling pathway level**
#------------------------------------------------------------------------------
#5. Compute communication at the pathway level
cellchat <- computeCommunProbPathway(cellchat)
#------------------------------------------------------------------------------
# **Calculate the aggregated cell-cell communication network**
#------------------------------------------------------------------------------
#6. Aggregate the inferred cell-cell communication network
cellchat <- aggregateNet(cellchat)
#7. Visualize the communication network
groupSize <- as.numeric(table(cellchat@idents))
netVisual_circle(cellchat@net$count, vertex.weight = groupSize, weight.scale = T, label.edge= F, title.name = "Number of interactions")
netVisual_circle(cellchat@net$weight, vertex.weight = groupSize, weight.scale = T, label.edge= F, title.name = "Interaction weights/strength")
mat <- cellchat@net$weight
par(mfrow = c(3,4), xpd=TRUE)
for (i in 1:nrow(mat)) {
mat2 <- matrix(0, nrow = nrow(mat), ncol = ncol(mat), dimnames = dimnames(mat))
mat2[i, ] <- mat[i, ]
netVisual_circle(mat2, vertex.weight = groupSize, weight.scale = T, edge.weight.max = max(mat), title.name = rownames(mat)[i])
}
NA
5. Visualization of cell-cell communication network
#------------------------------------------------------------------------------
# **Part III: Visualization of cell-cell communication network**
#------------------------------------------------------------------------------
#----------------------------------------------------------------------------------------------------
# **Visualize each signaling pathway using Hierarchy plot, Circle plot or Chord diagram**
#----------------------------------------------------------------------------------------------------
pathways.show <- c("IL2")
# Hierarchy plot
# Here we define `vertex.receive` so that the left portion of the hierarchy plot shows signaling to fibroblast and the right portion shows signaling to immune cells
vertex.receiver = seq(1,4) # a numeric vector.
netVisual_aggregate(cellchat, signaling = pathways.show, vertex.receiver = vertex.receiver)
# Circle plot
par(mfrow=c(1,1))
netVisual_aggregate(cellchat, signaling = pathways.show, layout = "circle")
# Chord diagram
par(mfrow=c(1,1))
netVisual_aggregate(cellchat, signaling = pathways.show, layout = "chord")
# Heatmap
par(mfrow=c(1,1))
netVisual_heatmap(cellchat, signaling = pathways.show, color.heatmap = "Reds")Do heatmap based on a single object 
#> Do heatmap based on a single object
# Get the levels of cell types
cell_types <- levels(cellchat@idents)
# Chord diagram
group.cellType <- c(rep("MK", 4), rep("IL2", 4), rep("TNF", 4))
group.cellType <- rep(NA, length(cell_types))
# Assign each cell type to a group
# You'll need to adjust this based on your actual data and biological knowledge
group.cellType[1:8] <- "MK"
group.cellType[9:16] <- "IL2"
group.cellType[17:30] <- "TNF"
# Name the vector with the cell types
names(group.cellType) <- cell_types
# Check the result
print(group.cellType) Cluster_0 Cluster_1 Cluster_10 Cluster_11 Cluster_12 Cluster_13 Cluster_14 Cluster_15 Cluster_16 Cluster_17 Cluster_18 Cluster_19
"MK" "MK" "MK" "MK" "MK" "MK" "MK" "MK" "IL2" "IL2" "IL2" "IL2"
Cluster_2 Cluster_20 Cluster_21 Cluster_22 Cluster_23 Cluster_24 Cluster_25 Cluster_26 Cluster_27 Cluster_28 Cluster_29 Cluster_3
"IL2" "IL2" "IL2" "IL2" "TNF" "TNF" "TNF" "TNF" "TNF" "TNF" "TNF" "TNF"
Cluster_4 Cluster_5 Cluster_6 Cluster_7 Cluster_8 Cluster_9
"TNF" "TNF" "TNF" "TNF" "TNF" "TNF" netVisual_chord_cell(cellchat, signaling = pathways.show, group = group.cellType, title.name = paste0(pathways.show, " signaling network"))Plot the aggregated cell-cell communication network at the signaling pathway level
#> Plot the aggregated cell-cell communication network at the signaling pathway level
# Compute the contribution of each ligand-receptor pair to the overall signaling pathway and visualize cell-cell
#communication mediated by a single ligand-receptor pair
netAnalysis_contribution(cellchat, signaling = pathways.show)
pairLR.CXCL <- extractEnrichedLR(cellchat, signaling = pathways.show, geneLR.return = FALSE)
LR.show <- pairLR.CXCL[1,] # show one ligand-receptor pair
# Hierarchy plot
vertex.receiver = seq(1,4) # a numeric vector
netVisual_individual(cellchat, signaling = pathways.show, pairLR.use = LR.show, vertex.receiver = vertex.receiver)[[1]]
#> [[1]]
# Circle plot
netVisual_individual(cellchat, signaling = pathways.show, pairLR.use = LR.show, layout = "circle")[[1]]
#> [[1]]
# Chord diagram
netVisual_individual(cellchat, signaling = pathways.show, pairLR.use = LR.show, layout = "chord")[[1]]
# Access all the signaling pathways showing significant communications
pathways.show.all <- cellchat@netP$pathways
# check the order of cell identity to set suitable vertex.receiver
levels(cellchat@idents) [1] "Cluster_0" "Cluster_1" "Cluster_10" "Cluster_11" "Cluster_12" "Cluster_13" "Cluster_14" "Cluster_15" "Cluster_16"
[10] "Cluster_17" "Cluster_18" "Cluster_19" "Cluster_2" "Cluster_20" "Cluster_21" "Cluster_22" "Cluster_23" "Cluster_24"
[19] "Cluster_25" "Cluster_26" "Cluster_27" "Cluster_28" "Cluster_29" "Cluster_3" "Cluster_4" "Cluster_5" "Cluster_6"
[28] "Cluster_7" "Cluster_8" "Cluster_9" vertex.receiver = seq(1,4)
for (i in 1:length(pathways.show.all)) {
# Visualize communication network associated with both signaling pathway and individual L-R pairs
netVisual(cellchat, signaling = pathways.show.all[i], vertex.receiver = vertex.receiver, layout = "hierarchy")
# Compute and visualize the contribution of each ligand-receptor pair to the overall signaling pathway
gg <- netAnalysis_contribution(cellchat, signaling = pathways.show.all[i])
ggsave(filename=paste0(pathways.show.all[i], "_L-R_contribution.pdf"), plot=gg, width = 3, height = 2, units = 'in', dpi = 300)
}
#----------------------------------------------------------------------------------------------------
# **Visualize cell-cell communication mediated by multiple ligand-receptors or signaling pathways**
#----------------------------------------------------------------------------------------------------
#(A) Bubble plot
# (1) show all the significant interactions (L-R pairs) from some cell groups (defined by 'sources.use') to other cell groups (defined by 'targets.use')
netVisual_bubble(cellchat, sources.use = 4, targets.use = c(5:11), remove.isolate = FALSE)Comparing communications on a single object 
#> Comparing communications on a single object
# (2) show all the significant interactions (L-R pairs) associated with certain signaling pathways
netVisual_bubble(cellchat, sources.use = 4, targets.use = c(5:11), signaling = c("CCL","CXCL"), remove.isolate = FALSE)Comparing communications on a single object 
#> Comparing communications on a single object
# (3) show all the significant interactions (L-R pairs) based on user's input (defined by `pairLR.use`)
pairLR.use <- extractEnrichedLR(cellchat, signaling = c("CCL","CXCL","FGF"))There is no significant communication of CXCL
There is no significant communication of FGFnetVisual_bubble(cellchat, sources.use = c(3,4), targets.use = c(5:8), pairLR.use = pairLR.use, remove.isolate = TRUE)Comparing communications on a single object 
#> Comparing communications on a single object
# (B) Chord diagram
# show all the significant interactions (L-R pairs) from some cell groups (defined by 'sources.use') to other cell groups (defined by 'targets.use')
# show all the interactions sending from Inflam.FIB
netVisual_chord_gene(cellchat, sources.use = 4, targets.use = c(5:11), lab.cex = 0.5,legend.pos.y = 30)
# show all the interactions received by Inflam.DC
netVisual_chord_gene(cellchat, sources.use = c(1,2,3,4), targets.use = 8, legend.pos.x = 15)
# show all the significant interactions (L-R pairs) associated with certain signaling pathways
netVisual_chord_gene(cellchat, sources.use = c(1,2,3,4), targets.use = c(5:11), signaling = c("CCL","CXCL"),legend.pos.x = 8)
# show all the significant signaling pathways from some cell groups (defined by 'sources.use') to other cell groups (defined by 'targets.use')
netVisual_chord_gene(cellchat, sources.use = c(1,2,3,4), targets.use = c(5:11), slot.name = "netP", legend.pos.x = 10)
#----------------------------------------------------------------------------------------------------
# **Plot the signaling gene expression distribution using violin/dot plot**
#----------------------------------------------------------------------------------------------------
plotGeneExpression(cellchat, signaling = "CXCL", enriched.only = TRUE, type = "violin")There is no significant communication of CXCL
Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.
plotGeneExpression(cellchat, signaling = "CXCL", enriched.only = FALSE)Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.Warning: Default search for "data" layer in "RNA" assay yielded no results; utilizing "counts" layer instead.Scale for y is already present.
Adding another scale for y, which will replace the existing scale.
6. Systems analysis of cell-cell communication network
#----------------------------------------------------------------------------------------------------
# **Part IV: Systems analysis of cell-cell communication network**
#----------------------------------------------------------------------------------------------------
# (A) Compute and visualize the network centrality scores
ptm = Sys.time()
# Compute the network centrality scores
cellchat <- netAnalysis_computeCentrality(cellchat, slot.name = "netP") # the slot 'netP' means the inferred intercellular communication network of signaling pathwaysWarning: UNRELIABLE VALUE: One of the ‘future.apply’ iterations (‘future_sapply-1’) unexpectedly generated random numbers without declaring so. There is a risk that those random numbers are not statistically sound and the overall results might be invalid. To fix this, specify 'future.seed=TRUE'. This ensures that proper, parallel-safe random numbers are produced via the L'Ecuyer-CMRG method. To disable this check, use 'future.seed = NULL', or set option 'future.rng.onMisuse' to "ignore".Warning: UNRELIABLE VALUE: One of the ‘future.apply’ iterations (‘future_sapply-2’) unexpectedly generated random numbers without declaring so. There is a risk that those random numbers are not statistically sound and the overall results might be invalid. To fix this, specify 'future.seed=TRUE'. This ensures that proper, parallel-safe random numbers are produced via the L'Ecuyer-CMRG method. To disable this check, use 'future.seed = NULL', or set option 'future.rng.onMisuse' to "ignore".Warning: UNRELIABLE VALUE: One of the ‘future.apply’ iterations (‘future_sapply-3’) unexpectedly generated random numbers without declaring so. There is a risk that those random numbers are not statistically sound and the overall results might be invalid. To fix this, specify 'future.seed=TRUE'. This ensures that proper, parallel-safe random numbers are produced via the L'Ecuyer-CMRG method. To disable this check, use 'future.seed = NULL', or set option 'future.rng.onMisuse' to "ignore".Warning: UNRELIABLE VALUE: One of the ‘future.apply’ iterations (‘future_sapply-4’) unexpectedly generated random numbers without declaring so. There is a risk that those random numbers are not statistically sound and the overall results might be invalid. To fix this, specify 'future.seed=TRUE'. This ensures that proper, parallel-safe random numbers are produced via the L'Ecuyer-CMRG method. To disable this check, use 'future.seed = NULL', or set option 'future.rng.onMisuse' to "ignore".# Visualize the computed centrality scores using heatmap, allowing ready identification of major signaling roles of cell groups
netAnalysis_signalingRole_network(cellchat, signaling = pathways.show, width = 8, height = 2.5, font.size = 10)
#(B) Visualize dominant senders (sources) and receivers (targets) in a 2D space
# Signaling role analysis on the aggregated cell-cell communication network from all signaling pathways
gg1 <- netAnalysis_signalingRole_scatter(cellchat)Signaling role analysis on the aggregated cell-cell communication network from all signaling pathways#> Signaling role analysis on the aggregated cell-cell communication network from all signaling pathways
# Signaling role analysis on the cell-cell communication networks of interest
gg2 <- netAnalysis_signalingRole_scatter(cellchat, signaling = c("CXCL", "CCL"))Signaling role analysis on the cell-cell communication network from user's input#> Signaling role analysis on the cell-cell communication network from user's input
gg1 + gg2
# (C) Identify signals contributing the most to outgoing or incoming signaling of certain cell groups
# Signaling role analysis on the aggregated cell-cell communication network from all signaling pathways
ht1 <- netAnalysis_signalingRole_heatmap(cellchat, pattern = "outgoing")
ht2 <- netAnalysis_signalingRole_heatmap(cellchat, pattern = "incoming")
ht1 + ht2Warning: Heatmap/annotation names are duplicated: Relative strength
# Signaling role analysis on the cell-cell communication networks of interest
ht <- netAnalysis_signalingRole_heatmap(cellchat, signaling = c("CXCL", "CCL"))
#---------------------------------------------------------------------------------------------------------------------------
# **Identify global communication patterns to explore how multiple cell types and signaling pathways coordinate together**
#----------------------------------------------------------------------------------------------------------------------------
# (A) Identify and visualize outgoing communication pattern of secreting cells
library(NMF)
#> Loading required package: registry
#> Loading required package: rngtools
#> Loading required package: cluster
#> NMF - BioConductor layer [OK] | Shared memory capabilities [NO: bigmemory] | Cores 2/2
#> To enable shared memory capabilities, try: install.extras('
#> NMF
#> ')
#>
#> Attaching package: 'NMF'
#> The following objects are masked from 'package:igraph':
#>
#> algorithm, compare
library(ggalluvial)
selectK(cellchat, pattern = "outgoing")Warning: Continuous limits supplied to discrete scale.
ℹ Did you mean `limits = factor(...)` or `scale_*_continuous()`?gc() used (Mb) gc trigger (Mb) max used (Mb)
Ncells 8931046 477.0 17088110 912.7 17088110 912.7
Vcells 2255385642 17207.3 3872088100 29541.7 3480046454 26550.7gc() used (Mb) gc trigger (Mb) max used (Mb)
Ncells 8930921 477.0 17088110 912.7 17088110 912.7
Vcells 2255385451 17207.3 3872088100 29541.7 3480046454 26550.7gc() used (Mb) gc trigger (Mb) max used (Mb)
Ncells 8930921 477.0 17088110 912.7 17088110 912.7
Vcells 2255385451 17207.3 3872088100 29541.7 3480046454 26550.7gc() used (Mb) gc trigger (Mb) max used (Mb)
Ncells 8930921 477.0 17088110 912.7 17088110 912.7
Vcells 2255385451 17207.3 3872088100 29541.7 3480046454 26550.7gc() used (Mb) gc trigger (Mb) max used (Mb)
Ncells 8930921 477.0 17088110 912.7 17088110 912.7
Vcells 2255385451 17207.3 3872088100 29541.7 3480046454 26550.7gc() used (Mb) gc trigger (Mb) max used (Mb)
Ncells 8930921 477.0 17088110 912.7 17088110 912.7
Vcells 2255385451 17207.3 3872088100 29541.7 3480046454 26550.7nPatterns = 6
cellchat <- identifyCommunicationPatterns(cellchat, pattern = "outgoing", k = nPatterns)Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
# river plot
netAnalysis_river(cellchat, pattern = "outgoing")Please make sure you have load `library(ggalluvial)` when running this function
#> Please make sure you have load `library(ggalluvial)` when running this function
# dot plot
netAnalysis_dot(cellchat, pattern = "outgoing")
setwd(".")
# (B) Identify and visualize incoming communication pattern of target cells
selectK(cellchat, pattern = "incoming")Warning: Continuous limits supplied to discrete scale.
ℹ Did you mean `limits = factor(...)` or `scale_*_continuous()`? nPatterns = 3
cellchat <- identifyCommunicationPatterns(cellchat, pattern = "incoming", k = nPatterns)Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
Warning: Recycling array of length 1 in array-vector arithmetic is deprecated.
Use c() or as.vector() instead.
# river plot
netAnalysis_river(cellchat, pattern = "incoming")Please make sure you have load `library(ggalluvial)` when running this function
#> Please make sure you have load `library(ggalluvial)` when running this function
# dot plot
netAnalysis_dot(cellchat, pattern = "incoming")
#---------------------------------------------------------------------------------------------------------------------------
# **Manifold and classification learning analysis of signaling networks**
#----------------------------------------------------------------------------------------------------------------------------
# Identify signaling groups based on their functional similarity
#8. Identify signaling roles and patterns
cellchat <- computeNetSimilarity(cellchat, type = "functional")
library(uwot)
cellchat <- netEmbedding(cellchat, type = "functional", umap.method = "uwot")Manifold learning of the signaling networks for a single dataset options(future.globals.maxSize = 16000 * 1024^2) # Set to 16 GB
cellchat <- netClustering(cellchat, type = "functional", do.parallel = FALSE)Classification learning of the signaling networks for a single dataset Warning: UNRELIABLE VALUE: One of the ‘future.apply’ iterations (‘future_sapply-1’) unexpectedly generated random numbers without declaring so. There is a risk that those random numbers are not statistically sound and the overall results might be invalid. To fix this, specify 'future.seed=TRUE'. This ensures that proper, parallel-safe random numbers are produced via the L'Ecuyer-CMRG method. To disable this check, use 'future.seed = NULL', or set option 'future.rng.onMisuse' to "ignore".Warning: UNRELIABLE VALUE: One of the ‘future.apply’ iterations (‘future_sapply-2’) unexpectedly generated random numbers without declaring so. There is a risk that those random numbers are not statistically sound and the overall results might be invalid. To fix this, specify 'future.seed=TRUE'. This ensures that proper, parallel-safe random numbers are produced via the L'Ecuyer-CMRG method. To disable this check, use 'future.seed = NULL', or set option 'future.rng.onMisuse' to "ignore".Warning: UNRELIABLE VALUE: One of the ‘future.apply’ iterations (‘future_sapply-3’) unexpectedly generated random numbers without declaring so. There is a risk that those random numbers are not statistically sound and the overall results might be invalid. To fix this, specify 'future.seed=TRUE'. This ensures that proper, parallel-safe random numbers are produced via the L'Ecuyer-CMRG method. To disable this check, use 'future.seed = NULL', or set option 'future.rng.onMisuse' to "ignore".Warning: UNRELIABLE VALUE: One of the ‘future.apply’ iterations (‘future_sapply-4’) unexpectedly generated random numbers without declaring so. There is a risk that those random numbers are not statistically sound and the overall results might be invalid. To fix this, specify 'future.seed=TRUE'. This ensures that proper, parallel-safe random numbers are produced via the L'Ecuyer-CMRG method. To disable this check, use 'future.seed = NULL', or set option 'future.rng.onMisuse' to "ignore".Warning: all eigenvalues are requested, eigen() is used insteadWarning: all eigenvalues are requested, eigen() is used insteadnetVisual_embedding(cellchat, type = "functional", label.size = 3.5)Warning: Using alpha for a discrete variable is not advised.
# Identify signaling groups based on structure similarity
#install umap-learn before running these commands
#pip install umap-learn
cellchat <- computeNetSimilarity(cellchat, type = "structural")
cellchat <- netEmbedding(cellchat, type = "structural", umap.method = "uwot")Manifold learning of the signaling networks for a single dataset #> Manifold learning of the signaling networks for a single dataset
options(future.globals.maxSize = 16000 * 1024^2) # Set to 16 GB
cellchat <- netClustering(cellchat, type = "structural", do.parallel = FALSE)Classification learning of the signaling networks for a single dataset Warning: UNRELIABLE VALUE: One of the ‘future.apply’ iterations (‘future_sapply-1’) unexpectedly generated random numbers without declaring so. There is a risk that those random numbers are not statistically sound and the overall results might be invalid. To fix this, specify 'future.seed=TRUE'. This ensures that proper, parallel-safe random numbers are produced via the L'Ecuyer-CMRG method. To disable this check, use 'future.seed = NULL', or set option 'future.rng.onMisuse' to "ignore".Warning: UNRELIABLE VALUE: One of the ‘future.apply’ iterations (‘future_sapply-2’) unexpectedly generated random numbers without declaring so. There is a risk that those random numbers are not statistically sound and the overall results might be invalid. To fix this, specify 'future.seed=TRUE'. This ensures that proper, parallel-safe random numbers are produced via the L'Ecuyer-CMRG method. To disable this check, use 'future.seed = NULL', or set option 'future.rng.onMisuse' to "ignore".Warning: UNRELIABLE VALUE: One of the ‘future.apply’ iterations (‘future_sapply-3’) unexpectedly generated random numbers without declaring so. There is a risk that those random numbers are not statistically sound and the overall results might be invalid. To fix this, specify 'future.seed=TRUE'. This ensures that proper, parallel-safe random numbers are produced via the L'Ecuyer-CMRG method. To disable this check, use 'future.seed = NULL', or set option 'future.rng.onMisuse' to "ignore".Warning: UNRELIABLE VALUE: One of the ‘future.apply’ iterations (‘future_sapply-4’) unexpectedly generated random numbers without declaring so. There is a risk that those random numbers are not statistically sound and the overall results might be invalid. To fix this, specify 'future.seed=TRUE'. This ensures that proper, parallel-safe random numbers are produced via the L'Ecuyer-CMRG method. To disable this check, use 'future.seed = NULL', or set option 'future.rng.onMisuse' to "ignore".Warning: all eigenvalues are requested, eigen() is used insteadWarning: all eigenvalues are requested, eigen() is used instead#> Classification learning of the signaling networks for a single dataset
# Visualization in 2D-space
netVisual_embedding(cellchat, type = "structural", label.size = 3.5)Warning: Using alpha for a discrete variable is not advised.
netVisual_embeddingZoomIn(cellchat, type = "structural", nCol = 2)
7. Save the CellChat object
save(cellchat, file = "cellCHAT_Analysis_Melange.Robj")LS0tCnRpdGxlOiAiSW5mZXJlbmNlIGFuZCBhbmFseXNpcyBvZiBjZWxsLWNlbGwgY29tbXVuaWNhdGlvbiB1c2luZyBDZWxsQ2hhdC1NZWxhbmdlIgphdXRob3I6IE5hc2lyIE1haG1vb2QgQWJiYXNpCmRhdGU6ICIyMDIzLTA5LTEyIgpvdXRwdXQ6CiAgaHRtbF9ub3RlYm9vazogCiAgICB0b2M6IHRydWUKICAgIHRvY19mbG9hdDogdHJ1ZQogICAgdG9jX2NvbGxhcHNlZDogdHJ1ZQogICAgdGhlbWU6IGRhcmtseQotLS0KCiMgMS4gbG9hZCBsaWJyYXJpZXMKYGBge3Igc2V0dXAsIGluY2x1ZGU9RkFMU0V9CgpsaWJyYXJ5KFNldXJhdCkKbGlicmFyeShkcGx5cikKbGlicmFyeShDZWxsQ2hhdCkKbGlicmFyeShwYXRjaHdvcmspCm9wdGlvbnMoc3RyaW5nc0FzRmFjdG9ycyA9IEZBTFNFKQojIHJldGljdWxhdGU6OnVzZV9weXRob24oIi9Vc2Vycy9zdW9xaW5qaW4vYW5hY29uZGEzL2Jpbi9weXRob24iLCByZXF1aXJlZD1UKSAKCmBgYAoKIyAyLiBMb2FkIFNldXJhdCBPYmplY3QKYGBge3IgbG9hZF9zZXVyYXRPQkp9CgojTG9hZCBTZXVyYXQgT2JqZWN0IEludGVncmF0ZWQgT2JqZWN0CmxvYWQoIk1lbGFuZ2VfdjMuUm9iaiIpCgpNZWxhbmdlCgoKRGVmYXVsdEFzc2F5KE1lbGFuZ2UpIDwtICJTQ1QiCgpgYGAKCgojIDMuIERhdGEgaW5wdXQgJiBwcm9jZXNzaW5nIGFuZCBpbml0aWFsaXphdGlvbiBvZiBDZWxsQ2hhdCBvYmplY3QKYGBge3IgREFUQUlOUFVULCBmaWcuaGVpZ2h0PTYsIGZpZy53aWR0aD0xMH0KCiMtLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KIyAqKlBhcnQgSTogRGF0YSBpbnB1dCAmIHByb2Nlc3NpbmcgYW5kIGluaXRpYWxpemF0aW9uIG9mIENlbGxDaGF0IG9iamVjdCoqCiMtLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCiMgTG9hZGluZyBhbmQgUHJlcGFyaW5nIERhdGEKTWVsYW5nZSA8LSBBZGRNZXRhRGF0YShvYmplY3QgPSBNZWxhbmdlLCAKICAgICAgICAgICAgICAgICAgICAgICAgICBtZXRhZGF0YSA9IHBhc3RlMCgiQ2x1c3Rlcl8iLCBNZWxhbmdlJHNldXJhdF9jbHVzdGVycyksIAogICAgICAgICAgICAgICAgICAgICAgICAgIGNvbC5uYW1lID0gJ0NsdXN0ZXJzJykKCiNNZWxhbmdlW1siUk5BIl1dIDwtIEpvaW5MYXllcnMoTWVsYW5nZVtbIlJOQSJdXSkKCgojIEV4dHJhY3QgdGhlIG5vcm1hbGl6ZWQgZGF0YQojZGF0YV9pbnB1dCA8LSBHZXRBc3NheURhdGEoTWVsYW5nZSwgYXNzYXkgPSAiUk5BIiwgc2xvdCA9ICJkYXRhIikKCgpkYXRhLmlucHV0IDwtIE1lbGFuZ2VAYXNzYXlzJFNDVCRkYXRhCgoKIzEuICBDcmVhdGUgQ2VsbENoYXQgb2JqZWN0CmNlbGxjaGF0IDwtIGNyZWF0ZUNlbGxDaGF0KG9iamVjdCA9IGRhdGEuaW5wdXQsIG1ldGEgPSBNZWxhbmdlQG1ldGEuZGF0YSwgZ3JvdXAuYnkgPSAiQ2x1c3RlcnMiKQoKCgoKIzIuICBTZXQgdGhlIGxpZ2FuZC1yZWNlcHRvciBpbnRlcmFjdGlvbiBkYXRhYmFzZQoKY2VsbGNoYXREQiA8LSBDZWxsQ2hhdERCLmh1bWFuICMgb3IgQ2VsbENoYXREQi5tb3VzZQpzaG93RGF0YWJhc2VDYXRlZ29yeShjZWxsY2hhdERCKQoKIyBTaG93IHRoZSBzdHJ1Y3R1cmUgb2YgdGhlIGRhdGFiYXNlCmRwbHlyOjpnbGltcHNlKGNlbGxjaGF0REIkaW50ZXJhY3Rpb24pCgojIHVzZSBhIHN1YnNldCBvZiBDZWxsQ2hhdERCIGZvciBjZWxsLWNlbGwgY29tbXVuaWNhdGlvbiBhbmFseXNpcwpjZWxsY2hhdERCLnVzZSA8LSBzdWJzZXREQihjZWxsY2hhdERCLCBzZWFyY2ggPSAiU2VjcmV0ZWQgU2lnbmFsaW5nIiwga2V5ID0gImFubm90YXRpb24iKSAjIHVzZSBTZWNyZXRlZCBTaWduYWxpbmcKCiMgc2V0IHRoZSB1c2VkIGRhdGFiYXNlIGluIHRoZSBvYmplY3QKY2VsbGNoYXRAREIgPC0gY2VsbGNoYXREQi51c2UKCiMtLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KIyAqKlByZXByb2Nlc3NpbmcgdGhlIGV4cHJlc3Npb24gZGF0YSBmb3IgY2VsbC1jZWxsIGNvbW11bmljYXRpb24gYW5hbHlzaXMqKgojLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCiMzLiAgUHJlcHJvY2VzcyB0aGUgZXhwcmVzc2lvbiBkYXRhCmNlbGxjaGF0IDwtIHN1YnNldERhdGEoY2VsbGNoYXQpCmNlbGxjaGF0IDwtIGlkZW50aWZ5T3ZlckV4cHJlc3NlZEdlbmVzKGNlbGxjaGF0KQpjZWxsY2hhdCA8LSBpZGVudGlmeU92ZXJFeHByZXNzZWRJbnRlcmFjdGlvbnMoY2VsbGNoYXQpCgoKCgpgYGAKCiMgNC4gSW5mZXJlbmNlIG9mIGNlbGwtY2VsbCBjb21tdW5pY2F0aW9uIG5ldHdvcmsKYGBge3IgSU5GRVJFTkNFLCBmaWcuaGVpZ2h0PTYsIGZpZy53aWR0aD0xMH0KIy0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQojICoqUGFydCBJSTogSW5mZXJlbmNlIG9mIGNlbGwtY2VsbCBjb21tdW5pY2F0aW9uIG5ldHdvcmsqKgojLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgojLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCiMgKipDb21wdXRlIHRoZSBjb21tdW5pY2F0aW9uIHByb2JhYmlsaXR5IGFuZCBpbmZlciBjZWxsdWxhciBjb21tdW5pY2F0aW9uIG5ldHdvcmsqKgojLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgojNC4gIENvbXB1dGUgY29tbXVuaWNhdGlvbiBwcm9iYWJpbGl0eSBhbmQgaW5mZXIgY2VsbHVsYXIgY29tbXVuaWNhdGlvbiBuZXR3b3JrCmNlbGxjaGF0IDwtIGNvbXB1dGVDb21tdW5Qcm9iKGNlbGxjaGF0KQpjZWxsY2hhdCA8LSBmaWx0ZXJDb21tdW5pY2F0aW9uKGNlbGxjaGF0LCBtaW4uY2VsbHMgPSAxMCkKCiMtLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KIyAqKkluZmVyIHRoZSBjZWxsLWNlbGwgY29tbXVuaWNhdGlvbiBhdCBhIHNpZ25hbGluZyBwYXRod2F5IGxldmVsKioKIy0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQojNS4gIENvbXB1dGUgY29tbXVuaWNhdGlvbiBhdCB0aGUgcGF0aHdheSBsZXZlbApjZWxsY2hhdCA8LSBjb21wdXRlQ29tbXVuUHJvYlBhdGh3YXkoY2VsbGNoYXQpCgojLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCiMgKipDYWxjdWxhdGUgdGhlIGFnZ3JlZ2F0ZWQgY2VsbC1jZWxsIGNvbW11bmljYXRpb24gbmV0d29yayoqCiMtLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KIzYuICBBZ2dyZWdhdGUgdGhlIGluZmVycmVkIGNlbGwtY2VsbCBjb21tdW5pY2F0aW9uIG5ldHdvcmsKY2VsbGNoYXQgPC0gYWdncmVnYXRlTmV0KGNlbGxjaGF0KQoKCiM3LiAgVmlzdWFsaXplIHRoZSBjb21tdW5pY2F0aW9uIG5ldHdvcmsKZ3JvdXBTaXplIDwtIGFzLm51bWVyaWModGFibGUoY2VsbGNoYXRAaWRlbnRzKSkKCm5ldFZpc3VhbF9jaXJjbGUoY2VsbGNoYXRAbmV0JGNvdW50LCB2ZXJ0ZXgud2VpZ2h0ID0gZ3JvdXBTaXplLCB3ZWlnaHQuc2NhbGUgPSBULCBsYWJlbC5lZGdlPSBGLCB0aXRsZS5uYW1lID0gIk51bWJlciBvZiBpbnRlcmFjdGlvbnMiKQpuZXRWaXN1YWxfY2lyY2xlKGNlbGxjaGF0QG5ldCR3ZWlnaHQsIHZlcnRleC53ZWlnaHQgPSBncm91cFNpemUsIHdlaWdodC5zY2FsZSA9IFQsIGxhYmVsLmVkZ2U9IEYsIHRpdGxlLm5hbWUgPSAiSW50ZXJhY3Rpb24gd2VpZ2h0cy9zdHJlbmd0aCIpCgptYXQgPC0gY2VsbGNoYXRAbmV0JHdlaWdodApwYXIobWZyb3cgPSBjKDMsNCksIHhwZD1UUlVFKQpmb3IgKGkgaW4gMTpucm93KG1hdCkpIHsKICBtYXQyIDwtIG1hdHJpeCgwLCBucm93ID0gbnJvdyhtYXQpLCBuY29sID0gbmNvbChtYXQpLCBkaW1uYW1lcyA9IGRpbW5hbWVzKG1hdCkpCiAgbWF0MltpLCBdIDwtIG1hdFtpLCBdCiAgbmV0VmlzdWFsX2NpcmNsZShtYXQyLCB2ZXJ0ZXgud2VpZ2h0ID0gZ3JvdXBTaXplLCB3ZWlnaHQuc2NhbGUgPSBULCBlZGdlLndlaWdodC5tYXggPSBtYXgobWF0KSwgdGl0bGUubmFtZSA9IHJvd25hbWVzKG1hdClbaV0pCn0KCmBgYAojIDUuIFZpc3VhbGl6YXRpb24gb2YgY2VsbC1jZWxsIGNvbW11bmljYXRpb24gbmV0d29yawpgYGB7ciBWSVNVQUxJWkFUSU9OLCBmaWcuaGVpZ2h0PTYsIGZpZy53aWR0aD0xMH0KIy0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQojICoqUGFydCBJSUk6IFZpc3VhbGl6YXRpb24gb2YgY2VsbC1jZWxsIGNvbW11bmljYXRpb24gbmV0d29yayoqCiMtLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KIAojLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQojICoqVmlzdWFsaXplIGVhY2ggc2lnbmFsaW5nIHBhdGh3YXkgdXNpbmcgSGllcmFyY2h5IHBsb3QsIENpcmNsZSBwbG90IG9yIENob3JkIGRpYWdyYW0qKgojLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQogIApwYXRod2F5cy5zaG93IDwtIGMoIklMMiIpIAojIEhpZXJhcmNoeSBwbG90CiMgSGVyZSB3ZSBkZWZpbmUgYHZlcnRleC5yZWNlaXZlYCBzbyB0aGF0IHRoZSBsZWZ0IHBvcnRpb24gb2YgdGhlIGhpZXJhcmNoeSBwbG90IHNob3dzIHNpZ25hbGluZyB0byBmaWJyb2JsYXN0IGFuZCB0aGUgcmlnaHQgcG9ydGlvbiBzaG93cyBzaWduYWxpbmcgdG8gaW1tdW5lIGNlbGxzIAp2ZXJ0ZXgucmVjZWl2ZXIgPSBzZXEoMSw0KSAjIGEgbnVtZXJpYyB2ZWN0b3IuIApuZXRWaXN1YWxfYWdncmVnYXRlKGNlbGxjaGF0LCBzaWduYWxpbmcgPSBwYXRod2F5cy5zaG93LCAgdmVydGV4LnJlY2VpdmVyID0gdmVydGV4LnJlY2VpdmVyKQojIENpcmNsZSBwbG90CnBhcihtZnJvdz1jKDEsMSkpCm5ldFZpc3VhbF9hZ2dyZWdhdGUoY2VsbGNoYXQsIHNpZ25hbGluZyA9IHBhdGh3YXlzLnNob3csIGxheW91dCA9ICJjaXJjbGUiKQoKCiMgQ2hvcmQgZGlhZ3JhbQpwYXIobWZyb3c9YygxLDEpKQpuZXRWaXN1YWxfYWdncmVnYXRlKGNlbGxjaGF0LCBzaWduYWxpbmcgPSBwYXRod2F5cy5zaG93LCBsYXlvdXQgPSAiY2hvcmQiKQoKIyBIZWF0bWFwCnBhcihtZnJvdz1jKDEsMSkpCm5ldFZpc3VhbF9oZWF0bWFwKGNlbGxjaGF0LCBzaWduYWxpbmcgPSBwYXRod2F5cy5zaG93LCBjb2xvci5oZWF0bWFwID0gIlJlZHMiKQojPiBEbyBoZWF0bWFwIGJhc2VkIG9uIGEgc2luZ2xlIG9iamVjdAoKCgojIEdldCB0aGUgbGV2ZWxzIG9mIGNlbGwgdHlwZXMKY2VsbF90eXBlcyA8LSBsZXZlbHMoY2VsbGNoYXRAaWRlbnRzKQoKIyBDaG9yZCBkaWFncmFtCmdyb3VwLmNlbGxUeXBlIDwtIGMocmVwKCJNSyIsIDQpLCByZXAoIklMMiIsIDQpLCByZXAoIlRORiIsIDQpKSAKCmdyb3VwLmNlbGxUeXBlIDwtIHJlcChOQSwgbGVuZ3RoKGNlbGxfdHlwZXMpKQoKIyBBc3NpZ24gZWFjaCBjZWxsIHR5cGUgdG8gYSBncm91cAojIFlvdSdsbCBuZWVkIHRvIGFkanVzdCB0aGlzIGJhc2VkIG9uIHlvdXIgYWN0dWFsIGRhdGEgYW5kIGJpb2xvZ2ljYWwga25vd2xlZGdlCmdyb3VwLmNlbGxUeXBlWzE6OF0gPC0gIk1LIgpncm91cC5jZWxsVHlwZVs5OjE2XSA8LSAiSUwyIgpncm91cC5jZWxsVHlwZVsxNzozMF0gPC0gIlRORiIKCiMgTmFtZSB0aGUgdmVjdG9yIHdpdGggdGhlIGNlbGwgdHlwZXMKbmFtZXMoZ3JvdXAuY2VsbFR5cGUpIDwtIGNlbGxfdHlwZXMKCiMgQ2hlY2sgdGhlIHJlc3VsdApwcmludChncm91cC5jZWxsVHlwZSkKCm5ldFZpc3VhbF9jaG9yZF9jZWxsKGNlbGxjaGF0LCBzaWduYWxpbmcgPSBwYXRod2F5cy5zaG93LCBncm91cCA9IGdyb3VwLmNlbGxUeXBlLCB0aXRsZS5uYW1lID0gcGFzdGUwKHBhdGh3YXlzLnNob3csICIgc2lnbmFsaW5nIG5ldHdvcmsiKSkKIz4gUGxvdCB0aGUgYWdncmVnYXRlZCBjZWxsLWNlbGwgY29tbXVuaWNhdGlvbiBuZXR3b3JrIGF0IHRoZSBzaWduYWxpbmcgcGF0aHdheSBsZXZlbAoKCiAgCiMgQ29tcHV0ZSB0aGUgY29udHJpYnV0aW9uIG9mIGVhY2ggbGlnYW5kLXJlY2VwdG9yIHBhaXIgdG8gdGhlIG92ZXJhbGwgc2lnbmFsaW5nIHBhdGh3YXkgYW5kIHZpc3VhbGl6ZSBjZWxsLWNlbGwgCiNjb21tdW5pY2F0aW9uIG1lZGlhdGVkIGJ5IGEgc2luZ2xlIGxpZ2FuZC1yZWNlcHRvciBwYWlyCgpuZXRBbmFseXNpc19jb250cmlidXRpb24oY2VsbGNoYXQsIHNpZ25hbGluZyA9IHBhdGh3YXlzLnNob3cpCgpwYWlyTFIuQ1hDTCA8LSBleHRyYWN0RW5yaWNoZWRMUihjZWxsY2hhdCwgc2lnbmFsaW5nID0gcGF0aHdheXMuc2hvdywgZ2VuZUxSLnJldHVybiA9IEZBTFNFKQpMUi5zaG93IDwtIHBhaXJMUi5DWENMWzEsXSAjIHNob3cgb25lIGxpZ2FuZC1yZWNlcHRvciBwYWlyCiMgSGllcmFyY2h5IHBsb3QKdmVydGV4LnJlY2VpdmVyID0gc2VxKDEsNCkgIyBhIG51bWVyaWMgdmVjdG9yCm5ldFZpc3VhbF9pbmRpdmlkdWFsKGNlbGxjaGF0LCBzaWduYWxpbmcgPSBwYXRod2F5cy5zaG93LCAgcGFpckxSLnVzZSA9IExSLnNob3csIHZlcnRleC5yZWNlaXZlciA9IHZlcnRleC5yZWNlaXZlcikKIz4gW1sxXV0KIyBDaXJjbGUgcGxvdApuZXRWaXN1YWxfaW5kaXZpZHVhbChjZWxsY2hhdCwgc2lnbmFsaW5nID0gcGF0aHdheXMuc2hvdywgcGFpckxSLnVzZSA9IExSLnNob3csIGxheW91dCA9ICJjaXJjbGUiKQoKCgojPiBbWzFdXQojIENob3JkIGRpYWdyYW0KbmV0VmlzdWFsX2luZGl2aWR1YWwoY2VsbGNoYXQsIHNpZ25hbGluZyA9IHBhdGh3YXlzLnNob3csIHBhaXJMUi51c2UgPSBMUi5zaG93LCBsYXlvdXQgPSAiY2hvcmQiKQoKIyBBY2Nlc3MgYWxsIHRoZSBzaWduYWxpbmcgcGF0aHdheXMgc2hvd2luZyBzaWduaWZpY2FudCBjb21tdW5pY2F0aW9ucwpwYXRod2F5cy5zaG93LmFsbCA8LSBjZWxsY2hhdEBuZXRQJHBhdGh3YXlzCiMgY2hlY2sgdGhlIG9yZGVyIG9mIGNlbGwgaWRlbnRpdHkgdG8gc2V0IHN1aXRhYmxlIHZlcnRleC5yZWNlaXZlcgpsZXZlbHMoY2VsbGNoYXRAaWRlbnRzKQp2ZXJ0ZXgucmVjZWl2ZXIgPSBzZXEoMSw0KQpmb3IgKGkgaW4gMTpsZW5ndGgocGF0aHdheXMuc2hvdy5hbGwpKSB7CiAgIyBWaXN1YWxpemUgY29tbXVuaWNhdGlvbiBuZXR3b3JrIGFzc29jaWF0ZWQgd2l0aCBib3RoIHNpZ25hbGluZyBwYXRod2F5IGFuZCBpbmRpdmlkdWFsIEwtUiBwYWlycwogIG5ldFZpc3VhbChjZWxsY2hhdCwgc2lnbmFsaW5nID0gcGF0aHdheXMuc2hvdy5hbGxbaV0sIHZlcnRleC5yZWNlaXZlciA9IHZlcnRleC5yZWNlaXZlciwgbGF5b3V0ID0gImhpZXJhcmNoeSIpCiAgIyBDb21wdXRlIGFuZCB2aXN1YWxpemUgdGhlIGNvbnRyaWJ1dGlvbiBvZiBlYWNoIGxpZ2FuZC1yZWNlcHRvciBwYWlyIHRvIHRoZSBvdmVyYWxsIHNpZ25hbGluZyBwYXRod2F5CiAgZ2cgPC0gbmV0QW5hbHlzaXNfY29udHJpYnV0aW9uKGNlbGxjaGF0LCBzaWduYWxpbmcgPSBwYXRod2F5cy5zaG93LmFsbFtpXSkKICBnZ3NhdmUoZmlsZW5hbWU9cGFzdGUwKHBhdGh3YXlzLnNob3cuYWxsW2ldLCAiX0wtUl9jb250cmlidXRpb24ucGRmIiksIHBsb3Q9Z2csIHdpZHRoID0gMywgaGVpZ2h0ID0gMiwgdW5pdHMgPSAnaW4nLCBkcGkgPSAzMDApCn0KCiMtLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCiMgKipWaXN1YWxpemUgY2VsbC1jZWxsIGNvbW11bmljYXRpb24gbWVkaWF0ZWQgYnkgbXVsdGlwbGUgbGlnYW5kLXJlY2VwdG9ycyBvciBzaWduYWxpbmcgcGF0aHdheXMqKgojLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQojKEEpIEJ1YmJsZSBwbG90CgojICgxKSBzaG93IGFsbCB0aGUgc2lnbmlmaWNhbnQgaW50ZXJhY3Rpb25zIChMLVIgcGFpcnMpIGZyb20gc29tZSBjZWxsIGdyb3VwcyAoZGVmaW5lZCBieSAnc291cmNlcy51c2UnKSB0byBvdGhlciBjZWxsIGdyb3VwcyAoZGVmaW5lZCBieSAndGFyZ2V0cy51c2UnKQpuZXRWaXN1YWxfYnViYmxlKGNlbGxjaGF0LCBzb3VyY2VzLnVzZSA9IDQsIHRhcmdldHMudXNlID0gYyg1OjExKSwgcmVtb3ZlLmlzb2xhdGUgPSBGQUxTRSkKIz4gQ29tcGFyaW5nIGNvbW11bmljYXRpb25zIG9uIGEgc2luZ2xlIG9iamVjdAoKCgojICgyKSBzaG93IGFsbCB0aGUgc2lnbmlmaWNhbnQgaW50ZXJhY3Rpb25zIChMLVIgcGFpcnMpIGFzc29jaWF0ZWQgd2l0aCBjZXJ0YWluIHNpZ25hbGluZyBwYXRod2F5cwpuZXRWaXN1YWxfYnViYmxlKGNlbGxjaGF0LCBzb3VyY2VzLnVzZSA9IDQsIHRhcmdldHMudXNlID0gYyg1OjExKSwgc2lnbmFsaW5nID0gYygiQ0NMIiwiQ1hDTCIpLCByZW1vdmUuaXNvbGF0ZSA9IEZBTFNFKQojPiBDb21wYXJpbmcgY29tbXVuaWNhdGlvbnMgb24gYSBzaW5nbGUgb2JqZWN0CgoKIyAoMykgc2hvdyBhbGwgdGhlIHNpZ25pZmljYW50IGludGVyYWN0aW9ucyAoTC1SIHBhaXJzKSBiYXNlZCBvbiB1c2VyJ3MgaW5wdXQgKGRlZmluZWQgYnkgYHBhaXJMUi51c2VgKQpwYWlyTFIudXNlIDwtIGV4dHJhY3RFbnJpY2hlZExSKGNlbGxjaGF0LCBzaWduYWxpbmcgPSBjKCJDQ0wiLCJDWENMIiwiRkdGIikpCm5ldFZpc3VhbF9idWJibGUoY2VsbGNoYXQsIHNvdXJjZXMudXNlID0gYygzLDQpLCB0YXJnZXRzLnVzZSA9IGMoNTo4KSwgcGFpckxSLnVzZSA9IHBhaXJMUi51c2UsIHJlbW92ZS5pc29sYXRlID0gVFJVRSkKIz4gQ29tcGFyaW5nIGNvbW11bmljYXRpb25zIG9uIGEgc2luZ2xlIG9iamVjdAoKCiMgKEIpIENob3JkIGRpYWdyYW0KCiMgc2hvdyBhbGwgdGhlIHNpZ25pZmljYW50IGludGVyYWN0aW9ucyAoTC1SIHBhaXJzKSBmcm9tIHNvbWUgY2VsbCBncm91cHMgKGRlZmluZWQgYnkgJ3NvdXJjZXMudXNlJykgdG8gb3RoZXIgY2VsbCBncm91cHMgKGRlZmluZWQgYnkgJ3RhcmdldHMudXNlJykKIyBzaG93IGFsbCB0aGUgaW50ZXJhY3Rpb25zIHNlbmRpbmcgZnJvbSBJbmZsYW0uRklCCm5ldFZpc3VhbF9jaG9yZF9nZW5lKGNlbGxjaGF0LCBzb3VyY2VzLnVzZSA9IDQsIHRhcmdldHMudXNlID0gYyg1OjExKSwgbGFiLmNleCA9IDAuNSxsZWdlbmQucG9zLnkgPSAzMCkKCgojIHNob3cgYWxsIHRoZSBpbnRlcmFjdGlvbnMgcmVjZWl2ZWQgYnkgSW5mbGFtLkRDCm5ldFZpc3VhbF9jaG9yZF9nZW5lKGNlbGxjaGF0LCBzb3VyY2VzLnVzZSA9IGMoMSwyLDMsNCksIHRhcmdldHMudXNlID0gOCwgbGVnZW5kLnBvcy54ID0gMTUpCgojIHNob3cgYWxsIHRoZSBzaWduaWZpY2FudCBpbnRlcmFjdGlvbnMgKEwtUiBwYWlycykgYXNzb2NpYXRlZCB3aXRoIGNlcnRhaW4gc2lnbmFsaW5nIHBhdGh3YXlzCm5ldFZpc3VhbF9jaG9yZF9nZW5lKGNlbGxjaGF0LCBzb3VyY2VzLnVzZSA9IGMoMSwyLDMsNCksIHRhcmdldHMudXNlID0gYyg1OjExKSwgc2lnbmFsaW5nID0gYygiQ0NMIiwiQ1hDTCIpLGxlZ2VuZC5wb3MueCA9IDgpCgojIHNob3cgYWxsIHRoZSBzaWduaWZpY2FudCBzaWduYWxpbmcgcGF0aHdheXMgZnJvbSBzb21lIGNlbGwgZ3JvdXBzIChkZWZpbmVkIGJ5ICdzb3VyY2VzLnVzZScpIHRvIG90aGVyIGNlbGwgZ3JvdXBzIChkZWZpbmVkIGJ5ICd0YXJnZXRzLnVzZScpCm5ldFZpc3VhbF9jaG9yZF9nZW5lKGNlbGxjaGF0LCBzb3VyY2VzLnVzZSA9IGMoMSwyLDMsNCksIHRhcmdldHMudXNlID0gYyg1OjExKSwgc2xvdC5uYW1lID0gIm5ldFAiLCBsZWdlbmQucG9zLnggPSAxMCkKCgojLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQojICoqUGxvdCB0aGUgc2lnbmFsaW5nIGdlbmUgZXhwcmVzc2lvbiBkaXN0cmlidXRpb24gdXNpbmcgdmlvbGluL2RvdCBwbG90KioKIy0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCnBsb3RHZW5lRXhwcmVzc2lvbihjZWxsY2hhdCwgc2lnbmFsaW5nID0gIkNYQ0wiLCBlbnJpY2hlZC5vbmx5ID0gVFJVRSwgdHlwZSA9ICJ2aW9saW4iKQoKCnBsb3RHZW5lRXhwcmVzc2lvbihjZWxsY2hhdCwgc2lnbmFsaW5nID0gIkNYQ0wiLCBlbnJpY2hlZC5vbmx5ID0gRkFMU0UpCgpgYGAKCiMgNi4gU3lzdGVtcyBhbmFseXNpcyBvZiBjZWxsLWNlbGwgY29tbXVuaWNhdGlvbiBuZXR3b3JrCmBgYHtyIFNZU1RFTVMtQU5BTFlTSVMsIGZpZy5oZWlnaHQ9NiwgZmlnLndpZHRoPTEwfQojLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQojICoqUGFydCBJVjogU3lzdGVtcyBhbmFseXNpcyBvZiBjZWxsLWNlbGwgY29tbXVuaWNhdGlvbiBuZXR3b3JrKioKIy0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KIyAoQSkgQ29tcHV0ZSBhbmQgdmlzdWFsaXplIHRoZSBuZXR3b3JrIGNlbnRyYWxpdHkgc2NvcmVzCgpwdG0gPSBTeXMudGltZSgpCiMgQ29tcHV0ZSB0aGUgbmV0d29yayBjZW50cmFsaXR5IHNjb3JlcwpjZWxsY2hhdCA8LSBuZXRBbmFseXNpc19jb21wdXRlQ2VudHJhbGl0eShjZWxsY2hhdCwgc2xvdC5uYW1lID0gIm5ldFAiKSAjIHRoZSBzbG90ICduZXRQJyBtZWFucyB0aGUgaW5mZXJyZWQgaW50ZXJjZWxsdWxhciBjb21tdW5pY2F0aW9uIG5ldHdvcmsgb2Ygc2lnbmFsaW5nIHBhdGh3YXlzCiMgVmlzdWFsaXplIHRoZSBjb21wdXRlZCBjZW50cmFsaXR5IHNjb3JlcyB1c2luZyBoZWF0bWFwLCBhbGxvd2luZyByZWFkeSBpZGVudGlmaWNhdGlvbiBvZiBtYWpvciBzaWduYWxpbmcgcm9sZXMgb2YgY2VsbCBncm91cHMKbmV0QW5hbHlzaXNfc2lnbmFsaW5nUm9sZV9uZXR3b3JrKGNlbGxjaGF0LCBzaWduYWxpbmcgPSBwYXRod2F5cy5zaG93LCB3aWR0aCA9IDgsIGhlaWdodCA9IDIuNSwgZm9udC5zaXplID0gMTApCgoKIyhCKSBWaXN1YWxpemUgZG9taW5hbnQgc2VuZGVycyAoc291cmNlcykgYW5kIHJlY2VpdmVycyAodGFyZ2V0cykgaW4gYSAyRCBzcGFjZQoKIyBTaWduYWxpbmcgcm9sZSBhbmFseXNpcyBvbiB0aGUgYWdncmVnYXRlZCBjZWxsLWNlbGwgY29tbXVuaWNhdGlvbiBuZXR3b3JrIGZyb20gYWxsIHNpZ25hbGluZyBwYXRod2F5cwpnZzEgPC0gbmV0QW5hbHlzaXNfc2lnbmFsaW5nUm9sZV9zY2F0dGVyKGNlbGxjaGF0KQojPiBTaWduYWxpbmcgcm9sZSBhbmFseXNpcyBvbiB0aGUgYWdncmVnYXRlZCBjZWxsLWNlbGwgY29tbXVuaWNhdGlvbiBuZXR3b3JrIGZyb20gYWxsIHNpZ25hbGluZyBwYXRod2F5cwojIFNpZ25hbGluZyByb2xlIGFuYWx5c2lzIG9uIHRoZSBjZWxsLWNlbGwgY29tbXVuaWNhdGlvbiBuZXR3b3JrcyBvZiBpbnRlcmVzdApnZzIgPC0gbmV0QW5hbHlzaXNfc2lnbmFsaW5nUm9sZV9zY2F0dGVyKGNlbGxjaGF0LCBzaWduYWxpbmcgPSBjKCJDWENMIiwgIkNDTCIpKQojPiBTaWduYWxpbmcgcm9sZSBhbmFseXNpcyBvbiB0aGUgY2VsbC1jZWxsIGNvbW11bmljYXRpb24gbmV0d29yayBmcm9tIHVzZXIncyBpbnB1dApnZzEgKyBnZzIKCgojIChDKSBJZGVudGlmeSBzaWduYWxzIGNvbnRyaWJ1dGluZyB0aGUgbW9zdCB0byBvdXRnb2luZyBvciBpbmNvbWluZyBzaWduYWxpbmcgb2YgY2VydGFpbiBjZWxsIGdyb3VwcwojIFNpZ25hbGluZyByb2xlIGFuYWx5c2lzIG9uIHRoZSBhZ2dyZWdhdGVkIGNlbGwtY2VsbCBjb21tdW5pY2F0aW9uIG5ldHdvcmsgZnJvbSBhbGwgc2lnbmFsaW5nIHBhdGh3YXlzCmh0MSA8LSBuZXRBbmFseXNpc19zaWduYWxpbmdSb2xlX2hlYXRtYXAoY2VsbGNoYXQsIHBhdHRlcm4gPSAib3V0Z29pbmciKQpodDIgPC0gbmV0QW5hbHlzaXNfc2lnbmFsaW5nUm9sZV9oZWF0bWFwKGNlbGxjaGF0LCBwYXR0ZXJuID0gImluY29taW5nIikKaHQxICsgaHQyCgojIFNpZ25hbGluZyByb2xlIGFuYWx5c2lzIG9uIHRoZSBjZWxsLWNlbGwgY29tbXVuaWNhdGlvbiBuZXR3b3JrcyBvZiBpbnRlcmVzdApodCA8LSBuZXRBbmFseXNpc19zaWduYWxpbmdSb2xlX2hlYXRtYXAoY2VsbGNoYXQsIHNpZ25hbGluZyA9IGMoIkNYQ0wiLCAiQ0NMIikpCgoKIy0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQojICoqSWRlbnRpZnkgZ2xvYmFsIGNvbW11bmljYXRpb24gcGF0dGVybnMgdG8gZXhwbG9yZSBob3cgbXVsdGlwbGUgY2VsbCB0eXBlcyBhbmQgc2lnbmFsaW5nIHBhdGh3YXlzIGNvb3JkaW5hdGUgdG9nZXRoZXIqKgojLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKIyAoQSkgSWRlbnRpZnkgYW5kIHZpc3VhbGl6ZSBvdXRnb2luZyBjb21tdW5pY2F0aW9uIHBhdHRlcm4gb2Ygc2VjcmV0aW5nIGNlbGxzCgpsaWJyYXJ5KE5NRikKIz4gTG9hZGluZyByZXF1aXJlZCBwYWNrYWdlOiByZWdpc3RyeQojPiBMb2FkaW5nIHJlcXVpcmVkIHBhY2thZ2U6IHJuZ3Rvb2xzCiM+IExvYWRpbmcgcmVxdWlyZWQgcGFja2FnZTogY2x1c3RlcgojPiBOTUYgLSBCaW9Db25kdWN0b3IgbGF5ZXIgW09LXSB8IFNoYXJlZCBtZW1vcnkgY2FwYWJpbGl0aWVzIFtOTzogYmlnbWVtb3J5XSB8IENvcmVzIDIvMgojPiAgIFRvIGVuYWJsZSBzaGFyZWQgbWVtb3J5IGNhcGFiaWxpdGllcywgdHJ5OiBpbnN0YWxsLmV4dHJhcygnCiM+IE5NRgojPiAnKQojPiAKIz4gQXR0YWNoaW5nIHBhY2thZ2U6ICdOTUYnCiM+IFRoZSBmb2xsb3dpbmcgb2JqZWN0cyBhcmUgbWFza2VkIGZyb20gJ3BhY2thZ2U6aWdyYXBoJzoKIz4gCiM+ICAgICBhbGdvcml0aG0sIGNvbXBhcmUKbGlicmFyeShnZ2FsbHV2aWFsKQoKc2VsZWN0SyhjZWxsY2hhdCwgcGF0dGVybiA9ICJvdXRnb2luZyIpCgoKblBhdHRlcm5zID0gNgoKY2VsbGNoYXQgPC0gaWRlbnRpZnlDb21tdW5pY2F0aW9uUGF0dGVybnMoY2VsbGNoYXQsIHBhdHRlcm4gPSAib3V0Z29pbmciLCBrID0gblBhdHRlcm5zKQoKCiMgcml2ZXIgcGxvdApuZXRBbmFseXNpc19yaXZlcihjZWxsY2hhdCwgcGF0dGVybiA9ICJvdXRnb2luZyIpCiM+IFBsZWFzZSBtYWtlIHN1cmUgeW91IGhhdmUgbG9hZCBgbGlicmFyeShnZ2FsbHV2aWFsKWAgd2hlbiBydW5uaW5nIHRoaXMgZnVuY3Rpb24KCgojIGRvdCBwbG90Cm5ldEFuYWx5c2lzX2RvdChjZWxsY2hhdCwgcGF0dGVybiA9ICJvdXRnb2luZyIpCgpzZXR3ZCgiLiIpCgojIChCKSBJZGVudGlmeSBhbmQgdmlzdWFsaXplIGluY29taW5nIGNvbW11bmljYXRpb24gcGF0dGVybiBvZiB0YXJnZXQgY2VsbHMKCnNlbGVjdEsoY2VsbGNoYXQsIHBhdHRlcm4gPSAiaW5jb21pbmciKQoKCgogblBhdHRlcm5zID0gMwpjZWxsY2hhdCA8LSBpZGVudGlmeUNvbW11bmljYXRpb25QYXR0ZXJucyhjZWxsY2hhdCwgcGF0dGVybiA9ICJpbmNvbWluZyIsIGsgPSBuUGF0dGVybnMpCgoKIyByaXZlciBwbG90Cm5ldEFuYWx5c2lzX3JpdmVyKGNlbGxjaGF0LCBwYXR0ZXJuID0gImluY29taW5nIikKIz4gUGxlYXNlIG1ha2Ugc3VyZSB5b3UgaGF2ZSBsb2FkIGBsaWJyYXJ5KGdnYWxsdXZpYWwpYCB3aGVuIHJ1bm5pbmcgdGhpcyBmdW5jdGlvbgoKCiMgZG90IHBsb3QKbmV0QW5hbHlzaXNfZG90KGNlbGxjaGF0LCBwYXR0ZXJuID0gImluY29taW5nIikKCgoKCiMtLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KIyAqKk1hbmlmb2xkIGFuZCBjbGFzc2lmaWNhdGlvbiBsZWFybmluZyBhbmFseXNpcyBvZiBzaWduYWxpbmcgbmV0d29ya3MqKgojLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKIyBJZGVudGlmeSBzaWduYWxpbmcgZ3JvdXBzIGJhc2VkIG9uIHRoZWlyIGZ1bmN0aW9uYWwgc2ltaWxhcml0eQojOC4gIElkZW50aWZ5IHNpZ25hbGluZyByb2xlcyBhbmQgcGF0dGVybnMKY2VsbGNoYXQgPC0gY29tcHV0ZU5ldFNpbWlsYXJpdHkoY2VsbGNoYXQsIHR5cGUgPSAiZnVuY3Rpb25hbCIpCmxpYnJhcnkodXdvdCkKY2VsbGNoYXQgPC0gbmV0RW1iZWRkaW5nKGNlbGxjaGF0LCB0eXBlID0gImZ1bmN0aW9uYWwiLCB1bWFwLm1ldGhvZCA9ICJ1d290IikKb3B0aW9ucyhmdXR1cmUuZ2xvYmFscy5tYXhTaXplID0gMTYwMDAgKiAxMDI0XjIpICAjIFNldCB0byAxNiBHQgpjZWxsY2hhdCA8LSBuZXRDbHVzdGVyaW5nKGNlbGxjaGF0LCB0eXBlID0gImZ1bmN0aW9uYWwiLCBkby5wYXJhbGxlbCA9IEZBTFNFKQpuZXRWaXN1YWxfZW1iZWRkaW5nKGNlbGxjaGF0LCB0eXBlID0gImZ1bmN0aW9uYWwiLCBsYWJlbC5zaXplID0gMy41KQoKIyBJZGVudGlmeSBzaWduYWxpbmcgZ3JvdXBzIGJhc2VkIG9uIHN0cnVjdHVyZSBzaW1pbGFyaXR5CiNpbnN0YWxsIHVtYXAtbGVhcm4gYmVmb3JlIHJ1bm5pbmcgdGhlc2UgY29tbWFuZHMKICNwaXAgaW5zdGFsbCB1bWFwLWxlYXJuCmNlbGxjaGF0IDwtIGNvbXB1dGVOZXRTaW1pbGFyaXR5KGNlbGxjaGF0LCB0eXBlID0gInN0cnVjdHVyYWwiKQpjZWxsY2hhdCA8LSBuZXRFbWJlZGRpbmcoY2VsbGNoYXQsIHR5cGUgPSAic3RydWN0dXJhbCIsIHVtYXAubWV0aG9kID0gInV3b3QiKQojPiBNYW5pZm9sZCBsZWFybmluZyBvZiB0aGUgc2lnbmFsaW5nIG5ldHdvcmtzIGZvciBhIHNpbmdsZSBkYXRhc2V0Cm9wdGlvbnMoZnV0dXJlLmdsb2JhbHMubWF4U2l6ZSA9IDE2MDAwICogMTAyNF4yKSAgIyBTZXQgdG8gMTYgR0IKY2VsbGNoYXQgPC0gbmV0Q2x1c3RlcmluZyhjZWxsY2hhdCwgdHlwZSA9ICJzdHJ1Y3R1cmFsIiwgZG8ucGFyYWxsZWwgPSBGQUxTRSkKIz4gQ2xhc3NpZmljYXRpb24gbGVhcm5pbmcgb2YgdGhlIHNpZ25hbGluZyBuZXR3b3JrcyBmb3IgYSBzaW5nbGUgZGF0YXNldAojIFZpc3VhbGl6YXRpb24gaW4gMkQtc3BhY2UKbmV0VmlzdWFsX2VtYmVkZGluZyhjZWxsY2hhdCwgdHlwZSA9ICJzdHJ1Y3R1cmFsIiwgbGFiZWwuc2l6ZSA9IDMuNSkKCgpuZXRWaXN1YWxfZW1iZWRkaW5nWm9vbUluKGNlbGxjaGF0LCB0eXBlID0gInN0cnVjdHVyYWwiLCBuQ29sID0gMikKCmBgYAojIDcuIFNhdmUgdGhlIENlbGxDaGF0IG9iamVjdApgYGB7ciBTQVZFLUNFTExDSEFULCBldmFsID0gRkFMU0V9CgoKc2F2ZShjZWxsY2hhdCwgZmlsZSA9ICJjZWxsQ0hBVF9BbmFseXNpc19NZWxhbmdlLlJvYmoiKQoKYGBgCg==