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 from PharmacomiR databaseSpidermiRprepare
: Preparing the dataSpidermiRanalyze
: : Analyze data from network dataSpidermiRanalyze_mirna_network
: Integration of microRNA-target interactions.SpidermiRanalyze_mirna_gene_complnet
: Integration of microRNA-target complete interactionsSpidermiRanalyze_mirnanet_pharm
: Integration of pharmacomiR in the networkSpidermiRanalyze_mirna_extra_cir
: Integration of Extracellular/Circulating miRNASpidermiRanalyze_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 interestSpidermiRanalyze_DEnetworkTCGA
: Integration with TCGA data in order to obtain a network of differentially expressed (DE) genes or miRNAs.SpidermiRvisualize
: To visualize the networkSpidermiRvisualize_mirnanet
: To Visualize the network.SpidermiRvisualize_BI
: To Visualize the network from a set of BI.SpidermiRvisualize_direction
: 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 barplotFeatures databases SpidermiR
:SpidermiR:Application Examples
for more information.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.
To install use the code below.
source("https://bioconductor.org/biocLite.R")
biocLite("SpidermiR")
SpidermiRquery
: Searching networkYou can easily search GeneMania data using the SpidermiRquery
function.
SpidermiRquery_species
: Searching by speciesThe user can query the species supported by GeneMania, using the function SpidermiRquery_species:
org<-SpidermiRquery_species(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 categoriesThe 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 )
net_type<-SpidermiRquery_networks_type(organismID=org[9,])
The list of network categories in Saccharomyces cerevisiae is shown below:
## [1] "Co-localization" "Predicted"
## [3] "Co-expression" "Physical Interactions"
## [5] "Genetic Interactions" "Shared protein domains"
## [7] "Other"
SpidermiRquery_spec_networks
: Searching by species, and network categoriesYou 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:
net_shar_prot<-SpidermiRquery_spec_networks(organismID = org[9,],
network = "SHpd")
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"
SpidermiRquery_disease
: Searching by miRNA-diseaseThe user can obtain a list of the diseases supported by SpidermiR, in order to focus only on miRNAs that have been already studied in a particular disease (retrieving data from miR2Disease (Jiang Q, Wang Y, Hao Y, Juan L, Teng M, Zhang X, Li M, Wang G, Liu Y. 2009)).
disease<-SpidermiRquery_disease(diseaseID)
Currently, the list of diseases supported is:
## [1] "kidney cancer"
## [2] "hepatocellular carcinoma (HCC)"
## [3] "lung cancer"
## [4] "non-small cell lung cancer (NSCLC)"
## [5] "ovarian cancer (OC)"
## [6] "colorectal cancer"
## [7] "head and neck squamous cell carcinoma (HNSCC)"
## [8] "malignant melanoma"
## [9] "Oral Squamous Cell Carcinoma (OSCC)"
## [10] "prostate cancer"
## [11] "Alzheimer's disease"
## [12] "breast cancer"
## [13] "cardiac hypertrophy"
## [14] "pancreatic ductal adenocarcinoma (PDAC)"
## [15] "uterine leiomyoma (ULM)"
## [16] "acute myeloid leukemia (AML)"
## [17] "Duchenne muscular dystrophy (DMD)"
## [18] "limb-girdle muscular dystrophies types 2A (LGMD2A)"
## [19] "miyoshi myopathy (MM)"
## [20] "Insulinoma"
## [21] "acute promyelocytic leukemia (APL)"
## [22] "cardiomyopathy"
## [23] "diffuse large B-cell lymphoma (DLBCL)"
## [24] "epithelial ovarian cancer (EOC)"
## [25] "pancreatic cancer"
## [26] "psoriasis"
## [27] "autism spectrum disorder (ASD)"
## [28] "acute lymphoblastic leukemia (ALL)"
## [29] "pituitary adenoma"
## [30] "glioblastoma multiforme (GBM)"
## [31] "chronic pancreatitis"
## [32] "neuroblastoma (NB)"
## [33] "serous ovarian cancer"
## [34] "vascular disease"
## [35] "anaplastic thyroid carcinoma (ATC)"
## [36] "Cerebellar neurodegeneration"
## [37] "heart failure"
## [38] "polycythemia vera (PV)"
## [39] "primary biliary cirrhosis (PBC)"
## [40] "nasopharyngeal carcinoma (NPC)"
## [41] "asthma"
## [42] "B-cell chronic lymphocytic leukemia"
## [43] "dermatomyositis (DM)"
## [44] "tongue squamous cell carcinoma"
## [45] "retinitis pigmentosa (RP)"
## [46] "Parkinson's disease"
## [47] "testicular germ cell tumor"
## [48] "nemaline myopathy (NM)"
## [49] "papillary thyroid carcinoma (PTC)"
## [50] "polymyositis (PM)"
## [51] "follicular lymphoma (FL)"
## [52] "Hodgkin's lymphoma"
## [53] "facioscapulohumeral muscular dystrophy (FSHD)"
## [54] "cholangiocarcinoma"
## [55] "rhabdomyosarcoma"
## [56] "myocardial injury"
## [57] "myocardial infarction"
## [58] "cervical cancer"
## [59] "cholesteatoma"
## [60] "chronic lymphocytic leukemia (CLL)"
## [61] "multiple sclerosis"
## [62] "schizophrenia"
## [63] "type 2 diabetes"
## [64] "gastric cancer (stomach cancer)"
## [65] "neuroblastoma"
## [66] "Intrahepatic cholangiocarcinoma (ICC)"
## [67] "bladder cancer"
## [68] "osteosarcoma"
## [69] "glioblastoma"
## [70] "glioma"
## [71] "medulloblastoma"
## [72] "Oral Carcinoma"
## [73] "non-alcoholic fatty liver disease (NAFLD)"
## [74] "esophageal cancer"
## [75] "head and neck cancer"
## [76] "laryngeal carcinoma"
## [77] "essential thrombocythemia (ET)"
## [78] "primary myelofibrosis"
## [79] "Polycystic Kidney Disease"
## [80] "sepsis"
## [81] "acute myelogeneous leukemia (AML)"
## [82] "multiple myeloma (MM)"
## [83] "mantle cell lymphoma (MCL)"
## [84] "MYC-rearranged lymphoma"
## [85] "malignant lymphoma"
## [86] "uveal melanoma"
## [87] "Down syndrome (DS)"
## [88] "metabolic disease"
## [89] "recurrent ovarian cancer"
## [90] "HCV infection"
## [91] "Huntington's disease (HD)"
## [92] "cancer"
## [93] "PFV-1 infection"
## [94] "Cowden Syndrome"
## [95] "T-cell leukemia"
## [96] "Becker muscular dystrophy (BMD)"
## [97] "lupus nephritis"
## [98] "neurodegeneration"
## [99] "Acute Promyelocytic Leukemia (APL)"
## [100] "thalassemia"
## [101] "Hepatocellular carcinoma (HCC)"
## [102] "endometriosis"
## [103] "medulloblastoma "
## [104] "renal clear cell carcinoma"
## [105] "pulmonary hypertension"
## [106] "adrenocortical carcinoma"
## [107] "squamous carcinoma"
## [108] "Malignant mesothelioma (MM)"
## [109] "myeloproliferative disorder"
## [110] "coronary artery disease"
## [111] "meningioma"
## [112] "prostate cance"
## [113] "prion disease"
## [114] "retinoblastoma"
## [115] "anxiety disorder"
## [116] "chronic myeloid leukemia (CML)"
## [117] "skin disease"
## [118] "homozygous sickle cell disease (HbSS)"
## [119] "congenital heart disease"
## [120] "Burkitt lymphoma"
## [121] "endometrial cancer"
## [122] "Inclusion body myositis (IBM)"
## [123] "hearing loss"
## [124] "melanoma"
## [125] "Waldenstrom Macroglobulinemia (WM)"
## [126] "diabetic nephropathy"
## [127] "thyroid cancer"
## [128] "lymphoproliferative disease"
## [129] "Obesity"
## [130] "Spinocerebellar ataxia 1"
## [131] "alcoholic liver disease (ALD)"
## [132] "Polycystic liver disease"
## [133] "teratocarcinoma"
## [134] "ulcerative colitis (UC)"
## [135] "vesicular stomatitis"
## [136] "Glomerulosclerosis"
## [137] "B-cell lymphoma"
## [138] "glomerular disease"
## [139] "follicular thyroid carcinoma (FTC)"
## [140] "hypertension"
## [141] "adenoma"
## [142] "hamartoma"
## [143] "lipoma"
## [144] "myoma"
## [145] "sarcoma"
## [146] "neutrophilia"
## [147] "diarrhea predominant irritable bowel syndrome (IBS-D)"
## [148] "HBV-related cirrhosis"
## [149] "frontotemporal dementia"
## [150] "rhabdomyosarcoma (RMS)"
## [151] "Head and neck cancer"
## [152] "tourette's syndrome"
SpidermiRdownload
: Downloading network dataThe user in this step can download the data, as previously queried.
SpidermiRdownload_net
: Download networkThe user can download the data (previously queried) with SpidermiRdownload_net
.
out_net<-SpidermiRdownload_net(net_shar_prot)
## [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 targetThe user can download the predicted miRNA-gene from 4 databases:DIANA, Miranda, PicTar and TargetScan
mirna<-c('hsa-miR-567','hsa-miR-566')
SpidermiRdownload_miRNAprediction(mirna_list=mirna)
SpidermiRdownload_miRNAextra_cir
:Download Extracellular Circulating microRNAsThe user can download extracellular circulating miRNAs from miRandola database
list<-SpidermiRdownload_miRNAextra_cir(miRNAextra_cir)
SpidermiRdownload_pharmacomir
: Download Pharmaco-miR Verified Sets from PharmacomiR databaseThe user can download Pharmaco-miR Verified Sets from PharmacomiR database (Rukov J, Wilentzik R, Jaffe I, Vinther J, Shomron N. 2013).
mir_pharmaco<-SpidermiRdownload_pharmacomir(pharmacomir=pharmacomir)
SpidermiRprepare
: Preparing the dataSpidermiRprepare_NET
: Prepare matrix of gene network with Ensembl Gene ID, and gene symbolsSpidermiRprepare_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.
geneSymb_net<-SpidermiRprepare_NET(organismID = org[9,],
data = out_net)
## [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 dataSpidermiRanalyze_mirna_network
: Integration of microRNA-target interactions.The user in this step obtained a network matrix with miRNA-target interactions starting from a specific network. The user can focus on miRNAs that have been already linked to a particular disease or take all miRNAs. miRNA-gene interactions include data from validated or predicted studies. You can filter the search by disease. The miRNA network can be filtered by disease using the name of the disease, as obtained from SpidermiRquery_disease
.
miRNA_NET<-SpidermiRanalyze_mirna_network(data=geneSymb_net,disease="prostate cancer",miR_trg="val")
The list of miRNA network is shown below:
## 'data.frame': 70 obs. of 2 variables:
## $ V1: chr "hsa-miR-127" "hsa-miR-141" "hsa-miR-210" "hsa-miR-222" ...
## $ V2: chr "XBP1" "SIP1" "PIM1" "MMP1" ...
SpidermiRanalyze_mirna_gene_complnet
: Integration of microRNA-target complete interactionsThe user in this step obtained a gene network matrix with miRNA-gene and gene-gene interaction. The user can focus on miRNAs that have been already linked to a particular disease or take all miRNAs. miRNA-gene interactions include data from validated or predicted studies. The miRNA network can be filtered by disease using the name of the disease, as obtained from SpidermiRquery_disease
.
miRNA_complNET<-SpidermiRanalyze_mirna_gene_complnet(data=geneSymb_net,disease="prostate cancer",miR_trg="val")
SpidermiRanalyze_mirnanet_pharm
: Integration of pharmacomiR in the networkThe user in this step can integrate the pharmacomiR database in order to link miRNA and drug effect in a specific network.
mir_pharmnet<-SpidermiRanalyze_mirnanet_pharm(mir_ph=mir_pharmaco,net=miRNA_NET)
SpidermiRanalyze_mirna_extra_cir
: Integration of Extracellular/Circulating miRNAThe user can select the extracellular/circulating miRNAs found in the obtained network. Extracellular/circulating miRNAs include data from mirandola database (Russo F, Di Bella S, Nigita G, Macca V, Lagana A, Giugno R, Pulvirenti A, Ferro A. 2012). The user using the following parameteres can specify the network type:
miRNA_NET_ext_circmT<-SpidermiRanalyze_mirna_extra_cir(data=miRNA_complNET,"mT")
miRNA_NET_ext_circmCT<-SpidermiRanalyze_mirna_extra_cir(data=miRNA_complNET,"mCT")
SpidermiRanalyze_direct_net
: Searching by biomarkers of interest with direct interactionStarting 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-miR-214","PTEN","FOXO1","hsa-miR-27a")
GIdirect_net<-SpidermiRanalyze_direct_net(data=miRNA_NET,BI=biomark_of_interest)
## [1] "hsa-miR-214 is not in the network or please check the correct name"
## [1] "PTEN is not in the network or please check the correct name"
## [1] "FOXO1 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': 1 obs. of 2 variables:
## $ V1: chr "hsa-miR-27a"
## $ V2: chr "PEX7"
SpidermiRanalyze_direct_subnetwork
: Network composed by only the nodes in a set of biomarkers of interestStarting from BI, this function finds sub-networks including all direct interactions involving only BI.
subnet<-SpidermiRanalyze_direct_subnetwork(data=miRNA_NET,BI=biomark_of_interest)
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.
GIdirect_net_neigh<-SpidermiRanalyze_subnetwork_neigh(data=miRNA_NET,BI=biomark_of_interest)
SpidermiRanalyze_degree_centrality
: Ranking degree centrality genesThis function finds the number of direct neighbours of a BI and allows the selection of those BI with a number of direct neighbours higher than a selected cut-off.
top10_cent_gene<-SpidermiRanalyze_degree_centrality(miRNA_NET,cut=10)
SpidermiRanalyze_Community_detection
: Find community detectionThis 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:
comm<- SpidermiRanalyze_Community_detection(data=miRNA_NET,type="FC")
SpidermiRanalyze_Community_detection_net
: Community detectionStarting 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).
cd_net<-SpidermiRanalyze_Community_detection_net(data=miRNA_NET,comm_det=comm,size=1)
SpidermiRanalyze_Community_detection_bi
: Community detection from a set of biomarkers of interestStarting 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.
gi=c("CF","ROCK1","KIT","CCND2")
mol<-SpidermiRanalyze_Community_detection_bi(data=comm,BI=gi)
SpidermiRanalyze_DEnetworkTCGA
: Integration with TCGA data in order to obtain a network of differentially expressed (DE) genes or miRNAs.SpidermiRanalyze_DEnetworkTCGA
integrates the information of differential analysis of TCGA data in the network. The final result will be a network with only DE genes or miRNAs depending whether the user chooses mRNA or miRNA TCGA data.
miRNA_cN <-data.frame(gA=c('IGFL3','GABRA1'),gB=c('IGFL2','KRT13'),stringsAsFactors=FALSE)
tumour<-c("TCGA-E9-A1RD-01A","TCGA-E9-A1RC-01A")
normal<-c("TCGA-BH-A18P-11A","TCGA-BH-A18L-11A")
de_int<-SpidermiRanalyze_DEnetworkTCGA(data=miRNA_cN,
TCGAmatrix=Data_CANCER_normUQ_filt,
tumour,
normal
)
## ----------------------- DEA -------------------------------
## there are Cond1 type Normal in 2 samples
## there are Cond2 type Tumor in 2 samples
## there are 15243 features as miRNA or genes
## I Need about 2 seconds for this DEA. [Processing 30k elements /s]
## ----------------------- END DEA -------------------------------
## [1] "The number of differentially expressed genes between the two classes is 1303"
## [1] "The network containing DEGs with direct interactions consistes of 2 nodes and 1 edges "
SpidermiRvisualize
: To visualize the networkSpidermiRvisualize_mirnanet
: To Visualize the network.The user can visualize a 3D representation of the network in different colours for miRNA, genes, and pharmaco. The user can manage the network directly moving the nodes and the edges, in order to interpret the results in the graphic way.
library(networkD3)
SpidermiRvisualize_mirnanet(data=mir_pharmnet[sample(nrow(mir_pharmnet), 100), ] )
SpidermiRvisualize_BI
: To Visualize the network from a set of BI.Starting from a graphical representation of a network, the user can highlight with a different color specific BI.
biomark_of_interest<-c("hsa-let-7b","MUC1","PEX7","hsa-miR-222")
SpidermiRvisualize_BI(data=mir_pharmnet[sample(nrow(mir_pharmnet), 100), ],BI=biomark_of_interest)
SpidermiRvisualize_direction
: To visualize the networklibrary(visNetwork)
SpidermiRvisualize_direction(data=mir_pharmnet[sample(nrow(mir_pharmnet), 100), ] )
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).
SpidermiRvisualize_plot_target(data=miRNA_NET[1:15,])
## NULL
SpidermiRvisualize_degree_dist
: plots the degree distribution of the networkThis function plots the cumulative frequency distribution of degree centrality of a community.
SpidermiRvisualize_degree_dist(data=miRNA_NET)
SpidermiRvisualize_adj_matrix
: plots the adjacency matrix of the networkIt plots the adjacency matrix of the community, representing the degree of connections among the nodes.
SpidermiRvisualize_adj_matrix(data=miRNA_NET[1:30,])
SpidermiRvisualize_3Dbarplot
: 3D barplotIt plots a summary representation of the networks with the number of edges, nodes and miRNAs.
SpidermiRvisualize_3Dbarplot(Edges_1net=1041003,Edges_2net=100016,Edges_3net=3008,Edges_4net=1493,Edges_5net=1598,NODES_1net=16502,NODES_2net=13338,NODES_3net=1429,NODES_4net=675,NODES_5net=712,nmiRNAs_1net=0,nmiRNAs_2net=74,nmiRNAs_3net=0,nmiRNAs_4net=0,nmiRNAs_5net=37)
Features databases SpidermiR
:Features of databases integrated in SpidermiR
are:
SpidermiR:Application Examples
for more information.We studied in this case study the relationships between miRNAs with higher centrality and their target proteins with higher centrality. Please see SpidermiR:Application Examples
for more information.
a2<-Case_Study2_loading_1_network(species)
b2<-Case_Study2_loading_2_network(data=a2)
c2<-Case_Study2_loading_3_network(sdas=a2,miRNA_NET=b2)
Session Information ******
sessionInfo()
## R version 3.4.1 (2017-06-30)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 16.04.3 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.5-bioc/R/lib/libRblas.so
## LAPACK: /home/biocbuild/bbs-3.5-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] visNetwork_2.0.1 networkD3_0.4 SpidermiR_1.7.4
## [4] miRNAtap_1.10.0 AnnotationDbi_1.38.2 IRanges_2.10.3
## [7] S4Vectors_0.14.4 Biobase_2.36.2 BiocGenerics_0.22.0
## [10] BiocStyle_2.4.1
##
## loaded via a namespace (and not attached):
## [1] backports_1.1.0 circlize_0.4.1
## [3] aroma.light_3.6.0 plyr_1.8.4
## [5] igraph_1.1.2 selectr_0.3-1
## [7] ConsensusClusterPlus_1.40.0 lazyeval_0.2.0
## [9] splines_3.4.1 BiocParallel_1.10.1
## [11] GenomeInfoDb_1.12.2 ggplot2_2.2.1
## [13] digest_0.6.12 foreach_1.4.3
## [15] htmltools_0.3.6 viridis_0.4.0
## [17] gdata_2.18.0 magrittr_1.5
## [19] memoise_1.1.0 cluster_2.0.6
## [21] doParallel_1.0.10 limma_3.32.7
## [23] ComplexHeatmap_1.14.0 Biostrings_2.44.2
## [25] readr_1.1.1 annotate_1.54.0
## [27] matrixStats_0.52.2 R.utils_2.5.0
## [29] colorspace_1.3-2 rvest_0.3.2
## [31] ggrepel_0.6.5 blob_1.1.0
## [33] dplyr_0.7.3 crayon_1.3.4
## [35] RCurl_1.95-4.8 jsonlite_1.5
## [37] roxygen2_6.0.1 genefilter_1.58.1
## [39] bindr_0.1 zoo_1.8-0
## [41] survival_2.41-3 miRNAtap.db_0.99.10
## [43] iterators_1.0.8 glue_1.1.1
## [45] survminer_0.4.0 gtable_0.2.0
## [47] zlibbioc_1.22.0 XVector_0.16.0
## [49] GetoptLong_0.1.6 DelayedArray_0.2.7
## [51] kernlab_0.9-25 shape_1.4.3
## [53] prabclus_2.2-6 DEoptimR_1.0-8
## [55] scales_0.5.0 DESeq_1.28.0
## [57] mvtnorm_1.0-6 edgeR_3.18.1
## [59] DBI_0.7 ggthemes_3.4.0
## [61] Rcpp_0.12.12 cmprsk_2.2-7
## [63] viridisLite_0.2.0 xtable_1.8-2
## [65] foreign_0.8-69 matlab_1.0.2
## [67] bit_1.1-12 mclust_5.3
## [69] km.ci_0.5-2 sqldf_0.4-11
## [71] htmlwidgets_0.9 httr_1.3.1
## [73] gplots_3.0.1 RColorBrewer_1.1-2
## [75] fpc_2.1-10 modeltools_0.2-21
## [77] pkgconfig_2.0.1 XML_3.98-1.9
## [79] R.methodsS3_1.7.1 flexmix_2.3-14
## [81] nnet_7.3-12 locfit_1.5-9.1
## [83] labeling_0.3 reshape2_1.4.2
## [85] rlang_0.1.2 munsell_0.4.3
## [87] tools_3.4.1 downloader_0.4
## [89] gsubfn_0.6-6 RSQLite_2.0
## [91] broom_0.4.2 devtools_1.13.3
## [93] evaluate_0.10.1 stringr_1.2.0
## [95] yaml_2.1.14 org.Hs.eg.db_3.4.1
## [97] knitr_1.17 bit64_0.9-7
## [99] robustbase_0.92-7 survMisc_0.5.4
## [101] caTools_1.17.1 purrr_0.2.3
## [103] dendextend_1.5.2 bindrcpp_0.2
## [105] TCGAbiolinks_2.5.9 nlme_3.1-131
## [107] EDASeq_2.10.0 whisker_0.3-2
## [109] R.oo_1.21.0 xml2_1.1.1
## [111] biomaRt_2.32.1 compiler_3.4.1
## [113] curl_2.8.1 testthat_1.0.2
## [115] tibble_1.3.4 geneplotter_1.54.0
## [117] stringi_1.1.5 highr_0.6
## [119] GenomicFeatures_1.28.5 lattice_0.20-35
## [121] trimcluster_0.1-2 Matrix_1.2-11
## [123] psych_1.7.8 commonmark_1.4
## [125] KMsurv_0.1-5 GlobalOptions_0.0.12
## [127] data.table_1.10.4 bitops_1.0-6
## [129] rtracklayer_1.36.4 GenomicRanges_1.28.5
## [131] R6_2.2.2 latticeExtra_0.6-28
## [133] hwriter_1.3.2 ShortRead_1.34.1
## [135] KernSmooth_2.23-15 gridExtra_2.3
## [137] codetools_0.2-15 MASS_7.3-47
## [139] gtools_3.5.0 assertthat_0.2.0
## [141] chron_2.3-50 SummarizedExperiment_1.6.4
## [143] proto_1.0.0 rprojroot_1.2
## [145] rjson_0.2.15 withr_2.0.0
## [147] mnormt_1.5-5 GenomicAlignments_1.12.2
## [149] Rsamtools_1.28.0 GenomeInfoDbData_0.99.0
## [151] diptest_0.75-7 hms_0.3
## [153] grid_3.4.1 tidyr_0.7.1
## [155] class_7.3-14 rmarkdown_1.6
## [157] ggpubr_0.1.5
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.”
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.”
Rukov J, Wilentzik R, Jaffe I, Vinther J, Shomron N. 2013. “Pharmaco MiR Linking MicroRNAs and Drug Effects.”
Russo F, Di Bella S, Nigita G, Macca V, Lagana A, Giugno R, Pulvirenti A, Ferro A. 2012. “MiRandola Extracellular Circulating MicroRNAs Database.”
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.”