This vignette illustrates how to easily use the CITMIC package. Here, we present a network-based computational approach, CTMIC, estimate cell infiltration by integrating gene transcriptomic data and individualized cell interaction network.
This network-based method major parts:
Constructing network and performing randomization. In our method, we fist constructed a Cell-GO bipartite network for each sample. The molecular function gene sets were downloaded from C5 Gene sets derived from the GO biological process ontology in the Molecular Signatures Database (MSigDB). We then defined an edge between a Cell and a GO term if they have at least a common gene, and give the weight of this edge that calculate by the Jaccard index and Gene Expression Profile(GEP). Next, we proposed to construct a Cell-Cell interaction network for each sample. The information from the Cell-GO bipartite graph was used to evaluate the extent of interaction between two Cells as the edge’s weight. We consider that the extent of interaction between two Cells will be stronger when there are more GO terms they target in common and the more related they are to these GO terms. Then, we used the eigenvector centrality measure to calculate how central each Cell is in this individualized Cell-Cell interaction network. Finally, cell infiltration in the patient dataset was obtained by normalizing the centrality scores of the cells.
This package provides the GetData_CITMIC
function to
return example data set and environment variables, such as the gene
expression profile and so on.
The function CITMIC
is used to calculate the Cell
infiltration(lnScore). This function requires users to input the gene
expression matrix and the damping is the restart probability of
random walk. These sample data are stored in this package.
The commands are as follows:
# Obtain the example data
GEP<-GetData_CITMIC("GEP")
# Calculate the zscore
lnScore<-CITMIC(GEP,damping=0.90)
head(lnScore)
## TCGA.EE.A2GJ.06A TCGA.EE.A2GI.06A TCGA.WE.A8ZM.06A
## Activated B cell 0.4155478 0.5319714 0.1549675
## Activated CD4+ T cells 0.6105840 0.6376140 0.3149371
## Activated CD8+ T cells 0.3987288 0.4678214 0.1503589
## aDCs 0.9157299 0.9146414 0.8568350
## Adipocytes 0.2190414 0.3170802 0.4874105
## Angiogenesis 0.5762410 0.5273987 0.6974899
## TCGA.DA.A1IA.06A TCGA.D3.A51H.06A
## Activated B cell 0.1314579 0.7365269
## Activated CD4+ T cells 0.4463615 0.7212100
## Activated CD8+ T cells 0.2170825 0.5600082
## aDCs 0.9070054 0.9060217
## Adipocytes 0.3171054 0.2778043
## Angiogenesis 0.5255911 0.5043280