Differential Network Analysis via Weighted Fused Conditional Gaussian Graphical Model

The development and prognosis of complex diseases usually involves changes in regulatory relationships among biomolecules. Understanding how the regulatory relationships change with genetic alterations can help to reveal the underlying biological mechanisms for complex diseases. Although several models have been proposed to estimate the differential network between two different states, they are not suitable to deal with situations where the molecules of interest are affected by other covariates. Nor can they make use of prior information that provides insights about the structures of biomolecular networks. In this study, we introduce a novel weighted fused conditional Gaussian graphical model to jointly estimate two state-specific biomolecular regulatory networks and their difference between two different states. Unlike previous differential network estimation methods, our model can take into account the related covariates and the prior network information when inferring differential networks. The effectiveness of our proposed model is first evaluated based on simulation studies. Experiment results demonstrate that our model outperforms other state-of-the-art differential networks estimation models in all cases. We then apply our model to identify the differential gene network between two subtypes of glioblastoma based on gene expression and miRNA expression data. Our model is able to discover known mechanisms of glioblastoma and provide interesting predictions.

[1]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[2]  Eloisa Arbustini,et al.  Cell density modulates SHC3 expression and survival of human glioblastoma cells through Fak activation , 2014, Journal of Neuro-Oncology.

[3]  Michael Sabel,et al.  Genomic and Expression Profiling of Glioblastoma Stem Cell–Like Spheroid Cultures Identifies Novel Tumor-Relevant Genes Associated with Survival , 2009, Clinical Cancer Research.

[4]  Marinka Zitnik,et al.  Gene network inference by fusing data from diverse distributions , 2015, Bioinform..

[5]  H. Zou,et al.  Sparse precision matrix estimation via lasso penalized D-trace loss , 2014 .

[6]  Mong-Lien Wang,et al.  MicroRNA-142-3p is involved in regulation of MGMT expression in glioblastoma cells , 2018, Cancer management and research.

[7]  Rajen Dinesh Shah,et al.  Variable selection with error control: another look at stability selection , 2011, 1105.5578.

[8]  Joshua M. Korn,et al.  Comprehensive genomic characterization defines human glioblastoma genes and core pathways , 2008, Nature.

[9]  Christian L. Müller,et al.  Sparse and Compositionally Robust Inference of Microbial Ecological Networks , 2014, PLoS Comput. Biol..

[10]  Xiaofeng Song,et al.  MicroRNA-137 inhibits growth of glioblastoma through EGFR suppression. , 2017, American journal of translational research.

[11]  Hongzhe Li,et al.  A SPARSE CONDITIONAL GAUSSIAN GRAPHICAL MODEL FOR ANALYSIS OF GENETICAL GENOMICS DATA. , 2011, The annals of applied statistics.

[12]  Le Ou-Yang,et al.  Identifying Gene Network Rewiring by Integrating Gene Expression and Gene Network Data , 2018, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[13]  Hiroyuki Ogata,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 1999, Nucleic Acids Res..

[14]  Xiaoyu Huang,et al.  miR-484/MAP2/c-Myc-positive regulatory loop in glioma promotes tumor-initiating properties through ERK1/2 signaling , 2018, Journal of Molecular Histology.

[15]  Jung Eun Shim,et al.  TRRUST: a reference database of human transcriptional regulatory interactions , 2015, Scientific Reports.

[16]  Yirong Wang,et al.  Analysis of hsa-miR-30a-5p Expression in Human Gliomas , 2013, Pathology & Oncology Research.

[17]  Feihu Huang,et al.  Learning Dynamic Conditional Gaussian Graphical Models , 2018, IEEE Transactions on Knowledge and Data Engineering.

[18]  K. Oguchi,et al.  Reduced PTEN expression and overexpression of miR-17-5p, -19a-3p, -19b-3p, -21-5p, -130b-3p, -221-3p and -222-3p by glioblastoma stem-like cells following irradiation. , 2015, Oncology letters.

[19]  Ling-Zhi Wang,et al.  Targetable BET proteins- and E2F1-dependent transcriptional program maintains the malignancy of glioblastoma , 2018, Proceedings of the National Academy of Sciences.

[20]  Quanquan Gu,et al.  Semiparametric Differential Graph Models , 2016, NIPS.

[21]  James U. Bowie,et al.  Network rewiring is an important mechanism of gene essentiality change , 2012, Scientific Reports.

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

