BACKGROUND
a major challenge in biomedical research is to understand the underlying mechanisms of human disease.
great effort has been spent on determining genes associated with human diseases.
however, most human diseases, and cancer in particular, cannot be attributed to single gene but arise due to complex interactions among multiple components of the cell, including genes, proteins, and mirnas
<cit> .
mirnas are a large family of gene regulators, found in both plants and animals, which impact gene activity by binding to the 3’utr of target mrnas leading to mrna degradation or translational inhibition
<cit> .
though mirnas are only 18- <dig> nucleotides, each can control the expression of hundreds of genes.
it is estimated that approximately half of the human genome is regulated by mirna-mediated gene control
<cit> .
mirnas play a key role in regulating diverse cellular functions, such as development, proliferation, apoptosis, and metabolism
<cit>  and are associated with a growing list of diseases including cancer
<cit> .
an increasing body of evidence suggests that mirnas impact gene expression in many cancer types including prostate cancer
<cit> .
several studies have investigated the role of mirnas in cancer using mrna and mirna expression profiling
<cit> .
better understanding the regulatory role of mirnas in cancer development and progression requires exploring their influence on other components of the cellular system they are a part of.
doing so, may lead to identifying predictive biomarkers and developing novel therapeutic strategies for cancer.
current major challenges in mirna research are prediction and experimental validation of mirna-target interactions, and determination of the functional role of mirnas.
computational prediction of mirnas is challenging in human genomes because of the imperfect pairing of the mirna with the corresponding target site
<cit> .
several factors can influence mirna-mediated gene control, like 3’utr length, number of mirna targets sites
<cit> , degree of complementary match
<cit> , amount of target mrna
<cit> , and the competition for targeted mrna
<cit> .
unfortunately, current sequence based predictions produce many false positive interactions and many of the predicted interactions may not be functional
<cit> , which means there may be no relationship between the expression levels of the mrna and the predicted targeting mirna.
several studies have tried to solve this by integrating gene expression data with sequence-based prediction to remove non-functional interactions and keep interactions that show negative correlation between mirna and their targets
<cit> .
thus, sequence-based methods provide a general view of the potential mirna targets but expression data or other cellular context information is required to more accurately predict mirna-target interactions.
determining the role of individual mirnas in cellular regulatory processes is still a major challenge.
the function of many mirnas remains unknown, and even for relatively well studied mirnas, only a handful of their targets have been characterized
<cit> .
delineating mirna function through knock-out and overexpression experiments in model organisms has had limited success, possibly because of functional redundancy among mirnas or among gene pathways regulated by mirnas
<cit> .
a mirna downregulates its targets, thus negative correlation in expression levels between a mirna and its direct targets indicates that the mirna is functional.
several studies have attempted to extract mirna-target modules based on the correlation between mirnas and targets
<cit>  and based on graph theory
<cit> .
however, these results are complicated by indirect effects - a single mirna may target many mrna targets that may influence other genes, thus negative correlation between mirna and targets does not indicate a direct interaction between mirna and target.
interactions between mirna and targets are not solely dependent on the 3’utr of the target, but depend on what other competing 3’utr targets are expressed in a given cellular context.
limited attempts have been made to investigate the impact of mirnas on protein interactors of the target.
it has been shown that protein-protein interaction  network topological features help to filter out false positive targets
<cit> , and help to prioritize mirnas in prostate cancer
<cit> .
recent evidence showed that some protein complexes are enriched with single mirna targets and some complexes are enriched with mirna cluster targets
<cit> .
for example, smad3-smad4-foxo <dig> complex is enriched with mir- <dig> targets, and mad1-sin3a-hdac <dig> complex is enriched with targets of the mir-510- <dig> and mir-1912- <dig> clusters.
other studies demonstrated that ppi context of mirna targets provides more representative information about mirna function compared to using only direct targets
<cit> .
direct targets of mirnas and their partners jointly showed higher modularity levels compared with mirna direct targets alone
<cit> .
analyzing properties of mirna targets is a promising approach to mirna function prediction.
mirpath
<cit>  is a computational tool developed to identify molecular pathways enriched in mirna targets set.
mirpath extracts mirna targets from other tools such as targetscan
<cit> , pita
<cit> , and then mirna function is predicted by assessing whether the predicted targets of a given mirna are enriched for particular functional annotations.
such enrichment based methods suffer from several limitations.
first, they solely depend on the mirna-target prediction algorithms that are noisy.
second, predicted mirna targets are usually large  and this leads to heterogeneous functional annotations that make it difficult to gain high confidence predictions.
integrating expression data is a promising approach to reduce noise in enrichment results.
the mirna body map
<cit>  is a web tool developed for mirna functional annotation in normal and diseased human tissues that integrates expression data to reduce heterogeneity in functional annotations.
fame
<cit>  is another tool with three main applications in the area of mirna functional analysis.
firstly, it infers mirna function directly using sets of genes sharing common annotations and secondly, infers mirna function indirectly using matched mrna/mirna expression data.
thirdly, fame predicts the function of genomic clusters of mirnas.
integrating the protein context of mirna targets is another promising dimension for mirna function prediction.
mirupnet
<cit>  is another mirna function prediction framework that predicts mirna function based on the upstream context of mirna and not downstream.
it infers the mirna function by functionally analyzing the context of its transcription factors in a protein-protein interaction network.
using information about tfs upstream of a mirna results in the discovery of additional biological processes not seen in mirna targets .
these observations shed light on the influence of mirnas on the ppi subnetwork involving the targets, and highlight the importance of considering target protein interactors when searching for functional mirna-target interactions.
in the post-genomics era, a crucial task in molecular biology is to understand mirna regulation in the context of biological networks.
since mirnas target proteins that are part of either protein complexes or signaling pathways, it is important to study the influence of mirnas on protein networks in disease progression.
characterizing the role of mirnas in the context of protein networks has emerged recently in several studies
<cit> .
by analyzing the interactions between mirnas and cellular signaling networks, mirnas were found to predominantly target proteins of the same signaling pathway and target highly connected scaffolds and most downstream network components such as signaling transcription factors.
mirnas were also found to less frequently target upstream components of the signaling pathways like membrane receptors and ligands
<cit> .
hsu et al <cit>  demonstrated that many mirna-targeted genes are hub proteins and bottleneck proteins in protein interaction networks  and thus have higher betweeness centrality.
when these hub or bottleneck proteins are repressed by individual or multiple mirnas, they may consequently influence large part of the interacting proteins and thus control key components of the ppin.
their analysis showed that the target proteins of individual mirnas tend to interact with more proteins than other non-mirna targets.
positive correlation between protein connectivity  and the number of mirnas targeting the corresponding protein has been observed by liang and li
<cit> .
this means that proteins with large numbers of partners in the ppin network need more mirnas to control their expression.
mirna induced influence can propagate in the regulatory network by targeting master transcription factors.
cui et al <cit>  found that 42% of  <dig> gene that are regulated by tfs, are mirna targets, and the average tf binding site count of mirna targets is significantly higher than that of non mirna targets.
this suggests that gene expression regulation by mirnas at the post-transcriptional level is coordinated with that of tfs at the transcriptional level and genes targeted by more mirnas have more tf binding sites.
in this work we introduce a new method to characterize mirna function based on its effect on the expression of the target and its neighbors in a functional interaction network.
unlike previous methods that weight mirna-target interactions based on sequence complementarity or gene expression correlation alone, we estimate the overall influence of a mirna on its target based on the target gene expression level and the gene expression levels of the interaction neighborhood of the target.
mirnas with high influence are validated using independent mirna expression datasets, and by analyzing the biological pathway enrichment of target protein modules.
we then used our mirna-target influence network to predict the overall influence of each mirna on individual prostate cancer patients to find those mirnas associated with aggressive cancer.
we show that mirnas with high influence on protein complexes and biological processes are likely involved in cancer progression and have potential prognostic significance.
methods
mirna targets
human mirna target predictions for mirna with conserved 3’utr were taken from targetscan  <dig>
<cit> , and experimentally validated mirna and their targets were taken from mirtarbase
<cit>  and mirecord
<cit> .
we used the union of mirtarbase and mirecord as a source of experimentally validated mirna-target interactions.
functional protein interaction  networks
we used combined undirected functional protein interactions  as described in
<cit> .
fpi includes annotated functional protein interactions from reactome, panther, cellmap, biocarta, kegg and tred, and includes interactions derived from physical protein interaction, co-expression data, domain-domain interaction data.
fpi was constructed using a naive bayes classifier  to distinguish high-likelihood fis from non-functional pairwise relationships.
we also used physical protein interactions from the hprd database
<cit> .
fpi functional interaction network includes hprd interactions, but the two networks have distinct topological features.
we also used another curated human signaling network from cui et al <cit> .
mirna and mrna expression data
we used mrna and mirna expression data from the mskcc prostate oncogenome project that is available at the gene expression omnibus
<cit> .
the data contains expression levels of  <dig> genes across  <dig> samples , and expression of  <dig> mirnas across  <dig> samples.
we used the expression data of  <dig> samples with both mrna and mirna data for our analysis.
to validate the mirna results we obtained using the taylor data, we used localized prostate cancer mirna expression data from independent prostate patient cohort  and prostate cell lines
<cit> .
mirna-target influence network construction
the initial mirna-target network is defined by a sequence based search ; seq =  <dig>  if t is a potential target for mirbased on the targetscan conserved sequence based prediction, otherwise seq =  <dig>  the relationship between the mirna expression level  and that of its  was computed using mutual information:
mimir,t=∑d∈mir∑r∈tplogppp
p is the joint probability density function of mir and t, and p and p are the marginal pdf’s of mirand t respectively.
we propose that the influence of a mirna  on its target depends on three variables.
first is the strength of the negative correlation between mirna and target expression profiles.
corrmir = mi, if seq =  <dig>  and corrmir =  <dig>  if seq =  <dig>  we only considered mir and t pairs with negative pearson’s correlation and seq =  <dig>  this step is needed to filter out mirna-target pairs with high mi due to positive correlation.
second is the direct impact of the mirna on the expression of the partners of the target.
we calculated mutual information  between the expression profiles of each target and its fpi partners , where corrfpi = mi  if ti is linked to tj in fpi, and corrfpi =  <dig> otherwise.
we used the maximum mi between the target and its partners to represent the direct influence of mirna on the target partners.
third is the indirect influence of mirna on the expression of the target through its partners.
the indirect impact of a mirna on its target through its partners is defined as w , where
w=∑kmi×mi where k is the partners of t in fpi.
we assess the influence of a mirna  on its potential target mrna  by integrating these three evidences.
figure
1a provides a schematic description of integrating gene expression and fpi to identify mirna influence on a protein target.
the mirti network was calculated by combining the three evidences of association between mirna and its targets as:
mirti=corrmirmir∗w∗max)
using the mirna-target influence network to measure the influence of mirnas on prostate cancer progression
we used the mirti network to predict mirnas with influence on gene expression profiles of prostate cancer samples .
we model the gene expression of prostate cancer samples as a linear combination of mirna effects on their targets
<cit> , as follows:
pcs=∑j=1mirmirtij∗βj+λpα
where
pα=∑j=1mir
is the elastic-net penalty.
pα is a compromise between the ridge regression penalty , and the lasso penalty .
in this model we set α to be  <dig>  as a middle value between  and  .
setting  will reduce sparsity achieved using lasso regression while still panelizing correlated predictors.
this penalty is particularly useful when there are many correlated predictor variables as in the case of mirnas.
we tried several values of α and we saw that when α changes from  <dig> to  <dig>  or  <dig>  it dramatically changes the minimum λ value and the number of non-zero values in the solution.
setting α to  <dig> produced large number of non-zero values in the solution  and α =  <dig> produced small number of non-zero values  that lead to a very sparse solution that might affect predictors of small effect on the response.
setting α to  <dig>  leads to a solution of medium sparsity with  <dig> non-zero elements .
however, changing α around  <dig>   does not change minimum λ value a lot  which does not impact number of non-zero elements.
we have selected α to be  <dig>  because when it is  <dig> , the minimum λ leads to the minimum mse.
β is the regression coefficient of each variable, which indicates how the expression level of each mirna can explain the gene expression profile of prostate cancer samples.
λ is a factor that determines the sparsity of the solution, as λ increases, the number of nonzero components of β decreases.
in this study, we selected  <dig> values of λ and used those that minimize the mean square error.
more details on λ optimization with respect to α is shown in additional file
1: figure s <dig>  elastic-net regression was fit using ten-fold cross validation.
we used glmnet package available at
http://www-stat.stanford.edu/ tibs/glmnet-matlab/ to solve the regression model.
for each patient we predicted the influence of the mirna set on the patient’s gene expression profile.
figure
1b describes the steps to construct the input of the model and its output.
the resulting mirna-patient influence profile was used to associate a mirna with a sample’s outcome.
detection of transcriptional activity centers in prostate cancer
several studies have shown that a functional interaction network provides information about the function of a gene of interest using the guilt by association concept
<cit> .
we define activity score for each gene based on the importance of the neighbor genes, similar to other studies
<cit> .
first, we computed the differential expression significance of each gene  in the prostate cancer gene expression data using the student’s t-test.
the bigger r is, the more significant the expression of the gene is.
we then used the prostate cancer corrfpi network to define the strength of the relationship between the gene of interest and its neighbors in the fpi network.
activityscore for gene is defined as:
activityscore=∗max∗r)
where n is the total number of neighbor genes.
using the mirna-target influence network to identify mirna influence on genes with high activity center scores
the activityscore measure reflects transcriptional activity of genes in modules rather than single genes.
here we model the expression activity of genes as a linear combination of the mirnas’ expression.
we used the mirti network that represents the influence of each mirna on each genes activity centers score to predict mirnas that explain the transcriptional activity center score using the same regression model as used above:
activityscore=∑j=1mirmirtij∗βj+λpα
the output of this model is a coefficient for each mirna that represents how each mirna explains the expression activity score of the genes.
RESULTS
global correlation between functional protein network topology and mirna regulation
to gain a global view of mirna regulation of the fpi network, we analyzed the relationship between the fpi network topological features and mirna regulation using two protein-protein interaction datasets  and two mirna-target networks  and mirecord .
we found a strong positive correlation between protein connectivity and the number of mirnas targeting the corresponding protein  and a negative correlation between protein clustering coefficient and number of targeting mirnas , across all networks we analyzed.
other network topological measures, like betweenness, did not show significant correlation.
this means that proteins with large numbers of partners in the fpi network need more mirnas to control their expression, and protein modules, such as complexes that are highly connected  need a smaller number of mirnas.
we also performed this analysis on randomly generated protein and mirna-target networks and did not observe any significant trend.
these results are in agreement with a recent studies by liang and li
<cit>  who found positive correlation between protein connectivity and average number of mirna target site types.
these findings motivated us to consider the protein network when estimating the influence of a mirna on its target.
spearman’s rank correlation of protein connectivity and clustering coefficient vs. number of targeting mirnas, significance of correlation is indicated by p-value in brackets, fdr <  <dig> .
prostate cancer mirnas target functionally associated genes
we manually collected a list of  <dig> mirnas that were experimentally validated to play a role in prostate cancer from several piblushed sources .
<dig> out of  <dig> have known targets as shown in table s <dig>  we also extracted  <dig> experimentally supported prostate mirnas from mir2disease
<cit>  and hmdd
<cit>  databases.
we next asked if the targets of the experimentally validated prostate mirnas are functionally related.
for each mirna, we predicted its targets using targetscan and then performed pathway enrichment analysis on the targets using the david online software .
enrichment analysis results  revealed that prostate related mirnas target proteins that are often in the same pathways and are functionally related.
figure
<dig> shows that the corresponding mirna targets are functionally related and connected in protein networks.
this motivated us to consider the target protein context to assess the influence of mirnas.
the fpi network helps reveal mirna-target modules that play a functional role in prostate cancer
we identified mirnas that are functionally active in prostate cancer, i.e.
they have a high influence on the protein partners of their targets using the mirna-target influence  network that was constructed based on prostate cancer expression data and fpi network.
we selected the top mirna-target interactions  with values greater than the upper quartile.
we observed that the targets of high influence mirnas are functionally related, which could explain, at a high-level, the mode of action of mirnas.
to investigate the functionality of the predicted mirna-target modules and further assess the robustness of our predictions, we followed two tracks.
first, the enriched pathways of mirna targets  were characterized.
enrichment analysis revealed that predicted targets are highly associated with focal adhesion and prostate cancer pathways .
as a control, we compared the predicted mirti network  that was generated using cytoscape
<cit>  with the corrmir network , which only represents correlation between mirnas and their target without considering the downstream regulation of mirnas on fpi target neighbors.
targets in corrmir also were enriched in focal adhesion and prostate cancer, but not as significantly .
we also randomly generated a functional mirna-target network by randomizing the fpi and the seq mirna-target networks and target did not show any significant enrichment.
in the second track, we extracted the expression profile of mirnas predicted to have high influence in two independent prostate cancer mirna expression datasets prostate cell line nci60
<cit> , and a prostate patient cohort .
here, it is worth mentioning that we do not extract differentially expressed genes from nci <dig> and gse <dig>  but rather we only extract the expression of our mirna list predicted to have high influence on targets.
in both datasets, mirnas predicted to have high influence on target interactors demonstrated diagnostic and prognostic significance.
mirnas predicted to have high influence based on our approach  are associated with prostate cancer better than those derived from the corrmir network and mirna-target network generated from the random fpi network based on multiple functional analysis strategies.
first, a large proportion of predicted high influence mirnas  are known to play a role in prostate cancer development .
several oncogene mirnas like mir- <dig>  mir <dig>  mir-125b, and mir-106b were identified in addition to several tumor suppressor mirnas like mir-34a, mir-20a, mir- <dig>  mir15a and mir-16
<cit> .
second, high influence mirnas are better able to accurately discriminate prostate cancer from normal samples.
mirnas with high-influence on protein complexes were able to classify patients in the taylor data
<cit>  into normal vs. cancer patients with 97% classification accuracy using a linear svm, better than prostate mirnas  that were extracted from the literature, which gave 92% classification accuracy.
randomly generated sets of mirnas of size  <dig> gave an average accuracy of 63%.
to further validate the robustness of the high influence mirnas, we extracted their expression profiles from two independent prostate mirna expression profiling studies  and showed that predicted high influence mirnas are accurate prostate cancer biomarkers.
both the high influence mirna and experimentally verified prostate mirna lists performed equally well on the nci <dig> data , but high-influence mirnas performed better than known prostate mirnas on gse <dig> .
randomly generated lists of the same size gave an average of 57% accuracy for nci <dig> and 52% for gse <dig>  third, predicted functional target modules of the high-impact mirnas are associated with multiple cancer pathways and prostate cancer related pathways, like tgf-b signaling pathway, and they are involved in several other cancers like glioma, melanoma and bladder .
also they were highly enriched with cell motion and cell migration go terms .
pathway enrichment map
<cit>  was used to show the map of enriched biological concepts.
this potentially indicates that these genes are important in prostate metastasis.
we further extracted  <dig> prognostic mirnas from the taylor data that are associated with aggressive prostate cancer by grouping prostate samples into aggressive prostate cancer  as one group and the other clusters as the other group.
the top  <dig> differentially expressed mirnas were extracted using sam.
<dig> out of our  <dig> high-influence mirnas were in common with aggressive mirnas from the taylor data .
<cit>  was used for network visualization.
since targetscan and fpi interaction data are noisy, we repeated the experiments using a highly curated human signaling network and curated mirna-target interactions from mirecord and mirtarbase.
the predicted interaction network  was found to be modular and partners of mirna targets are found in highly dense network regions.
<dig> mirnas were identified to have high influence on protein signaling network.
<dig> of them are in the list of  <dig> prostate mirnas  we collected from literature .
although the topological structure of the target modules are not very similar to modules in figure
<dig> due to difference in the protein networks and mirna-target interactions used to find them, the mirna targets’ partners are modular and form subnetworks of potential dysregulated proteins.
<dig> of the  <dig> mirnas predicted to have high influence on signaling network were in common with the  <dig> mirnas predicted to have high influence on the fpi network .
these results suggest that the completeness of protein interactions network plays a crucial role to identify high influence mirnas.
functional mirna-target modules are prognostic biomarkers that help identify patients with aggressive tumors
we investigated the relationship between the prostate cancer expression profile of mirnas and outcome using taylor expression data.
we used hierarchical clustering to group patients into two groups  based on the mirna expression profile across all samples.
our findings indicate that the predicted high influence mirnas are significantly associated with cancer recurrence .
we compared the prognostic effectiveness of the high impact mirnas with three mirna lists.
the first is mirnas with experimentally validated targets in prostate cancer  .
the second is differentially expressed mirnas between normal and primary prostate identified from taylor mirna expression data  , and the last list is aggressive prostate specific mirnas identified from taylor data  .
randomly generated lists of mirnas performed poorly to separate high vs. low risk patients .
the aggressive prostate specific mirnas are the most accurate mirnas to predict cancer recurrence; however, their mode of action is unclear.
the identified high influence mirnas, which have a predicted mode of action  based on our method, were effective at separating high vs. low risk patients .
we then used target expression from taylor data to find association between cancer recurrence and target expression, and we used target expression from the swedish prostate cohort  to find association between targets and cancer specific death.
results revealed that targets were associated with both recurrence  and cancer specific death .
unfortunately mirna expression is not available for the swedish cohort, so we could not conduct association between cancer specific death and mirnas.
altogether, the above results show that the mirna-target modules, which are the set of proteins targeted by one or more mirnas, are significant prognostic modules for prostate cancer.
<dig>  the expression of the functional mirnas was extracted from the taylor mirna data and hierarchical clustering was used to group patients into high and low risk categories.
results revealed that functional mirnas could be prognostic biomarkers.
patient specific mirna influence helps predict cancer recurrence
we next investigated if patient specific mirna influence has prognostic value for cancer recurrence.
using the mirti network and the gene expression profiles of patients as input to our regression model in eq., we predicted the influence of each mirna on the expression profile of each patient based on the influence of each mirna on its targets.
as a result, we obtained a patient specific mirna influence matrix.
we used this matrix  to identify clinically distinct patient groups.
we also used the seq network instead of mirti as input to the regression model for comparison.
we found that using mirti to predict patient specific mirna influence  is substantially more informative than using the binary seq network   in predicting cancer recurrence.
this further supports our results showing that considering protein interaction context to assess the clinical value of mirna is informative.
identification of mirnas of high-influence on principal regulators
we next used the mirti network to identify mirnas that can explain the activityscore of the target genes.
in the previous section, we used the mirti network to predict the gene expression profile of prostate cancer samples.
in this section, we propose another application for the mirti network to identify mirna influence on target gene expession.
integrating biological networks proved effective to help identify biomarkers that can explain most of the gene expression change
<cit> .
activityscore for each gene is calculated based on the observed effect  of its neighbors; genes that have high activityscore are defined as principal regulators.
we observed two sets of principal regulators; one with significant p and one with non-significant p .
the former set represents genes that are significantly differentially expressed  and their neighbors are also de.
these might be transcription factors that are affected in prostate cancer and they affect their target genes.
the latter set represents genes that are non−de but their neighbors are de.
these genes might be regulated at the post-translational level instead of at the transcriptional level.
for example, a change in the phosphorylation status of a gene might change its activity and thus affect the gene expression of its downstream genes.
the functional protein network of the latter set was enriched in zinc-finger proteins , which are known transcription factors.
figure
<dig> shows highly connected clusters of proteins of activityscore greater than the average score of all genes and have p-value greater than  <dig> .
these proteins are enriched in wnt and cadherin pathways and focal adhesion .
after calculating the activity score for each gene, we used the mirti network to predict the activityscore profile for all genes using our regression model in eq.. results revealed that mir- <dig>  mir- <dig>  mir- <dig>  mir-542-5p, mir- <dig>  mir- <dig>  and mir- <dig> are the mirnas that can positively explain the gene activity profile; this means that increasing expression level of mirnas will lead to increasing the transcription of activity centers.
mir- <dig> and mir- <dig> have been characterized as ongogenes
<cit>  and this supports the positive association between the two mirnas and activityscore.
mir- <dig> and mir-18b negatively explain the expression profile of the activityscore of genes which will have negative effect of activity centers.
this suggests that these mirnas might act as oncogenes or tumor suppressors.
since tumor heterogeneity affects the identification of robust cancer biomarkers, li et al.
<cit>  found that most cancer gene signatures are not robust and not reproducible.
thus they proposed a re-sampling based framework to identify robust cancer biomarkers.
in this work we asked the question whether re- sampling might have an effect on the activityscore profile.
to answer this question, we used significance analysis of microarray
<cit>  that is based on re-sampling to identify differentially expressed genes and then generate an activityscore profile using sam results.
we repeated sam analysis  <dig> times; each time we change the permutation number and generated an activityscore profile.
the resulting profiles demonstrated very significant correlation  which indicates that re-sampling does not affect the activityscore and that identified activity centers are robust and reproducible within our data set.
discussion
prostate cancer is one of the most commonly diagnosed malignant tumors in aged men in north america.
mirnas that are a family of regulatory molecules are significantly altered in prostate cancer
<cit> .
however, mirna’s mode of action and how the influence of prostate mirnas on target expression is involved in prostate cancer progression is not well known.
over- or under-expression of specific mirnas in different tumors makes them potential therapeutic targets and diagnostic or prognostic biomarkers; however, mirnas that are differentially expressed and influence their targets and target partners are important regulators and thus are more promising for diagnostics, prognostics or therapy.
in this work we use functional protein interactions to identify mirnas with high influence on targets and their partners.
we hypothesize that mirnas that influence a large number of interacting proteins are more important than those that only affect a few proteins.
we first showed that proteins that are highly connected have more regulating mirnas compared to those with low connectivity.
thus, identifying mirnas that regulate highly connected proteins is important to understand how to control propagation of gene expression changes via mirnas.
we showed that mirnas that have been experimentally verified to play a role in prostate cancer target functionally related genes.
this motivated us to investigate how mirnas that have high influence on protein partners of the target genes help us to better understand prostate cancer.
in this work we bridge a gap between systems biology and clinical biology by investigating the association between mirnas that have high influence on the system with the outcome of the system.
we built a mirna-target influence network  by following mirna influence of expression in prostate cancer of downstream genes in the fpi network and then proposed three applications of this network.
first, we used it to identify mirna target functional modules and complexes.
this revealed mirnas with high-influence on the target fpi neighborhood, which suggests that these mirna are important in prostate cancer.
the difference between high-influence mirnas and differentially expressed mirnas is that high-influence mirnas are differentially expressed and have differentially expressed targets and target interaction neighbors.
validating both mirna and targets in the functional modules against independent mirna expression datasets from prostate indicates that they are robust prostate cancer diagnostic biomarkers.
analyzing functional modules of mirna targets revealed several results.
first, target genes are enriched in prostate cancer and focal adhesion pathways, which may help explain the progression and metastasis process as our data includes metastatic samples.
functional modules are also of prognostic significance as they were associated with cancer recurrence and cancer specific death.
moreover, mirti network  revealed that some proteins like btbd <dig> ank <dig> col12a <dig> are highly repressed by several mirnas.
on the other hand, some mirnas  are highly influential on target partners as they regulate several connected proteins.
this suggests that mirnas have different mode of actions based on their influence on the expression of the target neighborhood.
this might help to define new regulatory classes of mirnas based on their mode of action.
the second application of mirti is to predict patient-specific mirna influence by using a regression model.
in this application we used the mirti network to predict the gene expression profile of the patients .
as a result of the regression model, we predict mirna-pcs network that shows how much each mirna explains the gene expression profile of a patient based on the weight with which it affects its targets.
we applied the regression model on all patients and generated a matrix that represents the influence of each mirna on each patient.
based on this mirna influence matrix we were able to group patients into aggressive and low risk cancer patients.
comparing the mirti with the seq network demonstrated that using mirna-target influence interactions gives more knowledge about mirna mode of action than using the binary seq weights that are based on only sequence predictions.
this result supports our initial conclusion that considering the downstream effect of mirna on protein partners of target is useful and has prognostic value.
we realized that both grouping patients based on mirna gene expression and based on patient-specific mirna influence from mirna-pcs network result in putting high risk patients in one group and low risk patients in the other group.
this indicates that the influence of each mirna on each patient is represented in the mrna expression of the patient.
the availability of differential mirna and mrna expression profiles from the same cancer samples enable functional analysis of mirnas in cancer, but there are few cancer cohorts that have expression levels of mirna and mrna from the same sample.
thus this result is very promising to predict the expression of mirnas in patients and predict their outcome without performing mirna expression profiling.
the third application of the mirti network is to predict mirnas with high-influence on genes with high activity center scores .
the activityscore profile of prostate cancer summarizes the activity of module proteins rather than the activity of single genes as in the second application.
here the mirti is used to predict the activityscore using the regression model.
the results emphasized the role of some mirnas already validated in prostate cancer , and identified novel mirnas like mir- <dig>  mir- <dig>  mir- <dig> and mir- <dig> that do not have a known mode of action in prostate cancer.
this means that these mirnas could be as important as the already validated mirnas, and could explain the summarized activity of the gene modules.
mirnas identified using the mirti and corrmir networks overlap; both networks identified mir- <dig> and mir- <dig> as important mirnas.
the advantage of using mirti over corrmir, seq and w  to identify mirna influence on target partners or on patient gene expression is that it produces two types of modules, unlike w  that favors the first type of modules and corrmir that favors the second type of module.
modules identified by our approach includes mirnas like mir- <dig> and mir- <dig> targeting highly interacting proteins, and mirnas like mir- <dig>  and mir- <dig> that target non-interacting complexes.
mirnas have been associated with clinical variables, prostate cancer recurrence and prostate cancer-specific death
<cit> .
however, the association between mirnas that target protein modules vs. clinical and survival data has not been well studied.
recent evidence showed that low mir- <dig> in human prostate tumors is associated with early disease recurrence
<cit> , and elevated levels of mir- <dig> is associated with high gleason score and higher risk of biochemical relapse
<cit> .
in this work we showed that mirnas identified using the mirti method are associated with cancer recurrence .
also, we showed that patient-specific mirna influences predicted using mirti are better prognostic biomarkers compared with binary, non-weighted mirna-target interactions.
this indicates that there is a link between the influence of mirna on target partners and its influence on outcome, but more analysis on larger cohorts and biological experiments are required to prove this result.
comparing the three applications of mirti revealed consistent results.
they all indicate the significant role of specific mirnas  in prostate cancer.
for instance, mir- <dig> and mir- <dig> are members of the same gene cluster and thus this supportes the effectiveness of integrating protein networks to identify mirnas with similar mode of action.
activityscore functional analysis indicates that zinc-finger proteins, zinc homeostasis, focal adhesion, and wnt signaling are enriched in genes with high activityscore .
evidence showed that zinc homeostasis is regulated by the mir-96-183- <dig> cluster.
this is in agreement with our results that demonstrate that mir- <dig> and mir- <dig> explain most of the genes activityscore that is significantly enriched in zinc homeostasis.
other predicted mirnas  may play a role in zinc homeostasis, focal adhesion, and cytoskeleton organization.
the large scale protein interactions and mirna target prediction data we used were useful to help elucidate the mechanistic role of mirnas in disease progression.
although the interaction datasets are far from complete and suffer from noise, our results were consistent across choice of ppi network.
using additional protein interaction networks, different mirna target prediction algorithms, and different expression data sets will likely reveal more mirnas with high-influence on cancer progression.
another future direction for this work is designing a systematic method to combine the three variables that determine the influence of mirnas on the target partners.
finally, this study on bridging the gap between clinical bioinformatics and network-based biomarkers provides clear evidence that protein interaction information is useful to identify diagnostic and prognostic cancer biomarkers, and to ameliorate the understanding of the functional mechanisms of mirnas.
CONCLUSIONS
we have developed a novel method to identify active mirna-target modules relevant to prostate cancer progression and outcome.
mirnas with high influence on protein networks are valuable biomarkers that can be used in clinical investigations for prostate cancer treatment.
combining the effects of mirnas on targets and target partners provides better understanding of mirnas function.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
ma, gb and ag designed the study.
ma wrote the initial manuscript.
gb edited the manuscript.
qm and ra proof-read the manuscript.
all authors contributed to the writing and approved the final manuscript.
supplementary material
additional file 1
experimentally validated prostate mirnas and supplementary figures.
this file lists mirnas that have been studied in prostate cancer and showed to play a role in prostate cancer progression.
it also shows some targets of these mirnas that have been validated in prostate cancer context.
click here for file
acknowledgements
m. alshalalfa and r. alhajj would like to thank nserc for funding.

 a major challenge in biomedical research is to understand the underlying mechanisms of human disease. great effort has been spent on determining genes associated with human diseases. however, most human diseases, and cancer in particular, cannot be attributed to single gene but arise due to complex interactions among multiple components of the cell, including genes, proteins, and mirnas. mirnas are a large family of gene regulators, found in both plants and animals, which impact gene activity by binding to the 3’utr of target mrnas leading to mrna degradation or translational inhibition. though mirnas are only 18- nucleotides, each can control the expression of hundreds of genes. it is estimated that approximately half of the human genome is regulated by mirna-mediated gene control. mirnas play a key role in regulating diverse cellular functions, such as development, proliferation, apoptosis, and metabolism and are associated with a growing list of diseases including cancer. an increasing body of evidence suggests that mirnas impact gene expression in many cancer types including prostate cancer. several studies have investigated the role of mirnas in cancer using mrna and mirna expression profiling. better understanding the regulatory role of mirnas in cancer development and progression requires exploring their influence on other components of the cellular system they are a part of. doing so, may lead to identifying predictive biomarkers and developing novel therapeutic strategies for cancer. current major challenges in mirna research are prediction and experimental validation of mirna-target interactions, and determination of the functional role of mirnas. computational prediction of mirnas is challenging in human genomes because of the imperfect pairing of the mirna with the corresponding target site. several factors can influence mirna-mediated gene control, like 3’utr length, number of mirna targets sites, degree of complementary match, amount of target mrna, and the competition for targeted mrna. unfortunately, current sequence based predictions produce many false positive interactions and many of the predicted interactions may not be functional, which means there may be no relationship between the expression levels of the mrna and the predicted targeting mirna. several studies have tried to solve this by integrating gene expression data with sequence-based prediction to remove non-functional interactions and keep interactions that show negative correlation between mirna and their targets. thus, sequence-based methods provide a general view of the potential mirna targets but expression data or other cellular context information is required to more accurately predict mirna-target interactions. determining the role of individual mirnas in cellular regulatory processes is still a major challenge. the function of many mirnas remains unknown, and even for relatively well studied mirnas, only a handful of their targets have been characterized. delineating mirna function through knock-out and overexpression experiments in model organisms has had limited success, possibly because of functional redundancy among mirnas or among gene pathways regulated by mirnas. a mirna downregulates its targets, thus negative correlation in expression levels between a mirna and its direct targets indicates that the mirna is functional. several studies have attempted to extract mirna-target modules based on the correlation between mirnas and targets and based on graph theory. however, these results are complicated by indirect effects - a single mirna may target many mrna targets that may influence other genes, thus negative correlation between mirna and targets does not indicate a direct interaction between mirna and target. interactions between mirna and targets are not solely dependent on the 3’utr of the target, but depend on what other competing 3’utr targets are expressed in a given cellular context. limited attempts have been made to investigate the impact of mirnas on protein interactors of the target. it has been shown that protein-protein interaction network topological features help to filter out false positive targets, and help to prioritize mirnas in prostate cancer. recent evidence showed that some protein complexes are enriched with single mirna targets and some complexes are enriched with mirna cluster targets. for example, smad3-smad4-foxo complex is enriched with mir- targets, and mad1-sin3a-hdac complex is enriched with targets of the mir-510- and mir-1912- clusters. other studies demonstrated that ppi context of mirna targets provides more representative information about mirna function compared to using only direct targets. direct targets of mirnas and their partners jointly showed higher modularity levels compared with mirna direct targets alone. analyzing properties of mirna targets is a promising approach to mirna function prediction. mirpath is a computational tool developed to identify molecular pathways enriched in mirna targets set. mirpath extracts mirna targets from other tools such as targetscan, pita, and then mirna function is predicted by assessing whether the predicted targets of a given mirna are enriched for particular functional annotations. such enrichment based methods suffer from several limitations. first, they solely depend on the mirna-target prediction algorithms that are noisy. second, predicted mirna targets are usually large and this leads to heterogeneous functional annotations that make it difficult to gain high confidence predictions. integrating expression data is a promising approach to reduce noise in enrichment results. the mirna body map is a web tool developed for mirna functional annotation in normal and diseased human tissues that integrates expression data to reduce heterogeneity in functional annotations. fame is another tool with three main applications in the area of mirna functional analysis. firstly, it infers mirna function directly using sets of genes sharing common annotations and secondly, infers mirna function indirectly using matched mrna/mirna expression data. thirdly, fame predicts the function of genomic clusters of mirnas. integrating the protein context of mirna targets is another promising dimension for mirna function prediction. mirupnet is another mirna function prediction framework that predicts mirna function based on the upstream context of mirna and not downstream. it infers the mirna function by functionally analyzing the context of its transcription factors in a protein-protein interaction network. using information about tfs upstream of a mirna results in the discovery of additional biological processes not seen in mirna targets. these observations shed light on the influence of mirnas on the ppi subnetwork involving the targets, and highlight the importance of considering target protein interactors when searching for functional mirna-target interactions. in the post-genomics era, a crucial task in molecular biology is to understand mirna regulation in the context of biological networks. since mirnas target proteins that are part of either protein complexes or signaling pathways, it is important to study the influence of mirnas on protein networks in disease progression. characterizing the role of mirnas in the context of protein networks has emerged recently in several studies. by analyzing the interactions between mirnas and cellular signaling networks, mirnas were found to predominantly target proteins of the same signaling pathway and target highly connected scaffolds and most downstream network components such as signaling transcription factors. mirnas were also found to less frequently target upstream components of the signaling pathways like membrane receptors and ligands. hsu et al demonstrated that many mirna-targeted genes are hub proteins and bottleneck proteins in protein interaction networks and thus have higher betweeness centrality. when these hub or bottleneck proteins are repressed by individual or multiple mirnas, they may consequently influence large part of the interacting proteins and thus control key components of the ppin. their analysis showed that the target proteins of individual mirnas tend to interact with more proteins than other non-mirna targets. positive correlation between protein connectivity and the number of mirnas targeting the corresponding protein has been observed by liang and li. this means that proteins with large numbers of partners in the ppin network need more mirnas to control their expression. mirna induced influence can propagate in the regulatory network by targeting master transcription factors. cui et al found that 42% of gene that are regulated by tfs, are mirna targets, and the average tf binding site count of mirna targets is significantly higher than that of non mirna targets. this suggests that gene expression regulation by mirnas at the post-transcriptional level is coordinated with that of tfs at the transcriptional level and genes targeted by more mirnas have more tf binding sites. in this work we introduce a new method to characterize mirna function based on its effect on the expression of the target and its neighbors in a functional interaction network. unlike previous methods that weight mirna-target interactions based on sequence complementarity or gene expression correlation alone, we estimate the overall influence of a mirna on its target based on the target gene expression level and the gene expression levels of the interaction neighborhood of the target. mirnas with high influence are validated using independent mirna expression datasets, and by analyzing the biological pathway enrichment of target protein modules. we then used our mirna-target influence network to predict the overall influence of each mirna on individual prostate cancer patients to find those mirnas associated with aggressive cancer. we show that mirnas with high influence on protein complexes and biological processes are likely involved in cancer progression and have potential prognostic significance. methods mirna targets human mirna target predictions for mirna with conserved 3’utr were taken from targetscan, and experimentally validated mirna and their targets were taken from mirtarbase and mirecord. we used the union of mirtarbase and mirecord as a source of experimentally validated mirna-target interactions. functional protein interaction networks we used combined undirected functional protein interactions as described in. fpi includes annotated functional protein interactions from reactome, panther, cellmap, biocarta, kegg and tred, and includes interactions derived from physical protein interaction, co-expression data, domain-domain interaction data. fpi was constructed using a naive bayes classifier to distinguish high-likelihood fis from non-functional pairwise relationships. we also used physical protein interactions from the hprd database. fpi functional interaction network includes hprd interactions, but the two networks have distinct topological features. we also used another curated human signaling network from cui et al. mirna and mrna expression data we used mrna and mirna expression data from the mskcc prostate oncogenome project that is available at the gene expression omnibus. the data contains expression levels of genes across samples, and expression of mirnas across samples. we used the expression data of samples with both mrna and mirna data for our analysis. to validate the mirna results we obtained using the taylor data, we used localized prostate cancer mirna expression data from independent prostate patient cohort and prostate cell lines. mirna-target influence network construction the initial mirna-target network is defined by a sequence based search; seq = if t is a potential target for mirbased on the targetscan conserved sequence based prediction, otherwise seq = the relationship between the mirna expression level and that of its was computed using mutual information: mimir,t=∑d∈mir∑r∈tplogppp p is the joint probability density function of mir and t, and p and p are the marginal pdf’s of mirand t respectively. we propose that the influence of a mirna on its target depends on three variables. first is the strength of the negative correlation between mirna and target expression profiles. corrmir = mi, if seq = and corrmir = if seq = we only considered mir and t pairs with negative pearson’s correlation and seq = this step is needed to filter out mirna-target pairs with high mi due to positive correlation. second is the direct impact of the mirna on the expression of the partners of the target. we calculated mutual information between the expression profiles of each target and its fpi partners, where corrfpi = mi if ti is linked to tj in fpi, and corrfpi = otherwise. we used the maximum mi between the target and its partners to represent the direct influence of mirna on the target partners. third is the indirect influence of mirna on the expression of the target through its partners. the indirect impact of a mirna on its target through its partners is defined as w, where w=∑kmi×mi where k is the partners of t in fpi. we assess the influence of a mirna on its potential target mrna by integrating these three evidences. figure 1a provides a schematic description of integrating gene expression and fpi to identify mirna influence on a protein target. the mirti network was calculated by combining the three evidences of association between mirna and its targets as: mirti=corrmirmir∗w∗max) using the mirna-target influence network to measure the influence of mirnas on prostate cancer progression we used the mirti network to predict mirnas with influence on gene expression profiles of prostate cancer samples. we model the gene expression of prostate cancer samples as a linear combination of mirna effects on their targets, as follows: pcs=∑j=1mirmirtij∗βj+λpα where pα=∑j=1mir is the elastic-net penalty. pα is a compromise between the ridge regression penalty, and the lasso penalty. in this model we set α to be as a middle value between and. setting will reduce sparsity achieved using lasso regression while still panelizing correlated predictors. this penalty is particularly useful when there are many correlated predictor variables as in the case of mirnas. we tried several values of α and we saw that when α changes from to or it dramatically changes the minimum λ value and the number of non-zero values in the solution. setting α to produced large number of non-zero values in the solution and α = produced small number of non-zero values that lead to a very sparse solution that might affect predictors of small effect on the response. setting α to leads to a solution of medium sparsity with non-zero elements. however, changing α around does not change minimum λ value a lot which does not impact number of non-zero elements. we have selected α to be because when it is, the minimum λ leads to the minimum mse. β is the regression coefficient of each variable, which indicates how the expression level of each mirna can explain the gene expression profile of prostate cancer samples. λ is a factor that determines the sparsity of the solution, as λ increases, the number of nonzero components of β decreases. in this study, we selected values of λ and used those that minimize the mean square error. more details on λ optimization with respect to α is shown in additional file 1: figure s elastic-net regression was fit using ten-fold cross validation. we used glmnet package available at http://www-stat.stanford.edu/ tibs/glmnet-matlab/ to solve the regression model. for each patient we predicted the influence of the mirna set on the patient’s gene expression profile. figure 1b describes the steps to construct the input of the model and its output. the resulting mirna-patient influence profile was used to associate a mirna with a sample’s outcome. detection of transcriptional activity centers in prostate cancer several studies have shown that a functional interaction network provides information about the function of a gene of interest using the guilt by association concept. we define activity score for each gene based on the importance of the neighbor genes, similar to other studies. first, we computed the differential expression significance of each gene in the prostate cancer gene expression data using the student’s t-test. the bigger r is, the more significant the expression of the gene is. we then used the prostate cancer corrfpi network to define the strength of the relationship between the gene of interest and its neighbors in the fpi network. activityscore for gene is defined as: activityscore=∗max∗r) where n is the total number of neighbor genes. using the mirna-target influence network to identify mirna influence on genes with high activity center scores the activityscore measure reflects transcriptional activity of genes in modules rather than single genes. here we model the expression activity of genes as a linear combination of the mirnas’ expression. we used the mirti network that represents the influence of each mirna on each genes activity centers score to predict mirnas that explain the transcriptional activity center score using the same regression model as used above: activityscore=∑j=1mirmirtij∗βj+λpα the output of this model is a coefficient for each mirna that represents how each mirna explains the expression activity score of the genes.  global correlation between functional protein network topology and mirna regulation to gain a global view of mirna regulation of the fpi network, we analyzed the relationship between the fpi network topological features and mirna regulation using two protein-protein interaction datasets and two mirna-target networks and mirecord. we found a strong positive correlation between protein connectivity and the number of mirnas targeting the corresponding protein and a negative correlation between protein clustering coefficient and number of targeting mirnas, across all networks we analyzed. other network topological measures, like betweenness, did not show significant correlation. this means that proteins with large numbers of partners in the fpi network need more mirnas to control their expression, and protein modules, such as complexes that are highly connected need a smaller number of mirnas. we also performed this analysis on randomly generated protein and mirna-target networks and did not observe any significant trend. these results are in agreement with a recent studies by liang and li who found positive correlation between protein connectivity and average number of mirna target site types. these findings motivated us to consider the protein network when estimating the influence of a mirna on its target. spearman’s rank correlation of protein connectivity and clustering coefficient vs. number of targeting mirnas, significance of correlation is indicated by p-value in brackets, fdr <. prostate cancer mirnas target functionally associated genes we manually collected a list of mirnas that were experimentally validated to play a role in prostate cancer from several piblushed sources. out of have known targets as shown in table s we also extracted experimentally supported prostate mirnas from mir2disease and hmdd databases. we next asked if the targets of the experimentally validated prostate mirnas are functionally related. for each mirna, we predicted its targets using targetscan and then performed pathway enrichment analysis on the targets using the david online software. enrichment analysis results revealed that prostate related mirnas target proteins that are often in the same pathways and are functionally related. figure shows that the corresponding mirna targets are functionally related and connected in protein networks. this motivated us to consider the target protein context to assess the influence of mirnas. the fpi network helps reveal mirna-target modules that play a functional role in prostate cancer we identified mirnas that are functionally active in prostate cancer, i.e. they have a high influence on the protein partners of their targets using the mirna-target influence network that was constructed based on prostate cancer expression data and fpi network. we selected the top mirna-target interactions with values greater than the upper quartile. we observed that the targets of high influence mirnas are functionally related, which could explain, at a high-level, the mode of action of mirnas. to investigate the functionality of the predicted mirna-target modules and further assess the robustness of our predictions, we followed two tracks. first, the enriched pathways of mirna targets were characterized. enrichment analysis revealed that predicted targets are highly associated with focal adhesion and prostate cancer pathways. as a control, we compared the predicted mirti network that was generated using cytoscape with the corrmir network, which only represents correlation between mirnas and their target without considering the downstream regulation of mirnas on fpi target neighbors. targets in corrmir also were enriched in focal adhesion and prostate cancer, but not as significantly. we also randomly generated a functional mirna-target network by randomizing the fpi and the seq mirna-target networks and target did not show any significant enrichment. in the second track, we extracted the expression profile of mirnas predicted to have high influence in two independent prostate cancer mirna expression datasets prostate cell line nci60, and a prostate patient cohort. here, it is worth mentioning that we do not extract differentially expressed genes from nci and gse but rather we only extract the expression of our mirna list predicted to have high influence on targets. in both datasets, mirnas predicted to have high influence on target interactors demonstrated diagnostic and prognostic significance. mirnas predicted to have high influence based on our approach are associated with prostate cancer better than those derived from the corrmir network and mirna-target network generated from the random fpi network based on multiple functional analysis strategies. first, a large proportion of predicted high influence mirnas are known to play a role in prostate cancer development. several oncogene mirnas like mir- mir mir-125b, and mir-106b were identified in addition to several tumor suppressor mirnas like mir-34a, mir-20a, mir- mir15a and mir-16. second, high influence mirnas are better able to accurately discriminate prostate cancer from normal samples. mirnas with high-influence on protein complexes were able to classify patients in the taylor data into normal vs. cancer patients with 97% classification accuracy using a linear svm, better than prostate mirnas that were extracted from the literature, which gave 92% classification accuracy. randomly generated sets of mirnas of size gave an average accuracy of 63%. to further validate the robustness of the high influence mirnas, we extracted their expression profiles from two independent prostate mirna expression profiling studies and showed that predicted high influence mirnas are accurate prostate cancer biomarkers. both the high influence mirna and experimentally verified prostate mirna lists performed equally well on the nci data, but high-influence mirnas performed better than known prostate mirnas on gse. randomly generated lists of the same size gave an average of 57% accuracy for nci and 52% for gse third, predicted functional target modules of the high-impact mirnas are associated with multiple cancer pathways and prostate cancer related pathways, like tgf-b signaling pathway, and they are involved in several other cancers like glioma, melanoma and bladder. also they were highly enriched with cell motion and cell migration go terms. pathway enrichment map was used to show the map of enriched biological concepts. this potentially indicates that these genes are important in prostate metastasis. we further extracted prognostic mirnas from the taylor data that are associated with aggressive prostate cancer by grouping prostate samples into aggressive prostate cancer as one group and the other clusters as the other group. the top differentially expressed mirnas were extracted using sam. out of our high-influence mirnas were in common with aggressive mirnas from the taylor data. was used for network visualization. since targetscan and fpi interaction data are noisy, we repeated the experiments using a highly curated human signaling network and curated mirna-target interactions from mirecord and mirtarbase. the predicted interaction network was found to be modular and partners of mirna targets are found in highly dense network regions. mirnas were identified to have high influence on protein signaling network. of them are in the list of prostate mirnas we collected from literature. although the topological structure of the target modules are not very similar to modules in figure due to difference in the protein networks and mirna-target interactions used to find them, the mirna targets’ partners are modular and form subnetworks of potential dysregulated proteins. of the mirnas predicted to have high influence on signaling network were in common with the mirnas predicted to have high influence on the fpi network. these results suggest that the completeness of protein interactions network plays a crucial role to identify high influence mirnas. functional mirna-target modules are prognostic biomarkers that help identify patients with aggressive tumors we investigated the relationship between the prostate cancer expression profile of mirnas and outcome using taylor expression data. we used hierarchical clustering to group patients into two groups based on the mirna expression profile across all samples. our findings indicate that the predicted high influence mirnas are significantly associated with cancer recurrence. we compared the prognostic effectiveness of the high impact mirnas with three mirna lists. the first is mirnas with experimentally validated targets in prostate cancer. the second is differentially expressed mirnas between normal and primary prostate identified from taylor mirna expression data, and the last list is aggressive prostate specific mirnas identified from taylor data. randomly generated lists of mirnas performed poorly to separate high vs. low risk patients. the aggressive prostate specific mirnas are the most accurate mirnas to predict cancer recurrence; however, their mode of action is unclear. the identified high influence mirnas, which have a predicted mode of action based on our method, were effective at separating high vs. low risk patients. we then used target expression from taylor data to find association between cancer recurrence and target expression, and we used target expression from the swedish prostate cohort to find association between targets and cancer specific death. results revealed that targets were associated with both recurrence and cancer specific death. unfortunately mirna expression is not available for the swedish cohort, so we could not conduct association between cancer specific death and mirnas. altogether, the above results show that the mirna-target modules, which are the set of proteins targeted by one or more mirnas, are significant prognostic modules for prostate cancer. the expression of the functional mirnas was extracted from the taylor mirna data and hierarchical clustering was used to group patients into high and low risk categories. results revealed that functional mirnas could be prognostic biomarkers. patient specific mirna influence helps predict cancer recurrence we next investigated if patient specific mirna influence has prognostic value for cancer recurrence. using the mirti network and the gene expression profiles of patients as input to our regression model in eq., we predicted the influence of each mirna on the expression profile of each patient based on the influence of each mirna on its targets. as a result, we obtained a patient specific mirna influence matrix. we used this matrix to identify clinically distinct patient groups. we also used the seq network instead of mirti as input to the regression model for comparison. we found that using mirti to predict patient specific mirna influence is substantially more informative than using the binary seq network in predicting cancer recurrence. this further supports our results showing that considering protein interaction context to assess the clinical value of mirna is informative. identification of mirnas of high-influence on principal regulators we next used the mirti network to identify mirnas that can explain the activityscore of the target genes. in the previous section, we used the mirti network to predict the gene expression profile of prostate cancer samples. in this section, we propose another application for the mirti network to identify mirna influence on target gene expession. integrating biological networks proved effective to help identify biomarkers that can explain most of the gene expression change. activityscore for each gene is calculated based on the observed effect of its neighbors; genes that have high activityscore are defined as principal regulators. we observed two sets of principal regulators; one with significant p and one with non-significant p. the former set represents genes that are significantly differentially expressed and their neighbors are also de. these might be transcription factors that are affected in prostate cancer and they affect their target genes. the latter set represents genes that are non−de but their neighbors are de. these genes might be regulated at the post-translational level instead of at the transcriptional level. for example, a change in the phosphorylation status of a gene might change its activity and thus affect the gene expression of its downstream genes. the functional protein network of the latter set was enriched in zinc-finger proteins, which are known transcription factors. figure shows highly connected clusters of proteins of activityscore greater than the average score of all genes and have p-value greater than. these proteins are enriched in wnt and cadherin pathways and focal adhesion. after calculating the activity score for each gene, we used the mirti network to predict the activityscore profile for all genes using our regression model in eq.. results revealed that mir- mir- mir- mir-542-5p, mir- mir- and mir- are the mirnas that can positively explain the gene activity profile; this means that increasing expression level of mirnas will lead to increasing the transcription of activity centers. mir- and mir- have been characterized as ongogenes and this supports the positive association between the two mirnas and activityscore. mir- and mir-18b negatively explain the expression profile of the activityscore of genes which will have negative effect of activity centers. this suggests that these mirnas might act as oncogenes or tumor suppressors. since tumor heterogeneity affects the identification of robust cancer biomarkers, li et al. found that most cancer gene signatures are not robust and not reproducible. thus they proposed a re-sampling based framework to identify robust cancer biomarkers. in this work we asked the question whether re- sampling might have an effect on the activityscore profile. to answer this question, we used significance analysis of microarray that is based on re-sampling to identify differentially expressed genes and then generate an activityscore profile using sam results. we repeated sam analysis times; each time we change the permutation number and generated an activityscore profile. the resulting profiles demonstrated very significant correlation which indicates that re-sampling does not affect the activityscore and that identified activity centers are robust and reproducible within our data set. discussion prostate cancer is one of the most commonly diagnosed malignant tumors in aged men in north america. mirnas that are a family of regulatory molecules are significantly altered in prostate cancer. however, mirna’s mode of action and how the influence of prostate mirnas on target expression is involved in prostate cancer progression is not well known. over- or under-expression of specific mirnas in different tumors makes them potential therapeutic targets and diagnostic or prognostic biomarkers; however, mirnas that are differentially expressed and influence their targets and target partners are important regulators and thus are more promising for diagnostics, prognostics or therapy. in this work we use functional protein interactions to identify mirnas with high influence on targets and their partners. we hypothesize that mirnas that influence a large number of interacting proteins are more important than those that only affect a few proteins. we first showed that proteins that are highly connected have more regulating mirnas compared to those with low connectivity. thus, identifying mirnas that regulate highly connected proteins is important to understand how to control propagation of gene expression changes via mirnas. we showed that mirnas that have been experimentally verified to play a role in prostate cancer target functionally related genes. this motivated us to investigate how mirnas that have high influence on protein partners of the target genes help us to better understand prostate cancer. in this work we bridge a gap between systems biology and clinical biology by investigating the association between mirnas that have high influence on the system with the outcome of the system. we built a mirna-target influence network by following mirna influence of expression in prostate cancer of downstream genes in the fpi network and then proposed three applications of this network. first, we used it to identify mirna target functional modules and complexes. this revealed mirnas with high-influence on the target fpi neighborhood, which suggests that these mirna are important in prostate cancer. the difference between high-influence mirnas and differentially expressed mirnas is that high-influence mirnas are differentially expressed and have differentially expressed targets and target interaction neighbors. validating both mirna and targets in the functional modules against independent mirna expression datasets from prostate indicates that they are robust prostate cancer diagnostic biomarkers. analyzing functional modules of mirna targets revealed several results. first, target genes are enriched in prostate cancer and focal adhesion pathways, which may help explain the progression and metastasis process as our data includes metastatic samples. functional modules are also of prognostic significance as they were associated with cancer recurrence and cancer specific death. moreover, mirti network revealed that some proteins like btbd ank col12a are highly repressed by several mirnas. on the other hand, some mirnas are highly influential on target partners as they regulate several connected proteins. this suggests that mirnas have different mode of actions based on their influence on the expression of the target neighborhood. this might help to define new regulatory classes of mirnas based on their mode of action. the second application of mirti is to predict patient-specific mirna influence by using a regression model. in this application we used the mirti network to predict the gene expression profile of the patients. as a result of the regression model, we predict mirna-pcs network that shows how much each mirna explains the gene expression profile of a patient based on the weight with which it affects its targets. we applied the regression model on all patients and generated a matrix that represents the influence of each mirna on each patient. based on this mirna influence matrix we were able to group patients into aggressive and low risk cancer patients. comparing the mirti with the seq network demonstrated that using mirna-target influence interactions gives more knowledge about mirna mode of action than using the binary seq weights that are based on only sequence predictions. this result supports our initial conclusion that considering the downstream effect of mirna on protein partners of target is useful and has prognostic value. we realized that both grouping patients based on mirna gene expression and based on patient-specific mirna influence from mirna-pcs network result in putting high risk patients in one group and low risk patients in the other group. this indicates that the influence of each mirna on each patient is represented in the mrna expression of the patient. the availability of differential mirna and mrna expression profiles from the same cancer samples enable functional analysis of mirnas in cancer, but there are few cancer cohorts that have expression levels of mirna and mrna from the same sample. thus this result is very promising to predict the expression of mirnas in patients and predict their outcome without performing mirna expression profiling. the third application of the mirti network is to predict mirnas with high-influence on genes with high activity center scores. the activityscore profile of prostate cancer summarizes the activity of module proteins rather than the activity of single genes as in the second application. here the mirti is used to predict the activityscore using the regression model. the results emphasized the role of some mirnas already validated in prostate cancer, and identified novel mirnas like mir- mir- mir- and mir- that do not have a known mode of action in prostate cancer. this means that these mirnas could be as important as the already validated mirnas, and could explain the summarized activity of the gene modules. mirnas identified using the mirti and corrmir networks overlap; both networks identified mir- and mir- as important mirnas. the advantage of using mirti over corrmir, seq and w to identify mirna influence on target partners or on patient gene expression is that it produces two types of modules, unlike w that favors the first type of modules and corrmir that favors the second type of module. modules identified by our approach includes mirnas like mir- and mir- targeting highly interacting proteins, and mirnas like mir- and mir- that target non-interacting complexes. mirnas have been associated with clinical variables, prostate cancer recurrence and prostate cancer-specific death. however, the association between mirnas that target protein modules vs. clinical and survival data has not been well studied. recent evidence showed that low mir- in human prostate tumors is associated with early disease recurrence, and elevated levels of mir- is associated with high gleason score and higher risk of biochemical relapse. in this work we showed that mirnas identified using the mirti method are associated with cancer recurrence. also, we showed that patient-specific mirna influences predicted using mirti are better prognostic biomarkers compared with binary, non-weighted mirna-target interactions. this indicates that there is a link between the influence of mirna on target partners and its influence on outcome, but more analysis on larger cohorts and biological experiments are required to prove this result. comparing the three applications of mirti revealed consistent results. they all indicate the significant role of specific mirnas in prostate cancer. for instance, mir- and mir- are members of the same gene cluster and thus this supportes the effectiveness of integrating protein networks to identify mirnas with similar mode of action. activityscore functional analysis indicates that zinc-finger proteins, zinc homeostasis, focal adhesion, and wnt signaling are enriched in genes with high activityscore. evidence showed that zinc homeostasis is regulated by the mir-96-183- cluster. this is in agreement with our results that demonstrate that mir- and mir- explain most of the genes activityscore that is significantly enriched in zinc homeostasis. other predicted mirnas may play a role in zinc homeostasis, focal adhesion, and cytoskeleton organization. the large scale protein interactions and mirna target prediction data we used were useful to help elucidate the mechanistic role of mirnas in disease progression. although the interaction datasets are far from complete and suffer from noise, our results were consistent across choice of ppi network. using additional protein interaction networks, different mirna target prediction algorithms, and different expression data sets will likely reveal more mirnas with high-influence on cancer progression. another future direction for this work is designing a systematic method to combine the three variables that determine the influence of mirnas on the target partners. finally, this study on bridging the gap between clinical bioinformatics and network-based biomarkers provides clear evidence that protein interaction information is useful to identify diagnostic and prognostic cancer biomarkers, and to ameliorate the understanding of the functional mechanisms of mirnas.  we have developed a novel method to identify active mirna-target modules relevant to prostate cancer progression and outcome. mirnas with high influence on protein networks are valuable biomarkers that can be used in clinical investigations for prostate cancer treatment. combining the effects of mirnas on targets and target partners provides better understanding of mirnas function. competing interests the authors declare that they have no competing interests. authors’ contributions ma, gb and ag designed the study. ma wrote the initial manuscript. gb edited the manuscript. qm and ra proof-read the manuscript. all authors contributed to the writing and approved the final manuscript. supplementary material additional file 1 experimentally validated prostate mirnas and supplementary figures. this file lists mirnas that have been studied in prostate cancer and showed to play a role in prostate cancer progression. it also shows some targets of these mirnas that have been validated in prostate cancer context. click here for file acknowledgements m. alshalalfa and r. alhajj would like to thank nserc for funding.
