Identification of Hub Genes in Liver Hepatocellular Carcinoma Based on Weighted Gene Co-expression Network Analysis

Abstract Liver hepatocellular carcinoma (LIHC) is a malignant cancer with high incidence and poor prognosis. To investigate the correlation between hub genes and progression of LIHC and to provided potential prognostic markers and therapy targets for LIHC. Our study mainly used The Cancer Genome Atlas (TCGA) LIHC database and the gene expression profiles of GSE54236 from the Gene Expression Omnibus (GEO) to explore the differential co-expression genes between LIHC and normal tissues. The differential co-expression genes were extracted by Weighted Gene Co-expression Network Analysis (WGCNA) and differential gene expression analysis methods. The Genetic Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were carried out to annotate the function of differential genes. Then the hub genes were validated using protein-protein interaction (PPI) network. And the expression level and prognostic analysis were performed. The probable associations between the expression of hub genes and both tumor purity and infiltration of immune cells were explored by TIMER. A total of 68 differential co-expression genes were extracted. These genes were mainly enriched in complement activation (biological process), collagen trimer (cellular component), carbohydrate binding and receptor ligand activity (molecular function) and cytokine − cytokine receptor interaction. Then we demonstrated that the 10 hub genes (CFP, CLEC1B, CLEC4G, CLEC4M, FCN2, FCN3, PAMR1 and TIMD4) were weakly expressed in LIHC tissues, the qRT-PCR results of clinical samples showed that six genes were significantly downregulated in LIHC patients compared with adjacent tissues. Worse overall survival (OS) and disease-free survival (DFS) in LIHC patients were associated with the lower expression of CFP, CLEC1B, FCN3 and TIMD4. Ten hub genes had positive association with tumor purity. CFP, CLEC1B, FCN3 and TIMD4 could serve as novel potential molecular targets for prognosis prediction in LIHC.