Hepatocellular carcinoma (HCC) is a globally prevailing cancer with a low 5-year survival rate. Little is known about its intricate gene expression profile. Single-cell RNA sequencing is an indispensable tool to explore the genetic characteristics of HCC at a more detailed level. In this study, we profiled the gene expression of single cells from human HCC tumor and para-tumor tissues using the Smart-seq 2 sequencing method. Based on differentially expressed genes, we identified heterogeneous subclones in HCC tissues, including five HCC and two hepatocyte subclones. We then carried out hub-gene co-network and functional annotations analysis followed pseudo-time analysis with regulated transcriptional factor co-networks to determine HCC cellular trajectory. We found that MLX interacting protein like (MLXIPL) was commonly upregulated in the single cells and tissues and associated with a poor survival rate in HCC. Mechanistically, MLXIPL activation is crucial for promoting cell proliferation and inhibits cell apoptosis by accelerating cell glycolysis. Taken together, our work identifies the heterogeneity of HCC subclones, and suggests MLXIPL might be a promising therapeutic target for HCC.