Case Study n.1: Role of miRNAs in shared protein domains network in Prostate Cancer
In this case study, we downloaded shared protein domains network in Homo Sapiens, using SpidermiRquery, SpidermiRprepare, and SpidermiRdownload.
Then, we focused on role of miRNAs in this network. We integrated miRNA information using SpidermiRanalyze. We obtained a big network, and in order to understand the underlying biological process of a set of biomarker of interest (e.g. from lab) we performed an analysis to identify their neighbor biomarkers in the shared protein domains network.
SpidermiRvisualize was used to see the results.
org<-SpidermiRquery_species(species)
net_shar_prot<-SpidermiRquery_spec_networks(organismID = org[6,],network = "SHpd")
out_net<-SpidermiRdownload_net(net_shar_prot)
geneSymb_net<-SpidermiRprepare_NET(organismID = org[6,],data = out_net)
miRNA_complNET<-SpidermiRanalyze_mirna_gene_complnet(data=geneSymb_net,disease="prostate cancer",miR_trg="val")
biomark_of_interest<-read.delim("C:/Users/UserInLab05/Google Drive/MIRNA AND GENEMANIA/1 case study/deg_prostate.txt",header=FALSE)
subnet<-SpidermiRanalyze_direct_subnetwork(data=miRNA_complNET,BI=biomark_of_interest$V1)
comm2<- SpidermiRanalyze_Community_detection(data=subnet,type="FC")
cd_net<-SpidermiRanalyze_Community_detection_net(data=subnet,comm_det=comm2,size=2)
SpidermiRvisualize_mirnanet(data=cd_net)
miRNA_NET<-SpidermiRanalyze_mirna_network(data=geneSymb_net,disease="prostate cancer")
cd_net_miRNA<-SpidermiRanalyze_Community_detection_net(data=miRNA_NET,comm_det=comm2,size=2)
SpidermiRvisualize_mirnanet(data=cd_net_miRNA)
Case Study n.2: miRNAs regulating degree centrality genes in physical interactions network in breast cancer
org<-SpidermiRquery_species(species)
net_PHint<-SpidermiRquery_spec_networks(organismID = org[6,],network = "PHint")
out_net<-SpidermiRdownload_net(net_PHint)
geneSymb_net<-SpidermiRprepare_NET(organismID = org[6,],data = out_net)
ds<-do.call("rbind", geneSymb_net)
data1<-as.data.frame(ds[!duplicated(ds), ])
sdas<-cbind(data1$gene_symbolA,data1$gene_symbolB)
sdas<-as.data.frame(sdas[!duplicated(sdas), ])
miRNA_NET<-SpidermiRanalyze_mirna_network(data=geneSymb_net,disease="breast cancer")
topwhol<-SpidermiRanalyze_degree_centrality(sdas)
top10_cent_gene<-SpidermiRanalyze_degree_centrality(miRNA_NET)
miRNA_degree<-top10_cent_gene[grep("hsa",top10_cent_gene$dfer),]
Session Information ******
## R version 3.3.1 (2016-06-21)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 14.04.4 LTS
##
## 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] visNetwork_1.0.2 networkD3_0.2.13 SpidermiR_1.2.5
## [4] miRNAtap_1.6.0 AnnotationDbi_1.34.4 IRanges_2.6.1
## [7] S4Vectors_0.10.3 Biobase_2.32.0 BiocGenerics_0.18.0
## [10] BiocStyle_2.0.3
##
## loaded via a namespace (and not attached):
## [1] circlize_0.3.9
## [2] aroma.light_3.2.0
## [3] plyr_1.8.4
## [4] igraph_1.0.1
## [5] ConsensusClusterPlus_1.36.0
## [6] splines_3.3.1
## [7] BiocParallel_1.6.6
## [8] GenomeInfoDb_1.8.7
## [9] ggplot2_2.1.0
## [10] TH.data_1.0-7
## [11] digest_0.6.10
## [12] foreach_1.4.3
## [13] BiocInstaller_1.22.3
## [14] htmltools_0.3.5
## [15] gdata_2.17.0
## [16] magrittr_1.5
## [17] memoise_1.0.0
## [18] doParallel_1.0.10
## [19] cluster_2.0.4
## [20] limma_3.28.21
## [21] ComplexHeatmap_1.10.2
## [22] Biostrings_2.40.2
## [23] readr_1.0.0
## [24] annotate_1.50.0
## [25] matrixStats_0.50.2
## [26] R.utils_2.4.0
## [27] sandwich_2.3-4
## [28] colorspace_1.2-6
## [29] rvest_0.3.2
## [30] ggrepel_0.5
## [31] dplyr_0.5.0
## [32] crayon_1.3.2
## [33] RCurl_1.95-4.8
## [34] jsonlite_1.1
## [35] hexbin_1.27.1
## [36] graph_1.50.0
## [37] TxDb.Hsapiens.UCSC.hg19.knownGene_3.2.2
## [38] roxygen2_5.0.1
## [39] genefilter_1.54.2
## [40] supraHex_1.10.0
## [41] iterators_1.0.8
## [42] survival_2.39-5
## [43] miRNAtap.db_0.99.7
## [44] zoo_1.7-13
## [45] ape_3.5
## [46] gtable_0.2.0
## [47] zlibbioc_1.18.0
## [48] XVector_0.12.1
## [49] GetoptLong_0.1.5
## [50] kernlab_0.9-25
## [51] Rgraphviz_2.16.0
## [52] shape_1.4.2
## [53] prabclus_2.2-6
## [54] DEoptimR_1.0-6
## [55] scales_0.4.0
## [56] DESeq_1.24.0
## [57] mvtnorm_1.0-5
## [58] edgeR_3.14.0
## [59] DBI_0.5-1
## [60] GGally_1.2.0
## [61] ggthemes_3.2.0
## [62] Rcpp_0.12.7
## [63] xtable_1.8-2
## [64] matlab_1.0.2
## [65] mclust_5.2
## [66] preprocessCore_1.34.0
## [67] sqldf_0.4-10
## [68] htmlwidgets_0.7
## [69] httr_1.2.1
## [70] gplots_3.0.1
## [71] RColorBrewer_1.1-2
## [72] fpc_2.1-10
## [73] modeltools_0.2-21
## [74] reshape_0.8.5
## [75] XML_3.98-1.4
## [76] R.methodsS3_1.7.1
## [77] flexmix_2.3-13
## [78] nnet_7.3-12
## [79] labeling_0.3
## [80] munsell_0.4.3
## [81] tools_3.3.1
## [82] downloader_0.4
## [83] gsubfn_0.6-6
## [84] RSQLite_1.0.0
## [85] devtools_1.12.0
## [86] evaluate_0.9
## [87] stringr_1.1.0
## [88] yaml_2.1.13
## [89] org.Hs.eg.db_3.3.0
## [90] knitr_1.14
## [91] robustbase_0.92-6
## [92] caTools_1.17.1
## [93] dendextend_1.3.0
## [94] coin_1.1-2
## [95] TCGAbiolinks_2.0.13
## [96] EDASeq_2.6.2
## [97] nlme_3.1-128
## [98] whisker_0.3-2
## [99] formatR_1.4
## [100] R.oo_1.20.0
## [101] xml2_1.0.0
## [102] biomaRt_2.28.0
## [103] curl_2.1
## [104] testthat_1.0.2
## [105] affyio_1.42.0
## [106] tibble_1.2
## [107] geneplotter_1.50.0
## [108] stringi_1.1.2
## [109] highr_0.6
## [110] GenomicFeatures_1.24.5
## [111] lattice_0.20-34
## [112] trimcluster_0.1-2
## [113] Matrix_1.2-7.1
## [114] GlobalOptions_0.0.10
## [115] parmigene_1.0.2
## [116] data.table_1.9.6
## [117] bitops_1.0-6
## [118] dnet_1.0.9
## [119] rtracklayer_1.32.2
## [120] GenomicRanges_1.24.3
## [121] R6_2.2.0
## [122] latticeExtra_0.6-28
## [123] affy_1.50.0
## [124] hwriter_1.3.2
## [125] ShortRead_1.30.0
## [126] KernSmooth_2.23-15
## [127] gridExtra_2.2.1
## [128] codetools_0.2-15
## [129] MASS_7.3-45
## [130] gtools_3.5.0
## [131] assertthat_0.1
## [132] chron_2.3-47
## [133] SummarizedExperiment_1.2.3
## [134] proto_0.3-10
## [135] rjson_0.2.15
## [136] withr_1.0.2
## [137] GenomicAlignments_1.8.4
## [138] Rsamtools_1.24.0
## [139] multcomp_1.4-6
## [140] diptest_0.75-7
## [141] grid_3.3.1
## [142] class_7.3-14
## [143] rmarkdown_1.0
