Recently, computational approaches integrating copy number aberrations (CNAs) and gene expression (GE) have been extensively studied to identify cancer-related genes and pathways. In this work, we integrate these two data sets with protein-protein interaction (PPI) information to find cancer-related functional modules. To integrate CNA and GE data, we first built a gene-gene relationship network from a set of seed genes by enumerating all types of pairwise correlations, e.g. GE-GE, CNA-GE, and CNA-CNA, over multiple patients. Next, we propose a voting-based cancer module identification algorithm by combining topological and data-driven properties (VToD algorithm) by using the gene-gene relationship network as a source of data-driven information, and the PPI data as topological information. We applied the VToD algorithm to 266 glioblastoma multiforme (GBM) and 96 ovarian carcinoma (OVC) samples that have both expression and copy number measurements, and identified 22 GBM modules and 23 OVC modules. Among 22 GBM modules, 15, 12, and 20 modules were significantly enriched with cancer-related KEGG, BioCarta pathways, and GO terms, respectively. Among 23 OVC modules, 19, 18, and 23 modules were significantly enriched with cancer-related KEGG, BioCarta pathways, and GO terms, respectively. Similarly, we also observed that 9 and 2 GBM modules and 15 and 18 OVC modules were enriched with cancer gene census (CGC) and specific cancer driver genes, respectively. Our proposed module-detection algorithm significantly outperformed other existing methods in terms of both functional and cancer gene set enrichments. Most of the cancer-related pathways from both cancer data sets found in our algorithm contained more than two types of gene-gene relationships, showing strong positive correlations between the number of different types of relationship and CGC enrichment -values (0.64 for GBM and 0.49 for OVC). This study suggests that identified modules containing both expression changes and CNAs can explain cancer-related activities with greater insights.
References
[1]
Hahn WC, Weinberg RA (2002) Modelling the molecular circuitry of cancer. Nat Rev Cancer 2: 331–341.
[2]
Vogelstein B, Kinzler KW (2004) Cancer genes and the pathways they control. Nat Med 10: 789–799.
[3]
Davies H, Bignell GR, Cox C, Stephens P, Edkins S, et al. (2002) Mutations of the BRAF gene in human cancer. Nature 417: 949–954.
[4]
Wan PTC, Garnett MJ, Roe SM, Lee S, Niculescu-Duvaz D, et al. (2004) Mechanism of Activation of the RAF-ERK Signaling Pathway by Oncogenic Mutations of B-RAF. Cell 116: 855–867.
[5]
Santarosa M, Ashworth A (2004) Haploinsufficiency for tumour suppressor genes: when you don't need to go all the way. Biochimica et Biophysica Acta (BBA) - Reviews on Cancer 1654: 105–122.
[6]
Chuang HY, Lee E, Liu YT, Lee D, Ideker T (2007) Network-based classification of breast cancer metastasis. Mol Syst Biol 3: 140.
[7]
Vaske CJ, Benz SC, Sanborn JZ, Earl D, Szeto C, et al. (2010) Inference of patient-specific pathway activities from multi-dimensional cancer genomics data using PARADIGM. Bioinformatics 26: i237–i245.
[8]
TCGA (2008) Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature 455: 1061–1068.
[9]
Cerami E, Demir E, Schultz N, Taylor BS, Sander C (2010) Automated Network Analysis Identifies Core Pathways in Glioblastoma. PLoS ONE 5: e8918.
[10]
Feuk L, Carson AR, Scherer SW (2006) Structural variation in the human genome. Nat Rev Genet 7: 85–97.
[11]
Beroukhim R, Getz G, Nghiemphu L, Barretina J, Hsueh T, et al. (2007) Assessing the significance of chromosomal aberrations in cancer: Methodology and application to glioma. Proceedings of the National Academy of Sciences 104: 20007–20012.
[12]
Taylor BS, Barretina J, Socci ND, DeCarolis P, Ladanyi M, et al. (2008) Functional Copy-Number Alterations in Cancer. PLoS ONE 3: e3179.
[13]
Hur Y, Lee H (2011) Wavelet-based identification of DNA focal genomic aberrations from single nucleotide polymorphism arrays. BMC Bioinformatics 12: 146.
Jornsten R, Abenius T, Kling T, Schmidt L, Johansson E, et al.. (2011) Network modeling of the transcriptional effects of copy number aberrations in glioblastoma. Mol Syst Biol 7.
[16]
Akavia UD, Litvin O, Kim J, Sanchez-Garcia F, Kotliar D, et al. (2010) An Integrated Approach to Uncover Drivers of Cancer. Cell 143: 1005–1017.
[17]
Witten DM, Tibshirani R, Hastie T (2009) A penalized matrix decomposition, with applications to sparse principal components and canonical correlation analysis. Biostatistics 10: 515–534.
[18]
Horvath S, Zhang B, Carlson M, Lu KV, Zhu S, et al. (2006) Analysis of oncogenic signaling networks in glioblastoma identifies ASPM as a molecular target. Proceedings of the National Academy of Sciences 103: 17402–17407.
[19]
Choi JK, Yu U, Yoo OJ, Kim S (2005) Differential coexpression analysis using microarray data and its application to human cancer. Bioinformatics 21: 4348–4355.
[20]
Mo W, Fu X, Han X, Yang G, Zhang J, et al. (2009) A stochastic model for identifying differential gene pair co-expression patterns in prostate cancer progression. BMC Genomics 10: 340.
[21]
Gorringe KL, George J, Anglesio MS, Ramakrishna M, Etemadmoghadam D, et al. (2010) Copy Number Analysis Identifies Novel Interactions Between Genomic Loci in Ovarian Cancer. PLoS ONE 5: e11408.
[22]
Klijn C, Bot J, Adams DJ, Reinders M, Wessels L, et al. (2010) Identification of Networks of Co-Occurring, Tumor-Related DNA Copy Number Changes Using a Genome-Wide Scoring Approach. PLoS Comput Biol 6: e1000631.
[23]
Cui Q (2010) A Network of Cancer Genes with Co-Occurring and Anti-Co-Occurring Mutations. PLoS ONE 5: e13180.
[24]
Masica DL, Karchin R (2011) Correlation of Somatic Mutation and Expression Identifies Genes Important in Human Glioblastoma Progression and Survival. Cancer Research 71: 4550–4561.
[25]
Lee H, Kong SW, Park PJ (2008) Integrative analysis reveals the direct and indirect interactions between DNA copy number aberrations and gene expression changes. Bioinformatics 24: 889–896.
[26]
Solvang H, Lingjaerde O, Frigessi A, Borresen-Dale AL, Kristensen V (2011) Linear and non-linear dependencies between copy number aberrations and mRNA expression reveal distinct molecular pathways in breast cancer. BMC Bioinformatics 12: 197.
[27]
Xu C, Liu Y, Wang P, Fan W, Rue T, et al. (2010) Integrative analysis of DNA copy number and gene expression in metastatic oral squamous cell carcinoma identifies genes associated with poor survival. Molecular Cancer 9: 143.
[28]
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M (2012) KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res 40: D109–14.
[29]
BioCarta. Available: http://www.biocarta.com/. Accessed 2011 Oct 19.
[30]
Ashburner M, Ball C, Blake J, Botstein D, Butler H, et al. (2000) Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 25: 25–9.
[31]
Futreal PA, Coin L, Marshall M, Down T, Hubbard T, et al. (2004) A census of human cancer genes. Nat Rev Cancer 4: 177–183.
[32]
Hwang W, Cho YR, Zhang A, Ramanathan M (2006) A novel functional module detection algorithm for protein-protein interaction networks. Algorithms for Molecular Biology 1: 24.
[33]
Parsons DW, Jones S, Zhang X, Lin JCH, Leary RJ, et al. (2008) An Integrated Genomic Analysis of Human Glioblastoma Multiforme. Science 321: 1807–1812.
[34]
Bast RC, Hennessy B, Mills GB (2009) The biology of ovarian cancer: new opportunities for translation. Nat Rev Cancer 9: 415–428.
[35]
Sarela A, Farmery S, Markham A, Guillou P (1999) The candidate tumour suppressor gene, ing1, is retained in colorectal carcinomas. European Journal of Cancer 35: 1264–1267.
[36]
Nakanishi A, Han X, Saito H, Taguchi K, Ohta Y, et al. (2007) Interference with brca2, which localizes to the centrosome during s and early m phase, leads to abnormal nuclear division. Biochemical and Biophysical Research Communications 355: 34–40.
[37]
Tallen G, Farhangi S, Tamannai M, Holtkamp N, Mangoldt D, et al. (2009) The inhibitor of growth 1 (ING1) proteins suppress angiogenesis and differentially regulate angiopoietin expression in glioblastoma cells. Oncol Res 18: 95–105.
[38]
Quiros S, Roos WP, Kaina B (2011) Rad51 and brca2 - new molecular targets for sensitizing glioma cells to alkylating anticancer drugs. PLoS ONE 6: e27183.
[39]
Li DM, Sun H (1998) Pten/mmac1/tep1 suppresses the tumorigenicity and induces g1 cell cycle arrest in human glioblastoma cells. Proceedings of the National Academy of Sciences 95: 15406–15411.
[40]
Liu Y, Shete S, Etzel CJ, Scheurer M, Alexiou G, et al. (2010) Polymorphisms of lig4, btbd2, hmga2, and rtel1 genes involved in the double-strand break repair pathway predict glioblastoma survival. Journal of Clinical Oncology 28: 2467–2474.
[41]
Watanabe A, Ogiwara H, Ehata S, Mukasa A, Ishikawa S, et al. (2011) Homozygously deleted gene dach1 regulates tumor-initiating activity of glioma cells. Proceedings of the National Academy of Sciences 108: 12384–12389.
[42]
del Mar Inda M, Fan X, Mu？oz J, Perot C, Fauvet D, et al. (2003) Chromosomal abnormalities in human glioblastomas: Gain in chromosome 7p correlating with loss in chromosome 10q. Molecular Carcinogenesis 36: 6–14.
[43]
Henson J, Schnitker B, Correa K, von Deimling A, Fassbender F, et al. (1994) The retinoblastoma gene is involved in malignant progression of astrocytomas. Ann Neurol 36: 714–21.
[44]
Carrero P, Okamoto K, Coumailleau P, O'Brien S, Tanaka H, et al. (2000) Redox-regulated recruitment of the transcriptional coactivators creb-binding protein and src-1 to hypoxia-inducible factor 1. Molecular and Cellular Biology 20: 402–415.
[45]
Donovan L, Potter N, Warr T, Pilkington G (2012) A Prominin-1-Rich Pediatric Glioblastoma: Biologic Behavior Is Determined by Oxygen Tension-Modulated CD133 Expression but Not Accompanied by Underlying Molecular Profiles. Transl Oncol 5: 141–54.
[46]
Irminger-Finger I, Leung WC, Li J, Dubois-Dauphin M, Harb J, et al. (2001) Identification of bard1 as mediator between proapoptotic stress and p53-dependent apoptosis. Molecular Cell 8: 1255–1266.
[47]
Sobottka S, Haase M, Fitze G, Hahn M, Schackert H, et al. (2000) Frequent loss of heterozygosity at the 19p13.3 locus without lkb1/stk11 mutations in human carcinoma metastases to the brain. Journal of Neuro-Oncology 49: 187–195.
[48]
Kanu OO, Hughes B, Di C, Lin N, Fu J, et al. (2009) Glioblastoma Multiforme Oncogenomics and Signaling Pathways. Clinical Medicine Insights: Oncology 3: 39.
[49]
Chen H, Ye D, Xie X, Chen B, Lu W (2004) Vegf, vegfrs expressions and activated stats in ovarian epithelial carcinoma. Gynecologic Oncology 94: 630–635.
[50]
Rosenthal LA, Winestock KD, Finbloom DS (1997) Il-2 and il-7 induce heterodimerization of stat5 isoforms in human peripheral blood t lymphoblasts. Cellular Immunology 181: 172–181.
[51]
Hennighausen L, Robinson G (2005) Information networks in the mammary gland. Nat Rev Mol Cell Biol 6: 715–25.
[52]
Bouilly J, Sonigo C, Auffret J, Gibori G, Binart N (2012) Prolactin signaling mechanisms in ovary. Molecular and Cellular Endocrinology 356: 80–87.
[53]
Grandis JR, Drenning SD, Chakraborty A, Zhou MY, Zeng Q, et al. (1998) Requirement of stat3 but not stat1 activation for epidermal growth factor receptor- mediated cell growth in vitro. The Journal of Clinical Investigation 102: 1385–1392.
[54]
Guren T, Abrahamsen H, Thoresen G, Babaie E, Berg T, et al. (1999) EGF-induced activation of Stat1, Stat3, and Stat5b is unrelated to the stimulation of DNA synthesis in cultured hepatocytes. Biochem Biophys Res Commun 258: 565–71.
[55]
Colomiere M, Ward A, Riley C, Trenerry M, Cameron-Smith D, et al. (2009) Cross talk of signals between EGFR and IL-6R through JAK2/STAT3 mediate epithelial-mesenchymal transition in ovarian carcinomas. Br J Cancer 100: 134–44.
[56]
Benito M, Valverde A, Lorenzo M (1996) IGF-I: a mitogen also involved in differentiation processes in mammalian cells. Int J Biochem Cell Biol 28: 499–510.
[57]
Ren Z, Schaefer T (2002) ErbB-2 activates Stat3 alpha in a Src- and JAK2-dependent manner. J Biol Chem 277: 38486–93.
[58]
Jones R, Gordus A, Krall J, MacBeath G (2006) A quantitative protein interaction network for the ErbB receptors using protein microarrays. Nature 439: 168–74.
[59]
Fujitani Y, Hibi M, Fukada T, Takahashi-Tezuka M, Yoshida H, et al. (1997) An alternative pathway for STAT activation that is mediated by the direct interaction between JAK and STAT. Oncogene 14: 751–61.
[60]
Barahmand-Pour F, Meinke A, Groner B, Decker T (1998) Jak2-Stat5 interactions analyzed in yeast. J Biol Chem 273: 12567–75.
[61]
Stortelers C, Souriau C, van Liempt E, van de Poll M, van Zoelen E (2002) Role of the N-terminus of epidermal growth factor in ErbB-2/ErbB-3 binding studied by phage display. Biochemistry 41: 8732–41.
[62]
Thornton C, Yaka R, Dinh S, Ron D (2003) H-Ras modulates N-methyl-D-aspartate receptor function via inhibition of Src tyrosine kinase activity. J Biol Chem 278: 23823–9.
[63]
Yip A, Horvath S (2007) Gene network interconnectedness and the generalized topological overlap measure. BMC Bioinformatics 8: 22.
[64]
Langfelder P, Zhang B, Horvath S (2008) Defining clusters from a hierarchical cluster tree: the Dynamic Tree Cut package for R. Bioinformatics. 24: 719–720.
[65]
Girvan M, Newman MEJ (2002) Community structure in social and biological networks. Proceedings of the National Academy of Sciences 99: 7821–7826.
[66]
Ulitsky I, Shamir R (2007) Identification of functional modules using network topology and highthroughput data. BMC Systems Biology 1: 8.
[67]
Li M, Wu X, Wang J, Pan Y (2012) Towards the identification of protein complexes and functional modules by integrating ppi network and gene expression data. BMC Bioinformatics 13: 109.
[68]
Nepusz T, Yu H, Paccanaro A (2012) Detecting overlapping protein complexes in protein-protein interaction networks. Nat Meth 9: 471–472.
[69]
Keshava Prasad TS, Goel R, Kandasamy K, Keerthikumar S, Kumar S, et al. (2009) Human protein reference database?2009 update. Nucleic Acids Research 37: D767–D772.
[70]
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, et al. (2003) Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Research 13: 2498–2504.