[23]  Tristan Mary-Huard,et al.  Structured regularization for conditional Gaussian graphical models , 2014, Statistics and Computing.

[24]  Feihu Huang,et al.  Joint Estimation of Multiple Conditional Gaussian Graphical Models , 2018, IEEE Transactions on Neural Networks and Learning Systems.

[25]  Kai Li,et al.  Epidermal growth factor receptor in glioblastoma , 2017, Oncology letters.

[26]  Xiao-Tong Yuan,et al.  Partial Gaussian Graphical Model Estimation , 2012, IEEE Transactions on Information Theory.

[27]  Jian Ma,et al.  Covariate Adjusted Precision Matrix Estimation via Nonconvex Optimization , 2018, ICML.

[28]  Thomas Thorne,et al.  NetDiff – Bayesian model selection for differential gene regulatory network inference , 2016, Scientific Reports.

[29]  Jing Xu,et al.  Expression of transcription factor E2F1 and telomerase in glioblastomas: mechanistic linkage and prognostic significance. , 2005, Journal of the National Cancer Institute.

[30]  J. Zico Kolter,et al.  Sparse Gaussian Conditional Random Fields: Algorithms, Theory, and Application to Energy Forecasting , 2013, ICML.

[31]  Quanquan Gu,et al.  Identifying gene regulatory network rewiring using latent differential graphical models , 2016, Nucleic acids research.

[32]  Hong Yan,et al.  Incorporating prior information into differential network analysis using non‐paranormal graphical models , 2017, Bioinform..

[33]  Hsien-Da Huang,et al.  miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions , 2017, Nucleic Acids Res..

[34]  Hong Yan,et al.  Differential network analysis from cross-platform gene expression data , 2016, Scientific Reports.

[35]  Kim-Anh Do,et al.  DINGO: differential network analysis in genomics , 2015, Bioinform..

[36]  Ulf Leser,et al.  Comparative assessment of differential network analysis methods , 2016, Briefings Bioinform..

[37]  Patrick Danaher,et al.  The joint graphical lasso for inverse covariance estimation across multiple classes , 2011, Journal of the Royal Statistical Society. Series B, Statistical methodology.

[38]  Kyung-Ah Sohn,et al.  Joint Estimation of Structured Sparsity and Output Structure in Multiple-Output Regression via Inverse-Covariance Regularization , 2012, AISTATS.

[39]  Ivan Smirnov,et al.  Clinical Stratification of Glioblastoma Based on Alterations in Retinoblastoma Tumor Suppressor Protein (RB1) and Association With the Proneural Subtype , 2012, Journal of neuropathology and experimental neurology.

[40]  Chun Luo,et al.  Knockdown of cyclin D1 inhibits proliferation, induces apoptosis, and attenuates the invasive capacity of human glioblastoma cells , 2012, Journal of Neuro-Oncology.

[41]  T. Cai,et al.  Direct estimation of differential networks. , 2014, Biometrika.

[42]  Larry A. Wasserman,et al.  Stability Approach to Regularization Selection (StARS) for High Dimensional Graphical Models , 2010, NIPS.

[43]  Erik P. Sulman,et al.  CDK4/6 inhibition is more active against the glioblastoma proneural subtype. , 2017, Oncotarget.

[44]  Le Ou-Yang,et al.  Identifying differential networks based on multi-platform gene expression data. , 2016, Molecular bioSystems.

[45]  S. Gabriel,et al.  Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. , 2010, Cancer cell.

[46]  Hongyu Zhao,et al.  Sparse Estimation of Conditional Graphical Models With Application to Gene Networks , 2012, Journal of the American Statistical Association.

[47]  Genevera I. Allen,et al.  A Local Poisson Graphical Model for Inferring Networks From Sequencing Data , 2013, IEEE Transactions on NanoBioscience.

[48]  Israel Steinfeld,et al.  miRNA-mRNA Integrated Analysis Reveals Roles for miRNAs in Primary Breast Tumors , 2011, PloS one.

[49]  M. Yuan,et al.  Model selection and estimation in the Gaussian graphical model , 2007 .

[50]  Min Wu,et al.  Node-based learning of differential networks from multi-platform gene expression data. , 2017, Methods.

[51]  R. Tibshirani,et al.  Sparse inverse covariance estimation with the graphical lasso. , 2008, Biostatistics.

[52]  Yurii Nesterov,et al.  Smooth minimization of non-smooth functions , 2005, Math. Program..