Hepatocellular carcinoma (HCC) is a leading cause of global cancer mortality. However, little is known about the precise molecular mechanisms involved in tumor formation and pathogenesis. The primary goal of this study was to elucidate genome-wide molecular networks involved in development of HCC with multiple etiologies by exploring high quality microarray data. We undertook a comparative network analysis across 264 human microarray profiles monitoring transcript changes in healthy liver, liver cirrhosis, and HCC with viral and alcoholic etiologies. Gene co-expression profiling was used to derive a consensus gene relevance network of HCC progression that consisted of 798 genes and 2,012 links. The HCC interactome was further confirmed to be phenotype-specific and non-random. Additionally, we confirmed that co-expressed genes are more likely to share biological function, but not sub-cellular localization. Analysis of individual HCC genes revealed that they are topologically central in a human protein-protein interaction network. We used quantitative RT-PCR in a cohort of normal liver tissue (n = 8), hepatitis C virus (HCV)-induced chronic liver disease (n = 9), and HCC (n = 7) to validate co-expressions of several well-connected genes, namely ASPM, CDKN3, NEK2, RACGAP1, and TOP2A. We show that HCC is a heterogeneous disorder, underpinned by complex cross talk between immune response, cell cycle, and mRNA translation pathways. Our work provides a systems-wide resource for deeper understanding of molecular mechanisms in HCC progression and may be used further to define novel targets for efficient treatment or diagnosis of this disease.
References
[1]
Parkin DM, Bray F, Ferlay J, Pisani P (2005) Global cancer statistics, 2002. CA Cancer J Clin 55: 74–108.
[2]
Altekruse SF, McGlynn KA, Reichman ME (2009) Hepatocellular carcinoma incidence, mortality, and survival trends in the United States from 1975 to 2005. Journal of clinical oncology: official journal of the American Society of Clinical Oncology 27: 1485–1491.
[3]
Schutte K, Bornschein J, Malfertheiner P (2009) Hepatocellular carcinoma–epidemiological trends and risk factors. Dig Dis 27: 80–92.
[4]
Moller H, Mellemgaard A, Lindvig K, Olsen JH (1994) Obesity and cancer risk: a Danish record-linkage study. Eur J Cancer 30A: 344–350.
[5]
Seitz HK, Stickel F (2006) Risk factors and mechanisms of hepatocarcinogenesis with special emphasis on alcohol and oxidative stress. Biol Chem 387: 349–360.
[6]
El-Serag HB, Hampel H, Javadi F (2006) The association between diabetes and hepatocellular carcinoma: a systematic review of epidemiologic evidence. Clin Gastroenterol Hepatol 4: 369–380.
[7]
Hu L, Lau SH, Tzang CH, Wen JM, Wang W, et al. (2004) Association of Vimentin overexpression and hepatocellular carcinoma metastasis. Oncogene 23: 298–302.
[8]
Su WH, Chao CC, Yeh SH, Chen DS, Chen PJ, et al. (2007) OncoDB.HCC: an integrated oncogenomic database of hepatocellular carcinoma revealed aberrant cancer target genes and loci. Nucleic Acids Res 35: D727–731.
[9]
Weng G, Bhalla US, Iyengar R (1999) Complexity in biological signaling systems. Science 284: 92–96.
[10]
Hartwell LH, Hopfield JJ, Leibler S, Murray AW (1999) From molecular to modular cell biology. Nature 402: C47–52.
[11]
Drozdov I, Tsoka S, Ouzounis CA, Shah AM (2010) Genome-wide expression patterns in physiological cardiac hypertrophy. BMC Genomics 11: 557.
[12]
Aggarwal A, Guo DL, Hoshida Y, Yuen ST, Chu KM, et al. (2006) Topological and functional discovery in a gene coexpression meta-network of gastric cancer. Cancer Res 66: 232–241.
[13]
Guimera R, Nunes Amaral LA (2005) Functional cartography of complex metabolic networks. Nature 433: 895–900.
[14]
Zhang Y, Wang S, Li D, Zhnag J, Gu D, et al. (2011) A systems biology-based classifier for hepatocellular carcinoma diagnosis. PLoS One 6: e22426.
[15]
He B, Zhang H, Shi T (2011) A comprehensive analysis of the dynamic biological networks in HCV induced hepatocarcinogenesis. PLoS One 6: e18516.
[16]
Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, et al. (2004) A gene atlas of the mouse and human protein-encoding transcriptomes. Proceedings of the National Academy of Sciences of the United States of America 101: 6062–6067.
[17]
Jemal A, Bray F, Center MM, Ferlay J, Ward E, et al. (2011) Global cancer statistics. CA: a cancer journal for clinicians 61: 69–90.
[18]
Elo LL, Jarvenpaa H, Oresic M, Lahesmaa R, Aittokallio T (2007) Systematic construction of gene coexpression networks with applications to human T helper cell differentiation process. Bioinformatics 23: 2096–2103.
[19]
Prieto C, Risueno A, Fontanillo C, De las Rivas J (2008) Human gene coexpression landscape: confident network derived from tissue transcriptomic profiles. PLoS One 3: e3911.
[20]
Schadt EE (2009) Molecular networks as sensors and drivers of common human diseases. Nature 461: 218–223.
[21]
Barabasi AL, Gulbahce N, Loscalzo J (2011) Network medicine: a network-based approach to human disease. Nat Rev Genet 12: 56–68.
[22]
Clauset A, Shalizi CR, Newman MEJ (2009) Power-law distributions in empirical data. SIAM Review 51: 661–703.
[23]
Maslov S, Sneppen K (2002) Specificity and stability in topology of protein networks. Science 296: 910–913.
[24]
Forbes SA, Bindal N, Bamford S, Cole C, Kok CY, et al. (2011) COSMIC: mining complete cancer genomes in the Catalogue of Somatic Mutations in Cancer. Nucleic acids research 39: D945–950.
[25]
Callaway DS, Newman ME, Strogatz SH, Watts DJ (2000) Network robustness and fragility: percolation on random graphs. Physical review letters 85: 5468–5471.
[26]
Buldyrev SV, Parshani R, Paul G, Stanley HE, Havlin S (2010) Catastrophic cascade of failures in interdependent networks. Nature 464: 1025–1028.
[27]
Yu H, Kim PM, Sprecher E, Trifonov V, Gerstein M (2007) The importance of bottlenecks in protein networks: correlation with gene essentiality and expression dynamics. PLoS Comput Biol 3: e59.
[28]
Rives AW, Galitski T (2003) Modular organization of cellular networks. Proceedings of the National Academy of Sciences of the United States of America 100: 1128–1133.
[29]
Jeong H, Mason SP, Barabasi AL, Oltvai ZN (2001) Lethality and centrality in protein networks. Nature 411: 41–42.
[30]
von Mering C, Huynen M, Jaeggi D, Schmidt S, Bork P, et al. (2003) STRING: a database of predicted functional associations between proteins. Nucleic Acids Res 31: 258–261.
[31]
Hermjakob H, Montecchi-Palazzi L, Lewington C, Mudali S, Kerrien S, et al. (2004) IntAct: an open source molecular interaction database. Nucleic Acids Res 32: D452–455.
[32]
de Souto MC, Costa IG, de Araujo DS, Ludermir TB, Schliep A (2008) Clustering cancer gene expression data: a comparative study. BMC Bioinformatics 9: 497.
[33]
Berg J, Lassig M (2006) Cross-species analysis of biological networks by Bayesian alignment. Proceedings of the National Academy of Sciences of the United States of America 103: 10967–10972.
[34]
Lee HK, Hsu AK, Sajdak J, Qin J, Pavlidis P (2004) Coexpression analysis of human genes across many microarray data sets. Genome Res 14: 1085–1094.
[35]
Fliri AF, Loging WT, Volkmann RA (2009) Drug effects viewed from a signal transduction network perspective. J Med Chem 52: 8038–8046.
[36]
Florez AF, Park D, Bhak J, Kim BC, Kuchinsky A, et al. (2010) Protein network prediction and topological analysis in Leishmania major as a tool for drug target selection. BMC Bioinformatics 11: 484.
[37]
Barrenas F, Chavali S, Holme P, Mobini R, Benson M (2009) Network properties of complex human disease genes identified through genome-wide association studies. PLoS One 4: e8090.
[38]
Hoshida Y, Toffanin S, Lachenmayer A, Villanueva A, Minguez B, et al. (2010) Molecular classification and novel targets in hepatocellular carcinoma: recent advancements. Seminars in liver disease 30: 35–51.
[39]
Xulvi-Brunet R, Li H (2010) Co-expression networks: graph properties and topological comparisons. Bioinformatics 26: 205–214.
[40]
Friedel CC, Zimmer R (2007) Influence of degree correlations on network structure and stability in protein-protein interaction networks. BMC Bioinformatics 8: 297.
[41]
Cheung ST, Chen X, Guan XY, Wong SY, Tai LS, et al. (2002) Identify metastasis-associated genes in hepatocellular carcinoma through clonality delineation for multinodular tumor. Cancer Res 62: 4711–4721.
[42]
Durkin ME, Ullmannova V, Guan M, Popescu NC (2007) Deleted in liver cancer 3 (DLC-3), a novel Rho GTPase-activating protein, is downregulated in cancer and inhibits tumor cell growth. Oncogene 26: 4580–4589.
[43]
Wong CM, Yam JW, Ching YP, Yau TO, Leung TH, et al. (2005) Rho GTPase-activating protein deleted in liver cancer suppresses cell proliferation and invasion in hepatocellular carcinoma. Cancer Res 65: 8861–8868.
[44]
Wurmbach E, Chen YB, Khitrov G, Zhang W, Roayaie S, et al. (2007) Genome-wide molecular profiles of HCV-induced dysplasia and hepatocellular carcinoma. Hepatology 45: 938–947.
[45]
Yu G, Li F, Qin Y, Bo X, Wu Y, et al. (2010) GOSemSim: an R package for measuring semantic similarity among GO terms and gene products. Bioinformatics 26: 976–978.
[46]
Shi Z, Derow CK, Zhang B (2010) Co-expression module analysis reveals biological processes, genomic gain, and regulatory mechanisms associated with breast cancer progression. BMC Syst Biol 4: 74.
[47]
Forsyth CB, Tang Y, Shaikh M, Zhang L, Keshavarzian A (2010) Alcohol stimulates activation of Snail, epidermal growth factor receptor signaling, and biomarkers of epithelial-mesenchymal transition in colon and breast cancer cells. Alcohol Clin Exp Res 34: 19–31.
[48]
Sun J, Zhao Z (2010) A comparative study of cancer proteins in the human protein-protein interaction network. BMC Genomics 11: Suppl 3S5.
[49]
Zotenko E, Mestre J, O'Leary DP, Przytycka TM (2008) Why do hubs in the yeast protein interaction network tend to be essential: reexamining the connection between the network topology and essentiality. PLoS Comput Biol 4: e1000140.
[50]
Stuart JM, Segal E, Koller D, Kim SK (2003) A gene-coexpression network for global discovery of conserved genetic modules. Science 302: 249–255.
[51]
Masaki T, Shiratori Y, Rengifo W, Igarashi K, Matsumoto K, et al. (2000) Hepatocellular carcinoma cell cycle: study of Long-Evans cinnamon rats. Hepatology 32: 711–720.
[52]
Skawran B, Steinemann D, Becker T, Buurman R, Flik J, et al. (2008) Loss of 13q is associated with genes involved in cell cycle and proliferation in dedifferentiated hepatocellular carcinoma. Mod Pathol 21: 1479–1489.
[53]
Abdel-Hamid NM (2009) Premalignant variations in extracellular matrix composition in chemically induced hepatocellular carcinoma in rats. J Membr Biol 230: 155–162.
[54]
Chen RX, Xia YH, Xue TC, Ye SL (2010) Osteopontin promotes hepatocellular carcinoma invasion by up-regulating MMP-2 and uPA expression. Mol Biol Rep.
[55]
Braconi C, Meng F, Swenson E, Khrapenko L, Huang N, et al. (2009) Candidate therapeutic agents for hepatocellular cancer can be identified from phenotype-associated gene expression signatures. Cancer 115: 3738–3748.
[56]
Shuda M, Kondoh N, Tanaka K, Ryo A, Wakatsuki T, et al. (2000) Enhanced expression of translation factor mRNAs in hepatocellular carcinoma. Anticancer Res 20: 2489–2494.
[57]
El-Serag HB (2011) Hepatocellular carcinoma. The New England journal of medicine 365: 1118–1127.
[58]
Lin SY, Pan HW, Liu SH, Jeng YM, Hu FC, et al. (2008) ASPM is a novel marker for vascular invasion, early recurrence, and poor prognosis of hepatocellular carcinoma. Clin Cancer Res 14: 4814–4820.
[59]
Wong N, Yeo W, Wong WL, Wong NL, Chan KY, et al. (2009) TOP2A overexpression in hepatocellular carcinoma correlates with early age onset, shorter patients survival and chemoresistance. International journal of cancer Journal international du cancer 124: 644–652.
[60]
Minguez B, Hoshida Y, Villanueva A, Toffanin S, Cabellos L, et al. (2011) Gene-expression signature of vascular invasion in hepatocellular carcinoma. Journal of hepatology 55: 1325–1331.
[61]
Lee SW, Reimer CL, Fang L, Iruela-Arispe ML, Aaronson SA (2000) Overexpression of kinase-associated phosphatase (KAP) in breast and prostate cancer and inhibition of the transformed phenotype by antisense KAP expression. Molecular and cellular biology 20: 1723–1732.
[62]
Sun Y, Wang L, Jiang M, Huang J, Liu Z, et al. (2010) Secreted phosphoprotein 1 upstream invasive network construction and analysis of lung adenocarcinoma compared with human normal adjacent tissues by integrative biocomputation. Cell Biochem Biophys 56: 59–71.
[63]
Bae I, Rih JK, Kim HJ, Kang HJ, Haddad B, et al. (2005) BRCA1 regulates gene expression for orderly mitotic progression. Cell Cycle 4: 1641–1666.
[64]
Stone R 2nd, Sabichi AL, Gill J, Lee IL, Adegboyega P, et al. (2010) Identification of genes correlated with early-stage bladder cancer progression. Cancer Prev Res (Phila) 3: 776–786.
[65]
Jones WM, Chao AT, Zavortink M, Saint R, Bejsovec A (2010) Cytokinesis proteins Tum and Pav have a nuclear role in Wnt regulation. J Cell Sci 123: 2179–2189.
[66]
Mas VR, Maluf DG, Archer KJ, Yanek K, Williams B, et al. (2006) Differentially expressed genes between early and advanced hepatocellular carcinoma (HCC) as a potential tool for selecting liver transplant recipients. Mol Med 12: 97–104.
[67]
Small M, Xu X, Zhou J, Zhang J, Sun J, et al. (2008) Scale-free networks which are highly assortative but not small world. Phys Rev E Stat Nonlin Soft Matter Phys 77: 066112.
[68]
Pau Ni IB, Zakaria Z, Muhammad R, Abdullah N, Ibrahim N, et al. (2010) Gene expression patterns distinguish breast carcinomas from normal breast tissues: the Malaysian context. Pathol Res Pract 206: 223–228.
[69]
Ancona N, Maglietta R, Piepoli A, D'Addabbo A, Cotugno R, et al. (2006) On the statistical assessment of classifiers using DNA microarray data. BMC Bioinformatics 7: 387.
[70]
Koziol JA, Feng AC, Jia Z, Wang Y, Goodison S, et al. (2009) The wisdom of the commons: ensemble tree classifiers for prostate cancer prognosis. Bioinformatics 25: 54–60.
[71]
Parkinson H, Kapushesky M, Kolesnikov N, Rustici G, Shojatalab M, et al. (2009) ArrayExpress update–from an archive of functional genomics experiments to the atlas of gene expression. Nucleic Acids Res 37: D868–872.
[72]
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, et al. (2003) Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 4: 249–264.
[73]
Birney E, Andrews TD, Bevan P, Caccamo M, Chen Y, et al. (2004) An overview of Ensembl. Genome Res 14: 925–928.
[74]
Carter SL, Brechbuhler CM, Griffin M, Bond AT (2004) Gene co-expression network topology provides a framework for molecular characterization of cellular state. Bioinformatics 20: 2242–2250.
[75]
Newman ME (2006) Modularity and community structure in networks. Proc Natl Acad Sci U S A 103: 8577–8582.
[76]
Drozdov I, Ouzounis CA, Shah AM, Tsoka S (2011) Functional Genomics Assistant (FUGA): a toolbox for the analysis of complex biological networks. BMC research notes 4: 462.
[77]
Blondel VD, Guillaume J-L, Lambiotte R, Lefebvre E (2008) Fast unfolding of communities in large network. J Stat Mech P10008:
[78]
Enright AJ, Van Dongen S, Ouzounis CA (2002) An efficient algorithm for large-scale detection of protein families. Nucleic Acids Res 30: 1575–1584.
[79]
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, et al. (2000) Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nature genetics 25: 25–29.
[80]
Huang da W, Sherman BT, Lempicki RA (2009) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4: 44–57.
[81]
Wex T, Treiber G, Lendeckel U, Malfertheiner P (2003) A two-step method for the extraction of high-quality RNA from endoscopic biopsies. Clin Chem Lab Med 41: 1033–1037.
[82]
Boyault S, Rickman DS, de Reynies A, Balabaud C, Rebouissou S, et al. (2007) Transcriptome classification of HCC is related to gene alterations and to new therapeutic targets. Hepatology 45: 42–52.
