Identification and validation of real hub genes in hepatocellular carcinoma based on weighted gene co-expression network analysis

Article type: Research Article

Authors: Qiao, Yu^a | Yuan, Fahu^{a; *} | Wang, Xin^b | Hu, Jun^b | Mao, Yurong^a | Zhao, Zhigang^{b; *}

Affiliations: [a] School of Medicine, Jianghan University, Wuhan, Hubei, China | [b] Department of Spine Surgery, Wuhan Fourth Hospital, Wuhan, Hubei, China

Correspondence: [*] Corresponding authors: Zhigang Zhao, Department of Spine Surgery, Wuhan Fourth Hospital, Gutian 3rd Road, Jiefang Avenue, Wuhan, Hubei 430035, China. Tel.: +86 2768831614; E-mail: qiaoyu@whsdsyy.wecom.work. Fahu Yuan, School of Medicine, Jianghan University, No. 8, Sanjiaohu Road, Wuhan Economic and Technological Development Zone, Wuhan, Hubei 430056, China. Tel.: +86 2784225149; E-mail: yuanfh@jhun.edu.cn.

Abstract: BACKGROUND: Hepatocellular carcinoma (HCC) is one of the most common liver malignancies in the world. With highly invasive biological characteristics and a lack of obvious clinical manifestations, hepatocellular carcinoma usually has a poor prognosis and ranks fourth in cancer mortality. The etiology and exact molecular mechanism of primary hepatocellular carcinoma are still unclear. OBJECTIVE: This work aims to help identify biomarkers of early HCC diagnosis or prognosis based on weighted gene co-expression network analysis (WGCNA). METHODS: Expression data and clinical information of HTSeq-Counts were downloaded from The Cancer Genome Atlas (TCGA) database, and gene expression map GSE121248 was downloaded from Gene Expression Omnibus (GEO). By differentially expressed genes (DEGs) and weighted gene co-expression network analysis (WGCNA) searched for modules in the two databases that had the same effect on the biological characteristics of HCC, and extracted the module genes with the highest positive correlation with HCC from two databases, and finally obtained overlapping genes. Then, we performed functional enrichment analysis on the overlapping genes to understand their potential biological functions. The top ten hub genes were screened according to MCC through the string database and Cytoscape software and then subjected to survival analysis. RESULTS: High expression of CDK1, CCNA2, CDC20, KIF11, DLGAP5, KIF20A, ASPM, CEP55, and TPX2 was associated with poorer overall survival (OS) of HCC patients. The DFS curve was plotted using the online website GEPIA2. Finally, based on the enrichment of these genes in the KEGG pathway, real hub genes were screened out, which were CDK1, CCNA2, and CDC20 respectively. CONCLUSIONS: High expression of these three genes was negatively correlated with survival time in HCC, and the expression of CDK1, CCNA2, and CDC20 were significantly higher in tumor tissues of HCC patients than in normal liver tissues as verified again by the HPA database. All in all, this provides a new feasible target for early and accurate diagnosis of HCC, clinical diagnosis, treatment, and prognosis.

Keywords: Bioinformatics analysis, gene expression omnibus (GEO), hepatocellular carcinoma, the cancer genome atlas (TCGA), weighted gene co-expression network analysis (WGCNA)

DOI: 10.3233/CBM-220151

Journal: Cancer Biomarkers, vol. 35, no. 2, pp. 227-243, 2022