This vignette demonstrates how a Protein-Protein interaction (PPI) graph may be constructed from the database stringDB.
Other networks can also be used with netSmooth. We mostly rely on networks from stringDB. StringDB has multiple species available such as human, mouse, zebrafish, C.elengas and D.melanogaster. It is also possible to prune the network differently. For our purposes we use the edges that have highest confidence score. Below, we are showing how to obtain and prune human network from stringDB. Specifically, we use the work flow below.
require(STRINGdb)
require(igraph)
require(biomaRt)
# 1. getSTRINGdb for human
string_db <- STRINGdb$new(species=9606)
human_graph <- string_db$get_graph()
# 2. get edges with high confidence score
edge.scores <- E(human_graph)$combined_score
ninetyth.percentile <- quantile(edge.scores, 0.9)
thresh <- data.frame(name='90th percentile',
val=ninetyth.percentile)
human_graph <- subgraph.edges(human_graph,
E(human_graph)[combined_score > ninetyth.percentile])
# 3. create adjacency matrix
adj_matrix <- as_adjacency_matrix(human_graph)
# 4. map gene ids to protein ids
### get gene/protein ids via Biomart
mart=useMart(host = 'grch37.ensembl.org',
biomart='ENSEMBL_MART_ENSEMBL',
dataset='hsapiens_gene_ensembl')
### extract protein ids from the human network
protein_ids <- sapply(strsplit(rownames(adj_matrix), '\\.'),
function(x) x[2])
### get protein to gene id mappings
mart_results <- getBM(attributes = c("ensembl_gene_id",
"ensembl_peptide_id"),
filters = "ensembl_peptide_id", values = protein_ids,
mart = mart)
### replace protein ids with gene ids
ix <- match(protein_ids, mart_results$ensembl_peptide_id)
ix <- ix[!is.na(ix)]
newnames <- protein_ids
newnames[match(mart_results[ix,'ensembl_peptide_id'], newnames)] <-
mart_results[ix, 'ensembl_gene_id']
rownames(adj_matrix) <- newnames
colnames(adj_matrix) <- newnames
ppi <- adj_matrix[!duplicated(newnames), !duplicated(newnames)]
nullrows <- Matrix::rowSums(ppi)==0
ppi <- ppi[!nullrows,!nullrows] ## ppi is the network with gene idssessionInfo()## R version 3.6.0 (2019-04-26)
## Platform: x86_64-w64-mingw32/x64 (64-bit)
## Running under: Windows Server 2012 R2 x64 (build 9600)
##
## Matrix products: default
##
## locale:
## [1] LC_COLLATE=C
## [2] LC_CTYPE=English_United States.1252
## [3] LC_MONETARY=English_United States.1252
## [4] LC_NUMERIC=C
## [5] LC_TIME=English_United States.1252
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] BiocStyle_2.12.0
##
## loaded via a namespace (and not attached):
## [1] BiocManager_1.30.4 compiler_3.6.0 magrittr_1.5
## [4] bookdown_0.9 htmltools_0.3.6 tools_3.6.0
## [7] yaml_2.2.0 Rcpp_1.0.1 stringi_1.4.3
## [10] rmarkdown_1.12 knitr_1.22 stringr_1.4.0
## [13] digest_0.6.18 xfun_0.6 evaluate_0.13