Identification of key genes in hepatocellular carcinoma and validation of the candidate gene, cdc25a, using gene set enrichment analysis, meta-analysis and cross-species comparison

The aim of the present study was to determine key pathways and genes involved in the pathogenesis of hepatocellular carcinoma (HCC) through bioinformatic analyses of HCC microarray data based on cross-species comparison. Microarray data of gene expression in HCC in different species were analyzed using gene set enrichment analysis (GSEA) and meta-analysis. Reverse transcription-quantitative polymerase chain reaction and western blotting were performed to determine the mRNA and protein expression levels of cdc25a, one of the identified candidate genes, in human, rat and tree shrew samples. The cell cycle pathway had the largest overlap between the GSEA and meta-analysis. Meta-analyses showed that 25 genes, including cdc25a, in the cell cycle pathway were differentially expressed. Cdc25a mRNA levels in HCC tissues were higher than those in normal liver tissues in humans, rats and tree shrews, and the expression level of cdc25a in HCC tissues was higher than in corresponding paraneoplastic tissues in humans and rats. In human HCC tissues, the cdc25a mRNA level was significantly correlated with clinical stage, portal vein tumor thrombosis and extrahepatic metastasis. Western blotting showed that, cdc25a protein levels were significantly upregulated in HCC tissues in humans, rats and tree shrews. In conclusion, GSEA and meta-analysis can be combined to identify key molecules and pathways involved in HCC. This study demonstrated that the cell cycle pathway and the cdc25a gene may be crucial in the pathogenesis and progression of HCC.