the objective of this study is to identify mirnas that have high influence on target protein complexes in prostate cancer as a case study. this could provide biomarkers or therapeutic targets relevant for prostate cancer treatment. mirnas with high influence on target protein complexes play a role in prostate cancer progression and are promising diagnostic or prognostic biomarkers. we uncovered several mirna-regulated protein modules which were enriched in focal adhesion and prostate cancer genes. several mirnas such as mir- mir- and mir- demonstrated high influence on their target protein complexes and could explain most of the gene expression changes in our analyzed prostate cancer data set. we describe a novel method to identify active mirna-target modules relevant to prostate cancer progression and outcome. mirnas with high influence on protein networks are valuable biomarkers that can be used in clinical investigations for prostate cancer treatment. to detect mirnas with high influence on target protein modules, we integrated mirna and mrna expression profiles with a sequence based mirna-target network and human functional and physical protein interactions. our findings demonstrate that a mirna’s functional role can be explained by its target protein connectivity within a physical and functional interaction network. the network of regulatory mirna-target and functional protein interactions is a rich source of information to elucidate the function and the prognostic value of mirnas in cancer.
 the use of biological molecular network information for diagnostic and prognostic purposes and elucidation of molecular disease mechanism is a key objective in systems biomedicine. the network of regulatory mirna-target and functional protein interactions is a rich source of information to elucidate the function and the prognostic value of mirnas in cancer. the objective of this study is to identify mirnas that have high influence on target protein complexes in prostate cancer as a case study. this could provide biomarkers or therapeutic targets relevant for prostate cancer treatment.  our findings demonstrate that a mirna’s functional role can be explained by its target protein connectivity within a physical and functional interaction network. to detect mirnas with high influence on target protein modules, we integrated mirna and mrna expression profiles with a sequence based mirna-target network and human functional and physical protein interactions. mirnas with high influence on target protein complexes play a role in prostate cancer progression and are promising diagnostic or prognostic biomarkers. we uncovered several mirna-regulated protein modules which were enriched in focal adhesion and prostate cancer genes. several mirnas such as mir- mir- and mir- demonstrated high influence on their target protein complexes and could explain most of the gene expression changes in our analyzed prostate cancer data set.  we describe a novel method to identify active mirna-target modules relevant to prostate cancer progression and outcome. mirnas with high influence on protein networks are valuable biomarkers that can be used in clinical investigations for prostate cancer treatment. mirnaprotein interactionssystems biologyhigh-influence mirna 