A powerful weighted statistic for detecting group differences of directed biological networks

Complex disease is largely determined by a number of biomolecules interwoven into networks, rather than a single biomolecule. Different physiological conditions such as cases and controls may manifest as different networks. Statistical comparison between biological networks can provide not only new insight into the disease mechanism but statistical guidance for drug development. However, the methods developed in previous studies are inadequate to capture the changes in both the nodes and edges, and often ignore the network structure. In this study, we present a powerful weighted statistical test for group differences of directed biological networks, which is independent of the network attributes and can capture the changes in both the nodes and edges, as well as simultaneously accounting for the network structure through putting more weights on the difference of nodes locating on relatively more important position. Simulation studies illustrate that this method had better performance than previous ones under various sample sizes and network structures. One application to GWAS of leprosy successfully identifies the specific gene interaction network contributing to leprosy. Another real data analysis significantly identifies a new biological network, which is related to acute myeloid leukemia. One potential network responsible for lung cancer has also been significantly detected. The source R code is available on our website.

[1]  P. Good,et al.  Permutation Tests: A Practical Guide to Resampling Methods for Testing Hypotheses , 1995 .

[2]  Huan Liu,et al.  A new chi-square approximation to the distribution of non-negative definite quadratic forms in non-central normal variables , 2009, Comput. Stat. Data Anal..

[3]  Eiliv Lund,et al.  Systems Epidemiology in Cancer , 2008, Cancer Epidemiology Biomarkers & Prevention.

[4]  David Warde-Farley,et al.  Dynamic modularity in protein interaction networks predicts breast cancer outcome , 2009, Nature Biotechnology.

[5]  T. Ottenhoff,et al.  Control of human host immunity to mycobacteria. , 2005, Tuberculosis.

[6]  Biao He,et al.  Activation of the Wnt pathway in non small cell lung cancer: evidence of dishevelled overexpression , 2003, Oncogene.

[7]  Y. Moreau,et al.  Finding the targets of a drug by integration of gene expression data with a protein interaction network. , 2013, Molecular bioSystems.

[8]  Steve Horvath,et al.  WGCNA: an R package for weighted correlation network analysis , 2008, BMC Bioinformatics.

[9]  A. Barabasi,et al.  Network medicine : a network-based approach to human disease , 2010 .

[10]  T. Ideker,et al.  Network-based classification of breast cancer metastasis , 2007, Molecular systems biology.

[11]  Joseph L. Fleiss,et al.  On the Distribution of a Linear Combination of Independent Chi Squares , 1971 .

[12]  Oliver Eickelberg,et al.  WNT signaling in lung disease: a failure or a regeneration signal? , 2010, American journal of respiratory cell and molecular biology.

[13]  Fernando M. A. Silva,et al.  Network comparison using directed graphlets , 2015, ArXiv.

[14]  Aziz Ghahary,et al.  Inflammatory Effects of Ex Vivo Human Th17 Cells Are Suppressed by Regulatory T Cells , 2010, The Journal of Immunology.

[15]  Nitai D. Mukhopadhyay,et al.  An inferential framework for biological network hypothesis tests , 2013, BMC Bioinformatics.

[16]  J. Cheverud,et al.  Lack of p21 expression links cell cycle control and appendage regeneration in mice , 2010, Proceedings of the National Academy of Sciences.

[17]  Daoxin Ma,et al.  Th17 and Treg Cells in Bone Related Diseases , 2013, Clinical & developmental immunology.

[18]  Emmanuel Barillot,et al.  Classification of microarray data using gene networks , 2007, BMC Bioinformatics.

[19]  Daoxin Ma,et al.  The Profile of T Helper Subsets in Bone Marrow Microenvironment Is Distinct for Different Stages of Acute Myeloid Leukemia Patients and Chemotherapy Partly Ameliorates These Variations , 2015, PloS one.

[20]  E. Schadt.  Molecular networks as sensors and drivers of common human diseases , 2009, Nature.

[21]  Zoran Levnajic,et al.  Revealing the Hidden Language of Complex Networks , 2014, Scientific Reports.

[22]  Abdul Salam Jarrah,et al.  The effect of negative feedback loops on the dynamics of boolean networks. , 2007, Biophysical journal.

[23]  J. Casanova,et al.  Genetic dissection of immunity to mycobacteria: the human model. , 2002, Annual review of immunology.

[24]  D. Philpott,et al.  The ubiquitin-editing enzyme A20 restricts nucleotide-binding oligomerization domain containing 2-triggered signals. , 2008, Immunity.

[25]  Tijana Milenkovic,et al.  Proper evaluation of alignment-free network comparison methods , 2015, Bioinform..

[26]  Hongzhe Li,et al.  In Response to Comment on "Network-constrained regularization and variable selection for analysis of genomic data" , 2008, Bioinform..

[27]  Robert Clarke,et al.  Differential dependency network analysis to identify condition-specific topological changes in biological networks , 2009, Bioinform..

[28]  R. Haring,et al.  Diving through the "-omics": the case for deep phenotyping and systems epidemiology. , 2012, Omics : a journal of integrative biology.

[29]  M. Netea,et al.  Genomewide association study of leprosy. , 2010, The New England journal of medicine.

[30]  Xiaoshuai Zhang,et al.  A powerful score-based statistical test for group difference in weighted biological networks , 2016, BMC Bioinformatics.

[31]  Adam M. Gustafson,et al.  Airway PI3K Pathway Activation Is an Early and Reversible Event in Lung Cancer Development , 2010, Science Translational Medicine.

[32]  Monica Chiogna,et al.  Gene set analysis exploiting the topology of a pathway , 2010, BMC Systems Biology.

[33]  Luonan Chen,et al.  Biomolecular Networks: Methods and Applications in Systems Biology , 2009 .

[34]  Michael A McGuckin,et al.  The interplay between endoplasmic reticulum stress and inflammation , 2012, Immunology and cell biology.

[35]  Xia Li,et al.  Network-based survival-associated module biomarker and its crosstalk with cell death genes in ovarian cancer , 2015, Scientific Reports.

[36]  Michael Q. Zhang,et al.  Network-based global inference of human disease genes , 2008, Molecular systems biology.

[37]  Kwang-Hyun Cho,et al.  Analysis of feedback loops and robustness in network evolution based on Boolean models , 2007, BMC Bioinformatics.

[38]  Jing Xu,et al.  Detection for pathway effect contributing to disease in systems epidemiology with a case–control design , 2015, BMJ Open.

[39]  D. Stewart,et al.  Review Wnt Signaling Pathway in Non–small Cell Lung Cancer Overview of the Wnt Canonical (β-catenin) and Noncanonical Signaling Pathways , 2022 .

[40]  Mayer Alvo,et al.  Testing for mean and correlation changes in microarray experiments: an application for pathway analysis , 2010, BMC Bioinformatics.

[41]  Yaling Yin,et al.  Network-Based Inference Framework for Identifying Cancer Genes from Gene Expression Data , 2013, BioMed research international.

[42]  S. Rosenberg,et al.  IL-2 administration increases CD4+ CD25(hi) Foxp3+ regulatory T cells in cancer patients. , 2006, Blood.

[43]  P. Sebastiani,et al.  Airway epithelial gene expression in the diagnostic evaluation of smokers with suspect lung cancer , 2007, Nature Medicine.