BMRF-MI: integrative identification of protein interaction network by modeling the gene dependency

Abstract

BackgroundIdentification of protein interaction network is a very important step for understanding the molecular mechanisms in cancer. Several methods have been developed to integrate protein-protein interaction (PPI) data with gene expression data for network identification. However, they often fail to model the dependency between genes in the network, which makes many important genes, especially the upstream genes, unidentified. It is necessary to develop a method to improve the network identification performance by incorporating the dependency between genes.ResultsWe proposed an approach for identifying protein interaction network by incorporating mutual information (MI) into a Markov random field (MRF) based framework to model the dependency between genes. MI is widely used in information theory to measure the uncertainty between random variables. Different from traditional Pearson correlation test, MI is capable of capturing both linear and non-linear relationship between random variables. Among all the existing MI estimators, we choose to use k-nearest neighbor MI (kNN-MI) estimator which is proved to have minimum bias. The estimated MI is integrated with an MRF framework to model the gene dependency in the context of network. The maximum a posterior (MAP) estimation is applied on the MRF-based model to estimate the network score. In order to reduce the computational complexity of finding the optimal network, a probabilistic searching algorithm is implemented. We further increase the robustness and reproducibility of the results by applying a non-parametric bootstrapping method to measure the confidence level of the identified genes. To evaluate the performance of the proposed method, we test the method on simulation data under different conditions. The experimental results show an improved accuracy in terms of subnetwork identification compared to existing methods. Furthermore, we applied our method onto real breast cancer patient data; the identified protein interaction network shows a close association with the recurrence of breast cancer, which is supported by functional annotation. We also show that the identified subnetworks can be used to predict the recurrence status of cancer patients by survival analysis.ConclusionsWe have developed an integrated approach for protein interaction network identification, which combines Markov random field framework and mutual information to model the gene dependency in PPI network. Improvements in subnetwork identification have been demonstrated with simulation datasets compared to existing methods. We then apply our method onto breast cancer patient data to identify recurrence related subnetworks. The experiment results show that the identified genes are enriched in the pathway and functional categories relevant to progression and recurrence of breast cancer. Finally, the survival analysis based on identified subnetworks achieves a good result of classifying the recurrence status of cancer patients.

Identification of cardiomyopathy-related core genes through human metabolic networks and expression data

“…Based on our previous study [ 44 ] and the calculation process in [ 45 ], the MRFms for module M incorporating Mutual Information (MI) was defined as where b is the number of edges, G 1 and G 2 are the set of DEGs and non-DEGs in the module; E u , E v and E z are the expression differences of genes u , v and z between different states (DCM_NF, ICM_NF, or DCM_ICM), and d v and d z are the degrees of genes v and z in the network, respectively. MI ( v , z ) is the mutual information of genes v and z .…”

Multi-omics analyses of drug repurposing reveal Acebutolol and Amiloride for osteoporosis treatment

“…First, the subnetwork searching started from every osteoporosis-related gene (seed gene) that formed an initial subnetwork with a single node, then randomly sampled a new gene in the background network to get a new subnetwork. Genes adjacent to the current network can be added, and the subnetwork gradually increased a predefined network score NetScore(G) , which was defined as a negative posterior potential function measured from an estimated discriminative Z-score of the gene between low and high BMD groups 32,34 . Simulated annealing stopped and outputted a final subnetwork when no significant improvement of the network score was achieved.…”

Identification of evolutionarily conserved DNA damage response genes that alter sensitivity to cisplatin

“…Indeed a growing number of components involved in DNA repair systems have been defined. Based on analysis of these components, it has become clear that for many vital cellular processes, changes in gene expression or function that affect phenotypes of interest can be dispersed throughout extended signaling networks [ 16 – 18 ]. In this context, even as genomic profiling has become more standard in the clinic [ 19 , 20 ], it remains a challenge to identify functionally essential components of the DDR response apparatus relevant to clinical resistance to cisplatin and other DNA damaging therapies.…”