Abstract
PURPOSE: The aim of this study was to develop a comprehensive differential gene profile for hepatocellular carcinoma (HCC) patients treated with sorafenib. METHODS: The RNA sequencing data and miRNA sequencing data of 114 HCC patients treated with sorafenib only and 326 HCC control patients treated without any chemotherapeutic drugs were studied using differential expression, functional enrichment, and protein-protein interaction analysis. RESULTS: Compared with HCC patients without any chemotherapy drugs, the sorafenib-treated patients develop 66 differentially expressed genes (DEGs), including 12 upregulated genes and 54 downregulated genes. Additionally, three differentially expressed miRNAs (DEMs) also show specific expression pattern. With further analysis, five primary genes including HTR2C, TRH, AGTR2, MCHR2, and SLC6A2 as well as three miRNAs (hsa-miR-4445, hsa-miR-466, and hsa-miR-2114) have been suggested as the potential targets for sorafenib. The specific gene expression of five genes has been validated in clinical HCC patients by ELISA method. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses indicate several significantly enriched biological processes (BPs), cellular components (CCs), and molecular functions (MFs). CONCLUSION: Since sorafenib is becoming increasingly important in HCC treatment clinically, this study will help us understand the potential targets and eliminate diverse existing side effects of it as well as explore several potential clinical biomarkers with comprehensive analysis of differential gene expression profile.