WMAXC: A Weighted Maximum Clique Method for Identifying Condition-Specific Sub-Network

Sub-networks can expose complex patterns in an entire bio-molecular network by extracting interactions that depend on temporal or condition-specific contexts. When genes interact with each other during cellular processes, they may form differential co-expression patterns with other genes across different cell states. The identification of condition-specific sub-networks is of great importance in investigating how a living cell adapts to environmental changes. In this work, we propose the weighted MAXimum clique (WMAXC) method to identify a condition-specific sub-network. WMAXC first proposes scoring functions that jointly measure condition-specific changes to both individual genes and gene-gene co-expressions. It then employs a weaker formula of a general maximum clique problem and relates the maximum scored clique of a weighted graph to the optimization of a quadratic objective function under sparsity constraints. We combine a continuous genetic algorithm and a projection procedure to obtain a single optimal sub-network that maximizes the objective function (scoring function) over the standard simplex (sparsity constraints). We applied the WMAXC method to both simulated data and real data sets of ovarian and prostate cancer. Compared with previous methods, WMAXC selected a large fraction of cancer-related genes, which were enriched in cancer-related pathways. The results demonstrated that our method efficiently captured a subset of genes relevant under the investigated condition.

[1]  Volkan Cevher,et al.  Sparse projections onto the simplex , 2012, ICML.

[2]  Daniela M Witten,et al.  Extensions of Sparse Canonical Correlation Analysis with Applications to Genomic Data , 2009, Statistical applications in genetics and molecular biology.

[3]  Natalie Wilson,et al.  Human Protein Reference Database , 2004, Nature Reviews Molecular Cell Biology.

[4]  I. Shih,et al.  Notch3 overexpression is related to the recurrence of ovarian cancer and confers resistance to carboplatin. , 2010, The American journal of pathology.

[5]  H. Hollema,et al.  The ErbB signalling pathway: protein expression and prognostic value in epithelial ovarian cancer , 2008, British Journal of Cancer.

[6]  W. Gerald,et al.  Targeting AKT/mTOR and ERK MAPK signaling inhibits hormone-refractory prostate cancer in a preclinical mouse model. , 2008, The Journal of clinical investigation.

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

[8]  M. T. Abreu-Martin,et al.  Mitogen-Activated Protein Kinase Kinase Kinase 1 Activates Androgen Receptor-Dependent Transcription and Apoptosis in Prostate Cancer , 1999, Molecular and Cellular Biology.

[9]  S. Murphy,et al.  The activated transforming growth factor‐beta signaling pathway in peritoneal metastases is a potential therapeutic target in ovarian cancer , 2012, International journal of cancer.

[10]  A. Sinclair,et al.  Dppa2 and Dppa4 Are Closely Linked SAP Motif Genes Restricted to Pluripotent Cells and the Germ Line , 2007, Stem cells.

[11]  R. Tibshirani,et al.  Gene expression profiling identifies clinically relevant subtypes of prostate cancer. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[12]  Richard M. Karp,et al.  Detecting Disease-Specific Dysregulated Pathways Via Analysis of Clinical Expression Profiles , 2008, RECOMB.

[13]  M. Burdick,et al.  L-selectin-mediated Lymphocyte-Cancer Cell Interactions under Low Fluid Shear Conditions* , 2008, Journal of Biological Chemistry.

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

[15]  R. Tibshirani,et al.  Significance analysis of microarrays applied to the ionizing radiation response , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[16]  Serban Nacu,et al.  Gene expression network analysis and applications to immunology , 2007, Bioinform..

[17]  M. Kattan,et al.  Preoperative plasma levels of transforming growth factor beta(1) (TGF-beta(1)) strongly predict progression in patients undergoing radical prostatectomy. , 2001, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[18]  Vladimir B. Bajic,et al.  Database for exploration of functional context of genes implicated in ovarian cancer , 2008, Nucleic Acids Res..

[19]  Hongyu Zhao,et al.  COSINE: COndition-SpecIfic sub-NEtwork identification using a global optimization method , 2011, Bioinform..

