- Open Access
A transcriptome profile in hepatocellular carcinomas based on integrated analysis of microarray studies
Diagnostic Pathology volume 12, Article number: 4 (2017)
Despite new treatment options for hepatocellular carcinomas (HCC) recently, 5-year survival remains poor, ranging from 50 to 70%, which may attribute to the lack of early diagnostic biomarkers. Thus, developing new biomarkers for early diagnosis of HCC, is extremely urgent, aiming to decrease HCC-related deaths.
In the study, we conducted a comprehensive characterization of gene expression data of HCC based on a bioinformatics method. The results were confirmed by real time polymerase chain reaction (RT-PCR) and TCGA database to prove the credibility of this integrated analysis.
After integrating analysis of seven HCC gene expression datasets, 1167 differential expressed genes (DEGs) were identified. These genes mainly participated in the process of cell cycle, oocyte meiosis, and oocyte maturation mediated by progesterone. The results of experiments and TCGA database validation in 10 genes was in full accordance with findings in integrated analysis, indicating the high credibility of our integrated analysis of different gene expression datasets. ASPM, CCT3, and NEK2 was showed to be significantly associated with overall survival of HCC patients in TCGA database.
This method of integrated analysis may be a useful tool to minish the heterogeneity of individual microarray, hopefully outputs more accurate HCC transcriptome profiles based on large sample size, and explores some potential biomarkers and therapy targets for HCC.
Hepatocellular carcinoma (HCC) is one of the most frequently occurring malignant tumors worldwide . Risk factors of HCC are well recognized including gender, infection by hepatitis B virus or hepatitis C virus, cirrhosis, metabolism diseases, toxins, excess alcohol consumption, and smoking. HCC varies with wide geography, and is more prevalent in Asia, Africa, and southern Europe. It has been well defined that experiencing surgery for early HCC patients could achieve a higher curative resection rate (80.5%) , and finally have a better survival rate. However, patients with early HCC frequently manifest non-typical symptoms, hence, most of patients are diagnosed with advanced HCC when seeing a doctor, resulting in a low 5-year survival rate, ranging from 50 and 70% . Therefore, developing biomarkers for early diagnosis is being emphasized to prolong survival in patients with HCC.
Over the last decades, large efforts have been made to promote the early diagnosis of HCC. Alpha-fetoprotein (AFP) has been the most commonly used tumor biomarker in the liver, testicles, and ovaries . Highly sensitive and specific biomarkers need to be developed in HCC diagnosis. Glypican-3 (GPC3), a membrane-associated heparan sulfate proteoglycan, is up-regulated in HCC. Additionally, GPC3 involved in hippo pathway to exert its function in HCC cell proliferation. GPC may be applied in clinical practice as a novel diagnostic biomarker .
Additionally, some researchers have attempted to employ prognostic markers for predicting HCC recurrence. Villa E et al. detected whole genome microarray expression profiling of 161 HCC samples, and revealed that five-gene signature (ANGPT2, NETO2, NR4A1, DLL4, ESM1) was able to predict fast growth and worst survival of HCC patients . The exploration of prognostic markers may facilitate individualized therapies.
Recently, detection of genome-wide gene transcripts expressed in a given tissue type is becoming more and more feasible with advent of high-throughput technologies, such as microarray and RNA-seq. The application of microarray-based gene expression profiling has produced tremendous information, and provided mechanistic insights into the oncogenic process of HCC . However, although many microarray studies of HCC have been performed [8–11], each of study holds a somewhat different view due to the heterogeneity caused by the variety in clinical samples, platform, analytical approach, etc. Toward this end, an integrated analysis of seven HCC gene expression datasets was conducted to identify differential expressed genes (DEGs) between tumor and normal tissues, revealing a common biological thread that linked the disparate microarray studies. Ten genes were selected for further real time polymerase chain reaction (RT-PCR) and TCGA database validation, to prove the credibility of this integrated analysis. We expected our study would be of some value for the future diagnosis and therapy of HCC in clinic.
Eligible HCC gene expression datasets
The raw gene expression datasets of HCC and control samples were selected and downloaded in the Gene Expression Omnibus (GEO) database. The datasets meeting the following criteria were included: i) the expression profile of whole genome; ii) data from the tumor and tumor-adjacent normal liver tissues from HCC patients in clinic; iii) raw data or standardized data. Cirrhotic liver tissue sets, non-human sets, and integrated analysis of gene expression profiles were excluded.
Identification of HCC gene expression profile
We selected the Z-score transformation  method to normalize raw data from different platforms. The MATrixLABoratory (MATLAB) software was applied to calculate differently expressed probe sets between tumor and tumor-adjacent normal tissue, using gene specific t-test. The genes with FDR ≤ 0.05 were selected as the significantly differentially expressed genes (DEGs). Heat map analysis was conducted using the “heatmap.2” function of the R/Bioconductor package “gplots” .
Gene ontology (GO) of differentially expressed genes
The GO and pathway enrichment was analyzed via the online software GENECODIS to facilitate the interpretation of biological roles of DEGs (http://genecodis.cnb.csic.es) . The GO functions of the DEGs were determined according to different categories including biological process, molecular functions, and cellular components. In addition, pathway enrichment analysis was based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database.
Protein-protein interaction (PPI) network construction
In order to find candidate genes involved in the oncogenesis and hepatic dysfunction of HCC, PPI networks of significantly DEGs were constructed according to the data from Biological General Repository for Interaction Datasets (BioGRID) (http://thebiogrid.org/). Among the candidate genes, the PPI networks of the top 20 most significantly dysregulated genes were visualized via Cytoscape .
RNA Isolation and RT-PCR validation
Tumor and matched adjacent normal liver tissues which were obtained from five HCC patients in the current study, were frozen immediately after surgery, and were stored at −135 °C for RNA extraction. Frozen sections were made and evaluated independently by senior pathologists. The study was approved by the First Affiliated Hospital of PLA General Hospital ethnics committee. The ethics committee approved the relating screening, inspection, and data collection of the patients, and all subjects signed a written informed consent form. All works were undertaken following the provisions of the Declaration of Helsinki.
The whole RNA of liver tissue for each sample was extracted using RNAeasy Mini Kit (Qiagen, Valencia, CA) according to the manufacture’s protocol. Ten genes were randomly selected from the 20 most significantly DEGs. Primers for the ten genes were designed using PrimerPlex 2.61 (PREMIER Biosoft, Palo Alto, CA) (Additional file 1: Table S1). Expression levels of genes were screened by SYBR (Applied Biosystems/Life Technologies, Carlsbad, CA) in ABI 7500 Real Time PCR System (Applied Biosystems, Carlsbad CA). Relative gene expression was calculated with Data Assist Software version 3.0 (Applied Biosystems/Life Technologies) and human actin gene was used as a reference. The expression level of each gene was determined according to the method of 2-△△ct.
TCGA database validation of selected genes in HCC patients
Through the online validation tools, the expression status of selected genes in HCC were determined in TCGA database (https://genome-cancer.ucsc.edu/), assessing their mRNA expression patterns in HCC patients (N = 423) . The selected genes were also evaluated for the overall survival time of HCC patients in correlation with their expression pattern (http://cbioportal.org) in the TCGA database (N = 442) .
Candidate genes involved in the occurrence of HCC
Seven microarray datasets of HCC were identified according to the including criteria. Among of them, GSE17548, GSE33006, GSE17856, and GSE1481 didn’t contain the gene expression data of tumor-adjacent normal liver tissues. 267 HCC samples and 67 control samples were enrolled in the integrated analysis. The information of each microarray dataset was shown in Table 1. Based on microarray datasets available for integrated analysis, a total of 1167 DEGs were identified, among which, 628 genes were up-regulated and 539 genes were down-regulated. The detailed information of the 20 most significantly up-regulated or down-regulated genes were shown in Additional file 1: Table S2. The top 50 most significantly DEGs were displayed in a heat map across different HCC microarray datasets (Fig. 1).
GO and KEGG analysis of the candidate genes
Mitotic cell cycle (GO: 0000278, 4.71E-36) and cell division (GO: 0051301, 7.83E-26) was significantly enriched upon the category of GO biological progress, and protein binding (GO: 0005515, 3.88E-85) and nucleotide binding (GO: 0000166, 2.71E-41) was significantly enriched upon the category of GO molecular function. While for the category of GO cellular component, cytoplasm (GO: 0005737, 1.77E-90) and nucleus (GO: 0005634, 5.20E-59) was significantly enriched (Table 2). Based on KEGG database, the 1167 DEGs were involved in 99 signal pathways, including cell cycle, oocyte meiosis, oocyte maturation mediated by progesterone, pathways in cancer, p53 signaling pathways, production of phagosome, metabolism of fatty acid, cytokines-cell factor receptor interactions, prion diseases, etc. (Table 3).
PPI Network Constructions
For PPI networks of the 20 most significantly dyregulated genes, they consisted of 377 edges and 503 nodes. Three hub proteins were identified in this network, including CCT3 (121°), NDC80 (98°), and ASPM (93°) (Fig. 2).
Experimental and TCGA database validation of selected genes in HCC patients
Ten genes (ASPM, CAP2, CCT3, NEK2, SNRPE, CLEC4M, DCN, ECM1, RND3 and SPINT2) were randomly retrieved from the 20 most significantly up-regulated or down-regulated genes, respectively. After performing RT-PCR, the expression levels of selected 10 genes in clinical samples were identical with the results of the integrated analysis. For the ten genes, the mRNA expression was statistically different between tumor and matched adjacent normal liver tissues (Fig. 3; Additional file 1: Table S3) (P < 0.01). Furthermore, results of TCGA database validation indicated that these genes showed similar expression trends to those obtained from the integrated analysis (Fig. 4). Among the ten genes, only the ASPM, CCT3, and NEK2 showed significant association with overall survival time of HCC patients in TCGA database (P < 0.05) (Fig. 5).
It is generally accepted that the altered gene expression pattern of a cancer tissue should be associated with the initiation and maintenance of the malignant phenotype. Previous studies have identified several HCC gene expression profiles [18–21]. However, there wasn’t a common pattern among disparate studies for HCC. While in this study, we integrated different microarray studies to identify a precise gene expression profile for HCC with more statistical power supported by large sample size. In the current study, an integrated analysis of seven HCC microarray datasets was conducted, and showed that 1167 DEGs were identified, among which 628 genes were up-regulated and 539 genes were down-regulated. These genes mainly participated in the process of cell cycle, oocyte meiosis, and oocyte maturation mediated by progesterone.
In the current study, further annotation and PPI network analysis of the 20 most significant DEGs were conducted. Most of the 20 genes were involved in the pathways of cell cycle, cytokines-cell factor receptor interactions, and intracellular signaling cascades, and their involvements in HCC have also been reported [22–26]. The functions of the 20 genes were in accordance with the results of GO and KEGG analysis. Three genes, including CCT3, NDC80, and ASPM were proved to be highly connected in the PPI network. CCT3, a subunit of CCT cluster, plays a role in assisting the folding of proteins involved in important biological processes. CCT3 was found to display a significantly different gene expression level in HCC compared to adjacent non-malignant liver tissues, arising from the occurrence of the amplicon 1q21-q22 , which is consistent with our result of RT-PCR validation. In addition, other genes’ expression status detected by RT-PCR was totally in accordance with the result of integrated analysis, suggesting that the bioinformatics method of integrated analysis was credible.
ASPM was highly expressed in fetal tissues but lowly in most adult tissues. Our result and previous evidences  found that ASPM and NEK2 mRNA was over-expressed in HCC. Moreover, we found that ASPM, NEK and CCT3 over-expression present significant association with overall survival of HCC patients based on TCGA validation, predicting enhanced invasive/metastatic potential of HCC and higher risk of early tumor recurrence. ASPM, NEK and CCT3 may be applied as potential prognostic biomarkers for HCC. CAP2 overexpression was also discovered in our study, and CAP2 has been suggested as a candidate biomarker of HCC owing to elevated level in the serum of HCC patients .
Among the 10 most significantly down-regulated genes, DCN, an extracellular matrix proteoglycan, has important biological functions in growth, development and diseases. Loss of the decorin gene, which are known to interfere with cellular events of tumorigenesis mainly by blocking various receptor tyrosine kinases such as EGFR, Met, IGF-IR, PDGFR and VEGFR2, is permissive for tumorigenic growth of HCC with decreasing levels of the cyclin-dependent kinase inhibitor p21WAF1/CIP1, suggesting potential utilization of DCN as an antitumor agent in HCC . RND3 down-regulation in HCC patients has been reported by several studies [26, 30, 31], and may be a metastasis suppressor gene in HCC.
However, the expression patterns of four genes among the 20 most significant DEGs in the current study were inconsistent with or ignored in the previous studies, including TBCE, SPINT2, ECM1, and KZAN. The function of KZAN was not identified, whereas the other three genes were all comprehensively studied. In the current study, the inconsistent results might inspire their roles in the oncogenesis and development of HCC with some novel views.
SPINT2 encodes a transmembrane protein with two extracellular Kunitz domains that inhibits a variety of serine proteases. The protein product of SPINT2 inhibits HGF activator, which prevents the formation of active hepatocyte growth factor, has been taken as a putative tumor suppressor . Previous studies mainly focus on the methylation of SPINT2 in HCC instead of its expression [33, 34]. Nevertheless, we have found that the expression level of SPINT2 was significantly suppressed in HCC expression profiles. The pattern was consistent with that in cell renal cell carcinoma , which might indicate its potential application as a novel HCC suppressor.
ECM1 encodes a soluble protein that is involved in endochondral bone formation, angiogenesis, and tumor biology. It interacts with a variety of extracellular and structural proteins, contributing to the maintenance of skin integrity and homeostasis . The expression of ECM1 is reported to be significantly up-regulated in HCC patients , however, the current analyses of expression profiles showed that expression of ECM1 was suppressed in HCC patients and were confirmed using RT-PCR. The discrepancy revealed the complicated functions of ECM1 in the oncogenesis and development of HCC.
In short, the current study gave an explicit elucidation of dysregulated genes in HCC by the integrated analysis of microarray datasets in GEO database, the biological function of these genes was significantly enriched in cell cycle. The results of RT-PCR and TCGA validation were consistent with that of integrated analysis, indicating the high credibility of this integrated analysis method. In addition, our study showed that some genes could be potentially valuable in the clinical diagnosis (such as ASPM, NEK2 and CCT3) and anticancer therapy (such as DCN, RND3) for HCC. Our study improved the understanding of the transcriptome status of HCC, and might shed a light on the further investigation on the mechanisms of HCC.
Biological General Repository for Interaction Datasets
Differential expressed genes
Gene Expression Omnibus
Kyoto Encyclopedia of Genes and Genomes
Real time polymerase chain reaction
Bosch FX, Ribes J, Borras J. Epidemiology of primary liver cancer. Semin Liver Dis. 1999;19:271–85.
Zhou XD, Tang ZY, Yang BH, Lin ZY, Ma ZC, Ye SL, Wu ZQ, Fan J, Qin LX, Zheng BH. Experience of 1000 patients who underwent hepatectomy for small hepatocellular carcinoma. Cancer. 2001;91:1479–86.
Bruix J, Llovet JM. Major achievements in hepatocellular carcinoma. Lancet. 2009;373:614–6.
Chaiteerakij R, Addissie BD, Roberts LR. Update on biomarkers of hepatocellular carcinoma. Clin Gastroenterol Hepatol. 2015;13:237–45.
Miao HL, Pan ZJ, Lei CJ, Wen JY, Li MY, Liu ZK, Qiu ZD, Lin MZ, Chen NP, Chen M. Knockdown of GPC3 inhibits the proliferation of Huh7 hepatocellular carcinoma cells through down-regulation of YAP. J Cell Biochem. 2013;114:625–31.
Villa E, Critelli R, Lei B, Marzocchi G, Camma C, Giannelli G, Pontisso P, Cabibbo G, Enea M, Colopi S, et al. Neoangiogenesis-related genes are hallmarks of fast-growing hepatocellular carcinomas and worst survival. Results from a prospective study. Gut. 2016;65:861–9.
Rhodes DR, Chinnaiyan AM. Integrative analysis of the cancer transcriptome. Nat Genet. 2005;37(Suppl):S31–7.
Kim BY, Lee JG, Park S, Ahn JY, Ju YJ, Chung JH, Han CJ, Jeong SH, Yeom YI, Kim S, et al. Feature genes of hepatitis B virus-positive hepatocellular carcinoma, established by its molecular discrimination approach using prediction analysis of microarray. Biochim Biophys Acta. 2004;1739:50–61.
Liao YL, Sun YM, Chau GY, Chau YP, Lai TC, Wang JL, Horng JT, Hsiao M, Tsou AP. 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:5578–89.
Hodo Y, Honda M, Tanaka A, Nomura Y, Arai K, Yamashita T, Sakai Y, Mizukoshi E, Sakai A, Sasaki M, et al. Association of interleukin-28B genotype and hepatocellular carcinoma recurrence in patients with chronic hepatitis C. Clin Cancer Res. 2013;19:1827–37.
Chung EJ, Sung YK, Farooq M, Kim Y, Im S, Tak WY, Hwang YJ, Kim YI, Han HS, Kim JC, Kim MK. Gene expression profile analysis in human hepatocellular carcinoma by cDNA microarray. Mol Cells. 2002;14:382–7.
Yang Z, Chen Y, Fu Y, Yang Y, Zhang Y, Li D. Meta-analysis of differentially expressed genes in osteosarcoma based on gene expression data. BMC Med Genet. 2014;15:80.
Reimers M, Carey VJ. Bioconductor: an open source framework for bioinformatics and computational biology. Methods Enzymol. 2006;411:119–34.
Tabas-Madrid D, Nogales-Cadenas R, Pascual-Montano A. GeneCodis3: a non-redundant and modular enrichment analysis tool for functional genomics. Nucleic Acids Res. 2012;40:W478–83.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Goldman M, Craft B, Swatloski T, Cline M, Morozova O, Diekhans M, Haussler D, Zhu J. The UCSC Cancer Genomics Browser: update 2015. Nucleic Acids Res. 2015;43:D812–7.
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2:401–4.
Zhou X, Zheng R, Zhang H, He T. Pathway crosstalk analysis of microarray gene expression profile in human hepatocellular carcinoma. Pathol Oncol Res. 2015;21:563–9.
Melis M, Diaz G, Kleiner DE, Zamboni F, Kabat J, Lai J, Mogavero G, Tice A, Engle RE, Becker S, et al. Viral expression and molecular profiling in liver tissue versus microdissected hepatocytes in hepatitis B virus-associated hepatocellular carcinoma. J Transl Med. 2014;12:230.
Smalling RL, Delker DA, Zhang Y, Nieto N, McGuiness MS, Liu S, Friedman SL, Hagedorn CH, Wang L. Genome-wide transcriptome analysis identifies novel gene signatures implicated in human chronic liver disease. Am J Physiol Gastrointest Liver Physiol. 2013;305:G364–74.
Tsuchiya M, Parker JS, Kono H, Matsuda M, Fujii H, Rusyn I. Gene expression in nontumoral liver tissue and recurrence-free survival in hepatitis C virus-positive hepatocellular carcinoma. Mol Cancer. 2010;9:74.
Uemura M, Fujimura Y, Ko S, Matsumoto M, Nakajima Y, Fukui H. Pivotal role of ADAMTS13 function in liver diseases. Int J Hematol. 2010;91:20–9.
Lin SY, Pan HW, Liu SH, Jeng YM, Hu FC, Peng SY, Lai PL, Hsu HC. ASPM is a novel marker for vascular invasion, early recurrence, and poor prognosis of hepatocellular carcinoma. Clin Cancer Res. 2008;14:4814–20.
Chen H, Jia WD, Li JS, Wang W, Xu GL, Ma JL, Ren WH, Ge YS, Yu JH, Liu WB, et al. Extracellular matrix protein 1, a novel prognostic factor, is associated with metastatic potential of hepatocellular carcinoma. Med Oncol. 2011;28 Suppl 1:S318–25.
Tahara D, Nakanishi T, Akazawa S, Yamaguchi Y, Yamamoto H, Akashi M, Chikuba N, Okuno S, Maeda Y, Kusumoto Y, et al. Lecithin-cholesterol acyltransferase and lipid transfer protein activities in liver disease. Metabolism. 1993;42:19–23.
Grise F, Sena S, Bidaud-Meynard A, Baud J, Hiriart JB, Makki K, Dugot-Senant N, Staedel C, Bioulac-Sage P, Zucman-Rossi J, et al. Rnd3/RhoE Is down-regulated in hepatocellular carcinoma and controls cellular invasion. Hepatology. 2012;55:1766–75.
Wong N, Chan A, Lee SW, Lam E, To KF, Lai PB, Li XN, Liew CT, Johnson PJ. Positional mapping for amplified DNA sequences on 1q21-q22 in hepatocellular carcinoma indicates candidate genes over-expression. J Hepatol. 2003;38:298–306.
Sakamoto M, Mori T, Masugi Y, Effendi K, Rie I, Du W. Candidate molecular markers for histological diagnosis of early hepatocellular carcinoma. Intervirology. 2008;51 Suppl 1:42–5.
Horvath Z, Kovalszky I, Fullar A, Kiss K, Schaff Z, Iozzo RV, Baghy K. Decorin deficiency promotes hepatic carcinogenesis. Matrix Biol. 2014;35:194–205.
Guo X, Cao C, Sun J, Zhang D, Liu L, Wu D. miR-128a is up-regulated in hepatocellular carcinoma and promotes tumor cell proliferation by targeting RND3. Nan Fang Yi Ke Da Xue Xue Bao. 2014;34:1408–13.
Luo H, Dong Z, Zou J, Zeng Q, Wu D, Liu L. Down-regulation of RhoE is associated with progression and poor prognosis in hepatocellular carcinoma. J Surg Oncol. 2012;105:699–704.
Morris MR, Gentle D, Abdulrahman M, Maina EN, Gupta K, Banks RE, Wiesener MS, Kishida T, Yao M, Teh B, et al. Tumor suppressor activity and epigenetic inactivation of hepatocyte growth factor activator inhibitor type 2/SPINT2 in papillary and clear cell renal cell carcinoma. Cancer Res. 2005;65:4598–606.
Moribe T, Iizuka N, Miura T, Kimura N, Tamatsukuri S, Ishitsuka H, Hamamoto Y, Sakamoto K, Tamesa T, Oka M. Methylation of multiple genes as molecular markers for diagnosis of a small, well-differentiated hepatocellular carcinoma. Int J Cancer. 2009;125:388–97.
Iizuka N, Oka M, Sakaida I, Moribe T, Miura T, Kimura N, Tamatsukuri S, Ishitsuka H, Uchida K, Terai S, et al. Efficient detection of hepatocellular carcinoma by a hybrid blood test of epigenetic and classical protein markers. Clin Chim Acta. 2011;412:152–8.
Smits P, Ni J, Feng P, Wauters J, Van Hul W, Boutaibi ME, Dillon PJ, Merregaert J. The human extracellular matrix gene 1 (ECM1): genomic structure, cDNA cloning, expression pattern, and chromosomal localization. Genomics. 1997;45:487–95.
Availability of data and materials
XWH conceived of the study. WFF, WRL and QXL carried out data analysis, participated in its design and coordination and drafted the manuscript. LQW, HYX, YJW and ZHX collected the transcriptome data. WQ, LGH, ZFY and ZH performed the analysis of raw data. ZX, PXM, and BY carried out Real time PCR. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study was approved by the First Affiliated Hospital of PLA General Hospital ethics committee (Reference number: 201403166). Written informed consents were obtained from the enrolled patients.
Table S1. Detail information of primers. Table S2. Information of the most significantly up-regulated or down-regulated DEGs in HCC. Table S3. The expression values of 10 genes on all 5 HCC cases. (DOC 89 kb)
About this article
Cite this article
Wang, F., Wang, R., Li, Q. et al. A transcriptome profile in hepatocellular carcinomas based on integrated analysis of microarray studies. Diagn Pathol 12, 4 (2017). https://doi.org/10.1186/s13000-016-0596-x
- Hepatocellular cancer
- Differentially expressed gene
- Integrated analysis
- Expression profile
- Real time polymerase chain reaction
- TCGA validation