Working with SpidermiR package
2020-01-13
Contents
- Introduction
- Installation
SpidermiRquery: Searching networkSpidermiRdownload: Downloading network dataSpidermiRdownload_net: Download networkSpidermiRdownload_miRNAprediction: Downloading miRNA predicted data targetSpidermiRdownload_miRNAvalidate: Downloading miRNA validated data targetSpidermiRdownload_miRNAextra_cir:Download Extracellular Circulating microRNAsSpidermiRdownload_pharmacomir: Download Pharmaco-miR Verified Sets
SpidermiRprepare: Preparing the dataSpidermiRanalyze: : Analyze data from network dataSpidermiRanalyze_direct_net: Searching by biomarkers of interest with direct interactionSpidermiRanalyze_direct_subnetwork: Network composed by only the nodes in a set of biomarkers of interestSpidermiRanalyze_subnetwork_neigh: Network composed by the nodes in the list of BI and all the edges among this brunch of nodes.SpidermiRanalyze_degree_centrality: Ranking degree centrality genesSpidermiRanalyze_Community_detection: Find community detectionSpidermiRanalyze_Community_detection_net: Community detectionSpidermiRanalyze_Community_detection_bi: Community detection from a set of biomarkers of interest
SpidermiRvisualize: To visualize the networkSpidermiRvisualize_plot_target: Visualize the plot with miRNAs and the number of their targets in the network.SpidermiRvisualize_degree_dist: plots the degree distribution of the networkSpidermiRvisualize_adj_matrix: plots the adjacency matrix of the networkSpidermiRvisualize_3Dbarplot: 3D barplot
Features databases SpidermiR:- References
Introduction
Biological systems are composed of multiple layers of dynamic interaction networks. These networks can be decomposed, for example, into: co-expression, physical, co-localization, genetic, pathway, and shared protein domains.
GeneMania provides us with an enormous collection of data sets for interaction network studies (Warde-Farley D, Donaldson S, Comes O, Zuberi K, Badrawi R, and others 2010). The data can be accessed and downloaded from different database, using a web portal. But currently, there is not a R-package to query and download these data.
An important regulatory mechanism of these network data involves microRNAs (miRNAs). miRNAs are involved in various cellular functions, such as differentiation, proliferation, and tumourigenesis. However, our understanding of the processes regulated by miRNAs is currently limited and the integration of miRNA data in these networks provides a comprehensive genome-scale analysis of miRNA regulatory networks.Actually, GeneMania doesn’t integrate the information of miRNAs and their interactions in the network.
SpidermiR allows the user to query, prepare, download network data (e.g. from GeneMania), and to integrate this information with miRNA data
with the possibility to analyze
these downloaded data directly in one single R package.
This techincal report gives a short overview of the essential SpidermiR methods and their application.
Installation
To install use the code below.
SpidermiRquery: Searching network
You can easily search GeneMania data using the SpidermiRquery function.
SpidermiRquery_species: Searching by species
The user can query the species supported by GeneMania, using the function SpidermiRquery_species:
The list of species is shown below:
| tabOrgd[, 2] | |
|---|---|
| 1 | Arabidopsis_thaliana |
| 2 | Caenorhabditis_elegans |
| 3 | Danio_rerio |
| 4 | Drosophila_melanogaster |
| 5 | Escherichia_coli |
| 6 | Homo_sapiens |
| 7 | Mus_musculus |
| 8 | Rattus_norvegicus |
| 9 | Saccharomyces_cerevisiae |
SpidermiRquery_networks_type: Searching by network categories
The user can query the network types supported by GeneMania for a specific specie, using the function SpidermiRquery_networks_type. The user can select a specific specie using an index obtained by the function SpidermiRquery_species (e.g. organismID=org[6,] is the input for Homo_sapiens,organismID=org[9,] is the input for Saccharomyces cerevisiae )
The list of network categories in Saccharomyces cerevisiae is shown below:
## [1] "Co-localization" "Predicted" "Co-expression"
## [4] "Physical Interactions" "Genetic Interactions" "Shared protein domains"
## [7] "Other"SpidermiRquery_spec_networks: Searching by species, and network categories
You can filter the search by species using organism ID (above reported), and the network category. The network category can be filtered using the following parameters:
- COexp Co-expression
- PHint Physical_interactions
- COloc Co-localization
- GENint Genetic_interactions
- PATH Pathway
- SHpd Shared_protein_domains
- pred predicted
The databases, which data are collected, are the output of this step. An example is shown below ( for Shared protein domains in Saccharomyces_cerevisiae data are collected in INTERPRO, and PFAM):
## [1] "http://genemania.org/data/current/Saccharomyces_cerevisiae/Shared_protein_domains.INTERPRO.txt"
## [2] "http://genemania.org/data/current/Saccharomyces_cerevisiae/Shared_protein_domains.PFAM.txt"SpidermiRdownload: Downloading network data
The user in this step can download the data, as previously queried.
SpidermiRdownload_net: Download network
The user can download the data (previously queried) with SpidermiRdownload_net.
## [1] "Downloading: http://genemania.org/data/current/Saccharomyces_cerevisiae/Shared_protein_domains.INTERPRO.txt ... reference n. 1 of 2"
## [1] "Downloading: http://genemania.org/data/current/Saccharomyces_cerevisiae/Shared_protein_domains.PFAM.txt ... reference n. 2 of 2"The list of SpidermiRdownload_net is shown below:
## List of 2
## $ :'data.frame': 47523 obs. of 3 variables:
## ..$ Gene_A: chr [1:47523] "Q0050" "Q0050" "Q0055" "Q0050" ...
## ..$ Gene_B: chr [1:47523] "Q0055" "Q0060" "Q0060" "Q0065" ...
## ..$ Weight: num [1:47523] 0.39 0.09 0.15 0.09 0.15 0.23 0.1 0.17 0.18 0.18 ...
## $ :'data.frame': 30228 obs. of 3 variables:
## ..$ Gene_A: chr [1:30228] "Q0050" "Q0055" "Q0055" "Q0060" ...
## ..$ Gene_B: chr [1:30228] "Q0055" "Q0060" "Q0065" "Q0065" ...
## ..$ Weight: num [1:30228] 0.39 0.12 0.12 0.22 0.14 0.14 0.14 0.15 0.15 0.34 ...SpidermiRdownload_miRNAprediction: Downloading miRNA predicted data target
The user can download the predicted miRNA-gene from 4 databases:DIANA, Miranda, PicTar and TargetScan using miRNAtap (Pajak M, Simpson TI 2019).
SpidermiRdownload_miRNAextra_cir:Download Extracellular Circulating microRNAs
The user can download extracellular circulating miRNAs from miRandola database
SpidermiRdownload_pharmacomir: Download Pharmaco-miR Verified Sets
The user can download Pharmaco-miR predicted Sets from a fisher test. This function will give an information about the significant interactions between the miRNAs and the drug of interest (in this case tamoxifen). p-value and FDR will be visualized. The fisher test will be calculated considering the interactions drug and gene derived by DGIdb and MATADOR database (Cotto KC, Wagner AH, Feng YY, Kiwala S, Coffman AC, Spies G, Wollam A, Spies NC, Griffith OL, Griffith M. 2018) (GUnther S, Kuhn M, Dunkel M, Campillos M, Senger C, Petsalaki E, Ahmed J, Urdiales EG, Gewiess A, Jensen LJ, Schneider R, Skoblo R, Russell RB, Bourne PE, Bork P, Preissner R. 2008).
list<-SpidermiRdownload_miRNAvalidate(validated)
drug="TAMOXIFEN"
drug_genetarget<-SpidermiRdownload_pharmacomir(list[1:100,],drug="TAMOXIFEN")## [1] "processing...... hsa-let-7a....n. 1 of......3"
## [1] "processing...... hsa-let-7a-3p....n. 2 of......3"
## [1] "processing...... hsa-let-7a-5p....n. 3 of......3"The output generated by SpidermiRdownload_pharmacomir will give a table with 6 columns: in the first column miRNAs associated with the drug of interest, in the second column the number of drug targets, the number of miRNA targets, the common genes between drug targets and miRNA targets, p-value and FDR.
| drug_target | miRNA_target | common_gene | Pvalue | FDR |
|---|
The list of drugs can be chosen between the drugs generated by the function SpidermiRdownload_drug_gene:
SpidermiRprepare: Preparing the data
SpidermiRprepare_NET: Prepare matrix of gene network with Ensembl Gene ID, and gene symbols
SpidermiRprepare_NET reads network data from SpidermiRdownload_net and enables user to prepare them for downstream analysis. In particular, it prepares matrix of gene network mapping Ensembl Gene ID to gene symbols. Gene symbols are needed to integrate miRNAdata.
## [1] "Preprocessing of the network n. 1 of 2"
## [1] "Preprocessing of the network n. 2 of 2"The network with gene symbols ID is shown below:
| Gene_A | Gene_B | Weight | gene_symbolA | gene_symbolB |
|---|---|---|---|---|
| Q0050 | Q0055 | 0.39 | AI1 | AI2 |
| Q0050 | Q0060 | 0.09 | AI1 | AI3 |
| Q0055 | Q0060 | 0.15 | AI2 | AI3 |
| Q0050 | Q0065 | 0.09 | AI1 | AI4 |
| Q0055 | Q0065 | 0.15 | AI2 | AI4 |
SpidermiRanalyze: : Analyze data from network data
SpidermiRanalyze_direct_net: Searching by biomarkers of interest with direct interaction
Starting from a set of biomarkers of interest (BI), genes, miRNA or both, given by the user, this function finds sub-networks including all direct interactions involving at least one of the BI.
biomark_of_interest<-c("hsa-let-7a","CDC34","hsa-miR-27a","PEX7","EPT1","FOX","hsa-miR-5a")
miRNA_NET <-data.frame(V1=c('hsa-let-7a','CASP3','BRCA','hsa-miR-7a','hsa-miR-5a','SMAD','SOX'),V2=c('CASP3','TAMOXIFEN','MYC','PTEN','FOX','HIF1','P53'),stringsAsFactors=FALSE)
GIdirect_net<-SpidermiRanalyze_direct_net(data=miRNA_NET,BI=biomark_of_interest)## [1] "CDC34 is not in the network or please check the correct name"
## [1] "hsa-miR-27a is not in the network or please check the correct name"
## [1] "PEX7 is not in the network or please check the correct name"
## [1] "EPT1 is not in the network or please check the correct name"The data frame of SpidermiRanalyze_direct_net, GIdirect_net, is shown below:
## 'data.frame': 2 obs. of 2 variables:
## $ V1: chr "hsa-let-7a" "hsa-miR-5a"
## $ V2: chr "CASP3" "FOX"SpidermiRanalyze_direct_subnetwork: Network composed by only the nodes in a set of biomarkers of interest
Starting from BI, this function finds sub-networks including all direct interactions involving only BI.
SpidermiRanalyze_subnetwork_neigh: Network composed by the nodes in the list of BI and all the edges among this brunch of nodes.
Starting from BI, this function finds sub-networks including all direct and indirect interactions involving at least one of BI.
SpidermiRanalyze_degree_centrality: Ranking degree centrality genes
This function finds the number of direct neighbours of a node in a network and allows the selection of those nodes with a number of direct neighbours higher than a selected cut-off.
SpidermiRanalyze_Community_detection: Find community detection
This function find the communities in the network, and describes them in terms of number of community elements (both genes and miRNAs). The function uses one of the algorithms currently implemented in (Csardi G, Nepusz T. 2006), selected by the user according to the user need.
The user can choose the algorithm in order to calculate the community structure:
- EB edge.betweenness.community
- FC fastgreedy.community
- WC walktrap.community
- SC spinglass.community
- LE leading.eigenvector.community
- LP label.propagation.community
SpidermiRanalyze_Community_detection_net: Community detection
Starting from one community to which some BI belong (the output of the previously described function) this function describes the community as network of elements (both genes and miRNAs).
SpidermiRanalyze_Community_detection_bi: Community detection from a set of biomarkers of interest
Starting from the community to which BI belong (the output of the previously described function), this function indicates if a set of BI is included within such community.
SpidermiRvisualize: To visualize the network
SpidermiRvisualize_plot_target: Visualize the plot with miRNAs and the number of their targets in the network.
For each BI of a community, the user can visualize a plot showing the number of direct neighbours of such BI (the degree centrality of such BI).
## NULLSpidermiRvisualize_degree_dist: plots the degree distribution of the network
This function plots the cumulative frequency distribution of degree centrality of a community.
SpidermiRvisualize_adj_matrix: plots the adjacency matrix of the network
It plots the adjacency matrix of the community, representing the degree of connections among the nodes.
SpidermiRvisualize_3Dbarplot: 3D barplot
It plots a summary representation of the networks with the number of edges, nodes and miRNAs.
Features databases SpidermiR:
Features of databases integrated in SpidermiR are:
Session Information ******
## R version 3.6.2 (2019-12-12)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 18.04.3 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.10-bioc/R/lib/libRblas.so
## LAPACK: /home/biocbuild/bbs-3.10-bioc/R/lib/libRlapack.so
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_US.UTF-8 LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## attached base packages:
## [1] parallel stats4 stats graphics grDevices utils datasets
## [8] methods base
##
## other attached packages:
## [1] SpidermiR_1.16.2 miRNAtap_1.20.0 AnnotationDbi_1.48.0
## [4] IRanges_2.20.2 S4Vectors_0.24.2 Biobase_2.46.0
## [7] BiocGenerics_0.32.0 BiocStyle_2.14.4
##
## loaded via a namespace (and not attached):
## [1] rappdirs_0.3.1 rtracklayer_1.46.0
## [3] ggthemes_4.2.0 R.methodsS3_1.7.1
## [5] tidyr_1.0.0 ggplot2_3.2.1
## [7] bit64_0.9-7 knitr_1.26
## [9] aroma.light_3.16.0 DelayedArray_0.12.2
## [11] R.utils_2.9.2 data.table_1.12.8
## [13] hwriter_1.3.2 KEGGREST_1.26.1
## [15] RCurl_1.95-4.12 doParallel_1.0.15
## [17] generics_0.0.2 GenomicFeatures_1.38.0
## [19] cowplot_1.0.0 callr_3.4.0
## [21] usethis_1.5.1 RSQLite_2.2.0
## [23] europepmc_0.3 chron_2.3-54
## [25] bit_1.1-15 enrichplot_1.6.1
## [27] xml2_1.2.2 ggsci_2.9
## [29] SummarizedExperiment_1.16.1 assertthat_0.2.1
## [31] viridis_0.5.1 xfun_0.12
## [33] hms_0.5.3 evaluate_0.14
## [35] fansi_0.4.1 progress_1.2.2
## [37] caTools_1.17.1.4 dbplyr_1.4.2
## [39] Rgraphviz_2.30.0 km.ci_0.5-2
## [41] igraph_1.2.4.2 DBI_1.1.0
## [43] geneplotter_1.64.0 htmlwidgets_1.5.1
## [45] EDASeq_2.20.0 matlab_1.0.2
## [47] purrr_0.3.3 ellipsis_0.3.0
## [49] selectr_0.4-2 dplyr_0.8.3
## [51] ggpubr_0.2.4 backports_1.1.5
## [53] bookdown_0.17 annotate_1.64.0
## [55] MAGeCKFlute_1.6.2 biomaRt_2.42.0
## [57] vctrs_0.2.1 remotes_2.1.0
## [59] pathview_1.26.0 withr_2.1.2
## [61] ggforce_0.3.1 triebeard_0.3.0
## [63] GenomicAlignments_1.22.1 prettyunits_1.1.0
## [65] cluster_2.1.0 DOSE_3.12.0
## [67] lazyeval_0.2.2 crayon_1.3.4
## [69] genefilter_1.68.0 labeling_0.3
## [71] edgeR_3.28.0 pkgconfig_2.0.3
## [73] tweenr_1.0.1 GenomeInfoDb_1.22.0
## [75] nlme_3.1-143 pkgload_1.0.2
## [77] devtools_2.2.1 rlang_0.4.2
## [79] lifecycle_0.1.0 downloader_0.4
## [81] BiocFileCache_1.10.2 rprojroot_1.3-2
## [83] polyclip_1.10-0 matrixStats_0.55.0
## [85] graph_1.64.0 urltools_1.7.3
## [87] Matrix_1.2-18 KMsurv_0.1-5
## [89] zoo_1.8-7 ggridges_0.5.2
## [91] GlobalOptions_0.1.1 processx_3.4.1
## [93] pheatmap_1.0.12 viridisLite_0.3.0
## [95] png_0.1-7 rjson_0.2.20
## [97] bitops_1.0-6 R.oo_1.23.0
## [99] ConsensusClusterPlus_1.50.0 KernSmooth_2.23-16
## [101] visNetwork_2.0.9 Biostrings_2.54.0
## [103] blob_1.2.0 TCGAbiolinks_2.14.0
## [105] shape_1.4.4 stringr_1.4.0
## [107] qvalue_2.18.0 ShortRead_1.44.1
## [109] gridGraphics_0.4-1 readr_1.3.1
## [111] jpeg_0.1-8.1 ggsignif_0.6.0
## [113] scales_1.1.0 memoise_1.1.0
## [115] magrittr_1.5 plyr_1.8.5
## [117] gplots_3.0.1.2 gdata_2.18.0
## [119] zlibbioc_1.32.0 compiler_3.6.2
## [121] RColorBrewer_1.1-2 clue_0.3-57
## [123] KEGGgraph_1.46.0 Rsamtools_2.2.1
## [125] miRNAtap.db_0.99.10 cli_2.0.1
## [127] XVector_0.26.0 ps_1.3.0
## [129] MASS_7.3-51.5 mgcv_1.8-31
## [131] tidyselect_0.2.5 stringi_1.4.5
## [133] highr_0.8 yaml_2.2.0
## [135] GOSemSim_2.12.0 askpass_1.1
## [137] locfit_1.5-9.1 latticeExtra_0.6-29
## [139] ggrepel_0.8.1 survMisc_0.5.5
## [141] grid_3.6.2 fastmatch_1.1-0
## [143] tools_3.6.2 circlize_0.4.8
## [145] rstudioapi_0.10 foreach_1.4.7
## [147] gridExtra_2.3 farver_2.0.2
## [149] ggraph_2.0.0 rvcheck_0.1.7
## [151] digest_0.6.23 BiocManager_1.30.10
## [153] proto_1.0.0 networkD3_0.4
## [155] Rcpp_1.0.3 GenomicRanges_1.38.0
## [157] broom_0.5.3 org.Hs.eg.db_3.10.0
## [159] httr_1.4.1 survminer_0.4.6
## [161] ComplexHeatmap_2.2.0 colorspace_1.4-1
## [163] rvest_0.3.5 XML_3.98-1.20
## [165] fs_1.3.1 splines_3.6.2
## [167] graphlayouts_0.5.0 ggplotify_0.0.4
## [169] sessioninfo_1.1.1 xtable_1.8-4
## [171] jsonlite_1.6 tidygraph_1.1.2
## [173] zeallot_0.1.0 testthat_2.3.1
## [175] R6_2.4.1 gsubfn_0.7
## [177] pillar_1.4.3 htmltools_0.4.0
## [179] glue_1.3.1 clusterProfiler_3.14.3
## [181] BiocParallel_1.20.1 DESeq_1.38.0
## [183] codetools_0.2-16 fgsea_1.12.0
## [185] pkgbuild_1.0.6 lattice_0.20-38
## [187] tibble_2.1.3 sva_3.34.0
## [189] sqldf_0.4-11 curl_4.3
## [191] gtools_3.8.1 GO.db_3.10.0
## [193] openssl_1.4.1 survival_3.1-8
## [195] limma_3.42.0 rmarkdown_2.0
## [197] desc_1.2.0 munsell_0.5.0
## [199] DO.db_2.9 GetoptLong_0.1.8
## [201] GenomeInfoDbData_1.2.2 iterators_1.0.12
## [203] reshape2_1.4.3 gtable_0.3.0References
Cotto KC, Wagner AH, Feng YY, Kiwala S, Coffman AC, Spies G, Wollam A, Spies NC, Griffith OL, Griffith M. 2018. “DGIdb 3.0 a Redesign and Expansion of the Drug Gene Interaction Database.”
Csardi G, Nepusz T. 2006. “The Igraph Software Package for Complex Network Research.”
Dweep H, Sticht C, Pandey P, Gretz N. 2011. “miRWalk - Database Prediction of Possible miRNA Binding Sites by ‘Walking’ the Genes of 3 Genomes.”
GUnther S, Kuhn M, Dunkel M, Campillos M, Senger C, Petsalaki E, Ahmed J, Urdiales EG, Gewiess A, Jensen LJ, Schneider R, Skoblo R, Russell RB, Bourne PE, Bork P, Preissner R. 2008. “SuperTarget and Matador Resources for Exploring Drug-Target Relationships.”
Jiang Q, Wang Y, Hao Y, Juan L, Teng M, Zhang X, Li M, Wang G, Liu Y. 2009. “miR2Disease a Manually Curated Database for microRNA Deregulation in Human Disease.”
Pajak M, Simpson TI. 2019. “miRNAtap microRNA Targets - Aggregated Predictions.”
Warde-Farley D, Donaldson S, Comes O, Zuberi K, Badrawi R, and others. 2010. “The Gene Mania Prediction Server Biological Network Integration for Gene Prioritization and Predicting Gene Function.”