Original ArticleGene expression profiling, pathway analysis and subtype classification reveal molecular heterogeneity in hepatocellular carcinoma and suggest subtype specific therapeutic targets

•mRNA expression data shows high molecular heterogeneity in HCC patients.•A total of 570 genes were differentially expressed (479 up- and 91 down-regulated) in HCC.•Four HCC subtypes predicted were unique in terms of pathway enrichment and median survival.

A very low 5-year survival rate among hepatocellular carcinoma (HCC) patients is mainly due to lack of early stage diagnosis, distant metastasis and high risk of postoperative recurrence. Hence ascertaining novel biomarkers for early diagnosis and patient specific therapeutics is crucial and urgent. Here, we have performed a comprehensive analysis of the expression data of 423 HCC patients (373 tumors and 50 controls) downloaded from The Cancer Genome Atlas (TCGA) followed by pathway enrichment by gene ontology annotations, subtype classification and overall survival analysis. The differential gene expression analysis using non-parametric Wilcoxon test revealed a total of 479 up-regulated and 91 down-regulated genes in HCC compared to controls. The list of top differentially expressed genes mainly consists of tumor/cancer associated genes, such as AFP, THBS4, LCN2, GPC3, NUF2, etc. The genes over-expressed in HCC were mainly associated with cell cycle pathways. In total, 59 kinases associated genes were found over-expressed in HCC, including TTK, MELK, BUB1, NEK2, BUB1B, AURKB, PLK1, CDK1, PKMYT1, PBK, etc. Overall four distinct HCC subtypes were predicted using consensus clustering method. Each subtype was unique in terms of gene expression, pathway enrichment and median survival. Conclusively, this study has exposed a number of interesting genes which can be exploited in future as potential markers of HCC, diagnostic as well as prognostic and subtype classification may guide for improved and specific therapy.