Hepatocellular carcinoma (HCC) is a dismal malignancy associated with multiple molecular changes. The purpose of this study was to identify the differentially expressed genes and analyze the biological processes related to HCC.
Methods and Study Design
Datasets of HCC were obtained from the NCBI Gene Expression Omnibus. Integrated analysis of differentially expressed genes was performed using the INMEX program. Then Gene Ontology enrichment analyses and pathway analysis were performed based on the Gene Ontology website and Kyoto Encyclopedia of Genes and Genomes. A protein-protein interaction network was constructed using the Cytoscape software; the netwerk served to find hub genes for HCC. Real-time RT-PCR was used to validate the microarray data for hub genes.
Results
We identified 273 genes that were differentially expressed in HCC. Gene Ontology enrichment analyses revealed response to cadmium ion, cellular response to cadmium ion, and cellular response to zinc ion for these genes. Pathway analysis showed that significant pathways included fatty acid metabolism, butanoate metabolism, and PPAR signaling pathway. The protein-protein interaction network indicated that CDH1, ECHS1, ACAA1, MT2A, and MYC were important genes which participated in many interactions. Experimental validation of the role of four upregulated genes (ECHS1, ACAA1, MT2A and MYC) in the progression of HCC was carried out.
Conclusions
Our study displayed genes that were consistently differentially expressed in HCC. The biological pathways and protein-protein interaction networks associated with those genes were also identified. We predicted that CDH1, ECHS1, ACAA1, MT2A, and MYC might be target genes for diagnosing HCC.
WhittakerS.MaraisR.ZhuA.X.The role of signaling pathways in the development and treatment of hepatocellular carcinoma. Oncogene.2010; 29(36): 4989–5005.
2.
AltekruseS.F.McGlynnK.A.ReichmanM.E.Hepatocellular carcinoma incidence, mortality, and survival trends in the United States from 1975 to 2005. J Clin Oncol.2009; 27(9): 1485–1491.
3.
El-SeragH.B.RudolphK.L.Hepatocellular carcinoma: epidemiology and molecular carcinogenesis. Gastroenterology.2007; 132(7): 2557–2576.
4.
StarleyB.Q.CalcagnoC.J.HarrisonS.A.Nonalcoholic fatty liver disease and hepatocellular carcinoma: a weighty connection. Hepatology.2010; 51(5): 1820–1832.
5.
El-SeragH.B.Epidemiology of viral hepatitis and hepatocellular carcinoma. Gastroenterology.2012; 142(6): 1264–1273, e1.
6.
HoshidaY.NijmanS.M.KobayashiM.Integrative transcriptome analysis reveals common molecular subclasses of human hepatocellular carcinoma. Cancer Res.2009; 69(18): 7385–7392.
7.
LiuW.PengY.TobinD.J.A new 12-gene diagnostic biomarker signature of melanoma revealed by integrated microarray analysis. Peer J.2013; 1: e49.
8.
KimT.H.ChoiS.J.LeeY.H.SongG.G.JiJ.D.Gene expression profile predicting the response to anti-TNF treatment in patients with rheumatoid arthritis; analysis of GEO datasets. Joint Bone Spine.2014; 81(4): 325–330.
9.
WrayC.J.KoT.C.TanF.K.Secondary use of existing public microarray data to predict outcome for hepatocellular carcinoma. J Surg Res.2014; 188(1): 137–142.
10.
DingB.XiY.GaoM.Gene expression profiles of entorhinal cortex in Alzheimer's disease. Am J Alzheimer's Dis Other Demen.2014; 29(6): 526–532.
11.
XuJ.ZhangC.SongW.Screening of differentially expressed genes associated with non-union skeletal fractures and analysis with a DNA microarray. Exp Ther Med.2014; 7(3): 609–614.
12.
WangH.B.JiangZ.B.LiM.Research on the Typical miRNA and Target Genes in Squamous Cell Carcinoma and Adenocarcinoma of Esophagus Cancer with DNA Microarray. Pathol Oncol Res.2014; 20(2): 245–252.
13.
YooB.K.EmdadL.SuZ.Z.Astrocyte elevated gene-1 regulates hepatocellular carcinoma development and progression. J Clin Invest.2009; 119(3): 465–477.
14.
BudhuA.ForguesM.YeQ-HPrediction of venous metastases, recurrence, and prognosis in hepatocellular carcinoma based on a unique immune response signature of the liver microenvironment. Cancer Cell.2006; 10(2): 99–111.
15.
JiaH-LYeQ-HQinL-XGene expression profiling reveals potential biomarkers of human hepatocellular carcinoma. Clin Cancer Res.2007; 13(4): 1133–1139.
BudhuA.RoesslerS.ZhaoX.Integrated metabolite and gene expression profiles identify lipid biomarkers associated with progression of hepatocellular carcinoma and patient outcomes. Gastroenterology.2013; 144(5): 1066–1075, e1.
18.
LiaoY.L.SunY.M.ChauG.Y.Identification of SOX4 target genes using phylogenetic footprinting-based prediction from expression microarrays suggests that overexpression of SOX4 potentiates metastasis in hepatocellular carcinoma. Oncogene.2008; 27(42): 5578–5589.
19.
HodoY.HondaM.TanakaA.Association of interleukin-28B genotype and hepatocellular carcinoma recurrence in patients with chronic hepatitis C. Clin Cancer Res.2013; 19(7): 1827–1837.
20.
XiaJ.FjellC.D.MayerM.L.PenaO.M.WishartD.S.HancockR.E.INMEX—a web-based tool for integrative meta-analysis of expression data. Nucleic Acids Res.2013; 41(Web Server issue): W63–70.
21.
IrizarryR.A.HobbsB.CollinF.Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics.2003; 4(2): 249–264.
22.
R Development Core Team. R: A language and environment for statistical computing.Vienna, Austria: the R Foundation for Statistical Computing; 2005.
23.
LiuX.FengR.DuL.The role of enoyl-CoA hydratase short chain 1 and peroxiredoxin 3 in PP2-induced apoptosis in human breast cancer MCF-7 cells. FEBS Lett.2010; 584(14): 3185–3192.
24.
BanerjeeS.GhoshalS.PorterT.D.Activation of AMP-kinase by policosanol requires peroxisomal metabolism. Lipids.2011; 46(4): 311–321.
25.
SongM.O.FreedmanJ.H.Expression of copper-responsive genes in HepG2 cells. Mol Cell Biochem.2005; 279(1-2): 141–147.
26.
ShenH.M.MichaelN.KimN.StorbU.The TATA binding protein, c-Myc and survivin genes are not somatically hypermutated, while Ig and BCL6 genes are hypermutated in human memory B cells. Int Immunol.2000; 12(7): 1085–1093.
27.
ChangY.WangS-XWangY-BECHS1 interacts with STAT3 and negatively regulates STAT3 signaling. FEBS Lett.2013; 587(6): 607–613.
28.
JanssenU.DavisE.M.Le BeauM.M.StoffelW.Human mitochondrial enoyl-CoA hydratase gene (ECHS1): structural organization and assignment to chromosome 10q26.2-q26.3. Genomics.1997; 40(3): 470–475.
29.
XiaoC-XYangX-NHuangQ-WECHS1 acts as a novel HBsAg-binding protein enhancing apoptosis through the mitochondrial pathway in HepG2 cells. Cancer Lett.2013; 330(1): 67–73.
30.
ZhuX.S.DaiY.C.ChenZ.X.XieJ.P.ZengW.LinY.Y.Enoyl-coenzyme A hydratase, short chain 1 involved in promoting the proliferation of hepatocellular carcinoma HepG2 cells. Tumor.2013; 33(5): 422–427.
31.
ZhuX.S.DaiY.C.ChenZ.X.Knockdown of ECHS1 protein expression inhibits hepatocellular carcinoma cell proliferation via suppression of Akt activity. Crit Rev Eukaryot Gene Expr.2013; 23(3): 275–282.
32.
CulkinD.J.GelderF.B.VickS.R.Evaluation of a new tumor marker for localized prostate cancer. Prostate.1992; 20(2): 117–122.
33.
PintoV.B.GelderF.B.MorrisD.M.Purification, partial characterization, and clinical evaluation of an adenocarcinoma-associated antigen. Cancer Res.1986; 46(12 Pt 1): 6520–6524.
34.
MarikarF.M.SunQ.M.HuaZ.C.Production of the polyclonal anti-human metallothionein 2A antibody with recombinant protein technology. Acta Biochim Biophys Sin (Shanghai).2006; 38(5): 305–309.
35.
CuiY.WangJ.ZhangX.ECRG2, a novel candidate of tumor suppressor gene in the esophageal carcinoma, interacts directly with metallothionein 2A and links to apoptosis. Biochem Biophys Res Commun.2003; 302(4): 904–915.
36.
ReineckeF.LevanetsO.OlivierY.Metallothionein isoform 2A expression is inducible and protects against ROS-mediated cell death in rotenone-treated HeLa cells. Biochem J.2006; 395(2): 405–415.
37.
WuK.J.GrandoriC.AmackerM.Direct activation of TERT transcription by c-MYC. Nat Genet.1999; 21(2): 220–224.
NesbitC.E.TersakJ.M.ProchownikE.V.MYC oncogenes and human neoplastic disease. Oncogene.1999; 18(19): 3004–3016.
40.
KimY.H.GirardL.GiacominiC.P.Combined microarray analysis of small cell lung cancer reveals altered apoptotic balance and distinct expression signatures of MYC family gene amplification. Oncogene.2006; 25(1): 130–138.
41.
HaywardD.G.FryA.M.Nek2 kinase in chromosome instability and cancer. Cancer Lett.2006; 237(2): 155–166.
42.
ZhangW.FletcherL.MuschelR.J.The role of Polo-like kinase 1 in the inhibition of centrosome separation after ionizing radiation. J Biol Chem.2005; 280(52): 42994–42999.
43.
FletcherL.CernigliaG.J.NiggE.A.YendT.J.MuschelR.J.Inhibition of centrosome separation after DNA damage: a role for Nek2. Radiat Res.2004; 162(2): 128–135.
44.
LiaoH.WinkfeinR.J.MackG.RattnerJ.B.YenT.J.CENP-F is a protein of the nuclear matrix that assembles onto kinetochores at late G2 and is rapidly degraded after mitosis. J Cell Biol.1995; 130(3): 507–518.
45.
MaL.ZhaoX.ZhuX.Mitosin/CENP-F in mitosis, transcriptional control, and differentiation. J Biomed Sci.2006; 13(2): 205–213.
46.
VergnolleM.A.TaylorS.S.Cenp-F links kinetochores to Ndel1/Nde1/Lis1/dynein microtubule motor complexes. Curr Biol.2007; 17(13): 1173–1179.
47.
O'BrienS.L.FaganA.FoxF.J.CENP-F expression is associated with poor prognosis and chromosomal instability in patients with primary breast cancer. Int J Cancer.2007; 120(7): 1434–1443.
48.
ShigeishiH.MizutaK.HigashikawaK.YonedaS.OnoS.KamataN.Correlation of CENP-F gene expression with tumor-proliferating activity in human salivary gland tumors. Oral Oncol.2005; 41(7): 716–722.
