BNrich: A Bayesian network approach to the pathway enrichment analysis

Abstract

MotivationOne of the most popular techniques in biological studies for analyzing high throughput data is pathway enrichment analysis (PEA). Many researchers apply the existing methods without considering the topology of pathways or at least they have overlooked a significant part of the structure, which may reduce the accuracy and generalizability of the results. Developing a new approach while considering gene expression data and topological features like causal relations regarding edge directions will help the investigators to achieve more accurate results. ResultsWe proposed a new pathway enrichment analysis based on Bayesian network (BNrich) as an approach in PEA. To this end, the cycles were eliminated in 187 KEGG human signaling pathways concerning intuitive biological rules and the Bayesian network structures were constructed. The constructed networks were simplified by the Least Absolute Shrinkage Selector Operator (LASSO), and their parameters were estimated using the gene expression data. We finally prioritize the impacted pathways by Fishers Exact Test on significant parameters. Our method integrates both edge and node related parameters to enrich modules in the affected signaling pathway network. In order to evaluate the proposed method, consistency, discrimination, false positive rate and empirical P-value criteria were calculated, and the results are compared to well-known enrichment methods such as signaling pathway impact analysis (SPIA), bi-level meta-analysis (BLMA) and topology-based pathway enrichment analysis (TPEA). AvailabilityThe R package is available on carn.