[20]  Y. Hamada,et al.  Autocrine expression of neurotrophins and their receptors in prostate cancer , 2001, International journal of urology : official journal of the Japanese Urological Association.

[21]  M. Kattan,et al.  Elevated levels of circulating interleukin-6 and transforming growth factor-beta1 in patients with metastatic prostatic carcinoma. , 1999, The Journal of urology.

[22]  A. Fiacco A Finite Algorithm for Finding the Projection of a Point onto the Canonical Simplex of R " , 2009 .

[23]  Benno Schwikowski,et al.  Discovering regulatory and signalling circuits in molecular interaction networks , 2002, ISMB.

[24]  Robert Clarke,et al.  Multilevel support vector regression analysis to identify condition-specific regulatory networks , 2010, Bioinform..

[25]  Pankaj Agarwal,et al.  Inferring pathways from gene lists using a literature-derived network of biological relationships , 2005, Bioinform..

[26]  E. van den Berg,et al.  Upregulation of the transcription factor TFEB in t(6;11)(p21;q13)-positive renal cell carcinomas due to promoter substitution. , 2003, Human molecular genetics.

[27]  Liang Chen,et al.  A statistical method for identifying differential gene-gene co-expression patterns , 2004, Bioinform..

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

[29]  C. Michelot A finite algorithm for finding the projection of a point onto the canonical simplex of ∝n , 1986 .

[30]  L. Moser,et al.  AN EXTREMAL PROBLEM IN GRAPH THEORY , 2001 .

[31]  Brad T. Sherman,et al.  Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources , 2008, Nature Protocols.

[32]  Aeilko H. Zwinderman,et al.  Sparse canonical correlation analysis for identifying, connecting and completing gene-expression networks , 2009, BMC Bioinformatics.

[33]  H. Kung,et al.  Signal Pathways in Up-regulation of Chemokines by Tyrosine Kinase MER/NYK in Prostate Cancer Cells , 2004, Cancer Research.

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

[35]  Boris Goldengorin,et al.  Handbook of combinatorial optimization , 2013 .

[36]  Rainer Breitling,et al.  Graph-based iterative Group Analysis enhances microarray interpretation , 2004, BMC Bioinformatics.

[37]  Christopher Thrasivoulou,et al.  A Novel Role for Wnt/Ca2+ Signaling in Actin Cytoskeleton Remodeling and Cell Motility in Prostate Cancer , 2010, PloS one.

[38]  Jing Zhu,et al.  Edge-based scoring and searching method for identifying condition-responsive protein-protein interaction sub-network , 2007, Bioinform..

[39]  T. Motzkin,et al.  Maxima for Graphs and a New Proof of a Theorem of Turán , 1965, Canadian Journal of Mathematics.

[40]  Pär Stattin,et al.  Transforming growth factor β1 is associated with angiogenesis, metastasis, and poor clinical outcome in prostate cancer , 1998 .

[41]  P. Erdös On an extremal problem in graph theory , 1970 .

[42]  H. Frierson,et al.  Activation of mitogen-activated protein kinase associated with prostate cancer progression. , 1999, Cancer research.

[43]  X-S Zhang,et al.  Identifying differentially expressed pathways via a mixed integer linear programming model. , 2009, IET systems biology.

[44]  O. Kallioniemi,et al.  FZD4 as a mediator of ERG oncogene-induced WNT signaling and epithelial-to-mesenchymal transition in human prostate cancer cells. , 2010, Cancer research.

[45]  Christina Backes,et al.  An integer linear programming approach for finding deregulated subgraphs in regulatory networks , 2011, Nucleic acids research.

[46]  Robert Clarke,et al.  Identifying protein interaction subnetworks by a bagging Markov random field-based method , 2012, Nucleic acids research.

[47]  Valesca Anschau,et al.  Dissecting Major Signaling Pathways throughout the Development of Prostate Cancer , 2013, Prostate cancer.

[48]  Natalie Wilson Human Protein Reference Database , 2004, Nature Reviews Genetics.

[49]  Vladimir B. Bajic,et al.  DDPC: Dragon Database of Genes associated with Prostate Cancer , 2010, Nucleic Acids Res..